Chapman–Enskog-metoden

På slutten av 1800-tallet er Boltzmann-ligningen som styrer dynamikken til det gassformige mediet i mikroskopisk skala og Euler- og Navier-Stokes-ligningene for det makroskopiske nivået kjent. Å gå fra den ene skalaen til den andre utgjør en del av Hilberts sjette problem. David Hilbert , forfatter av erklæringene om hovedproblemene som ble vurdert på slutten av 1800-tallet, legger grunnlaget for en metode i form av en utvikling som bærer hans navn (1912). Det vil ta noen år for Sydney Chapman og David Enskog å foreslå samtidig og uavhengig i 1916 og 1917 en løsning på dette problemet. [ 1 ][ 2 ]​ [ 3 ]​ Nylig har denne metoden blitt utvidet til tilfellet med en gass i termodynamisk ulikevekt , [ 4 ] er dette siste aspektet et veldig aktivt forskningsområde for tiden.

Chapman -Enskog- metoden er en forstyrrelsesmetode som består i å definere løsningen i form av en rekke distribusjonsfunksjoner i funksjoner av en "liten parameter" som kan sammenlignes med Knudsen-tallet . I rekkefølge null er Maxwell-Boltzmann-fordelingen og Euler-ligningene . Rekkefølge en lar oss kjenne uttrykket for varme- og momentumfluksene og det for transportkoeffisientene (diffusjonskoeffisienter på grunn av konsentrasjon, trykk- og temperaturgradienter, dynamiske og volumetriske viskositeter og konduktivitet). Av molekylære interaksjonspotensialer. Denne tilnærmingen gjør det mulig å finne Navier-Stokes-ligningene og rettferdiggjøre diffusjon ved termiske gradienter, ukjent på det tidspunktet Chapmans og Enskogs verk ble publisert. Denne metoden vil tillate at alle disse koeffisientene kan beregnes fra kunnskap om en av dem ved rekonstituering til et mål (vanligvis viskositeten) av et interaksjonspotensial som Lennard-Jones-potensialet .

