Stefan problem

Inom matematiken och dess tillämpningar, särskilt för fasövergångar i materia, är ett Stefan-problem en speciell typ av gränsvärdesproblem för ett system av partiella differentialekvationer (PDE), där gränsen mellan faserna kan röra sig med tiden. Det klassiska Stefan-problemet syftar till att beskriva utvecklingen av gränsen mellan två faser av ett material som genomgår en fasförändring , till exempel smältningen av ett fast ämne, såsom is till vatten . Detta åstadkoms genom att lösa värmeekvationer i båda regionerna, med förbehåll för givna gräns- och initialvillkor. Vid gränssnittet mellan faserna (i det klassiska problemet) ställs temperaturen in på fasändringstemperaturen. krävs ytterligare en ekvation, Stefan-villkoret . Detta är en energibalans som definierar positionen för det rörliga gränssnittet. Observera att denna utvecklande gräns är en okänd (hyper) yta ; därför är Stefan-problem exempel på problem med fria gränser .

Analoga problem uppstår till exempel i studiet av poröst medieflöde, matematisk ekonomi och kristalltillväxt från monomerlösningar.

Historisk anteckning

Problemet är uppkallat efter Josef Stefan (Jožef Stefan), den slovenske fysikern som introducerade den allmänna klassen av sådana problem runt 1890 i en serie om fyra artiklar om frysning av marken och bildandet av havsis . Men ungefär 60 år tidigare, 1831, hade ett motsvarande problem, angående bildandet av jordskorpan, studerats av Lamé och Clapeyron . Stefans problem medger en likhetslösning , detta kallas ofta för Neumann -lösningen, som påstås presenterades i en serie föreläsningar i början av 1860-talet.

En omfattande beskrivning av Stefans problems historia finns hos Rubinstein.

Premisser till den matematiska beskrivningen

Ur en matematisk synvinkel är faserna bara regioner där lösningarna för den underliggande PDE är kontinuerliga och differentierbara upp till PDE:s ordning. I fysikaliska problem representerar sådana lösningar egenskaper hos mediet för varje fas. De rörliga gränserna (eller gränssnitten ) är oändligt mycket tunna ytor som separerar intilliggande faser; därför kan lösningarna för den underliggande PDE och dess derivat drabbas av diskontinuiteter över gränssnitt.

De underliggande PDE:erna är inte giltiga vid fasändringsgränssnitten; därför behövs ytterligare ett villkor – Stefan-villkoret – för att erhålla stängning . Stefan-villkoret uttrycker den lokala hastigheten för en rörlig gräns, som en funktion av kvantiteter utvärderade på vardera sidan av fasgränsen, och härleds vanligtvis från en fysisk begränsning. I problem med värmeöverföring med fasförändring, till exempel, dikterar bevarande av energi att diskontinuiteten av värmeflödet vid gränsen måste förklaras av hastigheten för latent värmefrigöring (som är proportionell mot den lokala hastigheten för gränsytan).

Ekvationens regelbundenhet har studerats huvudsakligen av Luis Caffarelli och ytterligare förfinats genom arbete av Alessio Figalli , Xavier Ros-Oton och Joaquim Serra

Matematisk formulering

Det endimensionella enfasiga Stefan-problemet

Enfas Stefan-problemet bygger på ett antagande om att en av de materiella faserna kan försummas. Vanligtvis uppnås detta genom att anta att en fas är vid fasändringstemperaturen och följaktligen leder varje variation från detta till en fasändring. Detta är en matematiskt bekväm approximation, som förenklar analysen samtidigt som den visar de väsentliga idéerna bakom processen. En ytterligare standardförenkling är att arbeta i icke-dimensionellt format, så att temperaturen vid gränssnittet kan ställas in på noll och fjärrfältsvärden till eller .

