BoundPreserving Reconstruction of Tensor Quantities for Remap in ALE Fluid Dynamics
Abstract
We analyze several new and existing approaches for limiting tensor quantities in the context of deviatoric stress remapping in an ALE numerical simulation of elastic flow. Remapping and limiting of the tensor componentbycomponent is shown to violate radial symmetry of derived variables such as elastic energy or force. Therefore, we have extended the symmetrypreserving Vector Image Polygon algorithm, originally designed for limiting vector variables. This limiter constrains the vector (in our case a vector of independent tensor components) within the convex hull formed by the vectors from surrounding cells – an equivalent of the discrete maximum principle in scalar variables. We compare this method with a limiter designed specifically for deviatoric stress limiting which aims to constrain the J _{2} invariant that is proportional to the specific elastic energy and scale the tensor accordingly. We also propose a method which involves remapping and limiting the J _{2} invariant independently using known scalar techniques. The deviatoric stress tensor is then scaled to match this remapped invariant, which guarantees conservation in terms of elastic energy.
 Authors:
 Czech Technical Univ. in Prague, Praha (Czech Republic)
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Publication Date:
 Research Org.:
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Sponsoring Org.:
 USDOE Office of Science (SC). Advanced Scientific Computing Research (ASCR) (SC21)
 OSTI Identifier:
 1338785
 Report Number(s):
 LAUR1720068
 DOE Contract Number:
 AC5206NA25396
 Resource Type:
 Technical Report
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS
Citation Formats
Klima, Matej, Kucharik, MIlan, Shashkov, Mikhail Jurievich, and Velechovsky, Jan. BoundPreserving Reconstruction of Tensor Quantities for Remap in ALE Fluid Dynamics. United States: N. p., 2017.
Web. doi:10.2172/1338785.
Klima, Matej, Kucharik, MIlan, Shashkov, Mikhail Jurievich, & Velechovsky, Jan. BoundPreserving Reconstruction of Tensor Quantities for Remap in ALE Fluid Dynamics. United States. doi:10.2172/1338785.
Klima, Matej, Kucharik, MIlan, Shashkov, Mikhail Jurievich, and Velechovsky, Jan. Fri .
"BoundPreserving Reconstruction of Tensor Quantities for Remap in ALE Fluid Dynamics". United States.
doi:10.2172/1338785. https://www.osti.gov/servlets/purl/1338785.
@article{osti_1338785,
title = {BoundPreserving Reconstruction of Tensor Quantities for Remap in ALE Fluid Dynamics},
author = {Klima, Matej and Kucharik, MIlan and Shashkov, Mikhail Jurievich and Velechovsky, Jan},
abstractNote = {We analyze several new and existing approaches for limiting tensor quantities in the context of deviatoric stress remapping in an ALE numerical simulation of elastic flow. Remapping and limiting of the tensor componentbycomponent is shown to violate radial symmetry of derived variables such as elastic energy or force. Therefore, we have extended the symmetrypreserving Vector Image Polygon algorithm, originally designed for limiting vector variables. This limiter constrains the vector (in our case a vector of independent tensor components) within the convex hull formed by the vectors from surrounding cells – an equivalent of the discrete maximum principle in scalar variables. We compare this method with a limiter designed specifically for deviatoric stress limiting which aims to constrain the J2 invariant that is proportional to the specific elastic energy and scale the tensor accordingly. We also propose a method which involves remapping and limiting the J2 invariant independently using known scalar techniques. The deviatoric stress tensor is then scaled to match this remapped invariant, which guarantees conservation in terms of elastic energy.},
doi = {10.2172/1338785},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Fri Jan 06 00:00:00 EST 2017},
month = {Fri Jan 06 00:00:00 EST 2017}
}

Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme
From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian–Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense thatmore » 
An InvariantPreserving ALE Method for Solids under Extreme Conditions
We are proposing a fundamentally new approach to ALE methods for solids undergoing large deformation due to extreme loading conditions. Our approach is based on a physicallymotivated and mathematically rigorous construction of the underlying Lagrangian method, vector/tensor reconstruction, remapping, and interface reconstruction. It is transformational because it deviates dramatically from traditionally accepted ALE methods and provides the following set of unique attributes: (1) a threedimensional, finite volume, cellcentered ALE framework with advanced hypo/hyperelastoplastic constitutive theories for solids; (2) a new physically and mathematically consistent reconstruction method for vector/tensor fields; (3) advanced invariantpreserving remapping algorithm for vector/tensor quantities; (4) momentoffluid (MoF)more » 
SALE3D: a simplified ALE computer program for calculating threedimensional fluid flow
This report presents a simplified numerical fluiddynamics computing technique for calculating timedependent flows in three dimensions. An implicit treatment of the pressure equation permits calculation of flows far subsonic without stringent constraints on the time step. In addition, the grid vertices may be moved with the fluid in Lagrangian fashion or held fixed in an Eulerian manner, or moved in some prescribed manner to give a continuous rezoning capability. This report describes the combination of Implicit Continuousfluid Eulerian (ICE) and Arbitrary LagrangianEulerian (ALE) to form the ICEdALE technique in the framework of the SimplifiedALE (SALE3D) computer program, for which amore » 
SALE3D. ICEdALE Treatment of 3D Fluid Flow
SALE3D calculates threedimensional fluid flow at all speeds, from the incompressible limit to highly supersonic. An implicit treatment of the pressure calculation similar to that in the Implicit Continuousfluid Eulerian (ICE) technique provides this flow speed flexibility. In addition, the computing mesh may move with the fluid in a typical Lagrangian fashion, be held in an Eulerian manner, or move in some arbitrarily specified way to provide a continuous rezoning capability. This latitude results from use of an Arbitrary LagrangianEulerian (ALE) treatment of the mesh. The partial differential equations solved are the NavierStokes equations and the mass and internal energymore »