Harold Grad har foreslått en alternativ tilnærming som består i å søke etter løsningen ved hjelp av metodene for momenter av distribusjonsfunksjonen (1949). Boltzmann-ligningen multipliseres med ( er den mikroskopiske hastigheten til Boltzmann-ligningen og tensorprodusenten og integrert i hastighet. I denne typen metode viser likningen i n° av tid (n+1)°. Derfor må vi lag en hypotese for å "lukke" systemet. Grad antar løsningen uttrykt av en avkortet serie av hermitepolynomer. David Levermore har nylig (1996) foreslått en lukking ved bruk av en generell egenskap: løsningen maksimerer entropien til det lukkede systemet, som er mediets partikler [ 5 ] Av beregningskodene basert på disse metodene forblir de laboratoriets eiendom fordi de ikke gir en nevneverdig gevinst når det gjelder validitetsegenskap (i form av antall Knudsen ) sammenlignet med standardkodene som løser Navier-Stokes-ligningene, som har blitt gjenstand for betydelig utvikling .

Mikroskopiske og makroskopiske evolusjonsligninger

Mikroskopisk nivå

Vi skriver ned den statistiske fordelingsfunksjonen til hastigheten. . i øyeblikket ved punktet for partikkelen (atomet eller molekylet) som tilhører arten . Det sannsynlige antallet partikler ved volumhastigheter på dette tidspunktet er . Den statistiske fordelingen er derfor målt i s 3 ⋅m -6 .

Boltzmanns ligning er skrevet

hvor , er kollisjonsoperatoren (eller kjernen), er en kvadratisk integraloperator som er beskrevet nedenfor, og gir effekten av kollisjonene som vil bli antatt elastiske for å forenkle problemet: det er ingen utveksling mellom de interne og translasjonsgradene av frihet, det er ingen kjemisk reaksjon. Derfor er volumet av viskositet som følge av denne typen utveksling utelukket.

Det er like mange utbredelser som det finnes arter i miljøet. Til hver tilsvarer en Boltzmann-ligning koblet til de andre av de andre medlemmene som representerer homogene ( ) eller heterogene ( ) kollisjoner.

Den elastiske kollisjonen

Hastighetene før interaksjonen i en treghetsreferanseramme er og . Indeksene representerer utydelig den samme arten eller to forskjellige arter. Disse hastighetene holder og så blir vi plassert i et system sentrert i tyngdepunktet som har en konstant hastighet på grunn av bevaring av momentum. I dette systemet, som derfor er galileisk, er starthastigheten til partikkelen den relative hastigheten . Ved symmetri kan det bekreftes at banen vil være inneholdt i planet som inneholder origo og . Vi velger et referansepunkt som: (se figur). I denne referansen er avviket . I henhold til den relative hastighetspåvirkningsparameteren og vekselvirkningspotensialet antas kun å avhenge av avstanden i de to samvirkende partiklene. Hvis denne hypotesen er streng for interaksjonen mellom to atomer, kan den anses som brukbar for to molekyler: potensialet er da et statistisk gjennomsnittspotensial.

Interaksjonsutgangsretningen er definert av: . Slutthastighetene kan beregnes ut fra følgende betraktninger:

  • bevaring av momentum i samhandlingen innebærer
  • relativ hastighet har en konstant modul på grunn av bevaring av energi, så eller

Hastighetene etter interaksjonen er derfor:

Videre fører bevaring av vinkelmomentum under interaksjonen til .

Systemet som beskriver kollisjonen er reversibelt. Liouvilles teorem (Hamiltonsk mekanikk) lar oss skrive

Kjernekollisjonen

Det sannsynlige antallet partikler som passerer gjennom området per tidsenhet er . De samhandler med det sannsynlige antallet partikler i elementærvolumet. Antall partikler som forsvinner fra statistikken per tidsenhet er med;

På samme måte teller vi antall partikler som dukker opp

Gitt relasjonene ovenfor for kollisjonen. kollisjonsoperatøren er skrevet

Denne ligningen er navngitt som Wang Chang og Uhlenbeck- ligningen .

En ekvivalent formulering kan gis som introduserer differensialtverrsnittet og er definert:

Og dermed

Makroskopisk nivå

Variablene

Boltzmann-ligningen beskriver utviklingen av partikler på makroskopisk nivå. For å beskrive hver av artene på makroskopisk nivå, er dette definert:

- partikkeltetthet
- tettheten
- gjennomsnittshastighet - indre energi

Vi kan da definere verdier for settet med arter

- den totale hastigheten (barysentrisk hastighet for alle partikler)
- diffusjonshastigheten ( kinematisk mengde )

Noen hjelpevariabler ( det er Avogadros tall )

- Volumfraksjon Fortell meg hva
- massefraksjonen Fortell meg hva
- den molare massen valeur moyenne
Strømmene

Mengdeflyten er per definisjon mengden hvor

dette defineres ved å se på det produktdyadiske produktet

- massestrøm
- trykkstrammeren som representerer strømmen av momentum. Den er symmetrisk av konstruksjon.
- delvis Trykk er definert fra trykktensor - kurven .
- varmefluks representerer flyten av indre energi.

De globale strømmene oppnås ganske enkelt ved å legge sammen med trykket :. Vi kan da definere en temperatur fra tilstandsligningen

Evolusjonsligninger

Å multiplisere hver av Boltzann-ligningene suksessivt med hver av kollisjonsinvariantene og integrere over hastighetene og, om nødvendig, over arten, gir de makroskopiske evolusjonsligningene kalt de kontrakterte produktligningene .

Alle andre medlemmer er null på grunn av bevaringslovene

Dermed får vi et evolusjonssystem , og der flyter og gjenstår å bli avklart .

Dimensjonsløs Boltzmann-ligning

Det antas å være et homogent medium (bare en art er til stede).

For å estimere bidraget til hvert ledd i Boltzmann-ligningen er det nødvendig å modifisere dette. For dette er følgende referansemengder definert:

  • temperaturen som vi utleder en referansehastighet fra som er lydhastigheten for en perfekt gass
  • den makroskopiske referanselengden Relatert til problemet med å definere en referansetid Relatert til forplantningen av akustiske bølger
  • referansetrykk , som en referansetetthet og mengden
  • en mikroskopisk referansehastighet som tilsvarer medianhastigheten til en partikkel
  • en annen mikroskopisk referansehastighet som tilsvarer den gjennomsnittlige relative hastigheten mellom to partikler
  • en mikroskopisk referanselengde som tilsvarer den gjennomsnittlige frie banen hvor er differensialtverrsnittet.

Hvis vi nå definerer de reduserte variablene , , , og , skrives Boltzmann-ligningen

hvor:

  • Strouhal-nummer som veier den tidsmessige variasjonen av systemet,
  • Knudsen tall hvis inverse vekt er kollisjonsleddet, så tendensen til å gå tilbake til likevekt. I praksis, for at likningene til kontinuumet skal være gyldige, er det nødvendig at , typisk.

Formulering av løsningen

Vi skriver løsningen som en serie med en parameter i samme størrelsesorden som Knidsen-tallet

Med respekt for bevaringslovene, må hver av utviklingsvilkårene også respektere dem. Hvorfra er begrensningene for løsningen.

vi bruker denne tilnærmingen i Boltzmann-ligningen og skiller begrepene som tilsvarer hver potens av .

For å bestille null

Det er ganske enkelt oppnådd

Denne ligningen er verifisert hvis alle leddene som utgjør den er null, spesielt

Hvilket det innebærer

enten

Det er en kollisjonsvariant. Derfor er det skrevet som en lineær kombinasjon av de kanoniske kollisjonsinvariantene

Ved å introdusere dette uttrykket i ligningene som definerer de mikroskopiske variablene, blir parametrene for denne utviklingen identifisert og Maxwells hastighetsfordelingslov er funnet.

med. . Diffusjonsflukser kansellerer ut, så vel som varmefluks . Trykktensoren reduseres til sitt spor hvor . De tilsvarende makroskopiske ligningene er Euler-ligningene .

For å bestille en

Rekkefølge én viser en Fredholm-integralligning for det ukjente

Den vanskelige oppløsningen til denne ligningen gjør at man kan gi de forskjellige ukjente mengdene av Enskog-ligningene som deretter kan assimileres til Navier–Stokes-ligningene .

Kringkastingsstrømmen

Det oppnås i form av et lineært system kalt Stefan-Maxwell- systemet

Der vi ser utseendet til den binære diffusjonskoeffisienten og "multikomponent termisk diffusjonskoeffisient", (faktisk et dimensjonsløst tall) (som ikke er diffusjonskoeffisienter og som kan være negative) ved:

For et medium med arter er rekkevidden til dette systemet siden . din løsning

Vi finner diffusjonsbegrepene etter konsentrasjon, trykk og temperaturgradient ( Soret-effekt ). Det er multikomponent diffusjonskoeffisienten, en løsning av et lineært system som involverer de binære koeffisientene. Dette systemet er også et rank-term system Løsningen er ikke unik og innebærer uavhengige vilkår. Vi velger vanligvis for symmetriens skyld. Dette valget er vilkårlig.

Det er flere omtrentlige løsninger av Stefan-Maxwell-systemet som lar et eksplisitt uttrykk for diffusjonsfluksen oppnås under en nærliggende form av Ficks lov , som ikke er nøyaktig for en binær blanding.

Trykktensoren

Kompresjonsstrammeren har en klassisk form

hvor er enhetstensoren Den viskøse spenningstensoren

Et tilleggsbegrep dukker opp i trykket når uelastiske interaksjoner tas i betraktning. Dens innflytelse er svak, til og med helt ubetydelig for gasser med lav tetthet. [ 4 ]

Det bemerkes at Stokes-hypotesen naturlig er begrunnet med denne tilnærmingen.

Strømmen av varme

er gitt av

Dette er den termiske ledningsevnen. Det siste leddet i ligningen er følgen av Soret-effekten og kalles Dufour-effekten.

Transportkoeffisienter

Transportkoeffisientene uttrykkes i form av lineære systemer som involverer mengder av typen , som er utviklet i Sonine-Laguerre polynomer. Ekspansjonskoeffisientene uttrykkes i kollisjonsintegralfunksjoner. I praksis er vi fornøyd med den første ordren for ekspansjonen og kollisjonsintegralene er funksjoner av temperatur som er tabellert av forskjellige forfattere. I tillegg er det omtrentlige løsninger av de lineære systemene som gir de ulike koeffisientene i eksplisitt form.

Distribusjonsfunksjonen

Fordelingsfunksjonen er

hvor det oppstår i diffusjonsstrømmen

Denne fordelingsfunksjonen er nødvendig for beregningen av Knudsen-laget som gir veggbetingelsene for Navier-Stokes-ligningen.

For å bestille to

Her får vi igjen en Fredholm-integral for det ukjente

David Burnett foreslo i 1935 en løsning på denne ligningen. Dette har ulempen ved ikke å respektere teorien [ 6 ]​ H. Det ser ut til at stigningen i rekkefølge utgjør en blindvei, alle variantene som er foreslått til dags dato løser ikke dette problemet. [ 7 ]

Se også

Referanser

  1. (i) Sydney Chapman og Thomas George Cowling , The Mathematical Theory of Non-uniform Gases , Cambridge University Press ,1991 ( ISBN  0-521-40844-X )
  2. (i) Joseph Oakland Hirschfelder , Charles Francis Curtiss og Robert Byron Bird , Molecular Theory of Gases and Liquids , John Wiley and Sons ,1966 ( ISBN  978-0-471-40065-3 )
  3. (i) Gilberto Medeiros Kremer , «  The Methods of Chapman-Enskog and Grad and Applications  » , RTO-EN-AVT 194 , 2011 [1]
  4. a b Raymond Brun , Introduksjon til dynamique des gaz réactifs , Cépaduès ,2015 ( ISBN  978-2-364-93190-9 )
  5. (i) Charles David Levermore , "  Moment Closure Hierarchies for Kinetic Theories  " , Journal of Statistical Physics , vol.  83,nitten nittiseks ( les på nettet )
  6. Byung Chan Eu, Nonequilibrium Statistical Mechanics. Ensemblemetode. , springer ,1998 ( ISBN  978-90-481-5007-6 )
  7. Cédric Villani , Limites hydrodynamiques de l'équation de Boltzmann , koll.  "Seminaire Bourbaki" ( nr. 893  [ 2] ),2001