A compatible, energy and symmetry preserving Lagrangian hydrodynamics algorithm in three-dimensional Cartesian geometry

This work presents a numerical algorithm for the solution of fluid dynamics problems with moderate to high speed flow in three dimensions. Cartesian geometry is chosen owing to the fact that in this coordinate system no curvature terms are present that break the conservation law structure of the fluid equations. Written in Lagrangian form, these equations are discretized utilizing compatible, control volume differencing with a staggered-grid placement of the spatial variables. The concept of "compatibility" means that the forces used in the momentum equation to advance velocity are also incorporated into the internal energy equation so that these equations together define the total energy as a quantity that is exactly conserved in time in discrete form. Multiple pressures are utilized in each zone; they produce forces that resist spurious vorticity generation. This difficulty can severely limit the utility of the Lagrangian formulation in two dimensions and make this representation otherwise virtually useless in three dimensions. An edge-centered artificial viscosity whose magnitude is regulated by local velocity gradients is used to capture shocks. The particular difficulty of exactly preserving one-dimensional spherical symmetry in three-dimensional geometry is solved. This problem has both practical and pedagogical significance. The algorithm is suitable for both structured and unstructured grids. Limitations that symmetry preservation imposes on the latter type of grids are delineated.