A volume penalization method for solving conjugate scalar transport with interfacial jump conditions

Technische Samenvatting: Een Volume-penalisatiemethode voor het Oplossen van Conjunctief Scalair Transport met Interfaciale Sprongcondities

Probleemstelling
Conjunctief scalair transport (bijv. warmte- of massatransport) met interfaciale sprongcondities komt veel voor in thermische en chemische processen, zoals warmteoverdracht in microkanalen, chemische reacties aan interfaces en faseveranderingssystemen. In deze scenario's bestaan er vaak discontinuïteiten in de scalaire waarde zelf (bijv. temperatuur- of concentratiesprongen) en/of in de flux van deze waarde over de vloeistof-vast interface. Hoewel body-fitted mesh (BFM) simulaties deze condities kunnen afhandelen door de afzonderlijke fase-domeinen expliciet te koppelen, vereisen zij omslachtige gridgeneratie voor complexe geometrieën. Immersed Boundary Methods (IBM's), in het bijzonder de Volume Penalization Method (VPM), bieden een manier om willekeurige geometrieën in Cartesiaanse grids in te bedden zonder complexe meshing. Bestaande VPM-schema's hebben echter moeite met Neumann-randvoorwaarden en algemene interfaciale sprongcondities. Conventionele benaderingen genereren vaak niet-fysische scalaire velden binnen het vaste gebied en slagen er niet in om fluxsprongen of scalaire discontinuïteiten nauwkeurig af te dwingen zonder aanvullende, complexe koppelingsprocedures.

Methodologie
De auteurs stellen een nieuwe behandeling voor binnen het VPM-raamwerk om conjunctief scalair transport met algemene interfaciale sprongcondities op te lossen. De methodologie verloopt in drie hoofdfasen:

1. Level-Set Representatie: Een willekeurige vloeistof-vast interface wordt gerepresenteerd door een level-set functie ($\phi_0$), die wordt omgezet in een fase-indicator ($\phi$) die geleidelijk overgaat van 0 (vloeistof) naar 1 (vast) over een dunne interfaciale regio van dikte $2\delta$.
2. Behandeling van Neumann Randvoorwaarden: Voor een enkelvoudig-fasig advectie-diffusieprobleem met een Neumann-conditie, leiden de auteurs een bronterm af in divergentievorm, $\nabla \cdot (q_w \phi \mathbf{n})$, die aan de sturende vergelijking moet worden toegevoegd. In tegenstelling tot eerdere schema's die afgeleiden afdwingen via een penalisatieterm door het gehele vaste volume, is deze nieuwe term alleen niet-nul binnen de interfaciale regio. Dit zorgt ervoor dat de voorgeschreven scalaire flux $q_w$ precies wordt toegepast, terwijl de moleculaire diffusiviteit nul blijft binnen het vaste deel, waardoor niet-fysische diffusie in de vaste fase wordt voorkomen.
3. Generalisatie naar Conjunctief Transport met Sprongen: Om zowel scalaire sprongen ($c_f = \alpha c_s$) als fluxsprongen aan te kunnen pakken, introduceren de auteurs een "equivalente scalaire" $h = \alpha c$. Deze transformatie verwijdert de scalaire discontinuïteit, waardoor het gebruik van een verenigde sturende vergelijking voor beide fasen mogelijk wordt. De fluxsprongconditie wordt geherformuleerd tot een equivalente fluxsprong $q^*_w$, die zowel de voorgeschreven fluxsprong als het effect van de scalaire sprong incorporeert. De verenigde sturende vergelijking voor de equivalente scalaire $h$ is:
   $$\frac{\partial h}{\partial t} + \mathbf{u} \cdot \nabla h = \nabla \cdot (D \nabla h) + \nabla \cdot (q^*_w \phi \mathbf{n})$$
   waarbij $D$ de fase-afhankelijke diffusiviteit is.

Het schema is geïmplementeerd in de open-source finite volume code OpenFOAM (versie 8) met behulp van tweede-orde lineaire discretisatie voor ruimtelijke termen en het Euler impliciete schema voor transiënte termen.

