Self-similar structure of non-isothermal variable-density mushy Stefan problems and an improved low-Mach enthalpy method

Technische Samenvatting: Zelfgelijkende structuur van niet-isotherme variabele-densiteit mushy Stefan-problemen en een verbeterde low-Mach enthalpische methode

Probleemstelling
De enthalpische methode is een breed toegepaste benadering voor het simuleren van faseveranderingsproblemen (smelten en stollingsprocessen) op vaste roosters zonder expliciete interface-tracking. Het modelleert doorgaans de faseovergang over een eindig temperatuurinterval, $\Delta T = T_{liq} - T_{sol}$, wat een "mushy" regio (slushige zone) creëert. Hoewel de methode over het algemeen wordt aangenomen te convergeren naar het klassieke scherpe-interface Stefan-probleem wanneer $\Delta T \to 0$, is deze aanname grotendeels gevalideerd onder de vereenvoudigende conditie van gelijke vaste en vloeibare dichtheden ($\rho_S = \rho_L$).

In praktische toepassingen met significante volumeveranderingen (bijv. koken, condensatie of specifieke metaalstollingsprocessen) kan de dichtheidsratio $\rho_L/\rho_S$ aanzienlijk afwijken van eenheid. In dergelijke gevallen induceert faseverandering vloeistofbeweging door dichtheidsvariaties, waardoor het thermische probleem wordt gekoppeld aan massa- en momentumbehoud. De bestaande literatuur mist een systematische investigatie naar de nauwkeurigheid en convergentie-eigenschappen van enthalpische methoden voor problemen met hoge dichtheidsratio's en een eindige $\Delta T$. Bovطendien kan eerdere validatie tegen de scherpe-interface limiet ontoereikend zijn voor het beoordelen van methoden die ontworpen zijn om dichtheidsverandering-geïnduceerde stroming te behandelen, aange het mushy Stefan-probleem (eindige $\Delta T$) fundamenteel verschilt van het klassieke Stefan-probleem (infinitesimale $\Delta T$).

Methodologie
De auteurs pakken deze hiaten aan via een combinatie van theoretische analyse en numerieke ontwikkeling:

1. Theoretische formulering van het mushy Stefan-probleem:
   De auteurs formuleren een niet-isotherm, variabele-densiteit mushy Stefan-probleem. Ze demonstreren dat de regerende vergelijkingen een zelfgelijkende structuur bezitten, waarbij de partiële differentiaalvergelijkingen worden gereduceerd tot een randwaardeprobleem met twee onbekende gelijktijdigheidsparameters ($\lambda_s$ en $\lambda_l$) die de posities van de solidus- en liquidusfronten definiëren.

   • Gelijk-densiteit geval: Een gesloten vorm analytische gelijktijdige oplossing wordt afgeleid.
   • Ongelijk-densiteit geval: Het probleem reduceert tot een gekoppeld nietlineair systeem dat de temperatuur- en snelheidvelden in de mushy- en vloeibare regio's regelt. Dit systeem wordt numeriek opgelost met een randwaardesolver (MATLAB's bvp5c) met een continuatiestrategie om extreme dichtheidsratio's (tot $\rho_L/\rho_S = 2700$) en kleine $\Delta T$ aan te kunnen.

2. Analyse van enthalpie-definities:
   De studie onderzoekt twee definities voor specifieke enthalpie in de mushy-regio. De "oorspronkelijke" definitie leidt tot een nietlineaire relatie tussen vloeistoffractie ($\phi$) en temperatuur ($T$) wanneer dichtheden verschillen, terwijl de "verbeterde" definitie een lineaire $\phi$-$T$ relatie oplevert. De auteurs tonen aan dat de oorspronkelijke formulering spuure instabiliteiten vertoont bij hoge dichtheidsratio's, terwijl de verbeterde formulering stabiel blijft.

