Multifluid Hydrodynamic Simulation of Metallic-Plate Collision Using the VOF Method

Deze studie presenteert een één-dimensionale numerieke simulatie van de botsing tussen lood- en stalen platen bij explosief lassen, waarbij met de Volume-of-Fluid-methode en een multifluid Godunov-algoritme de golfbewegingen en de interface tussen de onmengbare fasen worden gemodelleerd en geverifieerd tegen experimentele data.

De kern

De Titel: Een Digitale Slag tussen Metaalplaten

Stel je voor dat je twee verschillende metalen platen hebt: één van lood (zacht en zwaar) en één van staal (hard en stevig). In dit onderzoek laten de wetenschappers deze platen met enorme snelheid op elkaar botsen, alsof ze een explosieve lasmethode testen. Maar in plaats van echt metaal te gebruiken, bouwen ze een digitale simulatie op een computer.

Het doel? Om te begrijpen wat er precies gebeurt in die fracties van een seconde waarin de metalen elkaar raken, vervormen en samensmelten.

Het Probleem: Een Chaos van Golfjes

Wanneer de loodplaat met 500 meter per seconde (dat is sneller dan een kogel!) op de staalplaat botst, gebeurt er iets magisch en chaotisch:

1. De metalen worden zo heet en onder zo'n enorme druk dat ze zich tijdelijk gedragen als een vloeistof, zelfs al zijn ze normaal gesproken vast.
2. Er ontstaan schokgolven die door het metaal razen, net als rimpelingen in een vijver als je een steen erin gooit.
3. Uiteindelijk komen deze golven terug en botsen ze tegen elkaar, waardoor er soms zelfs "negatieve druk" ontstaat (een soort zuigkracht die het metaal probeert uit te rekken).

Het lastige is: hoe simuleer je dit op een computer zonder dat het programma vastloopt of onnauwkeurige resultaten geeft? Vooral omdat je twee verschillende materialen (lood en staal) hebt die niet door elkaar heen lopen, maar wel als één vloeistof moeten worden behandeld tijdens de botsing.

De Oplossing: De "Vloeibare Bak" Methode (VOF)

De auteurs gebruiken een slimme techniek genaamd VOF (Volume of Fluid).

De Analogie:
Stel je een bak met water voor, en je giet er een beetje olie in. De olie en het water mengen zich niet; ze blijven gescheiden.

• In de computerwereld is elke "bakje" (een cel in het rooster) een klein stukje ruimte.
• Soms zit er in zo'n bakje alleen water (lood), soms alleen olie (staal), en soms een mengsel van beide aan de rand.
• De VOF-methode houdt precies bij hoeveel procent van elk bakje "lood" is en hoeveel "staal". Het tekent de grens tussen de twee materialen als een scherpe lijn, in plaats van ze wazig te laten worden.

De Uitdaging: De "Moeilijke Buurman"

Een groot probleem bij het simuleren van botsende metalen is dat ze heel snel van volume veranderen.

• Staal is als een harde rubberen bal: hij duwt terug als je erop duwt, maar verandert niet veel van vorm.
• Lood is als een zachte spons: hij plakt en verandert makkelijk van vorm.

In de simulatie moeten deze twee "buurmannen" in hetzelfde bakje precies dezelfde druk voelen. Als de computer dit niet goed regelt, krijg je onrealistische resultaten. De auteurs hebben een nieuwe manier bedacht om dit te berekenen, gebaseerd op een idee van de beroemde wiskundige Godunov. Ze zorgen ervoor dat de computer rekening houdt met het feit dat het ene materiaal makkelijker comprimeert dan het andere, terwijl ze toch één gemeenschappelijke druk houden.

De Resultaten: Een Perfecte Voorspelling

De wetenschappers lieten hun computerprogramma de botsing simuleren en vergeleken het met echte experimenten en andere rekenmethodes.

