Volume of fluid-metoden

Volume of fluid-metoden (VOF) er innen numerisk væskedynamikk (engelsk computational fluid dynamics forkortet CFD) en teknikk for modellering av frie overflater på væsker. Det vil si en numerisk teknikk for sporing og lokalisering den frie overflaten, grensesnittet mellom to væsker eller grensesnittet mellom en væske og et bevegelig fast legeme.

En illustrasjon på resultater fra væskesimulering ved bruk av VOF-metoden. Det viser hva som skjer med væskeoverflaten, dersom en slipper ned en kule.

Metoden er en Eulersk metoder, som er kjennetegnet ved at elementnettet enten er i ro eller beveger seg for å tilfredsstille varierende grensebetingelser over tid. Metoden sikrer i teorien at massen bevares i systemet under strømningene (adveksjon).

Navier-Stokes-ligningene som beskriver bevegelsene i strømningen, må løses separat. Det samme gjelder for alle andre adveksjonsalgoritmer.

Oversikt rediger

Volume of fluid-metoden består av:[1]

  • Et opplegg for å lokalisere overflaten.
  • En beregningsmetode som finner overflaten, der skarpe grensesnitt beveger seg gjennom elementnettet.
  • En metode for å sette randbetingelser på overflatene.

Metoden er basert på såkalte fraksjonsfunksjoner. Det er en skalarfunksjon, definert som integral av en væskes karakteristiske funksjon over volumet til ett kontrollvolum. Volumandelen væske spores gjennom hvert element i elementnettet. Når et element er tomt (uten væske), er verdien av funksjonen null. Når elementet er fullt er funksjonen lik en. Når det er delvis full har en et væske-grensesnitt i elementet, og funksjonen er mellom null og en. Retningen av væske-grensesnittet finnes der verdien av funksjonen endrer seg raskest.

Med denne metoden blir den frie overflaten ikke skarpt definert i enkeltelementene. For å oppnå nøyaktige resultater er en nødt til å dele elementene der funksjonen er mellom null og en, i mindre lokale elementnett.

VOF-metoden er i teorien beregningsvennlig, introduserer bare en ekstra likning, og krever dermed lite lageringsplass. Metoden har sin styrke i å håndtere svært ulineære problemer der den frie overflaten møter skarpe endringer. Den største vanskeligheten er glattingen av den frie overflaten.

Ved to væsker behandler noen overflaten mellom væskene som en grenseflate, mens andre beregner væskestrømningen i begge væskene.[2] Noen svakheter med volume of fluid-formuleringene:[3]

  • Følsomheten av gassområdet for trykkforandringer er vanligvis mye større enn i væskeområder. Dette gjør det vanskelig å oppnå konvergens for trykk og hastigheter. Noen ganger gir teknikken meget store regnetider.
  • På grunn av forskjellig respons på trykk, gass- og væskehastigheter på en grenseflate, blir tangenthastigheten i væsken og gassen i grensesnitt, vanligvis forskjellige. Grensesnittet beveger seg med en gjennomsnittlig hastighet, men dette fører ofte til urealistisk bevegelse av grenseflatene.

For nøyaktig sporing av skarpe væske-gass-grenseflater er det nødvendig å faktisk behandle grensesnittet som en diskontinuitet.

Diskretisering rediger

 
En illustrasjon av en grenseflate med VOF-metoden. Det mørke området er væsken. Tallet i hvert element er skalarfunksjonen (mellom null og en), som forteller hvor mye væske det er i hvert element. Væskeoverflaten tilnærmes i VOF-metoden med de røde strekene. Vektoren normalt på flaten er merket med "n".

For å unngå glatting av den frie overflaten, må ligningen løses med minst mulig diffusjon (masse som forsvinner). Da er en svært avhengig av metoden som brukes for å bevare massen (adveksjon). Mens en første ordens motstrømsløsning glatter grensesnittet, vil en tilsvarende medstrømsløsning føre til uberegnelig oppførsel, når strømmen ikke går langs rutenettet. Videre er høyere ordens løsninger ustabile og skaper svingninger.

Gjennom årene er det utviklet en rekke metoder for å bevare massen. Metodene kan grovt deles inn i:

  • Giver-mottaker-formuleringer (engelsk donor-acceptor formulation) som er basert på avgrensingskriterier og tilgjengelighetskriterier. Avgrensingskriteriet slår fast at verdien av funksjonene er avgrenset mellom null og en. Tilgjengelighetskriteriet sikrer at mengden væske over en grenseflate i et tidstrinn er mindre enn eller lik den mengden som er tilgjengelig i elementet væsken kom fra.
  • Høyere ordens differensiallikninger formulerer den konvektive transportligningen til høyere ordener.
  • Linjeteknikker som omgår problemene, ved ikke å spore grensesnittet i hver element. I stedet fordeles væsken i elementer med grensesnitt, ved hjelp av fordelingen av væskemengden i naboelementer.

Bølgeslag på oljeplattformer rediger

Ved bruk av volume of fluid-metoden kan en analysere bølgekammer som slår inn i dekk på oljeplattformer, fyr og vindmøller. En løser Navier-Stokes-ligningene med et elementnett som forflytter seg med bølgen, og der skalarfunksjonen forteller om elementet er fylt med væske eller ikke. Metoden tillater også at væsken kolliderer med legemer. Metoden har problemer med

  • Å lage økonomiske modeller som er tilstrekkelig nøyaktige.[4]
  • Å lage korrekte randbetingelser (engelsk boundary conditions).[5]
  • Å håndtere viskøse grenselag (engelsk boundary layers").[6]
  • Metoden tillater at væskepartikler skiller lag med resten av væsken, men det kan føre til uønsket tap av masse.[7]
  • Konservering av masse.[8]
  • Håndterer som regel bare todimensjonal strømning.[9]
  • Antakelser om inkompressibel væske medfører at en lokalt kan få uendelig store trykk. En må da midle dette trykket på en måte som gir et fysisk rimelig resultat.

Referanser rediger

  1. ^ Flow Science: VOF – What’s in a Name? - http://www.flow3d.com/home/resources/cfd-101/modeling-techniques/vof-whats-in-a-name Arkivert 6. oktober 2015 hos Wayback Machine.. Besøkt 6.10.2015
  2. ^ Flow Science: VOF – What’s in a Name? - http://www.flow3d.com/home/resources/cfd-101/modeling-techniques/vof-whats-in-a-name Arkivert 6. oktober 2015 hos Wayback Machine.. Besøkt 6.10.2015
  3. ^ Flow Science: Free Surface Modeling Methods - http://www.flow3d.com/home/resources/cfd-101/general-cfd/free-surface-modeling-methods Arkivert 6. oktober 2015 hos Wayback Machine.. Besøkt 6.10.2015
  4. ^ Vestbøstad, Tone M. "A numerical study of wave-in-deck impact using a two-dimensional constrained interpolation profile method." (2009).
  5. ^ A.E.P. Veldman og R.H.H. Huijsmans: Extreme wave impact on offshore platforms and coastal structures, Marin, 2008.
  6. ^ A.E.P. Veldman og R.H.H. Huijsmans: Extreme wave impact on offshore platforms and coastal structures, Marin, 2008.
  7. ^ A.E.P. Veldman og R.H.H. Huijsmans: Extreme wave impact on offshore platforms and coastal structures, Marin, 2008.
  8. ^ Vestbøstad, Tone M. "A numerical study of wave-in-deck impact using a two-dimensional constrained interpolation profile method." (2009).
  9. ^ A.E.P. Veldman og R.H.H. Huijsmans: Extreme wave impact on offshore platforms and coastal structures, Marin, 2008.