Sliding and multi-fluid velocities in Staggered Mesh Eulerian and MMALE codes
|
|
- Aubrey Copeland
- 5 years ago
- Views:
Transcription
1 Sliding and multi-fluid velocities in Staggered Mesh Eulerian and MMALE codes Gabi Luttwak 1 1 Rafael, P.O. Box 2250, Haifa 31021, Israel
2 The Velocity in Eulerian and MMALE Codes Most Eulerian and MMALE codes assume a common velocity in multimaterial cells Is this right and what are the consequences? Let explore the price and the benefits of defining separate velocities for each fluid MMALE (Luttwak & Rabie 1985,Luttwak1989), had multifluid velocities. We review its scheme. We introduce the multifluid velocities into The SMG scheme (Luttwak & Falcovitz 2005) Show some results Discuss the advantages price & limitations and look ahead for further improvements
3 Multi-fluid velocities What happens at an interface? In 1D both pressure (or the normal stress) and velocity are continuous In 2D (3D) slip lines (planes) form at density discontinuities: Even in a single fluid And almost always at material interfaces In Lagrangian codes Sliding-Impacting techniques are used to handle such situations
4 Sliding and Impact in Lagrange and SMALE In Lagrangian and almost-lagrangian ALE or SMALE codes a sliding/impacting contact algorithm enables the simulation of sliding, impact and separation Without it: impact/separation cannot be described Due to sliding, a row of cells near the interfaces undergoes nonphysical shear. This leads to a mesh size dependent and un-physical friction. If the relative sliding is large, this mesh deformation stops the calculation.
5 In Eulerian and MMALE codes: The fluids flow through the cells The mesh does not deform due to sliding/impact The calculation can proceed But, with a common interface velocity: Sliding introduces an effective friction and viscosity Therefore the strain field will be wrong near the interface in elastic-plastic flow this affects failure The rate of growth of instabilities (like R-T,R-M) might be slower.
6 MM-Euler and MM-ALE codes using multi-fluid velocities Most MM-Euler and MM-ALE codes use a common vertex velocities, for all species. Remarkable exceptions are: Luttwak & Rabie-1985 One of the first MMALE codes Walker & Anderson 1994 noticed the importance of this issue for elastic plastic flow and introduced cell centered multifluid velocities into the MM-Eulerian code CTH. Vitali & Benson 2006 use Belytschko X-FEM method to locally add multifluid velocity using an extended FEM method for staggered mesh MMALE formulations. Our purpose is to incorporate multi-fluid velocities into the SMG scheme for ALE hydrodynamics (Luttwak & Falcovitz 2005)
7 The old D MMALE code Like the SALE code, it used a staggered mesh in space and time: pressure, stress, density, internal energy and Q are cell centered The velocity is defined at the vertices at T n+1/2 These variables (except the Q) were also defined separately for each fluid. Unlike to SALE and to most other codes, each fluid had a separate vertex velocity and vertex mass. Both the cell centered and vertex material variables used a LIFO stack logic so that memory was used only at cells were a fluid is present Each fluid was accelerated separately: ai p 1 a i i A common normal to interface velocity was enforced by momentum conservation at the vertex.
8 The (1985) 2D MMALE code The normal to the interface and the wet volume in cell uniquely defines the wet polygon for 2 materials. For more than 2 materials the result depend on the order: The wet polygon of the first fluid is found as before. The wet volume for the next fluid defined by its normal and wet volume inside the polygon not filled by the previous fluids
9 MMALE (1984) The Advection Phase The advection phase was based on a remap: The vertices of the wet polygon in cell c filled with a fluid were moved with the fluid velocity there. The mesh was advanced with the grid velocity. i The overlap volume was found V c, nand it is the flux of fluid i from cell c to a nearby cell n. (Luttwak 1994,Luttwak&Cowler 1999) The densities in the flux volume were obtained by linear interpolation between donor and acceptor cell values with the flux values limited to prevent the generation of new extreme values. 1 c n 3
10 Vertex Variables Advection in the old MMALE-2D code Vertex variable update used SALE s method A cell centered mass and momentum was first advected For single material cells, the influx of mass and momentum to a cell were equally distributed to its vertices For multi-material cells, the fluxes were split according to the relative wet volume fraction of each fluid in the corner zone. The method was excessively diffusive and the distribution of the mass to the vertices for multi-material vertices caused problems
11 The Staggered Mesh Godunov (SMG/Q) Scheme SMG schemes for 3D ALE hydrodynamics were presented at Oxford 2005,APS 2003,2005. They use Riemann problem (RP) solutions instead of pseudo-viscosity to capture shocks. The vertex defined velocities may have jumps mid-way, on the faces separating the corner zones, near the cell center. These are simpler, collision RP. (with constant there). p, They are solved in the normal to shock direction, i.e. along the velocity difference. u L u R
12 The Staggered Mesh Godunov (SMG/Q) Scheme Let p* be the RP solution pressure. We look at q p* p cell as a uni-axial tensor pseudo-viscosity acting on the corner zone internal face along the shock direction Its impulse is exchanged between the two neighboring vertices. Its work, which must be dissipative and is added to the zone internal energy. The above procedure is similar to Caramana et al. compatible hydroscheme. The SMG scheme has some similarities with Christensen-1990 split-q and Caramana 1997 et al. edge viscosity.
13 The Staggered Mesh Godunov (SMG/Q) Scheme The velocities on both sides of the face, serving as data for the RP are evaluated from the vertex values using a, cell centered velocity gradient limited, to preserve a smooth velocity distribution. This serves as a shock detector. The resulting scheme captures sharp and smooth shocks, with an additional strong and natural mesh stabilizing effect. u L u L Limited slope u R C R Original slope Velocity jump for RP
14 Multi-material velocities The fluids have distinct vertex defined velocities and masses On a material interface a common normal to interface velocity component is enforced. When appropriate (fluid flowing away from the interface, negative pressure or normal stress, large void fraction) these normal velocities are left unchanged enabling interface separation The fluid vertex velocities are only stored for vertices sharing at least one cell containing that fluid.
15 The Set-up of the SMG-MMALE code Geometry Update (not necessary for pure Euler) compute face areas, volumes and volume changes. Lagrangian Phase The SMG Scheme Definition of pressure and stress in a MM zone Advection Phase (not present for pure Lagrange) Second order advection for cell-centered variables Staggered Mesh Momentum Advection HIS (Benson) Interface tracking: VOF (Luttwak 1984,2002) similar to Youngs scheme Grid Motion Algorithm Almost all parts need some change.
16 Lagrangian Phase Strain, velocity gradient and volume changes are specific for each fluid in the cell: To find v the velocity is integrated over the full cell. v, v, V If a fluid has a very small vertex mass, an extrapolated velocity is used obtained from the fluid velocity average in the cell. This prevents unphysical large gradients when the wet volume is small. A common vertex acceleration is computed integrating the pressure (without Q) over the (dual) cell using Wilkins scheme. The specific fluid accelerations are found using: a i 1 p a i i
17 SMG The collision RP is solved separately for each fluid The RP data uses the fluid velocities and densities The forces act only on the particular fluid momentum The associated work is added to that fluid internal energy
18 Interface Tracking For each fluid we find the normal to the interface: The normal is found by integrating the relative partial volumes over the cell faces. nˆ ; V i i V The normal and the wet volume V i determine the interface position and the wet part of the cell. Unlike to the 1985 MMALE 2D, we are ignoring possible overlap when several materials present in a cell. The normal at a vertex is found by integrating passing through the vertex (like for acceleration) i i i i over the faces
19 The Advection Phase The volume flux between two face neighbor cells is defined by the wet part of the volume enclosed by the common face old and new positions. The new face is obtained by moving its vertices with the fluid velocity relative to the grid. This is different for each fluid. The wet part is defined by cutting it with the interface plane found by the tracking algorithm advanced to T n+1. The remap algorithm used in the older 2D code is superior as it took into account diagonal advection but is more expensive L R
20 Vertex Variable Advection Similar to Benson HIS m L m, Treating I,J,K pair of faces separately, let L m R be the mass influx from left and right into a cell, ml m R is advected from the left face vertices (1,2,6,5 in the picture) to the corresponding vertices on the right. Each vertex takes an equal share. (A quarter for the hexa cell shown). Momentum advection is done separately for each pair of vertices, the velocity of the fluxed mass is evaluated using the donor vertex velocity and the cell limited velocity gradient. We keep the same scheme for the fluid fluxes. Thus each vertex shares an equal part of the cell fluid mass. However, until the vertex is wet (it is in the cell wet part) its velocity will be modified only by advection
21 LIFO stack memory set-up A structure (type in Fortran 90 ) is defined with all the fluid specific vertex variables. An array of this structure is allocated. The unused elements are connected in a linked list. When a fluid leaves all the cells surrounding the vertex its element is freed and put on the top of the stack to be first used. Each vertex points to an element of the array storing the data for the first fluid. The element itself points to the element containing the next fluid data or to a Null element (All empty vertices also point to the Null element). When there are no more free elements, another array is allocated with the elements linked in a stack. A similar logic is used to store the vertex variable fluxes. This is done in a separate structure as it must exist whenever there is a fluid in any vertex around.
22 Preliminary Tests Some results showing the capabilities of the method are presented. 2D problems run in the 3D code with one layer of cells and plane symmetry. The first is a trivial frictionless sliding of a copper block on steel. It exhibits some difficulties, which were mainly overcome as shown in the following picture The second are two test cases of normal plate penetration: 2m/ms impact of copper on steel 2m/ms impact of aluminum on steel
23 A (not so) trivial test Cu Cu 2m/ms 2mm/ s T=0 T=9.6 Cu Steel Steel Steel Steel
24 A (not so) trivial test of sliding The interface tracking produces some corner rounding. Taking a common normal velocity component there, starts an un-physical interaction. Checking if the moment of the impact is within the time step helps here, but makes the normal impact less smooth. The local interface location is found from the tracking algorithm. When this check is restricted to cells including more than 2 materials (including void) both problems run well. A better interface normal definition for corner cells sharing more than 2 materials will help. Perhaps the interface direction should be determined by the neighboring cells which share only 2 materials
25 2m/ms Normal Impact of Copper or Aluminum on Steel The impact was run both with the standard common vertex velocity and the new multifluid vertex velocity versions. A Steinberg EOS was used with the material parameters for Copper (Cu-OFHC), Steel (SS ) aluminum(al-2024-t4), taken from the AUTODYN library. The 3D code was again run in 2D mode using one layer of cells
26 2m/ms Normal Impact of Copper on Steel T=0 T=5
27 2m/ms Normal Impact of Copper on Steel T=15
28 Comparing T=5, 2mm/ms impact Cu on steel a) Multifluid velocities b) Common vertex velocities
29 Comparing T=7.5, 2mm/ms impact Cu on steel a) Multifluid velocities b) Common vertex velocities
30 Comparing T=10, 2mm/ms impact Cu on steel a) Multifluid velocities b) Common vertex velocities
31 Comparing T=15, 2mm/ms impact Cu on steel a) Multifluid velocities b) Common vertex velocities
32 Comparing T=5, 2mm/ms impact Al on steel
33 Comparing T=7.5, 2mm/ms impact Al on steel a) Multifluid velocities b) Common vertex velocities
34 Comparing T=10, 2mm/ms impact Al on steel a) Multifluid velocities b) Common vertex velocities
35 Comparing T=12.5, 2mm/ms impact Al on steel a) Multifluid velocities b) Common vertex velocities
36 The results of the normal impact tests The early time results are similar. At later stages the multi-fluid velocities allows separation of the interface which is a feature we look for. However, the multifluid-velocity version seems to produce too much fragmentation; with the common vertex velocity we get smoother results.
37 Concluding Remarks We have set-up a multi-fluid velocities version of the SMG scheme. The present multi-fluid velocity version of the code is not yet stable. Further improvements and more debugging is needed A better normal definition for cells sharing more than 2 materials The re-use of some of the features of the old MMALE 2D code should help: Remap type advection near the interfaces Multi-fluid velocities better describe the physics of the interface dynamics: at low speed, where strength and failure are important at higher speed where the growth of instabilities might be relevant
38 References 1. G. Luttwak, R.L. Rabie, The Multimaterial Arbitrary Lagrangian Eulerian code MMALE and its application to some problems of penetration and impact,la-ur ,(1985) 2. Luttwak G., "Numerical Simulation of Jet Formation Using an Adaptive Grid in the Multi-Material ALE Code MMALE", p45-52 in Proc. of the 3d HDP (High dynamic Pressure) Symp., La Grande Motte, France, (1989) 3. J.D. Walker, C.E. Anderson, Multimaterial Velocities for Mixed Cells, p ,High Pressure Science and Technology- 1993,ed. Schmidt et al.,aip,(1994) 4. G. Luttwak, J. Falcovitz, Staggered Mesh Godunov (SMG) Schemes for ALE Hydrodynamics, presented at the Numerical Methods for Multi-material Flows, Oxford, Sept. (2005) 5. E. Vitali, D. Benson, Int. J. Numer. Meth. Engng. 67,p ,(2006)
39 More References Christensen R.B., "Godunov Methods on a Staggered Mesh. An Improved Artificial Viscosity", L.L.N.L report UCRL-JC , (1990). Caramana E.J., Shaskov M.J., Whalen P.P., J. Comp. Phys. 144, p70, (1998). Wilkins M. L., "Calculation of Elastic-Plastic Flow, Meth. Comp. Phys., 3, p211, B. Alder et al. eds., Academic Press (1964). Luttwak, G., "Comparing Lagrangian Godunov and Pseudo- Viscosity Schemes for Multi-Dimensional Impact Simulations, p , Shock Compression of Condensed Matter-2001, Furnish M. D. et al Eds, (2002), AIP, CP620. Luttwak, G. Staggered Mesh Godunov (SMG) Schemes for Lagrangian Hydrodynamics, presented at the APS Topical Conference on Shock Compression of Condensed matter,baltimore,md,july31-aug5,2005
40 More References Luttwak G. "Second Order Discrete Rezoning", in High-Pressure Science and Technology -1993, Schmidt S.C. et al (editors), pp AIP-Press (1994) Luttwak G., Cowler M. S., "Remapping Techniques for Three Dimensional Meshes", presented at the 1999 Int. Workshop for New Models and Hydro-codes for Shock Wave Processes, held at the Univ. of Maryland, College Park, MD,(1999) Luttwak G., "Interface Tracking in Eulerian and MMALE Calculations",p , Shock Compression of Condensed Matter- 2001,Furnish M. D., et al. eds., AIP, CP620, (2002) Luttwak G., Cowler M. S., Birnbaum N., "Virtual Memory Techniques in Eulerian Calculations", p , Shock Compression of Condensed Matter-1999, Furnish M. D., et al. eds.,aip,cp505,(2000)
Vector 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 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 informationComparing Lagrangian Godunov and Pseudo-Viscosity Schemes for Multi-Dimensional Impact Simulations. Gabi Luttwak
Comparing Lagrangian Godunov and Pseudo-Viscosity Schemes for Multi-Dimensional Impact Simulations Gabi Luttwak Rafael, P.O.Box 2250, Haifa 31021,Israel Godunov or Pseudo-viscosity Modern Eulerian codes
More informationLIST OF PUBLICATIONS. Oved Y., Luttwak G., Rosenberg Z., "Shock Wave Propagation in Layered Composites", J. Comp. Mat.
LIST OF PUBLICATIONS Oved Y., Luttwak G., Rosenberg Z., "Shock Wave Propagation in Layered Composites", J. Comp. Mat., 12, p84, (1978) Luttwak G., Kivity Y., Betser A., "Effects of a Hypervelocity Jet
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 informationVIP (Vector Image Polygon) multi-dimensional slope limiters for scalar variables
VIP (Vector Image Polygon) multi-dimensional slope limiters for scalar variables Gabi Luttwak a, Joseph Falcovitz b a Rafael, Box 2250, Haifa 31021, Israel b Institute of Mathematics, The Hebrew University
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 informationInterface Tracking in Eulerian and MMALE Calculations
Interface Tracking in Eulerian and MMALE Calculation Gabi Luttwak Rafael P.O.Box 2250, Haifa 31021,Irael Interface Tracking in Eulerian and MMALE Calculation 3D Volume of Fluid (VOF) baed recontruction
More informationFluid-Structure Interaction in LS-DYNA: Industrial Applications
4 th European LS-DYNA Users Conference Aerospace / Fluid-Struct. Inter. Fluid-Structure Interaction in LS-DYNA: Industrial Applications M hamed Souli Universite des Sciences et Technologie de Lille Laboratoire
More informationALE Adaptive Mesh Refinement in LS-DYNA
12 th International LS-DYNA Users Conference FSI/ALE(2) ALE Adaptive Mesh Refinement in LS-DYNA Nicolas AQUELET Livermore Software Technology Corp. 7374 Las Positas Rd Livermore CA94550 aquelet@lstc.com
More informationSPH: Why and what for?
SPH: Why and what for? 4 th SPHERIC training day David Le Touzé, Fluid Mechanics Laboratory, Ecole Centrale de Nantes / CNRS SPH What for and why? How it works? Why not for everything? Duality of SPH SPH
More informationBest Practices Workshop: Overset Meshing
Best Practices Workshop: Overset Meshing Overview Introduction to Overset Meshes Range of Application Workflow Demonstrations and Best Practices What are Overset Meshes? Overset meshes are also known as
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 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 informationALE METHODS FOR DETERMINING STATIONARY SOLUTIONS OF METAL FORMING PROCESSES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 c ECCOMAS ALE METHODS FOR DETERMINING STATIONARY SOLUTIONS OF METAL FORMING PROCESSES
More informationSIMULATION OF A DETONATION CHAMBER TEST CASE
SIMULATION OF A DETONATION CHAMBER TEST CASE Daniel Hilding Engineering Research Nordic AB Garnisonen I4, Byggnad 5 SE-582 10 Linköping www.erab.se daniel.hilding@erab.se Abstract The purpose of a detonation
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 informationJ. Vira, M. Sofiev SILAM winter school, February 2013, FMI
Numerical aspects of the advection-diffusion equation J. Vira, M. Sofiev SILAM winter school, February 2013, FMI Outline Intro Some common requirements for numerical transport schemes Lagrangian approach
More informationCTH: A Software Family for Multi-Dimensional Shock Physics Analysis E. S. Hertel, Jr., R. L. Bell, M. G. Elrick, A. V. Farnsworth, G. I. Kerley, J.
CTH: A Software Family for Multi-Dimensional Shock Physics Analysis E. S. Hertel, Jr., R. L. Bell, M. G. Elrick, A. V. Farnsworth, G. I. Kerley, J. M. McGlaun, S. V. Petney, S. A. Silling, P. A. Taylor,
More informationModelling of Impact on a Fuel Tank Using Smoothed Particle Hydrodynamics
Modelling of Impact on a Fuel Tank Using Smoothed Particle Hydrodynamics R. Vignjevic a, T. De Vuyst a, J. Campbell a, N. Bourne b School of Engineering, Cranfield University, Bedfordshire, MK43 0AL, UK.
More informationA cell by cell anisotropic adaptive mesh Arbitrary Lagrangian Eulerian method for the numerical solution of the Euler equations. J. M.
THE UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS A cell by cell anisotropic adaptive mesh Arbitrary Lagrangian Eulerian method for the numerical solution of the Euler equations. J. M. Morrell Thesis
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 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 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 informationCloth Hair. and. soft bodies
Cloth Hair Lesson 11 and soft bodies Lesson 08 Outline Problem definition and motivations Modeling deformable solids with mass-spring model Position based dynamics Modeling cloths with mass-spring model
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 informationHigh performance computing using AUTODYN-3D
High performance computing using AUTODYN-3D M. S. Cowler', O. La'adan\ T. Ohta^ ' Century Dynamics Incorporated, USA. Hebrew University ofjerusalem, Israel. * CRC Research Institute Incorporated, Japan.
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 informationSupport for Multi physics in Chrono
Support for Multi physics in Chrono The Story Ahead Overview of multi physics strategy in Chrono Summary of handling rigid/flexible body dynamics using Lagrangian approach Summary of handling fluid, and
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 information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationFluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching
Fluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching Paul Groenenboom ESI Group Delft, Netherlands Martin Siemann German Aerospace Center (DLR) Stuttgart, Germany
More informationALE and Fluid-Structure Interaction in LS-DYNA March 2004
ALE and Fluid-Structure Interaction in LS-DYNA March 2004 Workshop Models 1. Taylor bar impact 2. One-dimensional advection test 3. Channel 4. Underwater explosion 5. Bar impacting water surface 6. Sloshing
More informationShape of Things to Come: Next-Gen Physics Deep Dive
Shape of Things to Come: Next-Gen Physics Deep Dive Jean Pierre Bordes NVIDIA Corporation Free PhysX on CUDA PhysX by NVIDIA since March 2008 PhysX on CUDA available: August 2008 GPU PhysX in Games Physical
More informationAcknowledgements. Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn. SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar
Philipp Hahn Acknowledgements Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar 2 Outline Motivation Lumped Mass Model Model properties Simulation
More informationANSYS/LS-Dyna. Modeling. Prepared by. M Senior Engineering Manager
ANSYS/LS-Dyna Multi-Material M t i l ALE Modeling Prepared by Steven Hale, M.S.M.E. M Senior Engineering Manager Multi-Material ALE Basics What is the Multi-Material ALE method (MMALE) MMALE stands for
More informationEssay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions
Table of Contents Essay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions and Scaling Analysis 1 About dimensional analysis 4 1.1 Thegoalandtheplan... 4 1.2 Aboutthisessay... 5 2 Models
More informationOVERVIEW OF ALE METHOD IN LS-DYNA. Ian Do, LSTC Jim Day, LSTC
1 OVERVIEW OF ALE METHOD IN LS-DYNA Ian Do, LSTC Jim Day, LSTC ELEMENT FORMULATIONS REVIEW A physical process involving fluids or fluid-like behavior may often be modeled in more than one way, using more
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 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 informationA METHOD TO INCREASE THE TIP VELOCITY OF A SHAPED CHARGE JET. William P. Walters and Daniel R. Scheffler
23 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16-20 APRIL 2007 A METHOD TO INCREASE THE TIP VELOCITY OF A SHAPED CHARGE JET William P. Walters and Daniel R. Scheffler US Army Research Laboratory,
More informationNUMERICAL SIMULATION OF STRUCTURAL DEFORMATION UNDER SHOCK AND IMPACT LOADS USING A COUPLED MULTI-SOLVER APPROACH
NUMERICAL SIMULATION OF STRUCTURAL DEFORMATION UNDER SHOCK AND IMPACT LOADS USING A COUPLED MULTI-SOLVER APPROACH X. Quan, N. K. Birnbaum, M.S. Cowler, B. I. Gerber Century Dynamics, Inc., 1001 Galaxy
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 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 informationTHE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS
March 18-20, 2013 THE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS Authors: M.R. Chiarelli, M. Ciabattari, M. Cagnoni, G. Lombardi Speaker:
More informationFINITE ELEMENT MODELING OF THICK PLATE PENETRATIONS IN SMALL CAL MUNITIONS
FINITE ELEMENT MODELING OF THICK PLATE PENETRATIONS IN SMALL CAL MUNITIONS Raymond Chaplin RDAR-MEM-I Picatinny Arsenal, NJ raymond.c.chaplin@us.army.mil 973-724-8562 Why Finite Element Modeling? Reduced
More informationGALAXY ADVANCED ENGINEERING, INC. P.O. BOX 614 BURLINGAME, CALIFORNIA Tel: (650) Fax: (650)
GALAXY ADVANCED ENGINEERING, INC. P.O. BOX 614 BURLINGAME, CALIFORNIA 94011 Tel: (650) 740-3244 Fax: (650) 347-4234 E-mail: bahmanz@aol.com PUFF-TFT/PC A Material Response Computer Code for PC Computer
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 informationRealtime Water Simulation on GPU. Nuttapong Chentanez NVIDIA Research
1 Realtime Water Simulation on GPU Nuttapong Chentanez NVIDIA Research 2 3 Overview Approaches to realtime water simulation Hybrid shallow water solver + particles Hybrid 3D tall cell water solver + particles
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
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 informationRealistic Animation of Fluids
1 Realistic Animation of Fluids Nick Foster and Dimitris Metaxas Presented by Alex Liberman April 19, 2005 2 Previous Work Used non physics-based methods (mostly in 2D) Hard to simulate effects that rely
More informationRotating Moving Boundary Analysis Using ANSYS 5.7
Abstract Rotating Moving Boundary Analysis Using ANSYS 5.7 Qin Yin Fan CYBERNET SYSTEMS CO., LTD. Rich Lange ANSYS Inc. As subroutines in commercial software, APDL (ANSYS Parametric Design Language) provides
More informationTable of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.
Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Preface page 2 Part I Fluid Dynamics and Waves 7 1 Elements of fluid dynamics 9 1.1
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 informationARBITRARY LAGRANGIAN-EULERIAN (ALE) METHODS IN COMPRESSIBLE FLUID DYNAMICS. Milan Kuchařík, Richard Liska, Pavel Váchal, Mikhail Shashkov
ARBITRARY LAGRANGIAN-EULERIAN (ALE) METHODS IN COMPRESSIBLE FLUID DYNAMICS Milan Kuchařík, Richard Liska, Pavel Váchal, Mikhail Shashkov Abstract The aim of this paper is to present an Arbitrary Lagrangian-Eulerian
More informationUniversity of Quebec at Chicoutimi (UQAC) Friction Stir Welding Research Group. A CUDA Fortran. Story
University of Quebec at Chicoutimi (UQAC) Friction Stir Welding Research Group A CUDA Fortran Story Summary At the University of Quebec at Chicoutimi (UQAC) friction stir welding (FSW) research group,
More informationTHREE-DIMENSIONAL ANALYSIS OF THE EXPLOSIVE INITIATION THRESHOLD FOR SIDE IMPACT ON A SHAPED CHARGE WARHEAD
Insensitive Munitions and Energetics Technology Symposium Tampa, Florida, USA, October 1997 THREE-DIMENSIONAL ANALYSIS OF THE EXPLOSIVE INITIATION THRESHOLD FOR SIDE IMPACT ON A SHAPED CHARGE WARHEAD David
More informationAPPLICATIONS OF NUMERICAL ANALYSIS IN SIMULATION ENGINEERING AND MANUFACTURING TECHNOLOGIES
APPLICATIONS OF NUMERICAL ANALYSIS IN SIMULATION ENGINEERING AND MANUFACTURING TECHNOLOGIES Haidar Amer 1 1 Stefan Cel Mare University of Suceava, amerhaidar85@gmail.com Abstract: This article summarizes
More information3D simulations of concrete penetration using SPH formulation and the RHT material model
3D simulations of concrete penetration using SPH formulation and the RHT material model H. Hansson Weapons and Protection, Swedish Defence Research Agency (FOI), Sweden Abstract This paper describes work
More informationFluid Simulation. [Thürey 10] [Pfaff 10] [Chentanez 11]
Fluid Simulation [Thürey 10] [Pfaff 10] [Chentanez 11] 1 Computational Fluid Dynamics 3 Graphics Why don t we just take existing models from CFD for Computer Graphics applications? 4 Graphics Why don t
More informationThe Immersed Smoothed Particle Galerkin Method in LS-DYNA for Material Failure Analysis of Fiber-Reinforced Solid Structures
The Immersed Smoothed Particle Galerkin Method in LS-DYNA for Material Failure Analysis of Fiber-Reinforced Solid Structures Wei Hu, C. T Wu Livermore Software Technology Corporation, 7374 Las Positas
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 informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationTechnical Report TR
Technical Report TR-2015-09 Boundary condition enforcing methods for smoothed particle hydrodynamics Arman Pazouki 1, Baofang Song 2, Dan Negrut 1 1 University of Wisconsin-Madison, Madison, WI, 53706-1572,
More informationAurélien Thinat Stéphane Cordier 1, François Cany
SimHydro 2012:New trends in simulation - Hydroinformatics and 3D modeling, 12-14 September 2012, Nice Aurélien Thinat, Stéphane Cordier, François Cany Application of OpenFOAM to the study of wave loads
More informationInteraction of Fluid Simulation Based on PhysX Physics Engine. Huibai Wang, Jianfei Wan, Fengquan Zhang
4th International Conference on Sensors, Measurement and Intelligent Materials (ICSMIM 2015) Interaction of Fluid Simulation Based on PhysX Physics Engine Huibai Wang, Jianfei Wan, Fengquan Zhang College
More informationAnalysis of Fluid-Structure Interaction Effects of Liquid-Filled Container under Drop Testing
Kasetsart J. (Nat. Sci.) 42 : 165-176 (2008) Analysis of Fluid-Structure Interaction Effects of Liquid-Filled Container under Drop Testing Chakrit Suvanjumrat*, Tumrong Puttapitukporn and Satjarthip Thusneyapan
More informationA meshfree weak-strong form method
A meshfree weak-strong form method G. R. & Y. T. GU' 'centre for Advanced Computations in Engineering Science (ACES) Dept. of Mechanical Engineering, National University of Singapore 2~~~ Fellow, Singapore-MIT
More informationComputational Fluid Dynamic Hydraulic Characterization: G3 Cube vs. Dolos Armour Unit. IS le Roux, WJS van der Merwe & CL de Wet
Computational Fluid Dynamic Hydraulic Characterization: G3 Cube vs. Dolos Armour Unit IS le Roux, WJS van der Merwe & CL de Wet Presentation Outline Scope. Assumptions and boundary values. Numerical mesh.
More informationSimulating Sinkage & Trim for Planing Boat Hulls. A Fluent Dynamic Mesh 6DOF Tutorial
Simulating Sinkage & Trim for Planing Boat Hulls A Fluent Dynamic Mesh 6DOF Tutorial 1 Introduction Workshop Description This workshop describes how to perform a transient 2DOF simulation of a planing
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 informationComputing & Verifying Compressible Fluid Dynamics:! The Good, The Bad and The Ugly!
: LA- UR 11-05852 Computing & Verifying Compressible Fluid Dynamics:! The Good, The Bad and The Ugly! Tariq Aslam! Los Alamos National Laboratory! WX-9: Shock and Detonation Physics! Slide 1 Background:!
More informationAREPO: a moving-mesh code for cosmological hydrodynamical simulations
AREPO: a moving-mesh code for cosmological hydrodynamical simulations E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh Springel, 2010 arxiv:0901.4107 Rubens
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 informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationSIMPLIFIED VOLUME OF FLUID METHOD (SVOF) FOR TWO-PHASE FLOWS
TASK QUARTERLY 12 No 3, 255 265 SIMPLIFIED VOLUME OF FLUID METHOD (SVOF) FOR TWO-PHASE FLOWS MACIEJ MAREK, WOJCIECH ANISZEWSKI AND ANDRZEJ BOGUSŁAWSKI Institute of Thermal Machinery, Czestochowa University
More informationSimulation of liquid cube fracture with SPH
Journal of Physics: Conference Series PAPER OPEN ACCESS Simulation of liquid cube fracture with SPH To cite this article: M N Davydov 2016 J. Phys.: Conf. Ser. 754 062001 View the article online for updates
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationCHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING
113 CHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING 6.1 INTRODUCTION Superplastic properties are exhibited only under a narrow range of strain rates. Hence, it
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationHigh-Order Finite Difference Schemes for computational MHD
High-Order 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 of Turin,,
More informationNumerical Simulations of Granular Materials Flow around Obstacles: The role of the interstitial gas
Numerical Simulations of Granular Materials Flow around Obstacles: The role of the interstitial gas Avi Levy, Dept. Mech. Eng., Ben Gurion University, Beer Sheva, Israel. Mohamed Sayed, CHC, National Research
More informationEffect of Soil Material Models on SPH Simulations for Soil- Structure Interaction
Effect of Soil Material Models on SPH Simulations for Soil- Structure Interaction Ronald F. Kulak 1 Len Schwer 2 1) RFK Engineering Mechanics Consultants 307 Warwick Drive Naperville, IL 60565, USA rfkemc@aol.com
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationLagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion
Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution Woods Hole, MA, 02543 July 31, 2006 Summary: This essay
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 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 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 informationSPH: Towards the simulation of wave-body interactions in extreme seas
SPH: Towards the simulation of wave-body interactions in extreme seas Guillaume Oger, Mathieu Doring, Bertrand Alessandrini, and Pierre Ferrant Fluid Mechanics Laboratory (CNRS UMR6598) Ecole Centrale
More informationSimulation of Fuel Sloshing Comparative Study
3. LS-DYNA Anwenderforum, Bamberg 2004 Netfreie Verfahren Simulation of Fuel Sloshing Comparative Study Matej Vesenjak 1, Heiner Müllerschön 2, Alexander Hummel 3, Zoran Ren 1 1 University of Maribor,
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 informationNavier-Stokes & Flow Simulation
Last Time? Navier-Stokes & Flow Simulation Optional Reading for Last Time: Spring-Mass Systems Numerical Integration (Euler, Midpoint, Runge-Kutta) Modeling string, hair, & cloth HW2: Cloth & Fluid Simulation
More informationParticle-based Fluid Simulation
Simulation in Computer Graphics Particle-based Fluid Simulation Matthias Teschner Computer Science Department University of Freiburg Application (with Pixar) 10 million fluid + 4 million rigid particles,
More information2-D Tank Sloshing Using the Coupled Eulerian- LaGrangian (CEL) Capability of Abaqus/Explicit
2-D Tank Sloshing Using the Coupled Eulerian- LaGrangian (CEL) Capability of Abaqus/Explicit Jeff D. Tippmann, Sharat C. Prasad 2, and Parthiv N. Shah ATA Engineering, Inc. San Diego, CA 923 2 Dassault
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 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 informationCode_Aster. SSNV209 Interface in contact rubbing with X-FEM
Titre : SSNV209 - Interface en contact frottant avec X-FEM Date : 21/07/2015 Page : 1/34 SSNV209 Interface in contact rubbing with X-FEM Summary: This problem corresponds to a quasi-static analysis of
More informationNumerical Methods. (Additional Notes from Talks by PFH)
Numerical Methods (Additional Notes from Talks by PFH) SPH Challenge: POPULAR METHODS FOR HYDRODYNAMICS HAVE PROBLEMS Lucy 77, Gingold & Monaghan 77 Reviews by: Springel 11, Price 12 Smoothed-Particle
More information