Betrakta ett semi-oändligt endimensionellt isblock initialt vid smälttemperatur för . Den mest välkända formen av Stefan-problem involverar smältning via en pålagd konstant temperatur vid den vänstra gränsen, vilket lämnar en region upptagen av vatten. Det smälta djupet, betecknat med , är en okänd funktion av tiden. Stefan-problemet definieras av

  • Värmeekvationen:
  • En fast temperatur, över smälttemperaturen, på den vänstra gränsen:
  • Gränssnittet vid smälttemperaturen är inställt på
  • Stefan-villkoret: beta Stefan-talet, förhållandet mellan latent till specifik förnuftig värme (där specifik anger att den är dividerad med massan). Observera att denna definition följer naturligt av icke-dimensionaliseringen och används i många texter, men den kan också definieras som motsatsen till detta .
  • Den initiala temperaturfördelningen:
  • Det inledande djupet för det smälta isblocket:
Neumann-lösningen, erhållen genom att använda självliknande variabler, indikerar att gränsens position ges av där uppfyller den transcendentala ekvationen
Temperaturen i vätskan ges då av

Ansökningar

Förutom modellering av smältning av fasta ämnen, används Stefan-problemet också som en modell för det asymptotiska beteendet (i tiden) av mer komplexa problem. Till exempel använder Pego matchade asymptotiska expansioner för att bevisa att Cahn-Hilliards lösningar för fasseparationsproblem beter sig som lösningar på ett icke-linjärt Stefan-problem på en mellanliggande tidsskala. Dessutom är lösningen av Cahn–Hilliard-ekvationen för en binär blandning rimligen jämförbar med lösningen av ett Stefan-problem. I denna jämförelse löstes Stefan-problemet med en front-tracking, moving-mesh-metod med homogena Neumann-randvillkor vid den yttre gränsen. Stefan-problem kan också användas för att beskriva andra fastransformationer än fast-fluid eller fluid-fluid.

Stefan-problemet har också en rik omvänd teori; i sådana problem är mätningsdjupet (eller kurvan eller hyperytan ) s det kända datumet och problemet är att hitta u eller f .

Avancerade former av Stefanproblem

Det klassiska Stefan-problemet handlar om stationära material med konstanta termofysiska egenskaper (vanligtvis oberoende av fas), en konstant fasändringstemperatur och, i exemplet ovan, en momentan växling från den initiala temperaturen till ett distinkt värde vid gränsen. I praktiken kan de termiska egenskaperna variera och gör det specifikt alltid när fasen ändras. Hoppet i densitet vid fasförändring inducerar en flytande rörelse: den resulterande kinetiska energin figurerar inte i standardenergibalansen. Med en momentan temperaturomkopplare är den initiala vätskehastigheten oändlig, vilket resulterar i en initial oändlig kinetisk energi. I själva verket är vätskeskiktet ofta i rörelse, vilket kräver advektion eller konvektionstermer i värmeekvationen . Smälttemperaturen kan variera med storlek, krökning eller hastighet på gränsytan. Det är omöjligt att omedelbart växla temperaturer och då svårt att hålla en exakt fast gränstemperatur. Vidare, på nanoskala kanske temperaturen inte ens följer Fouriers lag.

Ett antal av dessa problem har angripits under de senaste åren för en mängd olika fysiska tillämpningar. Vid stelning av underkylda smältor kan en analys där fasändringstemperaturen beror på gränssnittshastigheten hittas i Font et al . Nanoskala stelning, med variabel fasförändringstemperatur och energi/densitetseffekter modelleras in. Stelning med flöde i en kanal har studerats, i samband med lava och mikrokanaler, eller med en fri yta i samband med att vatten fryser över ett islager . En generell modell som inkluderar olika egenskaper i varje fas, variabel fasförändringstemperatur och värmeekvationer baserade på antingen Fouriers lag eller Guyer-Krumhansl-ekvationen analyseras i.

Se även

Anteckningar

Historiska referenser

Vetenskapliga och allmänna referenser

externa länkar