A subzone reconstruction algorithm for efficient staggered compatible remapping

Staggered-grid Lagrangian hydrodynamics algorithms frequently make use of subzonal discretization of state variables for the purposes of improved numerical accuracy, generality to unstructured meshes, and exact conservation of mass, momentum, and energy. For Arbitrary Lagrangian-Eulerian (ALE) methods using a geometric overlay, it is difficult to remap subzonal variables in an accurate and efficient manner due to the number of subzone-subzone intersections that must be computed. This becomes prohibitive in the case of 3D, unstructured, polyhedral meshes. A new procedure is outlined in this paper to avoid direct subzonal remapping. The new algorithm reconstructs the spatial profile of a subzonal variable using remapped zonal and nodal representations of the data. The reconstruction procedure is cast as an under-constrained optimization problem. Enforcing conservation at each zone and node on the remapped mesh provides the set of equality constraints; the objective function corresponds to a quadratic variation per subzone between the values to be reconstructed and a set of target reference values. Numerical results for various pure-remapping and hydrodynamics tests are provided. Ideas for extending the algorithm to staggered-grid radiation-hydrodynamics are discussed as well as ideas for generalizing the algorithm to include inequality constraints.

[1]  Mikhail Shashkov,et al.  Extension of efficient, swept-integration-based conservative remapping method for meshes with changing connectivity , 2008 .

[2]  Marvin L. Adams,et al.  A piecewise linear finite element discretization of the diffusion equation for arbitrary polyhedral grids , 2008, J. Comput. Phys..

[3]  Raphaël Loubère,et al.  The repair paradigm: New algorithms and applications to compressible flow , 2006 .

[4]  Pavel Váchal,et al.  Optimization-based synchronized flux-corrected conservative interpolation (remapping) of mass and momentum for arbitrary Lagrangian-Eulerian methods , 2010, J. Comput. Phys..

[5]  Jacek Gondzio,et al.  Multiple centrality corrections in a primal-dual method for linear programming , 1996, Comput. Optim. Appl..

[6]  Marvin L. Adams,et al.  Subcell balance methods for radiative transfer on arbitrary grids , 1997 .

[7]  J. L. Taylor.  XXXVII. An exact solution of the spherical blast wave problem , 1955 .

[8]  L. Sedov.  Similarity and Dimensional Methods in Mechanics , 1960 .

[9]  Jeffrey Grandy,et al.  Conservative Remapping and Region Overlays by Intersecting Arbitrary Polyhedra , 1999 .

[10]  J. Michael Owen,et al.  Arbitrary Lagrangian Eulerian remap treatments consistent with staggered compatible total energy conserving Lagrangian methods , 2013, J. Comput. Phys..

[11]  M. Shashkov,et al.  The Construction of Compatible Hydrodynamics Algorithms Utilizing Conservation of Total Energy , 1998 .

[12]  T. Belytschko,et al.  A uniform strain hexahedron and quadrilateral with orthogonal hourglass control , 1981 .

[13]  Dimitri P. Bertsekas,et al.  Nonlinear Programming , 1997 .

[14]  Todd S. Palmer.  Discretizing the diffusion equation on unstructured polygonal meshes in two dimensions , 2001 .

[15]  Jérôme Breil,et al.  A swept‐intersection‐based remapping method in a ReALE framework , 2013 .

[16]  Sanjay Mehrotra,et al.  On the Implementation of a Primal-Dual Interior Point Method , 1992, SIAM J. Optim..

[17]  P Anninos.  Kull ALE: II. Grid Motion on Unstructured Arbitrary Polyhedral Meshes , 2002 .

[18]  Len G. Margolin,et al.  Remapping, recovery and repair on a staggered grid , 2004 .

[19]  B. Després,et al.  Convergence and Sensitivity Analysis of Repair Algorithms in 1D , 2006 .

[20]  Raphaël Loubère,et al.  A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods , 2005 .

[21]  W. F. Noh.  Errors for calculations of strong shocks using an artificial viscosity and artificial heat flux , 1985 .

[22]  Raphaël Loubère,et al.  ReALE: A reconnection-based arbitrary-Lagrangian-Eulerian method , 2010, J. Comput. Phys..