The internal consistency, stability, and accuracy of the discrete, compatible formulation of Lagrangian hydrodynamics.

*(English)*Zbl 1161.76538Summary: This work explores the somewhat subtle meaning and consequences of the salient properties of the discrete, compatible formulation of Lagrangian hydrodynamics. In particular, since this formulation preserves total energy to roundoff error, the amount of error in the conservation of total energy cannot be used to gauge the internal consistency of calculations, as is often done with the older forms of this algorithm. However, the compatible formulation utilizes two definitions of zone volume: the first is the usual definition whereby the volume of a zone is defined as some prescribed function of the coordinates of the points that define it; the second is given as the integration in time of the continuity equation for zone volume as expressed in Lagrangian form. It is the use of this latter volume in the specific internal energy equation that enables total energy to be exactly conserved. These two volume definitions are generally not precisely equal. It is the analysis of this difference that forms the first part of this study. It is shown that this difference in zone volumes can be used to construct a practical internal consistency measure that not only takes the place of the lack of total energy conservation of the older forms of Lagrangian hydrodynamics, but is more general in that it can be defined on a single zone basis. It can also be used to ascertain the underlying spatial and temporal order of accuracy of any given set of calculations. The difference in these two definitions of zone volume may be interpreted as a type of entropy error. However, this entropy error is found to be significant only when a given calculation becomes numerically unstable, otherwise it remains at or far beneath truncation error levels. In fact, it can be utilized to provide an upper bound on the size of the spatial truncation error for a stable computation. It is also shown how this volume difference can be used as an indicator of numerical difficulties, since exact local conservation of total energy does not guarantee numerical stability or the quality of any numerical calculation. The discrete, compatible formulation of Lagrangian hydrodynamics utilizes a two level predictor/corrector-type of time integration scheme; a stability analysis, both analytical and numerical, is given. This analysis reveals a novel stability diagram that has not been heretofore published, and gives definitive information as to how the stabilizing corrector step should be centered in time.

##### MSC:

76M25 | Other numerical methods (fluid mechanics) (MSC2010) |

76N15 | Gas dynamics (general theory) |

##### Keywords:

Lagrangian scheme; stability; compressible fluids; compatible discretization; support operator
PDF
BibTeX
XML
Cite

\textit{A. L. Bauer} et al., J. Comput. Phys. 218, No. 2, 572--593 (2006; Zbl 1161.76538)

Full Text:
DOI

**OpenURL**

##### References:

[1] | H.A. Bethe, R. Christy, R.R. Davis, R.P. Feynman, S. Frankel, E.J. Konopinski, N. Livesay, N. Metropolis, E. Nelson, IBM calculations of implosion hydrodynamics, Los Alamos National Laboratory report, LA-94, June 20, 1944, unclassified, cover only. |

[2] | Richtmyer, R.D.; Morton, K.W., Difference methods for initial-value problems, (1967), Interscience New York, (Chapter 12) · Zbl 0155.47502 |

[3] | B. Alder, S. Fernbach, M. Rotenberg (Eds.), Methods in Computational Physics, vol. 3, 1964, various articles. |

[4] | VonNeumann, J.; Richtmyer, R.D., A method for the calculation of hydrodynamic shocks, J. appl. phys., 21, 232-237, (1950) · Zbl 0037.12002 |

[5] | W.B. Goad, WAT: a numerical method for two-dimensional unsteady fluid flow, Los Alamos National Laboratory Report, LAMS-2365, 1960. |

[6] | Caramana, E.J.; Whalen, P.P., Numerical preservation of symmetry properties of continuum problems, J. comput. phys., 141, 174-198, (1998) · Zbl 0933.76066 |

[7] | Favorskii, A.P., Variational-discrete models of hydrodynamics equations, Diff. equations, 16, 1308-1321, (1980) |

[8] | L.G. Margolin, T.F. Adams, Spatial differencing for finite difference codes, Los Alamos National Laboratory Report, LA-10249, 1985. |

[9] | Samarskii, A.A.; Tishkin, V.F.; Favorskii, A.P.; Shashkov, M.J., Operational finite difference schemes, Diff. equations, 17, 854-862, (1981) · Zbl 0485.65060 |

[10] | L.G. Margolin, A.E. Tarwater, A diffusion operator for Lagrangian meshes, Lawrence Livermore National Laboratory Report, UCRL-95652, 1986. |

[11] | Burton, D.E., Exact conservation of energy and momentum in staggered-grid hydrodynamics with arbitrary connectivity, Advances in the free-Lagrange method, (1990), Springer New York |

[12] | D.E. Burton, Multidimensional discretizations of conservation laws for unstructured polyhedral grids, Lawrence Livermore National Laboratory Report, UCRL-JC-118306, 1994. |

[13] | Caramana, E.J.; Burton, D.E.; Shashkov, M.J.; Whalen, P.P., The construction of compatible hydrodynamics algorithms utilizing conservation of total energy, J. comput. phys., 146, 227-262, (1998) · Zbl 0931.76080 |

[14] | R. Loubère, E.J. Caramana, The force/work differencing of exceptional points in the discrete, compatible formulation of Lagrangian hydrodynamics, LAUR-04-8906, J. Comput. Phys., in press, doi:10.1016/j.jcp.2005.11.022. |

[15] | Caramana, E.J.; Shashkov, M.J.; Whalen, P.P., Formulations of artificial viscosity for multi-dimensional shock wave computations, J. comput. phys., 144, 70-97, (1998) · Zbl 1392.76041 |

[16] | Caramana, E.J.; Shashkov, M.J., Elimination of artificial grid distortion and hourglass-type motions by means of Lagrangian subzonal masses and pressures, J. comput. phys., 142, 521-561, (1998) · Zbl 0932.76068 |

[17] | Caramana, E.J.; Loubère, R., Curl-q: a vorticity damping artificial viscosity for Lagrangian hydrodynamics calculations, LAUR-05-6618, J. comput. phys., 215, 385-391, (2006) · Zbl 1173.76380 |

[18] | Caramana, E.J.; Rousculp, C.L.; Burton, D.E., A compatible, energy and symmetry preserving Lagrangian hydrodynamics algorithm in three-dimensional Cartesian geometry, J. comput. phys., 157, 89-119, (2000) · Zbl 0961.76049 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.