1. Scherp beeld: Hun methode gaf een veel scherper beeld van de grens tussen lood en staal dan oudere methodes. Het was alsof ze van een wazige foto overstapten naar een 4K-foto.
2. Negatieve druk: Een van de coolste dingen is dat hun methode "negatieve druk" (zuigkracht) kon berekenen zonder te crashen. Normaal gesproken haken computers af als ze negatieve druk zien, omdat ze denken dat er een fout is. Dit model begrijpt echter dat dit een normaal fysiek verschijnsel is bij het uitrekken van metaal.
3. Tijdsynchronisatie: De tijd die het duurde voordat de "ontlastingsgolf" (de golf die de druk verlaagt) de grens tussen de platen bereikte, kwam exact overeen met wat in het echte leven en in andere simulaties werd gezien.

Waarom is dit belangrijk?

Dit onderzoek is niet alleen leuk wiskundig puzzelen. Het helpt ingenieurs om explosief lassen beter te begrijpen. Dit is een techniek waarbij metalen met een explosie aan elkaar worden gelast (bijvoorbeeld in de luchtvaart of scheepsbouw).

Door te weten hoe de golven zich gedragen en hoe de metalen zich gedragen als vloeistof tijdens de botsing, kunnen mensen in de toekomst:

• Beter lasverbindingen maken.
• Minder materiaal verspillen.
• Veiligere constructies bouwen.

Kortom: De auteurs hebben een digitale "slimme camera" gebouwd die kan kijken hoe twee verschillende metalen platen in een flits van tijd samensmelten, en ze hebben bewezen dat hun nieuwe rekenmethode dit veel nauwkeuriger doet dan de oude methodes.

Technische samenvatting

Titel: Multifluid Hydrodynamische Simulatie van de Klap van Metalen Platen met de VOF-methode

Auteurs: Fedor Belolutskiy, Elena Oparina, Svetlana Fortova
Publicatie: Mathematical Models and Computer Simulations (2025)

---

1. Probleemstelling

Het artikel richt zich op de numerieke simulatie van een explosieve lasproces (explosieve welding), specifiek de botsing van een loodplaat en een stalen plaat. Dit is een complex dynamisch proces waarbij metalen bij hoge snelheid botsen en tijdelijk in een "pseudo-vloeibare" condensatietoestand verkeren.

• Situatie: Een loodplaat (dikte 2 mm, dichtheid 11.300 kg/m³) beweegt met 500 m/s richting een stalen plaat (dikte 3 mm, dichtheid 7.900 kg/m³). Beide platen zijn omgeven door lucht.
• Uitdaging: Het modelleren van de interactie tussen onmengbare fasen (lood, staal, lucht) die elk hun eigen toestandsvergelijking (Equation of State - EoS) hebben. Een specifiek probleem is het optreden van trekspanningen (negatieve druk) in het metaal na het passeren van een ontlastingsgolf, wat numerieke instabiliteit kan veroorzaken in traditionele methoden.
• Doel: Het ontwikkelen en valideren van een numeriek model dat de golfvoortplanting, de beweging van het contactoppervlak en de vorming van negatieve drukken nauwkeurig kan simuleren zonder extra procedures voor negatieve drukken.

2. Methodologie

De auteurs gebruiken een multifluid hydrodynamische aanpak gebaseerd op de Euler-vergelijkingen in een Euleriaanse raamwerk.

• Wiskundig Model:

  • Het model behandelt de mengcellen als een "effectieve fase" met gemeenschappelijke druk en snelheid, maar behoudt individuele eigenschappen (dichtheid, inwendige energie) voor elke fase.
  • Er wordt uitgegaan van mechanisch evenwicht (gelijke druk en snelheid) en isentropische compressie.
  • De vergelijkingen omvatten behoudswetten voor massa, impuls en energie, aangevuld met transportvergelijkingen voor de volumefracties ($f^{(\alpha)}$) die rekening houden met de compressibiliteit van de verschillende fasen.
  • Voor de toestandsvergelijkingen wordt een stiffened-gas model gebruikt voor alle fasen (lood, staal, lucht).