3. Verbeterde low-Mach enthalpische methode:
   Op basis van de theoretische inzichten stellen de auteurs een verbeterde low-Mach enthalpische methode voor. Belangrijke kenmerken zijn:

   • Een herziene definitie van specifieke enthalpie die een lineaire vloeistoffractie–temperatuur relatie waarborgt.
   • Een kinematische snelheid-divergentie beperking afgeleid direct uit behoud van massa om dichtheidsverandering-geïnduceerde stroming te accounten.
   • Implementatie binnen de IBAMR-bibliotheek met gebruikmakend van een gestaffeld Cartesiaans rooster, Newton-iteraties voor de nietlineaire enthalpievergelijking, en een projectiemethode met een preconditioner om de stijve Carman–Kozeny dragkracht in de vaste regio te behandelen.

Belangrijkste Resultaten

• Zelfgelijkenis en Convergentie: De analyse bevestigt dat de enthalpische methode het mushy Stefan-probleem oplost, en niet het klassieke Stefan-probleem. De breedte van de mushy-regio groeit proportioneel met $\sqrt{t}$, terwijl het klassieke Stefan-front infinitesimaal dun blijft.
• Gevoeligheid voor Dichtheidsratio:
  • Voor gelijke en lage dichtheidsratio's ($\rho_L/\rho_S \approx 1$ tot $2$) benadert de mushy Stefan-oplossing de scherpe-interface oplossing nauwkeurig over een breed bereik van $\Delta T$-waarden. Bijgevolg zijn gelijke-densiteit benchmarks niet streng genoeg om methoden voor dichtheidsverandering-geïnduceerde stroming te valideren.
  • Voor grote dichtheidsratio's ($\rho_L/\rho_S \gg 1$) is de positie van het fasefront zeer gevoelig voor $\Delta T$. De mushy-oplossing benadert de scherpe-interface limiet alleen wanneer $\Delta T$ voldoende klein is.
• Prestaties van de Methode:
  • De verbeterde low-Mach enthalpische methode vertoont een uitstekende overeenstemming met de afgeleide mushy Stefan gelijktijdige oplossingen over een breed bereik van $\Delta T$-waarden en dichtheidsratio's (tot $\rho_L/\rho_S \approx 540$ in verificatiegevallen).
  • In contrast hiermee produceert de oorspronkelijke formulering interface-instabiliteiten en onfysische vortale structuren in de vloeibare regio bij grote dichtheidsratio's.
  • Convergentiesnelheden: De verbeterde methode bereikt ruimtelijke-temporele convergentiesnelheden tussen eerste en tweede orde voor zowel de fasefront-locatie als de temperatuurvelden, zelfs bij dichtheidsratio's van $O(3)$.

Betekenis en Claims
Het artikel claimt een langdurig gat in de enthalpische methode-literatuur te hebben gevuld door een rigoureuze theoretische benchmark (de gelijktijdige oplossing van het niet-isotherme variabele-densiteit mushy Stefan-probleem) te bieden voor de validatie van numerieke methoden bij een eindige $\Delta T$.

De auteurs stellen dat:

1. Validiteit van de Benchmark: Grote dichtheidsratio-gevallen bieden een aanzienlijk striktere validatiebenchmark voor moderne enthalpische methoden dan gevallen met gelijke dichtheden, aange de laatstgenoemde de gevoeligheid van de oplossing voor het faseveranderingsinterval niet kunnen vangen.
2. Methodologische Verbetering: De voorgestelde verbeteringen aan de low-Mach enthalpische methode (specifiek de herziene enthalpie-definitie en de lineaire $\phi$-$T$ relatie) essentieel zijn voor stabiliteit en nauwkeurigheid bij het simuleren van faseverandering met substantiële dichtheidscontrasten.
3. Verlengde Toepasbaarheid: Deze theoretische en numerieke ontwikkelingen breiden de toepasbaarheid van enthalpische methoden uit voorbij traditionele metaalstollingsprocessen (waar dichtheidsratio's nabij de eenheid liggen) naar processen die veel grotere dichtheidsratio's impliceren, zoals koken en condensatie.

Het werk stelt geen nieuwe experimentele toepassingen voor, maar vestigt een robuust numeriek kader en een theoretische basis voor het uitbreiden van enthalpie-gebaseerde simulaties naar een bredere klasse van variabele-densiteit faseveranderingsproblemen.