Belangrijkste Bijdragen

• Nieuwe Brontermformulering: Het artikel introduceert een bronterm in divergentievorm voor VPM die interfaciale fluxsprongen precies afdwingt zonder niet-fysische oplossingen in het vaste gebied te genereren.
• Verenigde Sturende Vergelijking: Door een equivalente scalaire en flux te definiëren, verenigt de methode de sturende vergelijkingen voor de vloeistof- en vaste fasen, waardoor de noodzaak voor expliciete interfacekoppeling of aparte domeinoplossingen wordt geëlimineerd.
• Uitbreiding naar Algemene Sprongcondities: De methode slaagt erin VPM succesvol uit te breiden om simultane sprongen in zowel de scalaire waarde als de flux te behandelen, een capaciteit die voorheen ontbrak in finite-volume-gebaseerde immersed boundary-schema's.
• Foutanalyse: De auteurs leveren een analytische foutanalyse die aantoont dat het schema een eerste-orde nauwkeurigheid bereikt ten opzichte van de gridgrootte, consistent met conventionele VPM's, zelfs in aanwezigheid van interfaciale kromming.

Resultaten
De voorgestelde methode werd gevalideerd via drie verschillende testgevallen:

1. Eendimensionale Diffusie: Een transiënt diffusieprobleem met een Neumann-randvoorwaarde werd opgelost. De resultaten vertoonden een goede overeenstemming met analytische oplossingen, met een gemiddelde relatieve fout van ongeveer 1,25%. Cruciaal is dat, in tegenstelling tot bestaande VPM-schema's, de huidige methode een nul scalaire veld binnen het vaste deel behield, waardoor niet-fysische diffusie werd vermeden. De fout bleef stabiel bij variërende diffusiviteitswaarden, terwijl bestaande schema's toenemende fouten vertoonden bij hogere diffusiviteit.
2. Vloeistof-Vast Gekoppelde Diffusie: Simulaties van een cilinder in een vierkant domein met zowel scalaire als fluxsprongen werden uitgevoerd. De VPM-resultaten werden vergeleken met body-fitted mesh (BFM) simulaties met behulp van chtMultiRegionFoam. De relatieve afwijkingen tussen de VPM- en BFM-resultaten lagen binnen de 3,0% (gemiddeld 1,46%) over vier testgevallen met variërende diffusiviteitsratio's en sprongmagnitudes.
3. Vloeistof-Vast Gekoppelde Advectie-Diffusie: Een stationaire laminaire stroming ($Re=40$) langs een cilinder met interfaciale sprongen werd gesimuleerd. De VPM-resultaten vertoonden opnieuw een sterke overeenstemming met BFM-benchmarks, met een relatieve afwijking van 1,59%. De methode reproduceerde de scalaire distributie en fluxsprongen nabij de interface accuraat.

Betekenis en Claims
Het artikel claimt dat de voorgestelde VPM een robuust en efficiënt alternatief biedt voor body-fitted meshes voor conjunctieve scalaire transportproblemen met complexe interfaciale geometrieën en sprongcondities. De primaire betekenis ligt in het vermogen om gekoppelde vloeistof-vast problemen op te lossen met een enkele verenigde sturende vergelijking terwijl lokale sprongcondities accuraat worden afgedwongen. De auteurs merken op dat deze aanpak de niet-fysische artefacten vermijdt die gebruikelijk zijn bij traditionele VPM-implementaties voor Neumann-condities.

De auteurs positioneren dit werk bescheiden als een fundamentele stap, waarbij zij suggereren dat het vermogen van de methode om complexe geometrieën te hanteren, het geschikt maakt voor toekomstige toepassingen in ruwe wanden en complexe warmteoverdrachtsoppervlakken. Ze benadrukken ook het potentieel om deze VPM te combineren met adjoint-methoden voor vorm- en topologieoptimalisatie in radiatieve, turbulente en conjunctieve warmteoverdrachtsproblemen, hoewel deze specifieke toepassingen worden geïdentificeerd als richtingen voor toekomstig werk in plaats van directe resultaten van de huidige studie.