Moment of fluid method for multimaterial flows
|
|
- Norma Stephens
- 5 years ago
- Views:
Transcription
1 Moment of fluid method for multimaterial flows M. Jemison, M. Sussman, M. Shashkov Department of Mathematics, Florida State University October 1, 2012
2 Recent Motivating Work Volume preserving Moment-of-Fluid reconstruction for multiple materials. Moment of Fluid method for deforming boundary problems (Ahn and Shashkov, JCP 2007; Ahn and Shashkov, JCP 2009). variable density cell centered pressure projection (Kwatra, Su, Gretarsson, Fedkiw 2009)
3 Laminar jet - density ratio 819:1 X. Li (UTRC), M. Arienti (UTRC), M. Soteriou (UTRC), M. Sussman (FSU).
4 Turbulent jet spray - density ratio 819:1 X. Li (UTRC), M. Arienti (UTRC), M. Soteriou (UTRC), M. Sussman (FSU).
5 Impinging jets - density ratio 819:1 Arienti (UTRC), Soteriou (UTRC), Sussman (FSU).
6 Flow through diesel nozzle with moving valve
7 Flow through diesel nozzle with moving valve
8 Simulation of Flow through diesel nozzle with moving valve
9 AMR grid
10 Six hole diesel nozzle Pressure (bar)! Figure: Snapshot of a simulation of liquid injection in a Bosch six-hole vertical diesel nozzle. Effective fine grid resolution 640x640x256.
11 Microfluidics Sussman (FSU), Jemison (FSU), Duffy (FSU), Roper (FSU)
12 Vortex Rings Sussman (FSU), Ohta (Tokushima)
13 Vortex Rings Sussman (FSU), Ohta (Tokushima)
14 Background Shock capturing (ENO) - Shu, Osher, JCP 1989 Front capturing (LS) - Sussman, Smereka, Osher, JCP 1994 Conservative LS - Olsson, Kreiss, Zahedi, JCP 2007 Front tracking - Unverdi, Tryggvason, JCP 1992 VOF - Kothe, Brackbill, Zemach, JCP 1992 CLSVOF - Sussman, Puckett; Sussman, Smith, Hussaini, et al. CLSVOF - Stern Particle LS - Enright, Fedkiw Refined LS - Herrmann, Pitsch
15 Grid
16 We propose Semi-Lagrangian technique for multimaterial problems, in contrast to a finite volume technique which is (in 1D, u is constant): ρ i = 1 xi+1/2 ρ(x)dx x x i 1/2 d ρ i (t) dt = f (ρ i+1/2, ρ+ i+1/2 ) f (ρ i 1/2, ρ+ i 1/2 ) x where d/dt discretized using high order TVD preserving RK and ρ i+1/2 derived from a high order ENO or WENO reconstruction.
17 Results motivating SL instead of Finite volume
18 Results motivating MOF/VOF instead of ENO/WENO reconstruction
19 Results motivating MOF/VOF instead of ENO/WENO reconstruction
20 Governing Equations ρ(φ) Du = p + (2µ(φ)D) σκ(φ) H(φ) + ρ(φ)gẑ Dt u = 0 Dφ Dt = 0 { 1 φ 0 H(φ) = 0 φ < 0 ρ(φ) = ρ L H(φ) + ρ G (1 H(φ)), µ(φ) = µ L H(φ) + µ G (1 H(φ)), κ(φ) = φ φ, D = u + ( u)t, 2
21 Splitting of advection terms from the pressure terms Kwatra, Su, Gretarsson (JCP, 2009) ρ t + (ρu) = 0 (ρu) t + (ρuu) = p (ρe) t + (ρue) = (up)
22 Splitting of advection terms from the pressure terms Kwatra, Su, Gretarsson (JCP, 2009) solve the following equations to get ρ n+1, u, and E. These equations are solved using a directionally split method in which ρ at material boundaries is derived from the multimaterial MOF reconstruction. ρ t + (ρu) = 0 (ρu) t + (ρuu) = 0 (ρe) t + (ρue) = 0
23 Splitting of advection terms from the pressure terms p n+1 p = ρc 2 u t u n+1 = u t p ρ n+1 p n+1 ρc 2 t 2 pn+1 ρ n+1 = p ρc 2 t 2 1 t u
24 Splitting of advection terms from the pressure terms u n+1 i+1/2 = u i+1/2 u n+1 i pn+1 i+1 t pn+1 i ρ i+1/2 x = u i t pn+1 i+1/2 pn+1 i 1/2 ρ i x (ρe) n+1 = (ρe) t (u n+1 p n+1 )
25 CLSMOF (illustration of Piecewise linear reconstruction)
26 Multimaterial MOF reconstruction For cell i 1, slope of Red (solid) cut determined first, then slope of green (fluid 1) cut determined second. Blue material occupies remaining unfilled region. For cell i, slope of blue (fluid 2) cut determined first. Green material occupies remaining unfilled region.
27 Multimaterial MOF reconstruction Slope of Red (solid) cut determined first, then slope of green (fluid 1) cut determined second. Blue material occupies remaining unfilled region.
28 MOF optimization problem ˆn points into the material being reconstructed. The starting guess is ˆn = x ref x unfilled x ref x unfilled. {x R 3 ˆn (x x i,j,k ) + b = 0} F ref (ˆn, b) F A (ˆn, b) = 0 E MOF = x ref x A (ˆn, b) 2 E MOF (Φ, Θ ) = f (Φ, Θ ) 2 = min (Φ,Θ) f (Φ, Θ) 2, f : R 2 R 3, f (Φ, Θ) = (x ref x A (Φ, Θ)) Constrained optimization problem solved using the Gauss-Newton method.
29 Gauss-Newton minimization algorithm 0. initial angles: (Φ 0, Θ 0 ). Tolerance: tol = 10 8 x. while not converged 1. find b k (Φ k, Θ k ). 2. find the centroid x k (b k, Φ k, Θ k ) 3. find the Jacobian matrix J k of f evaluated in (Φ k, Θ k ) and f k = f (Φ k, Θ k ) 4. stop if one of the following three conditions is fulfilled: J T k f k tol 10 2 x fk < tol k = 11 else continue 5. solve the linear least squares problem: find s k R 2 such that J k s k + f k 2 = min s R 2 J ks + f k 2 by means of the normal equations. (J T k J ks k = J T k f k) 6. update the angles: (Φ k+1, Θ k+1 ) = (Φ k, Θ k ) + s k 7. k := k + 1
30 Property of MOF reconstruction x ref x A Figure: The reference centroid, x ref, does not coincide with the exact centroid, x A, for a parabolic interface cutting the cell. The difference between the two centroids is proportional to the curvature.
31 Backwards Tracing A = x i 1/2 u i 1/2 t B = x i 1/2 C = x i+1/2 u i+1/2 t
32 Backwards Tracing (before and after) A = x i 1/2 u i 1/2 t B = x i 1/2 C = x i+1/2 u i+1/2 t D = x i 1/2 E = J i (x i 1/2 ) F = x i+1/2
33 Backwards Tracing - conserved variables AKA Eulerian Implicit ρ i,j,k = 1 x y z xi+1/2 yj+1/2 zk+1/2 x i 1/2 WLOG, drop j and k subscripts. ρ n+1 i = y j 1/2 z k 1/2 i+1 i =i 1 Ω D i Ω i ρn i (x)dx x ρ(x, y, z)dzdydx
34 Backwards Tracing - volume fraction AKA Eulerian Implicit (Scardovelli and Zaleski, Le Chenedac and Pitsch) 1 xi+1/2 yj+1/2 zk+1/2 F i,j,k = H(φ(x, y, z))dzdydx x y z x i 1/2 y j 1/2 z k 1/2 H(φ) = 1 if φ 0 (i.e. (x, y, z) inside the material being advected). φ n+1 i = φ n i (J 1 i (x), y, z) F n+1 i = i+1 i =i 1 J i (Ω Di Ω H(φn+1 i ) i x y z (x, y, z))dxdydz
35 Backwards Tracing - centroid AKA Eulerian Implicit x act i,j,k = 1 F i,j,k x y z xi+1/2 yj+1/2 zk+1/2 x i 1/2 y j 1/2 z k 1/2 H(φ) = 1 if φ 0 (i.e. (x, y, z) inside the material being advected). φ n+1 i = φ n i (J 1 i (x), y, z) xh(φ(x, y, z))dzdydx x n+1 i = i+1 i =i 1 J i (Ω Di Ω xh(φn+1 i ) i x y z F n+1 i (x, y, z))dxdydz
36 Backwards Tracing - mapping function Ω i = [x i 1/2, x i+1/2 ] Ω D i = [x i 1/2 u i 1/2 t, x i+1/2 u i+1/2 t] Ω T i = [x i 1/2, x i+1/2 ]
37 Backwards Tracing - mapping function Let J i (x) = Ax + B be a function which maps from cell i departure region to cell i target region. For backwards tracing, J i (x i 1/2 u i 1/2 t) = A(x i 1/2 u i 1/2 t) + B = x i 1/2 J i (x i+1/2 u i+1/2 t) = A(x i+1/2 u i+1/2 t) + B = x i+1/2 A = 1 1 (u i+1/2 u i 1/2 ) t x B = x i 1/2 A(x i 1/2 u i 1/2 t) As long as u t < x/2, the mapping has an inverse.
38 Backwards Tracing - inverse mapping function Let J 1 i (x) be the inverse map which is a function which maps from cell i target region to cell i departure region. If J i (x) = Ax + B, then J 1 i (x) = 1 A x B A Let Ω 1 = J 1 i (Ω) be the inverse map which finds the region Ω 1 which maps to Ω under the definition of J i. In otherwords, Ω 1 = J 1 i (Ω) = [J 1 i Ω = [a, b] (a), J 1 i (b)]
39 Backwards Tracing - mapping function to find target region Likewise, J i (Ω) is defined as the target region that the region Ω maps to: Ω = [a, b] J i (Ω) = [J i (a), J i (b)]
40 Backwards Tracing - interface advection defined by mapping Given a local level set function, After advection, we have φ n i (x, y, z) = n (x x i ) + α, φ n+1 i (x, y, z) = φ n i (J 1 i (x), y, z) = n 1 ( x A B A x i ) + n 2(y y j ) + n 3 (z z k ) = n 1 A (x J i(x i )) + n 2 (y y j ) + n 3 (z z k ) + α
41 Forwards Tracing A = J 1 i 1 (x i 1/2) B = x i 1/2 C = J 1 i (x i+1/2 ) J i 1 and J i define the mappings of cells i 1 and i contents forward, respectively.
42 Forwards Tracing (before and after) A = J 1 i 1 (x i 1/2) B = x i 1/2 C = J 1 i (x i+1/2 ) D = x i 1/2 E = J i (x i 1/2 ) F = x i+1/2
43 Forwards Tracing - conserved variables AKA Lagrangian Explicit ρ n+1 i = ρ i = 1 x i+1 i =i 1 xi+1/2 x i 1/2 ρ(x)dx J 1 i (Ω T i Ω i ) ρn i (x)dx x
44 Forwards Tracing - volume fraction AKA Lagrangian Explicit (Scardovelli and Zaleski, Le Chenedac and Pitsch) F i,j,k = 1 x y z xi+1/2 yj+1/2 zk+1/2 x i 1/2 y j 1/2 z k 1/2 H(φ(x, y, z))dzdydx H(φ) = 1 if φ 0 (i.e. (x, y, z) inside the material being advected). φ n+1 i (x, y, z) = φ n i (J 1 i (x), y, z) F n+1 i = i+1 i =i 1 Ω i Ω Ti H(φn+1 i x y z (x, y, z))dxdydz
45 Forwards Tracing - centroid AKA Lagrangian Explicit x act i,j,k = 1 F i,j,k x y z xi+1/2 yj+1/2 zk+1/2 x i 1/2 y j 1/2 z k 1/2 H(φ) = 1 if φ 0 (i.e. (x, y, z) inside the material being advected). φ n+1 i (x, y, z) = φ n i (J 1 i (x), y, z) xh(φ(x, y, z))dzdydx x n+1 i = i+1 i =i 1 Ω i Ω Ti xh(φn+1 i F n+1 i x y z (x, y, z))dxdydz
46 Forwards Tracing - mapping function Ω i = [x i 1/2, x i+1/2 ] Ω T i = [x i 1/2 + u i 1/2 t, x i+1/2 + u i+1/2 t] Ω D i = [x i 1/2, x i+1/2 ]
47 Forwards Tracing - mapping function Let J i (x) = Ax + B be a function which maps from cell i departure region to cell i target region. For forwards tracing, J i (x i 1/2 ) = Ax i 1/2 + B = x i 1/2 + u i 1/2 t J i (x i+1/2 ) = Ax i+1/2 + B = x i+1/2 + u i+1/2 t A = 1 + (u i+1/2 u i 1/2 ) t x B = x i 1/2 + u i 1/2 t Ax i 1/2 As long as u t < x/2, A is positive.
48 Forwards Tracing - inverse mapping function Let J 1 i (x) be the inverse map which is a function which maps from cell i target region to cell i departure region. If J i (x) = Ax + B, then J 1 i (x) = 1 A x B A Let Ω 1 = J 1 i (Ω) be the inverse map which finds the region Ω 1 which maps to Ω under the definition of J i. In otherwords, Ω 1 = J 1 i (Ω) = [J 1 i Ω = [a, b] (a), J 1 i (b)]
49 Forwards Tracing - mapping function to find target region Likewise, J i (Ω) is defined as the target region that the region Ω maps to: Ω = [a, b] J i (Ω) = [J i (a), J i (b)]
50 Forwards Tracing - advection of interface reconstruction by defined mapping Given a local level set function, φ n i (x, y, z) = n (x x i ) + α, After forward advection, we have φ n+1 i (x, y, z) = φ n i (J 1 i (x), y, z) = n 1 ( x A B A x i ) + n 2(y y j ) + n 3 (z z k ) = n 1 A (x J i (x i )) + n 2(y y j ) + n 3 (z z k ) + α
51 DS Semi-Lagrangian discretization - Strang Splitting 1. x direction; backwards projection. (Eulerian Implicit) 2. y direction; forwards projection. (Lagrangian Explicit) 3. z direction; backwards projection. (Eulerian Implicit) 4. z direction; forwards projection. (Lagrangian Explicit) 5. y direction; backwards projection. (Eulerian Implicit) 6. x direction; forwards projection. (Lagrangian Explicit)
52 2D volume conserved exactly
53 reversible single vortex
54 2D Unsplit MOF, reversible vortex T=1/2 RK2, backwards tracing: Figure: single vortex problem (T = 1/2) midway through the simulation and at the end. Shape deforms back to a circle. Two materials. At t = T = 1/2, symmetric difference error is 7.4E 5 (operator split symmetric difference error: 14.3E 5)
55 2D Unsplit multimaterial MOF, reversible vortex T=1/2 RK2, backwards tracing: Figure: single vortex problem (T = 1/2) midway through the simulation and at the end. Shape deforms back to a circle. Three materials. Initially, Material 1 (red) is the left half of the circle, Material 2 (green) is the right half, and Material 3 (blue) is the space outside the circle. At t = T = 1/2, symmetric difference error is 9.6E 5 for material 1, 10.1E 5 for material 2, and 14.8E 5 for material 3.
56 3D Unsplit MOF, reversible vortex T=3 RK2, backwards tracing: Figure: 3D single vortex problem (T = 3) midway through the simulation and at the end. Shape deforms back to a sphere. Two materials. At t = T = 3, symmetric difference error is (operator split symmetric difference error: )
57 3D Unsplit multimaterial MOF, reversible vortex T=3 RK2, backwards tracing: Figure: single vortex problem (T = 3) midway through the simulation and at the end. Shape deforms back to a sphere. Three materials. Initially, Material 1 (red) is the left half of the sphere, and Material 2 (green) is the right half. At t = T = 3, symmetric difference error is for material 1, for material 2, and for material 3.
58 Rotating Notched disk (Zalesak s problem)
59 Rotating Letter A
60 Rising gas bubble in liquid Figure: Left: condition 1 Right: condition 3 Figure: Left: condition 4 Right: condition 6
61 Surface tension driven vibrations of a drop
62 Oscillating Cylinder
63 Cylinder falling into a pool of liquid
64 Impinging jets - same material
65 Impinging jets - different materials
66 Future work The multimaterial MOF representation enables more accurate simulation of: multimaterial flows with minimal volume fluctuation - capturing corners and filaments. surface tension and contact line effects. simulating transport on deforming surfaces. predicting mass transfer on deforming surfaces. predicting boundary layer effects on underresolved grids. Improvements are still being made to the multimaterial MOF scheme in accelerating the multimaterial MOF reconstruction and unsplit multimaterial MOF advection algorithm. Adding the capability to robustly simulate compressible multiphase/multimaterial flows is one priority now.
Overview of Traditional Surface Tracking Methods
Liquid Simulation With Mesh-Based Surface Tracking Overview of Traditional Surface Tracking Methods Matthias Müller Introduction Research lead of NVIDIA PhysX team PhysX GPU acc. Game physics engine www.nvidia.com\physx
More informationConsiderations about level-set methods: Accuracy, interpolation and reinitialization
Considerations about level-set methods: Accuracy, interpolation and reinitialization Gustavo C. Buscaglia Coll: E. Dari, R. Ausas, L. Ruspini Instituto de Ciências Matemáticas e de Computação Universidade
More informationFLORIDA STATE UNIVERSITY COLLEGE OF ARTS AND SCIENCES AN ASYMPTOTICALLY PRESERVING METHOD FOR MULTIPHASE FLOW MATTHEW JEMISON
FLORIDA STATE UNIVERSITY COLLEGE OF ARTS AND SCIENCES AN ASYMPTOTICALLY PRESERVING METHOD FOR MULTIPHASE FLOW By MATTHEW JEMISON A Dissertation submitted to the Department of Mathematics in partial fulfillment
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More informationDrop Impact Simulation with a Velocity-Dependent Contact Angle
ILASS Americas 19th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Drop Impact Simulation with a Velocity-Dependent Contact Angle S. Afkhami and M. Bussmann Department
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationThe Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection
The Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection Jonathan Robey 1 September 12, 2016 1 Funded in part by a GAANN Fellowship 1 Motivation Mantle
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationGradient Free Design of Microfluidic Structures on a GPU Cluster
Gradient Free Design of Microfluidic Structures on a GPU Cluster Austen Duffy - Florida State University SIAM Conference on Computational Science and Engineering March 2, 2011 Acknowledgements This work
More informationA CONSERVATIVE LEVEL SET METHOD FOR SHARP INTERFACE MULTIPHASE FLOW SIMULATION
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal,4-7 June 2 A CONSERVATIVE LEVEL SET METHOD FOR SHARP INTERFACE MULTIPHASE FLOW
More informationConservative Level Set/Ghost Fluid Method for Simulating Primary Atomization
ILASS Americas 20th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 2007 Conservative Level Set/Ghost Fluid Method for Simulating Primary Atomization O. Desjardins, V. Moureau,
More informationDroplet collisions using a Level Set method: comparisons between simulation and experiments
Computational Methods in Multiphase Flow III 63 Droplet collisions using a Level Set method: comparisons between simulation and experiments S. Tanguy, T. Ménard & A. Berlemont CNRS-UMR6614-CORIA, Rouen
More informationLevel set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field
Level Set Methods Overview Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Convection Upwind ddifferencingi
More informationA 3D VOF model in cylindrical coordinates
A 3D VOF model in cylindrical coordinates Marmar Mehrabadi and Markus Bussmann Department of Mechanical and Industrial Engineering, University of Toronto Recently, volume of fluid (VOF) methods have improved
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 58 63 (2007) Published online 22 May 2007 in Wiley InterScience (www.interscience.wiley.com).151 An explicit and conservative remapping strategy for semi-lagrangian
More informationPrim ary break-up: DNS of liquid jet to improve atomization m odelling
Computational Methods in Multiphase Flow III 343 Prim ary break-up: DNS of liquid jet to improve atomization m odelling T. Ménard, P. A. Beau, S. Tanguy, F. X. Demoulin & A. Berlemont CNRS-UMR6614-CORIA,
More informationAn added mass partitioned algorithm for rigid bodies and incompressible flows
An added mass partitioned algorithm for rigid bodies and incompressible flows Jeff Banks Rensselaer Polytechnic Institute Overset Grid Symposium Mukilteo, WA October 19, 216 Collaborators Bill Henshaw,
More informationCS-184: Computer Graphics Lecture #21: Fluid Simulation II
CS-184: Computer Graphics Lecture #21: Fluid Simulation II Rahul Narain University of California, Berkeley Nov. 18 19, 2013 Grid-based fluid simulation Recap: Eulerian viewpoint Grid is fixed, fluid moves
More informationA mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates
A mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates Peter Hjort Lauritzen Atmospheric Modeling & Predictability Section National Center for Atmospheric
More informationCurvature computation in Volume of Fluid. 1. Institut Jean Le Rond d Alembert, CNRS & Université Pierre et Marie Curie -- UPMC Paris 6
Curvature computation in Volume of Fluid Stéphane Zaleski, Ruben Scardovelli 2 & Yue Ling. Institut Jean Le Rond d Alembert, CNRS & Université Pierre et Marie Curie -- UPMC Paris 6 2. DIENCA, Laboratorio
More informationImagery for 3D geometry design: application to fluid flows.
Imagery for 3D geometry design: application to fluid flows. C. Galusinski, C. Nguyen IMATH, Université du Sud Toulon Var, Supported by ANR Carpeinter May 14, 2010 Toolbox Ginzburg-Landau. Skeleton 3D extension
More informationInvestigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow
Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training
More informationMETHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS
XVIII International Conference on Water Resources CMWR 2010 J. Carrera (Ed) CIMNE, Barcelona 2010 METHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS Poorya A. Ferdowsi
More informationLeast Square Curvature Calculation Method for VOF Schemes
ILASS Americas 28th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2016 Least Square Curvature Calculation Method for VOF Schemes E. Cauble, M. Owkes Department of Mechanical
More informationA Hybrid Particle Level Set Method for Improved Interface Capturing
Journal of Computational Physics 183, 83 116 (2002) doi:10.1006/jcph.2002.7166 A Hybrid Particle Level Set Method for Improved Interface Capturing Douglas Enright,,1,2 Ronald Fedkiw,,1,2 Joel Ferziger,,2
More informationParallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes
Parallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes L. Jofre, O. Lehmkuhl, R. Borrell, J. Castro and A. Oliva Corresponding author: cttc@cttc.upc.edu Centre Tecnològic
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationImplicit surface tension model for stimulation of interfacial flows
surface tension model for stimulation of interfacial flows Vinh The University of Massachusetts Dartmouth March 3rd, 2011 About my project Project Advisor Dr. Mehdi Raessi Department of Mechanical Engineering
More informationPreliminary Spray Cooling Simulations Using a Full-Cone Water Spray
39th Dayton-Cincinnati Aerospace Sciences Symposium Preliminary Spray Cooling Simulations Using a Full-Cone Water Spray Murat Dinc Prof. Donald D. Gray (advisor), Prof. John M. Kuhlman, Nicholas L. Hillen,
More informationAdvective and conservative semi-lagrangian schemes on uniform and non-uniform grids
Advective and conservative semi-lagrangian schemes on uniform and non-uniform grids M. Mehrenberger Université de Strasbourg and Max-Planck Institut für Plasmaphysik 5 September 2013 M. Mehrenberger (UDS
More information(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more
Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon
More information2. Give an example of a non-constant function f(x, y) such that the average value of f over is 0.
Midterm 3 Review Short Answer 2. Give an example of a non-constant function f(x, y) such that the average value of f over is 0. 3. Compute the Riemann sum for the double integral where for the given grid
More informationMath 690N - Final Report
Math 690N - Final Report Yuanhong Li May 05, 008 Accurate tracking of a discontinuous, thin and evolving turbulent flame front has been a challenging subject in modelling a premixed turbulent combustion.
More informationUSE OF THE PARTICLE LEVEL SET METHOD FOR ENHANCED RESOLUTION OF FREE SURFACE FLOWS
USE OF THE PARTICLE LEVEL SET METHOD FOR ENHANCED RESOLUTION OF FREE SURFACE FLOWS a dissertation submitted to the program in scientific computing and computational mathematics and the committee on graduate
More informationFEM techniques for interfacial flows
FEM techniques for interfacial flows How to avoid the explicit reconstruction of interfaces Stefan Turek, Shu-Ren Hysing (ture@featflow.de) Institute for Applied Mathematics University of Dortmund Int.
More informationTwo-Phase flows on massively parallel multi-gpu clusters
Two-Phase flows on massively parallel multi-gpu clusters Peter Zaspel Michael Griebel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität Bonn Workshop Programming of Heterogeneous
More informationConservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow. Abstract
Conservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow Jing-Mei Qiu 1 and Chi-Wang Shu Abstract In this paper, we propose a semi-lagrangian finite difference
More informationVolume Tracking on Adaptively Refined Grids with Curvature Based Refinement
Volume Tracking on Adaptively Refined Grids with Curvature Based Refinement Mayank Malik, Markus Bussmann Department of Mechanical & Industrial Engineering, University of Toronto mayank.malik@utoronto.ca,
More informationA MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES
A MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES Khalil Bensassi, and Herman Deconinck Von Karman Institute for Fluid Dynamics Aeronautics & Aerospace
More informationDr. Allen Back. Nov. 21, 2014
Dr. Allen Back of Nov. 21, 2014 The most important thing you should know (e.g. for exams and homework) is how to setup (and perhaps compute if not too hard) surface integrals, triple integrals, etc. But
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More information>> 1 because the resulting discretization matrix is poorly conditioned. Consider the 1D discretization of eqn. 1.1
AN IMPROVED VARIABLE DENSITY PRESSURE PROJECTION SOLVER FOR ADAPTIVE MESHES AUSTEN DUFFY, ALAN KUHNLE, AND MARK SUSSMAN Abstract Tatebe [7] introduced the multigrid preconditioned conjugate gradient method
More informationWater. Notes. Free surface. Boundary conditions. This week: extend our 3D flow solver to full 3D water We need to add two things:
Notes Added a 2D cross-section viewer for assignment 6 Not great, but an alternative if the full 3d viewer isn t working for you Warning about the formulas in Fedkiw, Stam, and Jensen - maybe not right
More informationExample Simulations in OpenFOAM
Example Simulations in OpenFOAM Hrvoje Jasak h.jasak@wikki.co.uk Wikki Ltd, United Kingdom FSB, University of Zagreb, Croatia 18/Nov/2005 Example Simulations in OpenFOAM p.1/26 Outline Objective Present
More informationMetafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction
ALE simulations ua sus using Metafor eao 1. Introduction 2. Operator split 3. Convection schemes 4. Rezoning methods 5. Contact with friction 1 Introduction EULERIAN FORMALISM Undistorted mesh Ideal for
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationALE and AMR Mesh Refinement Techniques for Multi-material Hydrodynamics Problems
ALE and AMR Mesh Refinement Techniques for Multi-material Hydrodynamics Problems A. J. Barlow, AWE. ICFD Workshop on Mesh Refinement Techniques 7th December 2005 Acknowledgements Thanks to Chris Powell,
More informationComputational Astrophysics 5 Higher-order and AMR schemes
Computational Astrophysics 5 Higher-order and AMR schemes Oscar Agertz Outline - The Godunov Method - Second-order scheme with MUSCL - Slope limiters and TVD schemes - Characteristics tracing and 2D slopes.
More informationCFD-1. Introduction: What is CFD? T. J. Craft. Msc CFD-1. CFD: Computational Fluid Dynamics
School of Mechanical Aerospace and Civil Engineering CFD-1 T. J. Craft George Begg Building, C41 Msc CFD-1 Reading: J. Ferziger, M. Peric, Computational Methods for Fluid Dynamics H.K. Versteeg, W. Malalasekara,
More informationMultiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM3
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol., pp.491-497 (11 ARTICLE Multiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM3 Akio IKEBATA 1,, Yoshihiko MURAOKA 1 and
More informationMoving Interface Problems: Methods & Applications Tutorial Lecture II
Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,
More informationAdarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs)
OBJECTIVE FLUID SIMULATIONS Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs) The basic objective of the project is the implementation of the paper Stable Fluids (Jos Stam, SIGGRAPH 99). The final
More informationA TSTT integrated FronTier code and its applications in computational fluid physics
Institute of Physics Publishing Journal of Physics: Conference Series 16 (2005) 471 475 doi:10.1088/1742-6596/16/1/064 SciDAC 2005 A TSTT integrated FronTier code and its applications in computational
More informationAn adaptive discretization of incompressible flow using a multitude of moving Cartesian grids
An adaptive discretization of incompressible flow using a multitude of moving Cartesian grids R. Elliot English, Linhai Qiu,YueYu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationImaging of flow in porous media - from optimal transport to prediction
Imaging of flow in porous media - from optimal transport to prediction Eldad Haber Dept of EOS and Mathematics, UBC October 15, 2013 With Rowan Lars Jenn Cocket Ruthotto Fohring Outline Prediction is very
More informationAn adaptive discretization of compressible flow using a multitude of moving Cartesian grids
An adaptive discretization of compressible flow using a multitude of moving Cartesian grids Linhai Qiu, Wenlong Lu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationAn adaptive discretization of incompressible flow using a multitude of moving Cartesian grids
An adaptive discretization of incompressible flow using a multitude of moving Cartesian grids R. Elliot English, Linhai Qiu, Yue Yu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford,
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationChallenges and recent progress in developing numerical methods for multi-material ALE Hydrocodes
Challenges and recent progress in developing numerical methods for multi-material ALE Hydrocodes ICFD 25 year Anniversary Conference 15-16 th September 2008 Oxford University Andrew Barlow AWE Introduction
More informationNumerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1
Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1 Overview These notes are meant to be an overview of the various memory access patterns
More informationSurface Tension Approximation in Semi-Lagrangian Level Set Based Fluid Simulations for Computer Graphics
Surface Tension Approximation in Semi-Lagrangian Level Set Based Fluid Simulations for Computer Graphics Israel Pineda and Oubong Gwun Chonbuk National University israel_pineda_arias@yahoo.com, obgwun@jbnu.ac.kr
More informationVector Image Polygon (VIP) limiters in ALE Hydrodynamics
New Models and Hydrocodes, Paris, France,24-28 May 2010 Vector Image Polygon (VIP) limiters in ALE Hydrodynamics 6 4 8 3 5 2 79 1 2 1 1 1 Gabi Luttwak 1 and Joseph Falcovitz 2 1 Rafael, P.O. Box 2250,
More informationDr. Allen Back. Nov. 19, 2014
Why of Dr. Allen Back Nov. 19, 2014 Graph Picture of T u, T v for a Lat/Long Param. of the Sphere. Why of Graph Basic Picture Why of Graph Why Φ(u, v) = (x(u, v), y(u, v), z(u, v)) Tangents T u = (x u,
More informationWeno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach
Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach Florence Hubert and Rémi Tesson Abstract In this paper we develop a DDFV approach for WENO scheme on unstructred grids for
More informationFLOWING FLUIDS AND PRESSURE VARIATION
Chapter 4 Pressure differences are (often) the forces that move fluids FLOWING FLUIDS AND PRESSURE VARIATION Fluid Mechanics, Spring Term 2011 e.g., pressure is low at the center of a hurricane. For your
More informationBubble Dynamics using Free Surfaces in a VOF framework
Bubble Dynamics using Free Surfaces in a VOF framework Sevilla Meeting: Numerical Challenges in two-phase flows 27 October 2014 Léon MALAN Supervisor: Prof. Stéphane ZALESKI 1 Contents Introduction: Simulation
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationA Survey of Interface Tracking Methods in Multi-phase Fluid Visualization
A Survey of Interface Tracking Methods in Multi-phase Fluid Visualization Fang Chen 1 and Hans Hagen 2 1 University of Kaiserslautern AG Computer Graphics and HCI Group, TU Kaiserslautern, Germany chen@cs.uni-kl.de
More informationStrömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
More informationA brief description of the particle finite element method (PFEM2). Extensions to free surface
A brief description of the particle finite element method (PFEM2). Extensions to free surface flows. Juan M. Gimenez, L.M. González, CIMEC Universidad Nacional del Litoral (UNL) Santa Fe, Argentina Universidad
More informationDirections: 1) Delete this text box 2) Insert desired picture here
Directions: 1) Delete this text box 2) Insert desired picture here Multi-Disciplinary Applications using Overset Grid Technology in STAR-CCM+ CD-adapco Dmitry Pinaev, Frank Schäfer, Eberhard Schreck Outline
More informationNew formulations of the semi-lagrangian method for Vlasov-type equations
New formulations of the semi-lagrangian method for Vlasov-type equations Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 17 September 2008 In collaboration
More informationFree Surface Flow Simulations
Free Surface Flow Simulations Hrvoje Jasak h.jasak@wikki.co.uk Wikki Ltd. United Kingdom 11/Jan/2005 Free Surface Flow Simulations p.1/26 Outline Objective Present two numerical modelling approaches for
More informationHigh-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method
High-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method Per-Olof Persson University of California, Berkeley, Berkeley, CA 9472-384, U.S.A. We study some of the properties
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationSection 20.1: Choosing a General Multiphase Model. Section 20.2: Volume of Fluid (VOF) Model. Section 20.6: Setting Up a General Multiphase Problem
Chapter 20. Models General Multiphase This chapter discusses the general multiphase models that are available in FLUENT. Chapter 18 provides a brief introduction to multiphase modeling, Chapter 19 discusses
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationKEY STAR TECHNOLOGIES: DISPERSED MULTIPHASE FLOW AND LIQUID FILM MODELLING DAVID GOSMAN EXEC VP TECHNOLOGY, CD-adapco
KEY STAR TECHNOLOGIES: DISPERSED MULTIPHASE FLOW AND LIQUID FILM MODELLING DAVID GOSMAN EXEC VP TECHNOLOGY, CD-adapco INTRODUCTION KEY METHODOLOGIES AVAILABLE IN STAR-CCM+ AND STAR-CD 1. Lagrangian modelling
More informationConservative Semi-Lagrangian solvers on mapped meshes
Conservative Semi-Lagrangian solvers on mapped meshes M. Mehrenberger, M. Bergot, V. Grandgirard, G. Latu, H. Sellama, E. Sonnendrücker Université de Strasbourg, INRIA Grand Est, Projet CALVI, ANR GYPSI
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
More informationNumerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods
Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods Decheng Wan 1 and Stefan Turek 2 Institute of Applied Mathematics LS III, University of Dortmund,
More informationImmersed Boundary Method and Chimera Method applied to Fluid-
The numericsacademy Fixed Colloquium IBM on Moving Immersed IBM Boundary Applications Methods : Conclusion Current Status and Future Research Directions 15-17 June 2009, Academy Building, Amsterdam, the
More informationHigh-order, conservative, finite difference schemes for computational MHD
High-order, conservative, finite difference schemes for computational MHD A. Mignone 1, P. Tzeferacos 1 and G. Bodo 2 [1] Dipartimento di Fisica Generale, Turin University, ITALY [2] INAF Astronomic Observatory
More informationDevelopment of Hybrid Fluid Jet / Float Polishing Process
COMSOL Conference - Tokyo 2013 Development of Hybrid Fluid Jet / Float Polishing Process A. Beaucamp, Y. Namba Dept. of Mechanical Engineering, Chubu University, Japan Zeeko LTD, United Kingdom Research
More informationLevel Set Calculations for Incompressible Two-Phase Flows on a Dynamically Adaptive Grid
Journal of Scientific Computing, Vol. 31, Nos. 1/2, May 2007 ( 2006) DOI: 10.1007/s10915-006-9119-3 Level Set Calculations for Incompressible Two-Phase Flows on a Dynamically Adaptive Grid Yana Di, 1 Ruo
More informationIMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL
More informationweighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces.
weighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces. joint work with (S. Osher, R. Fedkiw and M. Kang) Desired properties for surface reconstruction:
More informationA Novel Approach to High Speed Collision
A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.
More informationA New Incompressibility Discretization for a Hybrid Particle MAC Grid Representation with Surface Tension
A New Incompressibility Discretization for a Hybrid Particle MAC Grid Representation with Surface Tension Wen Zheng, Bo Zhu, Byungmoon Kim, Ronald Fedkiw Stanford University, 353 Serra Mall Room 207, Stanford,
More informationSpatially Adaptive Techniques for Level Set Methods and Incompressible Flow
Spatially Adaptive Techniques for Level Set Methods and Incompressible Flow Frank Losasso Ronald Fedkiw Stanley Osher May 3, 2005 Abstract Since the seminal work of [92] on coupling the level set method
More informationSemi-Conservative Schemes for Conservation Laws
Semi-Conservative Schemes for Conservation Laws Rosa Maria Pidatella 1 Gabriella Puppo, 2 Giovanni Russo, 1 Pietro Santagati 1 1 Università degli Studi di Catania, Catania, Italy 2 Università dell Insubria,
More informationCloth Simulation. COMP 768 Presentation Zhen Wei
Cloth Simulation COMP 768 Presentation Zhen Wei Outline Motivation and Application Cloth Simulation Methods Physically-based Cloth Simulation Overview Development References 2 Motivation Movies Games VR
More informationFully conservative leak-proof treatment of thin solid structures immersed in compressible uids
Fully conservative leak-proof treatment of thin solid structures immersed in compressible uids Jón Tómas Grétarsson, Ron Fedkiw tanford University, 353 erra Mall Room 207, tanford, CA 94305 Abstract We
More informationSolver Settings. Introductory FLUENT Training ANSYS, Inc. All rights reserved. ANSYS, Inc. Proprietary
Solver Settings Introductory FLUENT Training 2006 ANSYS, Inc. All rights reserved. 2006 ANSYS, Inc. All rights reserved. 5-2 Outline Using the Solver Setting Solver Parameters Convergence Definition Monitoring
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More information