[en] A virtual test facility (VTF) for studying the three-dimensional dynamic response of solid materials subject to strong shock and detonation waves has been constructed as part of the research program of the Center for Simulating the Dynamic Response of Materials at the California Institute of Technology. The compressible fluid flow is simulated with a Cartesian finite volume method and treating the solid as an embedded moving body, while a Lagrangian finite element scheme is employed to describe the structural response to the hydrodynamic pressure loading. A temporal splitting method is applied to update the position and velocity of the boundary between time steps. The boundary is represented implicitly in the fluid solver with a level set function that is constructed on-the-fly from the unstructured solid surface mesh. Block-structured mesh adaptation with time step refinement in the fluid allows for the efficient consideration of disparate fluid and solid time scales. We detail the design of the employed object-oriented mesh refinement framework AMROC and outline its effective extension for fluid-structure interaction problems. Further, we describe the parallelization of the most important algorithmic components for distributed memory machines and discuss the applied partitioning strategies. As computational examples for typical VTF applications, we present the dynamic deformation of a tantalum cylinder due to the detonation of an interior solid explosive and the impact of an explosion-induced shock wave on a multi-material soft tissue body.
Disciplines :
Mechanical engineering
Author, co-author :
Deiterding, Ralf; California Institut of Technology - CALTECH
Radovitzky, Raúl; Massachusetts Institute of Technology - MIT > Aeronautics & Astronautics
Mauch, Sean; California Institut of Technology - CALTECH
Noels, Ludovic ; Université de Liège - ULiège > Département d'aérospatiale et mécanique > LTAS - Milieux continus et thermomécanique
Cummings, Julian; California Institut of Technology - CALTECH
Meiron, Dan; California Institut of Technology - CALTECH
Language :
English
Title :
A virtual test facility for the efficient simulation of solid material response under strong shock and detonation wave loading
Aivazis M, Goddard W, Meiron D, Ortiz M, Pool J, Shepherd J (2000) A virtual test facility for simulating the dynamic response of materials. Comput Sci Eng 2(2):42-53
Cummings J, Aivazis M, Samtaney R, Radovitzky R, Mauch S, Meiron D (2002) A virtual test facility for the simulation of dynamic response in materials. J Supercomput 23:39-50
Mauch S, Meiron D, Radovitzky R, Samtaney R (2003) Coupled Eulerian-Lagrangian simulations using a level set method. In: Bathe K (ed) 2nd M.I.T. conference on computational fluid and solid mechanics, Cambridge, MA, June 17-20
Arienti M, Hung P, Morano E, Shepherd J (2003) A level set approach to Eulerian-Lagrangian coupling. J Comput Phys 185:213-251
Cirak F, Radovitzky R (2003) A general algorithm for coupling Lagrangian-shell with Eulerian-fluid formulations. In: Proceedings of the IUTAM symposium on integrated modeling of fully coupled fluid-structure interactions using analysis, computations and experiments, New Brunswick, NJ, June 1-6, 2003. International Union of Theoretical and Applied Mechanics
Löhner R, Baum J, Charman C, Pelessone D (2003) Fluid-structure interaction simulations using parallel computers. In: Lecture notes in computer science, vol 2565. Springer, Berlin Heidelberg New York, pp 3-23
Cirak F, Radovitzky R (2005) A Lagrangian-Eulerian shell-fluid coupling algorithm based on level sets. Comput Struct 83:491-498
Löhner R, Cebral J, Yang C, Baum J, Mestreau E, Charman C, Pelessone D, (2004) Large-scale fluid-structure interaction simulations. Comput Sci Eng 6(3):27-37
Hu H, Patankar N, Zhu M (2001) Direct numerical simulations of fluid-solid systems using the arbitrary Lagrangian-Eulerian technique. J Comput Phys 169(2):427-462
Mittal R, Iaccarino G (2005) Immersed boundary methods. Annu Rev Fluid Mech 37:239-261
Quirk J (1994) An alternative to unstructured grids for computing gas dynamics flows around arbitrarily complex two-dimensional bodies. Comput Fluids 23:125-142
Berger M, Helzel C (2002) Grid aligned H-box methods for conservation laws in complex geometries. In: Proceedings of the 3rd international symposium on finite volumes for complex applications, Porquerolles
Fedkiw R (2002) Coupling an Eulerian fluid calculation to a Lagrangian solid calculation with the ghost fluid method. J Comput Phys 175:200-224
Tam D, Radovitzky R, Samtaney R (2005) An algorithm for modelling the interaction of a flexible rod with a two-dimensional high-speed flow. Int J Numer Methods Eng 64(8):1057-1077
Pantano C, Deiterding R, Hill D, Pullin D (2005) A low-numerical dissipation patch-based adaptive mesh refinement method for large-eddy simulation of compressible flows. In: Proceedings of the Cyprus international symposium on complex effects in large eddy simulations, University of Cyprus, Nicosia
Deiterding R (2005a) Detonation structure simulation with AMROC. In: Lecture notes in computer science, vol 3726. Springer, Berlin Heidelberg New York, pp 916-927
Deiterding R, Bader G (2005) High-resolution simulation of detonations detailed chemistry. In: Warnecke G (ed) Analysis and numerics of conservation laws. Springer, Berlin Heidelberg New York, pp 69-91
Deiterding R (2005b) A high-resolution method for realistic detonation structure simulation. In: Asakura F et al (eds) Proceegins of the 10th international conference on hyperbolic problems: theory, numerics, applications, Yokohama Publishers
Deiterding R (2003) Parallel adaptive simulation of multidimensional detonation structures. PhD Thesis, Brandenburgische Technische Universität Cottbus
Deiterding R (2002) Efficient simulation of multi-dimensional detonation phenomena. In: Handlikova A et al (eds) Proceedings of the ALGORITMY 2002, 16th conference on scientific computing. Slovak University of Technology, Bratislava, pp 94-101
Fickett W, Davis W (1979) Detonation. University of California Press, Berkeley, Los Angeles
Clarke J, Karni S, Quirk J, Roe P, Simmonds L, Toro E (1993) Numerical computation of two-dimensional unsteady detonation waves in high energy solids. J Comput Phys 106:215-233
Toro E (1999) Riemann solvers and numerical methods for fluid dynamics. Springer, Berlin Heidelberg New York
Einfeldt B, Munz C, Roe P, Sjögreen B (1991) On Godunov-type methods near low densities. J Comput Phys 92:273-295
Fedkiw R, Aslam T, Merriman B, Osher S (1999) A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method). J Comput Phys 152:457-492
Osher S, Fedkiw R (2003) Level set methods and dynamic implicit surfaces. Springer, Berlin Heidelberg New York
Berger M, Colella P (1988) Local adaptive mesh refinement for shock hydrodynamics. J Comput Phys 82:64-84
Bell J, Berger M, Saltzman J, Welcome M (1994) Three-dimensional adaptive mesh refinement for hyperbolic conservation laws. SIAM J Sci Comp 15(1):127-138
Rendleman C, Beckner V, Lijewski M, Crutchfield W, Bell J (2000) Parallelization of structured, hierarchical adaptive mesh refinement algorithms. Comput Vis Sci 3:147-157
Parashar M, Browne J (1997) System engineering for high performance computing software: the HDDA/DAGH infrastructure for implementation of parallel structured adaptive mesh refinement. In: Structured adaptive mesh refinement grid methods, IMA volumes in mathematics and its applications. Springer, Berlin Heidelberg New York
Deiterding R (2005c) AMROC-blockstructured adaptive mesh refinement in object-oriented C++. http://www.am-roc.sourceforge.net
Deiterding R (2005d) Construction and application of an AMR algorithm for distributed memory computers. In: Plewa T et al (eds) Lecture notes in computational science and engineering, vol 41. Springer, Berlin Heidelberg New York, pp 361-372
Parashar M, Browne J (1996) On partitioning dynamic adaptive grid hierarchies. In: Proceedings of the 29th annual Hawaii international conference on system sciences
MacNeice P, Olson K, Mobarry C, deFainchtein R, Packer C (2000) PARAMESH: a parallel adaptive mesh refinement community toolkit. Comput Phys Commun 126:330-354
Booch G, Rumbaugh J, Jacobsen I (1999) The unified modeling language user guide. Addison-Wesley, Reading
Marsden J, Hughes T (1993) Mathematical foundations of elasticity. Dover Publications, New York
Lew A, Radovitzky R, Ortiz M (2002) An artificial-viscosity method for the Lagrangian analysis of shocks in solids with strength on unstructured, arbitrary-order tetrahedral meshes. J Comput-Aided Mater Des 8:213-231
Von Neumann J, Richtmyer R (1950) A method for the numerical calculation of hydrodynamic shocks. J Appl Phys 21:232-243
Lubliner J (1972) On the thermodynamic foundations of non-linear solid mechanics. Int J Non-Linear Mech 7:237-254
Radovitzky R, Ortiz M (1999) Error estimation and adaptive meshing in strongly nonlinear dynamic problems. Comput Methods Appl Mech Eng 172:203-240
Ortiz M, Stainier L (1999) The variational formulation of viscoplastic constitutive updates. Comput Methods Appl Mech Eng 171:419-444
Cuitino A, Ortiz M (1992) Material-independent method for extending stress update algorithms from small-strain plasticity to finite plasticity with multiplicative kinematics. Eng Comput 9(4):437-451
Stainier L, Cuitino A, Ortiz M (2002) A micromechanical model of hardening, rate sensitivity and thermal softening in BCC single crystals. J Mech Phys Solids 50:1511-1545
Kambouchev N, Fernandez J, Radovitzky R (2005) Poly-convex model for materials with cubic anisotropy. http://www.arxiv.org/abs/cond-mat/0505178
Benson D (1992) Computational methods in Lagrangian and Eulerian hydrocodes. Comput Methods Appl Mech Eng 99:235-394
Ciarlet P (1976) Numerical analysis of the finite element method. Les Presses de L'Universite de Montreal, Quebec
Belytschko T (1983) An overview of semidiscretization and time integration procedures. In: Belytschko T, Hughes T (eds) Computational methods for transient analysis. North-Holland, Amsterdam, pp 1-65
Hughes T (1983) Analysis of transient algorithms with particular reference to stability behavior. In: Belytschko T, Hughes T (eds) Computational methods for transient analysis. North-Holland, Amsterdam, pp 67-155
Kane C, Marsden J, Ortiz M, West M (2000) Variational integrators and the Newmark algorithm for conservative and dissipative mechanical systems. Int J Numer Methods Eng 49:1295-1325
Deiterding R, Cirak F, Mauch S, Meiron D (2006) A virtual test facility for simulating detonation-induced fracture of thin flexible shells. In: Alexandrov V et al (eds) Lecture notes in computer science, vol 3992. Springer, Berlin Heidelberg New York, pp 122-130
Karypis G (2001) METIS-family of multilevel partitioning algorithms. http://www.users.cs.umn.edu/~karypis/metis
Liu A, Joe B (1996) Quality local refinement of tetrahedral meshes based on 8-subtetrahedron subdivision. Math Comput 65:1183-1200
Stainier L, Radovitzky R, Ortiz M (1999) Constitutive modeling of Tantalum crystals and polycrystals. In: Croitoro E (ed) Proceedings of the 1st Canadian conference on nonlinear solid mechanics, Victoria, BC, pp 203-213
Radovitzky R, Cuitino A (2003) Direct numerical simulation of polycrystals. In: Proceedings of the 44th AIAA/ASME/ASCE/AHS structures, structural dynamics, and materials conference, Norfolk, AIAA
Zhao Z, Radovitzky R, Kuchnicki S, Cuitino A (2003) Dynamical evolution of texture in polycrystalline metals under high strain-rate conditions. In: Proceedings of the international symposium on plasticity and its current applications, Quebec City
Cuitino A, Radovitzky R (eds) (2004) Multiscale material modeling and simulation. In: Sandia NM Proceedings of the symposia on multiscale material modeling and simulation, Cambridge, MA. Modelling Simul Mater Sci Eng, vol 12, issue 4, Institute of Physics Publishing
Zhao Z, Radovitzky R, Cuitino A (2004) A study of surface roughening in fcc metals using direct numerical simulation. Acta Mater 52:5791-5804
Kuchnicki S, Cuitino A, Radovitzky R (2005) Efficient and robust constitutive integrators for single-crystal plasticity modeling. Int J Plast 22:1988-2011
Sethian J (1999) Level set methods and fast marching methods. Cambridge University Press, Cambridge
Johnson D, Cohen E (1998) In: Proceedings of the IEEE international conference on robotics and automation, Leuven, pp 3678-3684
Mauch S (2003) Efficient algorithms for solving static Hamilton-Jacobi equations. PhD Thesis, California Institute of Technology
Foley J, van Dam A, Feiner S, Hughes J (1996) Computer graphics: principles and practice. Addison-Wesley, Reading
Watt A (1993) 3D computer graphics. Addison-Wesley, Reading
Park K, Felippa C, Deruntz J (1997) Stabilization of staggered solution procedures for fluid-structure interaction analysis. In: Belytschko T, Geers T (eds) Computational methods for fluid-structure interaction problems, New York, pp 94-124
Zhang Q, Hisada T (2004) Studies of the strong coupling and weak coupling methods in FSI analysis. Int J Numer Methods Eng 60:2013-2029
Taylor G (1950) The formation of a blast wave by a very intense explosion, I theorical discussion. Proc R Soc Ser A 201(1065):159-174
Brode H (1955) Numerical solution of spherical blast waves. J Appl Phys 26(6):766-775
Cohen R, Gülseren O (2001) Thermal equation of state of tantalum. Phys Rev B 63:3363
Steinberg D, Cochran S, Guinan M (1980) A constitutive model for metals applicale at high-strain rate. J Appl Phys 51:1498-1504
Mader C (1979) Numerical modeling of detonations. University of California Press, Berkeley, Los Angeles
