Available online at ScienceDirect. Procedia IUTAM 15 (2015 )
|
|
- Claude Lawrence
- 6 years ago
- Views:
Transcription
1 Available online at ScienceDirect Procedia IUTAM 15 (2015 ) IUTAM Symposium on Multiphase flows with phase change: challenges and opportunities, Hyderabad, India (December 08 December 11, 2014) Recent Numerical and Algorithmic Advances within the Volume Tracking Framework for Modeling Interfacial Flows Marianne M. François Los Alamos National Laboratory,Los Alamos, NM 87545, USA LA-UR Abstract A review of recent advances made in numerical methods and algorithms within the volume tracking framework is presented. The volume tracking method, also known as the volume-of-fluid method has become an established numerical approach to model and simulate interfacial flows. Its advantage is its strict mass conservation. However, because the interface is not explicitly tracked but captured via the material volume fraction on a fixed mesh, accurate estimation of the interface position, its geometric properties and modeling of interfacial physics in the volume tracking framework remain difficult. Several improvements have been made over the last decade to address these challenges. In this paper, the multimaterial interface reconstruction method via power diagram, curvature estimation via heights and mean values and the balanced-force algorithm for surface tension are highlighted The The Authors. Authors. Published Published by by Elsevier Elsevier B.V. B.V. This is an open access article under the CC BY-NC-ND license ( Peer-review Peer-review under under responsibility responsibility of Indian of Indian Institute Institute of Technology, of Technology, Hyderabad. Hyderabad. Keywords: volume tracking, volume-of-fluid, curvature, interface reconstruction, surface tension 1. Introduction Interfacial flows are multi-material flows comprised of two or more immiscible materials demarcated by interfaces. They are ubiquitous in nature and are encountered in many industrial applications. A prototypical example of such flows with phase change occurs in materials processing (casting, welding, additive manufacturing, etc.), where liquid-solid phase change (solidification/melting) takes place. An illustration of a simulation of laser spot welding using the Truchas code 1,2 is shown in Figure 1. High-fidelity simulations of interfacial flows require development of physical models and development of accurate numerical methods and robust algorithms. The position of the interface between the fluids and the interface physics needs to be predicted as part of the solution of The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license ( Peer-review under responsibility of Indian Institute of Technology, Hyderabad. doi: /j.piutam
2 Marianne M. François / Procedia IUTAM 15 ( 2015 ) the flow equations. In addition, the discontinuity in fluid properties at the interface and interfacial physics (such as surface tension, phase change) present additional numerical challenges. There is a range of numerical techniques to model interfacial flows at the continuum level. They can be classified into three categories: (1) Lagrangian method, (2) mixed Eulerian-Lagrangian and (3) Eulerian method. In Lagrangian methods, the mesh on which the governing equations are solved adapts and moves with the interface. Lagrangian methods are particularly well suited and accurate for problems with small deformation (e.g. for fluidstructure interaction problem). In mixed Eulerian-Lagrangian methods, the governing equations are solved on a fixed background mesh and the interface is represented by marker points and moved in a Lagrangian way. Front-tracking methods 3 and immersed boundary methods 4 belong to the mixed Eulerian-Lagrangian class of method. These methods allow accurate representation of the interface, however they rely on ad-hoc algorithms to handle change in interface topology and ensure mass conservation. In Eulerian methods, the governing equations are solved on a fixed mesh and the interface is represented through an indicator function that is defined on each computational cell of the mesh and is advected along with the fluid velocity through an advection equation. Level-set 5,6 and volume tracking 7,8,9 (also known as Volume-Of-Fluid, VOF) methods belong to this category. These methods can handle large deformation and account for change in interface topology. However, because the interface is not explicitly tracked but captured, accurate estimation of the interface geometric properties, its advection and modeling of interfacial physics such as surface tension are challenging. In the level-set method, the indicator function is the distance function from the interface. In volume tracking method, the indicator function is the material volume fractions. The values of the volume fractions are between 0 and 1, as they represent the amount of a given material volume with respect to the computational cell volume. By its nature/design, volume tracking methods are conservative, whereas level-set methods are not and require correction algorithm. In this paper, we focus on the presentation of recent improvements that we have proposed in the numerical methods and algorithms within the context of the volume tracking method. Specifically, we highlight multimaterial interface reconstruction via power diagram, curvature estimation via heights and mean values and surface tension modelling via balanced-force algorithm. (a) (b) Fig. 1. Simulation of laser spot welding of 304 stainless steel with Truchas 1. (a) Temperature contours of a moving heat source on stainless steel piece; (b) Weld pool with flow circulation shown by velocity vector. Nomenclature Subscripts 1 fluid 1 2 fluid 2 c cell centered f face centered k material/fluid index x derivative with respect to x (horizontal direction) Superscripts n time level (N) normal component
3 272 Marianne M. François / Procedia IUTAM 15 ( 2015 ) (T) δ γ κ μ ρ σ A f, F n P t T Y tangential component delta function surface tension gradient interface curvature fluid dynamic viscosity fluid density surface tension computational cell face area volume fraction surface force gravitational acceleration unit interface normal computational cell face normal pressure time temperature fluid velocity vector height function in vertical direction 2. Multimaterial interface reconstruction via power diagram In volume tracking method, the interface position is not known, since the interface is not explicitly tracked by marker points, but is captured through the volume fraction. Knowing the interface position accurately is of primary importance in volume tracking as it is used to estimate the material fluxes geometrically to avoid numerical diffusion. These material fluxes are then used to estimate consistent mass and momentum advection terms that is of primary importance for simulating large density ratio 8. Finding the interface position given the volume fraction is commonly named interface reconstruction. Given the volume fractions in a computational cell, the goal is to find the interface position that matches the volume fractions assuming the shape of the interface. Early volume tracking method 10,11 used a simple interface reconstruction by assuming that the interface is aligned with one of the coordinate axis. This is referred to as simple line interface calculation (SLIC). In 1982, Youngs 12 extended the method to allow the interface to have an orientation within the cell based on the interface normal. The interface normal is calculated as the gradient of the volume fraction using finite-difference formulas, Green Gauss formula or least-square technique. This interface representation is referred to as piecewise linear interface calculation (PLIC). A parabolic interface reconstruction was proposed by Renardy and Renardy 13, but is much more complex to implement. Several studies have been conducted in PLIC area to improve the interface normal calculation 14, to generalize the method to unstructured meshes 7, to seek analytical solution on arbitrary mesh 15 and to develop smoothing algorithm to avoid discontinuity of interface planes across cells 16,20. While many improvements have been made in PLIC method, they have mostly focused on single interface reconstruction, i.e. interface between two materials. Reconstructing multiple interfaces within a computational cell (i.e. where more than two materials are present) has received little attention. The most common approach for multimaterial interface reconstruction has been to treat material one-by-one leading to a reconstruction that is strongly dependent on the order in which the materials are treated. The onionskin approach 12,17 consists of considering each material in a specific order, by cumulating the volume fractions in order to reconstruct the material interface sequentially. For the first interface, the interface is reconstructed using the volume fraction from the first material. Then, the volume fraction of the first material is added to the volume fraction of the second material and the interface between the first two materials and the remainder ones is reconstructed. This procedure is repeated until all interface planes have been reconstructed, which corresponds to
4 Marianne M. François / Procedia IUTAM 15 ( 2015 ) the number of material present in the cell minus one. The onion-skin technique works well for layer structures (for interface planes that are close of being parallel) but can result in overlapping layers for other cases. A more general and correct approach is the nested dissection method 18, where each material is separated from the others in a specified order. In the method, a pure polygon (or polyhedron) representing the first material is marked out from the cell, leaving a mixed polygon for the remaining materials. Then, a polygon representing the second material is marked out from the mixed polygon and the process continues until the last material is treated. With the correct material ordering, the interface reconstructed by one of the above methods is close to the correct configuration. However an incorrect ordering results in substantial degradation of the interface during advection. Motivated by the need to develop a material-order independent multimaterial interface reconstruction, Schofield et al 19,20 proposed a unique method that employs power diagram (weighted Voronoi diagram) to generate the multiple interfaces. In this method, the positions of materials in the computational cells are first approximated (either using particles or a slope-limited linear reconstruction of volume fraction). These approximate positions then serve as the generators to build the weighted Voronoi diagram. The weights are then adjusted iteratively in order to match the volume fractions of each material. The power diagram then results in a unique partition of the computational cells into multi-material regions. All materials are treated simultaneously and, as such, have no material order dependency. An illustration of a triple point configuration and its reconstruction using power diagram is shown in Figure 2. (a) (b) Fig. 2. Triple point configuration. (a) Exact interface; (b) Material independent interface reconstruction via power diagram. 3. Curvature estimation via heights and mean values Another challenge in the volume tracking framework, is the calculation of interface curvature. Estimating curvatures, involve taking spatial second-order derivatives of abruptly varying volume fractions. Various methods are presented and compared in Cummins et al 21. In this section, we review how to calculate interface curvature via heights and mean values. This is also referred to as the height function method. The height function technique 6,9,21 consists of computing curvature to second-order accuracy given integral mean values (or columns heights) associated with three adjacent columns on uniform or nonuniform orthogonal grid and to fourth-order accuracy given mean values associated with five adjacent columns 22. In multi-dimension, the mean values are computed in the dominant direction of the interface normal component. Once the mean values are computed, the curvature is calculated by differentiating the mean values (Y) as: (1) In two-dimension, the stencil used to compute the height function is 3x7, and in three-dimension, the stencil is 3x3x7. An illustration of the two-dimensional case is shown in Figure 3. Recent studies on the height function technique have focused on improving the estimation of the height function itself by choosing the best local rectangular set of mesh squares 23,24,25,26.
5 274 Marianne M. François / Procedia IUTAM 15 ( 2015 ) Fig. 3. Illustration of the height function method in three vertical columns with volume fractions shown in each computational cell. 4. Balanced-Force Algorithm for Modeling Surface Tension Modeling surface tension force, an interfacial force, in the volume tracking framework represents some challenges, since the interface is embedded within the computational mesh and not explicitly tracked. In this section, we present the balanced forced-algorithm of Francois et al 9 to model surface tension driven flow for the case of two fluids. We first review the governing equations for the fluid flow and volume tracking, and then present the formulation for surface tension force and the balanced-force algorithm Governing Equations The governing equations consist of the volume tracking equation, also known as volume fraction advection equation and of the mass and momentum conservation equations. The flow is assumed to be incompressible and governed by a single set of mass and momentum conservation equations on a fixed grid. The fluids are assumed to be Newtonian and immiscible and denoted by subscripts 1 and 2. In computational cells occupied with fluid 1, the volume fraction, f is unity, and in cells occupied with fluid 2, f is zero. For cells containing the interface bounding fluid 1 and 2, f lies between zero and unity. The volume fractions of fluid 1 and 2 sum to unity in every computational cell. The governing equations are the equation for the advection of the volume fraction f : f t + u f = 0 (2) and the mass and momentum conservation equations: u = 0 (3)
6 ( ) ρu t + ρu u Marianne M. François / Procedia IUTAM 15 ( 2015 ) ( ) = P + μ ( ( u + T u ))+ ρg + F S (4) where ρ and μ are the fluid density and viscosity, respectively, defined as: ρ = ρ 1 f + ρ 2 1 μ = μ 1 f + μ 2 1 ( f ) ( f ) (5) (6) The expression for density in Eq. (5) results from mass conservation, whereas the expression for mixture viscosity in Eq. (6) is a model. The system of equations (Eq. 3-4) is solved using a pressure correction projection method as described in Francois et al 9 and the volume fraction advection equation (Eq. 2) is solved by first geometrically reconstructing the interface by piecewise linear interface planes (PLIC) and secondly by geometrically calculating the advection fluxes to keep a sharp representation of the interface Surface Tension Models The Continuum Surface Force (CSF) method 27 is widely used in VOF, level-set, front-tracking methods. It consists of transforming the interfacial force to volume force via delta function. This method allows a smooth and continuous representation of the interfacial jump condition. Another method that has mainly been used in the context of level-set function but was extended by Francois et al 9 for volume tracking method is the Ghost Fluid Method (GFM) 28,29. In this method, the interfacial jump conditions due to surface tension are imposed via ghost values. This method allows a sharp representation of the interfacial jumps. Both methods are known to generate unphysical flow near the interface spurious currents particularly when surface tension forces are dominant. We have identified the source for these spurious currents (numerical errors) to several main factors. The first factor is the coupling of the surface tension force with flow algorithm. The formulation should seek proper balance of forces in the discretized momentum equation. If the flow algorithm ensures proper balance of forces, there should be no spurious current, as shown in Francois et al 9 when exact interface curvature is prescribed. The second factor is accurate estimations of curvature from volume fractions. Using the second-order accurate height function method in conjunction with the balanced-force algorithm resulted in drastic reduction of spurious currents. We now detail the continuum surface force (CSF) model of Brackbill et al 27 and its implementation using the balanced-force algorithm of Francois et al 9. The surface tension force is given as a volumetric force: F S = ( σκ ˆn + s σ )δ where σ is the surface tension, κ is the interface curvature, n ˆ is the interface unit normal vector, s is the surface (N) gradient, and δ is the delta function. The surface tension force is decomposed into a normal F S and tangential F S (T) volumetric force: where F S = F S (N) + F S (T) F S (N) = σκ ˆ n δ = σκ f F (T) S = ( s σ )δ = ( σ ˆn ( ˆn. )σ )δ (10) (7) (8) (9)
7 276 Marianne M. François / Procedia IUTAM 15 ( 2015 ) If the surface tension coefficient varies with temperature, i.e. σ = σ ref + γ ( T T ref ) with γ = σ T (11) the tangential component will be non-zero and it will generate thermo-capillary forces along the interface also known as Marangoni effects. In our balanced-force algorithm 9, both the normal and tangential surface tension force terms appear in the predictor step of the flow algorithm: ρ n+1 uc * = ρ n uc n Δt f u c n k ρ k f k u n f. n n+1 ( f )A f Δtρ f P c + Δt μ f A f nf u + T u + ΔtF (T ) Sc f ( ) ρ f (N ) F Sf ρ f n f c (12) The normal component of surface tension needs to balance the pressure jump across the interface. Therefore it is computed on the computational cell faces since the pressure gradient is evaluated at faces. This is to ensure proper balance of forces. Similarly, the tangential component of surface tension needs to balance the viscous stresses. They are both computed at cell centers. Verification and validation test cases of our surface tension force implementation can be found in Francois et al 9 for the normal component of surface tension and in Francois et al 30 for the implementation of the tangential surface tension component. Our surface tension model is implemented within the flow solver of the Truchas code. 4. Summary Over the last decade, substantial improvements have been made in the volume tracking method. In this paper, we have highlighted three recent advances: (1) a material-order independent multimaterial interface reconstruction, (2) the height function method to calculate interface curvature from volume fractions information and (3) the balancedforce algorithm for surface tension. These improvements have resulted in substantial reduction of numerical errors (drastic reduction of spurious currents in surface tension dominated flow) and reduction in the number of numerical input parameters from the user (material ordering and smoothing length to compute curvature). More accurate efficient, robust computational methods are always sought as well as modeling additional physics. Acknowledgements At LANL, I have had the privilege to work with outstanding colleagues and develop collaborations. In particular, I would like to thank Sharen Cummins, Edward Dendy, Rao Garimella, Douglas Kothe, Kin Lam, Raphael Loubere, Lucie Parietti, Sam Schofield, Jim Sicilian, Blair Swartz, Matthew Williams. I would also like to thank the organizers of the IUTAM symposium for their invitation. This work was performed under the auspices of the National Nuclear Security Administration of the US Department of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA References 1. Parietti L, Lam K, Validation of a multiphysics software package for simulations of pulsed laser welding at Los Alamos National Laboratory, Technical Report LA-UR , Los Alamos National Laboratory, The Telluride Team, Truchas: Physics and Algorithms, Technical Report LA-UR , Los Alamos National Laboratory 2003.
8 Marianne M. François / Procedia IUTAM 15 ( 2015 ) Tryggvason G, Bunner B, Esmaeeli A, Juric D, Al-Rawahi N, Tauber W, Han J, Nas S, Jan YJ, A Front-Tracking Method for the Computations of Multiphase Flow, J Comput Phys 2001;169: Francois M, Shyy W, Computations of drop dynamics with the immersed boundary method, part 1: numerical algorithm and buoyancyinduced effect, Numerical Heat Transfer: Part B: Fundamentals 2003;44: Osher S, Fedkiw RP, Level set methods: an overview and some recent results, J Comput Phys 2001;169: Sussman M, A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles, J Comput Phys 2003;187: Rider WJ, Kothe DB. Reconstructing volume tracking. J Comput Phys 1998;141: Bussmann M, Kothe DB, Sicilian JM, Modeling high density ratio incompressible interfacial flows, Technical Report FEDSM , in Proceedings of the Sixth ASME/JSME Joint Fluids Engineering Conference, July 14-18, Montreal, Canada, Francois MM, Cummins SJ, Dendy ED, Kothe DB, Sicilian JM, Williams MW., A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework, J Comput Phys 2006;213: Hirt CW, Nichols BD. Volume of fluid (VOF) method for the dynamcs of free boundaries. J Comput Phys 1981;39: Noh WF, Woodward P. SLIC (Simple line interface calculation). In Lecture Notes in Physics, vol. 59, van de Vooren AI, Zandbergen PJ (eds). Springer: Berlin, 1976; Youngs DL. Time dependent multi-material flow with large fluid distortion. In Numerical Methods for Fluid Dynamics, Morton KW, Baines MJ (eds). Academic Press: New York, 1982; Renardy Y, Renardy M, PROST: A parabolic reconstruction of surface tension for the volume-of-fluid method, J Comput Phys 2002;183: Pilliod Jr. JE, Puckett EG, Second-order accurate volume-of-fluid algorithms for tracking material interfaces, J Comput Phys 2004;199: Diot S, François MM, Dendy ED, An interface reconstruction method based on analytical formulae for 2D planar and axisymmetric arbitrary convex cells, J Comput Phys 2014;275: Swartz BK, The second-order sharpening of blurred smooth borders, Mathematics of Computation 1989;52: Choi BY, Bussmann M. A piecewise linear approach to volume tracking a triple point, Int J Numer Meth Fluids 2007; 53: Dyadechko V, Shashkov MJ. Multi-material interface reconstruction from moment data, Technical Report LA-UR , Los Alamos National Laboratory, Schofield SP, Garimella RV, Francois MM, Loubere R, Material order-independent interface reconstruction using power diagrams, Int J Numer Meth Fluids 2008;56: Schofield SP, Garimella RV, Francois MM, Loubere R, A second-order accurate material-order-independent interface reconstruction technique for multi-material flow simulations, J Comput Phys 2009;228: Cummins SJ, Francois MM, Kothe DB, Estimating curvature from volume fractions, Computers and Structures 2005;83: Francois MM, Swartz BK, Interface curvature via volume fractions, heights, and mean values on nonuniform rectangular grids, J Comput Phys 2010;229: Veldman AEP, Gerrits J, Luppes R, Helder JA, Vreeburg JPB, The numerical simulation of liquid sloshing on board spacecraft, J Comput Phys 2007;224: Afkhami S, Bussmann M, Height functions for applying contact angles to 2D VOF simulations, Int J Numer Meth Fluids 2008;57: Popinet S, An accurate adaptive solver for surface-tension-driven interfacial flows, J Comput Phys 2009;228: Hernandez J, Lopez J, Gomez P, Zanzi C, Faura F, A new volume of fluid method in three dimensions Part I: Multidimensional advection method with face-matched flux polyhedra, Int J Numer Meth Fluids 2008;58: Brackbill JU, Kothe DB, and Zemach C. A continuum method for modeling surface tension, J Comput Phys 1992;100: Kang M., Fedkiw RP, Liu XD, A boundary condition capturing method for multiphase incompressible flow, J Sci Comput 2000;15: Liu XD, Fedkiw RP, Kang M, A boundary condition capturing method for Poisson s equation on irregular domains, J Comput Phys 2000;160: Francois MM, Sicilian JM, Kothe DB, Modeling of thermocapillary forces within a volume tracking algorithm. In Modeling of Casting, Welding and Advanced Solidification processes XI. Charles-Andre Gandin and Michel Bellet, TMS, 2006.
A 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 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 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 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 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 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 informationSIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH
Fifth International Conference on CFD in the Process Industries CSIRO, Melbourne, Australia 3-5 December 006 SIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH Manasa Ranjan Behera and Murali.
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 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 informationHeight functions for applying contact angles to 3D VOF simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2009; 61:827 847 Published online 8 December 2008 in Wiley InterScience (www.interscience.wiley.com)..1974 Height functions
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 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 informationThis tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:
Tutorial 22. Modeling Solidification Introduction This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following: Define a
More informationAvailable online at ScienceDirect. Procedia Engineering 136 (2016 ) Dynamic analysis of fuel tank
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 136 (2016 ) 45 49 The 20 th International Conference: Machine Modeling and Simulations, MMS 2015 Dynamic analysis of fuel tank
More informationAccurate interface normal and curvature estimates on three-dimensional unstructured non-convex polyhedral meshes
Center for Turbulence Research Annual Research Briefs 2014 53 Accurate interface rmal and curvature estimates on three-dimensional unstructured n-convex polyhedral meshes By C. B. Ivey AND P. Moin 1. Motivation
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 informationOverview 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 informationMoment of fluid method for multimaterial flows
Moment of fluid method for multimaterial flows M. Jemison, M. Sussman, M. Shashkov Department of Mathematics, Florida State University October 1, 2012 Recent Motivating Work Volume preserving Moment-of-Fluid
More informationIMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING
Image Anal Stereol 2001;20:143-148 Original Research Paper IMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING DAVID GARCIA 1, GUY COURBEBAISSE 2 AND MICHEL JOURLIN 3 1 European Polymer Institute (PEP),
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 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 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 informationSimulation of Transient and Three-Dimensional Coating Flows Using a Volume-of-Fluid Technique
Simulation of Transient and Three-Dimensional Coating Flows Using a Volume-of-Fluid Technique C.W. Hirt and J.E. Richardson, Flow Science, Inc., 1257 40th, Los Alamos, NM 87544 and Ken S. Chen, Sandia
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 informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
More informationAvailable online at ScienceDirect. Procedia Engineering 111 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 111 (2015 ) 115 120 XXV R-S-P seminar, Theoretical Foundation of Civil Engineering (24RSP) (TFoCE 2015) Numerical simulation
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 informationMarker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 1, pp. 1-41 Commun. Comput. Phys. January 2009 REVIEW ARTICLE Marker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
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 informationRobust Simulation of Sparsely Sampled Thin Features in SPH-Based Free Surface Flows
Copyright of figures and other materials in the paper belong to original authors. Robust Simulation of Sparsely Sampled Thin Features in SPH-Based Free Surface Flows Xiaowei He et al. ACM SIGGRAPH 2015
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationTwo-dimensional Numerical Simulation of Multiphase Flows in MATLAB
Two-dimensional Numerical Simulation of Multiphase Flows in MATLAB Saadi Daftari A Thesis in The Department of Mechanical and Industrial Engineering Presented in Partial Fulfillment of the Requirements
More informationOn the thickness of discontinuities computed by THINC and RK schemes
The 9th Computational Fluid Dynamics Symposium B7- On the thickness of discontinuities computed by THINC and RK schemes Taku Nonomura, ISAS, JAXA, Sagamihara, Kanagawa, Japan, E-mail:nonomura@flab.isas.jaxa.jp
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 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 informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationPROST: A Parabolic Reconstruction of Surface Tension for the Volume-of-Fluid Method
Journal of Computational Physics 183, 400 421 (2002) doi:10.1006/jcph.2002.7190 PROST: A Parabolic Reconstruction of Surface Tension for the Volume-of-Fluid Method Yuriko Renardy 1 and Michael Renardy
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 informationDeforming fluid domains within the finite element method
Deforming fluid domains within the finite element method Five mesh-based tracking methods in comparison S. Elgeti H. Sauerland Aachen Institute for Advanced Study in Computational Engineering Science Financial
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 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 informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationAnimation of Bubbles in Liquid
EUROGRAPHICS 2003 / P. Brunet and D. Fellner (Guest Editors) Volume 22 (2003), Number 3 Animation of Bubbles in Liquid Jeong-Mo Hong and Chang-Hun Kim Department of Computer Science, Korea University Abstract
More informationNumerical simulations of dam-break floods with MPM
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (206) 000 000 www.elsevier.com/locate/procedia st International Conference on the Material Point Method, MPM 207 Numerical
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 informationA Method for Capturing Sharp Fluid Interfaces on Arbitrary Meshes
Journal of Computational Physics 153, 26 50 (1999) Article ID jcph.1999.6276, available online at http://www.idealibrary.com on A Method for Capturing Sharp Fluid Interfaces on Arbitrary Meshes O. Ubbink
More informationVOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 67 69 (997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS MURRAY RUDMAN Commonwealth Scientific and Industrial Research Organisation
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationTutorial 17. Using the Mixture and Eulerian Multiphase Models
Tutorial 17. Using the Mixture and Eulerian Multiphase Models Introduction: This tutorial examines the flow of water and air in a tee junction. First you will solve the problem using the less computationally-intensive
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 informationMESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP
Vol. 12, Issue 1/2016, 63-68 DOI: 10.1515/cee-2016-0009 MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP Juraj MUŽÍK 1,* 1 Department of Geotechnics, Faculty of Civil Engineering, University
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 informationCharacteristic Aspects of SPH Solutions
Characteristic Aspects of SPH Solutions for Free Surface Problems: Source and Possible Treatment of High Frequency Numerical Oscillations of Local Loads. A. Colagrossi*, D. Le Touzé & G.Colicchio* *INSEAN
More informationCFD Simulation and Validation of Flow Pattern Transition Boundaries during Moderately Viscous Oil-Water Two-Phase Flow through Horizontal Pipeline
Vol:7, No:, 203 CFD Simulation and Validation of Flow Pattern Transition Boundaries during Moderately Viscous Oil-Water Two-Phase Flow through Horizontal Pipeline Anand B. Desamala, Anjali Dasari, Vinayak
More informationmidas NFX 2017R1 Release Note
Total Solution for True Analysis-driven Design midas NFX 2017R1 Release Note 1 midas NFX R E L E A S E N O T E 2 0 1 7 R 1 Major Improvements Midas NFX is an integrated finite element analysis program
More informationThermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA
Thermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA Jingxiao Xu 1, Jason Wang 2 1 LSTC 2 LSTC 1 Abstract Smooth particles hydrodynamics is a meshfree, Lagrangian particle method
More informationInterfaces and multicomponent fluids
Interfaces and multicomponent fluids Junseok Kim and John Lowengrub Department of Mathematics University of California, Irvine, USA E-mail: lowengrb@math.uci.edu 20 December 2004 Introduction Many important
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 information2D numerical simulation of ocean waves
2D numerical simulation of ocean waves Qingjie. Du,*, Y.C. Dennis. Leung Department of Mechanical Engineering, The University of Hong Kong, Hong Kong, China * Corresponding author. Tel: +852 51743593,
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 informationMicrowell Mixing with Surface Tension
Microwell Mixing with Surface Tension Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract For many applications in the pharmaceutical
More informationAn accurate adaptive solver for surface-tension-driven interfacial flows
An accurate adaptive solver for surface-tension-driven interfacial flows Stéphane Popinet To cite this version: Stéphane Popinet. An accurate adaptive solver for surface-tension-driven interfacial flows.
More informationLS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.
12 th International LS-DYNA Users Conference FSI/ALE(1) LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA Part 1 Facundo Del
More informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
More informationAn accurate numerical method for computing surface tension forces in CFD codes
1 An accurate numerical method for computing surface tension forces in CFD codes Numerical experiments with surface tension M. Coquerelle 1 S. Glockner 1 1 Numerical Fluid Mechanics (MFN), département
More informationFinite Volume Methodology for Contact Problems of Linear Elastic Solids
Finite Volume Methodology for Contact Problems of Linear Elastic Solids H. Jasak Computational Dynamics Ltd. Hythe House 200 Shepherds Bush Road London W6 7NY, England E-mail: h.jasak@cd.co.uk H.G. Weller
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
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 informationCFD Simulation for Stratified Oil-Water Two-Phase Flow in a Horizontal Pipe
CFD Simulation for Stratified Oil-Water Two-Phase Flow in a Horizontal Pipe Adib Zulhilmi Mohd Alias, a, Jaswar Koto, a,b,* and Yasser Mohamed Ahmed, a a) Department of Aeronautics, Automotive and Ocean
More informationLecture 1.1 Introduction to Fluid Dynamics
Lecture 1.1 Introduction to Fluid Dynamics 1 Introduction A thorough study of the laws of fluid mechanics is necessary to understand the fluid motion within the turbomachinery components. In this introductory
More informationThis is an author-deposited version published in : Eprints ID : 19938
Open Archive TOULOUSE Archive Ouverte (OATAO) OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. This is an author-deposited
More informationA Sharp Interface Cartesian Grid Method for Simulating Flows with Complex Moving Boundaries
Journal of Computational Physics 174, 345 380 (2001) doi:10.1006/jcph.2001.6916, available online at http://www.idealibrary.com on A Sharp Interface Cartesian Grid Method for Simulating Flows with Complex
More informationThree Dimensional Numerical Simulation of Turbulent Flow Over Spillways
Three Dimensional Numerical Simulation of Turbulent Flow Over Spillways Latif Bouhadji ASL-AQFlow Inc., Sidney, British Columbia, Canada Email: lbouhadji@aslenv.com ABSTRACT Turbulent flows over a spillway
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 informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More informationA Contact Angle Model for the Parallel Free Surface Lattice Boltzmann Method in walberla Stefan Donath (stefan.donath@informatik.uni-erlangen.de) Computer Science 10 (System Simulation) University of Erlangen-Nuremberg
More informationChapter 2 Water and Bubbles
Chapter 2 Water and Bubbles 2.1 Animation of Bubbles in Liquid Abstract This section introduces a new fluid animation technique in which liquid and gas interact with each other, using the example of bubbles
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
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 informationDynamic Computational Modeling of the Glass Container Forming Process
Dynamic Computational Modeling of the Glass Container Forming Process Matthew Hyre 1, Ryan Taylor, and Morgan Harris Virginia Military Institute, Lexington, Virginia, USA Abstract Recent advances in numerical
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationSliding and multi-fluid velocities in Staggered Mesh Eulerian and MMALE codes
Sliding and multi-fluid velocities in Staggered Mesh Eulerian and MMALE codes Gabi Luttwak 1 1 Rafael, P.O. Box 2250, Haifa 31021, Israel The Velocity in Eulerian and MMALE Codes Most Eulerian and MMALE
More informationINTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS
INTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS Yoshinobu Yamamoto Department of Nuclear Engineering, Kyoto University Yoshida Sakyo Kyoto, 66-85,
More informationNumerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
More informationModelling of Levitation Melting using a Fixed Mesh Method
International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, October 27-29, 2008 Modelling of Levitation Melting using a Fixed Mesh Method D. Hectors, E. Toorman, K. Van Reusel
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 13: The Lecture deals with:
The Lecture deals with: Some more Suggestions for Improvement of Discretization Schemes Some Non-Trivial Problems with Discretized Equations file:///d /chitra/nptel_phase2/mechanical/cfd/lecture13/13_1.htm[6/20/2012
More informationA Finite Particle Approach for the Simulation of Multiphase Flows
ILASS-Americas 29th Annual Conference on Liquid Atomization and Spray Systems, Atlanta, GA, May 2017 A Finite Particle Approach for the Simulation of Multiphase Flows E.A. Wenzel and S.C. Garrick Department
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationC. A. D. Fraga Filho 1,2, D. F. Pezzin 1 & J. T. A. Chacaltana 1. Abstract
Advanced Computational Methods and Experiments in Heat Transfer XIII 15 A numerical study of heat diffusion using the Lagrangian particle SPH method and the Eulerian Finite-Volume method: analysis of convergence,
More informationThermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA
Thermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA INTRODUCTION: Jingxiao Xu, Jason Wang LSTC Heat transfer is very important in many industrial and geophysical problems. Many
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More informationEVALUATION OF VOLUME TRACKING ALGORITHMS FOR GAS-LIQUID TWO-PHASE FLOWS
Proceedings of ASME FEDSM 03 4 th ASME_JSME Joint Fluids Engineering Conference Honolulu, Hawaii, USA, July 6-10, 2003 F E DSM200 3-4 5164 EVALUATION OF VOLUME TRACKING ALGORITHMS FOR GAS-LIQUID TWO-PHASE
More informationReconstructing Volume Tracking 1
JOURNAL OF COMPUTATIONAL PHYSICS 141, 112 152 (1998) ARTICLE NO. CP985906 Reconstructing Volume Tracking 1 William J. Rider and Douglas B. Kothe Applied Theoretical and Computational Physics Division,
More information