• Numerieke Methode:

  • VOF (Volume-of-Fluid): De interface tussen de onmengbare fasen wordt getrackt met de VOF-methode, die bekend staat om zijn vermogen om scherpe interfaces te behouden en massa te conserveren.
  • Oplossingsalgoritme: Een Godunov-type (finite-volume) algoritme wordt gebruikt.
    • De lokale Riemann-problemen aan de celgrenzen worden opgelost met de HLLC-solver (Harten-Lax-van Leer-Contact).
    • Voor hogere resolutie wordt MUSCL-Hancock reconstructie toegepast (tweede orde nauwkeurigheid).
    • Een druk-relaxatie stap wordt geïntegreerd om het mechanisch evenwicht in mengcellen te herstellen na elke tijdstap.
  • Behandeling van Negatieve Druk: Het model is ontworpen om trekspanningen (negatieve druk) natuurlijk te registreren als negatieve drukwaarden, zonder dat er extra "curing"-procedures nodig zijn, wat de robuustheid van het schema verhoogt.

3. Belangrijkste Bijdragen

1. Aangepaste CGF-methode: De auteurs passen een methode voor (oorspronkelijk ontwikkeld door Colella, Glaz en Ferguson voor gassen) toe op samengeperste media (condensed matter) met verschillende toestandsvergelijkingen.
2. Robuustheid bij Trekspanningen: Het numerieke schema kan succesvol simuleren wanneer de druk negatief wordt (trekspanning), een situatie die vaak leidt tot instabiliteit in andere methoden.
3. Lage Numerieke Diffusie: Door gebruik te maken van de HLLC-solver in combinatie met MUSCL-reconstructie, wordt de numerieke diffusie aan de interfaces sterk verminderd, wat leidt tot een hogere resolutie van de golven en interfaces.
4. Validatie: De resultaten worden vergeleken met experimentele data en eerdere simulaties (o.a. met het Baer-Nunziato-model), waarbij een uitstekende overeenkomst wordt aangetoond.

4. Resultaten

De simulatie werd uitgevoerd voor een tijdsduur van 2 $\mu$s. De belangrijkste bevindingen zijn:

• Golfvoortplanting:
  • De botsing genereert schokgolven die zich door de metalen voortplanten.
  • Wanneer deze golven de vrije oppervlakken (grens met lucht) bereiken, ontstaan er ontlastingsgolven (unloading waves).
  • De ontlastingsgolf uit het staal bereikt het contactoppervlak rond $t \approx 1,13$ $\mu$s, wat leidt tot een toename van de contactsnelheid van 191 m/s naar 425,2 m/s.
  • Rond $t \approx 1,4$ $\mu$s snijden de ontlastingsgolven elkaar in het lood, wat leidt tot negatieve drukken en een afname van de dichtheid.
• Vergelijking van Methoden:
  • Simulaties met MUSCL-reconstructie tonen aanzienlijk minder numerieke diffusie dan die zonder reconstructie.
  • De resultaten van het huidige model op een fijnere rooster (3.601 cellen) komen binnen 4,6% overeen met een zeer gedetailleerde simulatie (10.000 cellen).
  • De aankomsttijd van de ontlastingsgolf (1,13 $\mu$s) komt overeen met experimentele data en eerdere numerieke studies (1,05 $\mu$s).
• Negatieve Druk: Het model slaagt erin de vorming van negatieve drukken in het lood (na $t \approx 1,5$ $\mu$s) correct te modelleren, terwijl andere methoden hier vaak tegenaan lopen.

5. Betekenis en Conclusie

Dit onderzoek demonstreert dat een multifluid hydrodynamisch model, gebaseerd op de VOF-methode en de HLLC-solver, zeer geschikt is voor het simuleren van explosief lassen.

• Nauwkeurigheid: Het model biedt een hoge resolutie van contactinterfaces en golven met lage numerieke diffusie.
• Toepasbaarheid: Het vermogen om trekspanningen en negatieve drukken te hanteren maakt het model ideaal voor het bestuderen van complexe dynamische processen in metalen.
• Toekomst: De auteurs plannen om deze methode uit te breiden naar twee dimensies om instabiliteiten (zoals Rayleigh-Taylor) aan het contactoppervlak tijdens explosief lassen verder te onderzoeken.

Samenvattend biedt dit artikel een robuust en nauwkeurig numeriek kader voor het simuleren van botsingen tussen metalen platen, waarbij de beperkingen van eerdere methoden op het gebied van negatieve drukken en interface-resolutie worden overwonnen.