Differentiable Lagrangian Shock Hydrodynamics with Application to Stable Shock Acceleration of Density Interfaces
By: Kevin Korner , Brandon Talamini , Julian Andrej and more
Potential Business Impact:
Stops explosions from mixing and failing.
We develop a gradient based optimization approach for the equations of compressible, Lagrangian hydrodynamics and demonstrate how it can be employed to automatically uncover strategies to control hydrodynamic instabilities arising from shock acceleration of density interfaces. Strategies for controlling the Richtmyer-Meshkov instability (RMI) are of great benefit for inertial confinement fusion (ICF) where shock interactions with many small imperfections in the density interface lead to instabilities which rapidly grow over time. These instabilities lead to mixing which, in the case of laser driven ICF, quenches the runaway fusion process ruining the potential for positive energy return. We demonstrate that control of these instabilities can be achieved by optimization of initial conditions with (> 100) parameters. Optimizing over a large parameter space like this is not possible with gradient-free optimization strategies. This requires computation of the gradient of the outputs of a numerical solution to the equations of Lagrangian hydrodynamics with respect to the inputs. We show that the efficient computation of these gradients is made possible via a judicious application of (i) adjoint methods, the exact formal representation of sensitivities involving partial differential equations, and (ii) automatic differentiation (AD), the algorithmic calculation of derivatives of functions. Careful regularization of multiple operators including artificial viscosity and timestep control is required. We perform design optimization of > 100 parameter energy field driving the Richtmyer Meshkov instability showing significant suppression while simultaneously enhancing the acceleration of the interface relative to a nominal baseline case.
Similar Papers
A sharp-interface discontinuous Galerkin method for simulation of two-phase flow of real gases based on implicit shock tracking
Fluid Dynamics
Simulates explosions and gas flow very accurately.
Enhanced Diffuse Interface Method for Multiphase Flow Simulations Across All Mach Numbers
Numerical Analysis
Improves computer models for fast and slow liquids.
High-order staggered Lagrangian hydrodynamics (I): framework of the discretization scheme
Numerical Analysis
Makes computer simulations of explosions more accurate.