Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Rheology of dense fiber suspensions: Origin of yield stress, shear thinning, and normal stress differences

We use high-fidelity computational model of fiber suspension to investigate suspension rheology. We investigate the interplay between the hydrodynamic, noncontact attractive and repulsive, and interfiber contact interactions. The shear-thinning viscosity and finite yield stress obtained from the Immerse Boundary Method simulations align quantitatively with experimental findings from the literature. The study demonstrates that attractive interactions lead to both yield stress and shear thinning behavior in rigid fiber suspensions. This discovery holds significance as it contributes to the ongoing debate on the source of yield stress in fiber suspensions. The proposed model is used to quantify normal stresses in addition to shear thinning and yield stress.

Ardekani, Arezoo↗
DOE OSTI2014

Risk Evaluation for CO2 Geosequestration in the Knox Supergroup, Illinois Basin Final Report

This report describes a process and provides seed information for identifying and evaluating risks pertinent to a hypothetical carbon dioxide (CO2) capture and sequestration (CCS) project. In the envisioned project, the target sequestration reservoir rock is the Potosi Formation of the Knox Supergroup. The Potosi is identified as a potential target formation because (1) at least locally, it contains vuggy to cavernous layers that have very high porosity, and (2) it is present in areas where the deeper Mt. Simon Sandstone (a known potential reservoir unit) is absent or nonporous. The key report content is discussed in Section 3.3, which describes two lists of Features, Events, and Processes (FEPs) that should be considered during the design stage of such a project. These lists primarily highlight risk elements particular to the establishment of the Potosi as the target formation in general. The lists are consciously incomplete with respect to risk elements that would be relevant for essentially all CCS projects regardless of location or geology. In addition, other risk elements specific to a particular future project site would have to be identified. Sources for the FEPs and scenarios listed here include the iconic Quintessa FEPs list developed for the International Energy Agency Greenhouse Gas (IEAGHG) Programme; previous risk evaluation projects executed by Schlumberger Carbon Services; and new input solicited from experts currently working on aspects of CCS in the Knox geology. The projects used as sources of risk information are primarily those that have targeted carbonate reservoir rocks similar in age, stratigraphy, and mineralogy to the Knox-Potosi. Risks of using the Potosi Formation as the target sequestration reservoir for a CCS project include uncertainties about the levels of porosity and permeability of that rock unit; the lateral consistency and continuity of those properties; and the ability of the project team to identify suitable (i.e., persistently porous and permeable) injection depths within the overall formation. Less direct implications include the vertical position of the Potosi within the rock column and the absence of a laterally extensive shale caprock immediately overlying the Potosi. Based on modeling work done partly in association with this risk report, risks that should also be evaluated include the ability of available methods to predict and track the development of a CO2 plume as it migrates away from the injection point(s). The geologic and hydrodynamic uncertainties present risks that are compounded at the stage of acquiring necessary drilling and injection permits. It is anticipated that, in the future, a regional geologic study or CO2-emitter request may identify a small specific area as a prospective CCS project site. At that point, the FEPs lists provided in this report should be evaluated by experts for their relative levels of risk. A procedure for this evaluation is provided. The higher-risk FEPs should then be used to write project-specific scenarios that may themselves be evaluated for risk. Then, actions to reduce and to manage risk can be described and undertaken. The FEPs lists provided as Appendix 2 should not be considered complete, as potentially the most important risks are ones that have not yet been thought of. But these lists are intended to include the most important risk elements pertinent to a Potosi-target CCS project, and they provide a good starting point for diligent risk identification, evaluation, and management.

Hnottavange-Telleen, Ken↗
DOE OSTI2026

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗
DOE OSTI2016

Hydrodynamic effects on β-amyloid (16-22) peptide aggregation

Computer simulations based on simplified representations are routinely used to explore the early steps of amyloid aggregation. However, when protein models with implicit solvent are employed, these simulations miss the effect of solvent induced correlations on the aggregation kinetics and lifetimes of metastable states. In this work, we apply the multi-scale Lattice Boltzmann Molecular Dynamics technique (LBMD) to investigate the initial aggregation phases of the amyloid Aβ{sub 16−22} peptide. LBMD includes naturally hydrodynamic interactions (HIs) via a kinetic on-lattice representation of the fluid kinetics. The peptides are represented by the flexible OPEP coarse-grained force field. First, we have tuned the essential parameters that control the coupling between the molecular and fluid evolutions in order to reproduce the experimental diffusivity of elementary species. The method is then deployed to investigate the effect of HIs on the aggregation of 100 and 1000 Aβ{sub 16−22} peptides. We show that HIs clearly impact the aggregation process and the fluctuations of the oligomer sizes by favouring the fusion and exchange dynamics of oligomers between aggregates. HIs also guide the growth of the leading largest cluster. For the 100 Aβ{sub 16−22} peptide system, the simulation of ∼300 ns allowed us to observe the transition from ellipsoidal assemblies to an elongated and slightly twisted aggregate involving almost the totality of the peptides. For the 1000 Aβ{sub 16−22} peptides, a system of unprecedented size at quasi-atomistic resolution, we were able to explore a branched disordered fibril-like structure that has never been described by other computer simulations, but has been observed experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Numerical methods for the weakly compressible Generalized Langevin Model in Eulerian reference frame

A well established approach for the computation of turbulent flow without resolving all turbulent flow scales is to solve a filtered or averaged set of equations, and to model non-resolved scales by closures derived from transported probability density functions (PDF) for velocity fluctuations. Effective numerical methods for PDF transport employ the equivalence between the Fokker–Planck equation for the PDF and a Generalized Langevin Model (GLM), and compute the PDF by transporting a set of sampling particles by GLM (Pope (1985) [1]). The natural representation of GLM is a system of stochastic differential equations in a Lagrangian reference frame, typically solved by particle methods. A representation in a Eulerian reference frame, however, has the potential to significantly reduce computational effort and to allow for the seamless integration into a Eulerian-frame numerical flow solver. GLM in a Eulerian frame (GLMEF) formally corresponds to the nonlinear fluctuating hydrodynamic equations derived by Nakamura and Yoshimori (2009) [12]. Unlike the more common Landau–Lifshitz Navier–Stokes (LLNS) equations these equations are derived from the underdamped Langevin equation and are not based on a local equilibrium assumption. Similarly to LLNS equations the numerical solution of GLMEF requires special considerations. In this paper we investigate different numerical approaches to solving GLMEF with respect to the correct representation of stochastic properties of the solution. We find that a discretely conservative staggered finite-difference scheme, adapted from a scheme originally proposed for turbulent incompressible flow, in conjunction with a strongly stable (for non-stochastic PDE) Runge–Kutta method performs better for GLMEF than schemes adopted from those proposed previously for the LLNS. We show that equilibrium stochastic fluctuations are correctly reproduced.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Constrained optimization framework for interface-aware sub-scale dynamics models for voids closure in Lagrangian hydrodynamics

In hydrocodes, voids are used to represent vacuum and model free boundaries between vacuum and real materials. We give a systematic description of a new treatment of void closure in the framework of the multimaterial arbitrary Lagrangian–Eulerian (ALE) methods. This includes a new formulation of the interface-aware sub-scale-dynamics (IA-SSD) closure model for multimaterial cells with voids, which is used in the Lagrangian stage of our indirect ALE scheme. The results of the comprehensive testing of the new model are presented for one- and two-dimensional multimaterial calculations in the presence of voids. Finally, we also present a sneak peek of a realistic shaped charge calculation in the presence of voids and solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Shell stability and conditions analyzed using a new method of extracting shell areal density maps from spectrally resolved images of direct-drive inertial confinement fusion implosions

In warm target direct-drive ICF implosion experiments performed at the OMEGA laser facility, plastic microballoons doped with a titanium tracer layer in the shell and filled with deuterium gas were imploded using a low-adiabat shaped laser pulse. Continuum radiation emitted in the core is transmitted through the tracer layer and the resulting spectrum recorded with a gated multi-monochromatic x-ray imager (MMI). Titanium K-shell line absorption spectra observed in the data are due to transitions in L-shell titanium ions driven by the backlighting continuum. The MMI data consist of an array of spectrally resolved images of the implosion. These 2-D space-resolved titanium spectral features constrain the plasma conditions and areal density of the titanium doped region of the shell. The MMI data were processed to obtain narrow-band images and space resolved spectra of titanium spectral features. Shell areal density maps, ρL(x,y), extracted using a new method using both narrow-band images and space resolved spectra are confirmed to be consistent within uncertainties. We report plasma conditions in the titanium-doped region of electron temperature (T e ) = 400±28eV, electron number density (N e ) = 8.5x10 24 ±2.5x10 24 cm -3 , and average areal density <ρR> = 86±7mg/cm 2 . Fourier analysis of areal density maps reveals shell modulations caused by hydrodynamic instability growth near the fuel-shell interface in the deceleration phase. We observe significant structure in modes l = 2-9, dominated by l = 2. We extract a target breakup fraction of 7.1±1.5% from our Fourier analysis. A new method for estimating mix width is evaluated against existing literature and our target breakup fraction. We estimate a mix width of 10.5±1μm.

42 ENGINEERING↗
DOE OSTI2016

Observable consequences of event-by-event fluctuations of HBT radii

We explore the effects of event-by-event fluctuations of Hanbury Brown–Twiss (HBT) radii and show how they can be observed experimentally. The relation of measured HBT radii extracted from ensemble-averaged correlation functions to the mean of their event-by-event probability distribution is clarified. Here, we propose a method to experimentally determine the mean and variance of this distribution and test it on an ensemble of fluctuating events generated with the viscous hydrodynamic code VISH2+1. Using the same code, the sensitivity of the mean and variance of the HBT radii to the specific QGP shear viscosity $\eta$/s is studied. We report sensitivity of the mean pion HBT radii and their variances to the temperature dependence of \eta/s near the quark-hadron transition at a level similar (10-20%) to that which was previously observed for elliptic and quadrangular flow of charged hadrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

iPTF15dtg: a double-peaked Type Ic supernova from a massive progenitor

Type Ic supernovae (SNe Ic) arise from the core-collapse of H- (and He-) poor stars, which could either be single Wolf-Rayet (WR) stars or lower-mass stars stripped of their envelope by a companion. Their light curves are radioactively powered and usually show a fast rise to peak (~10-15 d), without any early (in the first few days) emission bumps (with the exception of broad-lined SNe Ic) as sometimes seen for other types of stripped-envelope SNe (e.g., Type IIb SN 1993J and Type Ib SN 2008D). Here, we have studied iPTF15dtg, a spectroscopically normal SN Ic with an early excess in the optical light curves followed by a long (~30 d) rise to the main peak. It is the first spectroscopically-normal double-peaked SN Ic to be observed. Our aim is to determine the properties of this explosion and of its progenitor star. Methods. Optical photometry and spectroscopy of iPTF15dtg was obtained with multiple telescopes. The resulting light curves and spectral sequence are analyzed and modeled with hydrodynamical and analytical models, with particular focus on the early emission. iPTF15dtg is a slow rising SN Ic, similar to SN 2011bm. Hydrodynamical modeling of the bolometric properties reveals a large ejecta mass (~10 M ⊙ ) and strong 56 Ni mixing. The luminous early emission can be reproduced if we account for the presence of an extended (≳500 R ⊙ ), low-mass (≳0.045 M ⊙ ) envelope around the progenitor star. Alternative scenarios for the early peak, such as the interaction with a companion, a shock-breakout (SBO) cooling tail from the progenitor surface, or a magnetar-driven SBO are not favored. In conclusion, the large ejecta mass and the presence of H- and He-free extended material around the star suggest that the progenitor of iPTF15dtg was a massive (≳35 M ⊙ ) WR star that experienced strong mass loss.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Managing photon flux in a miniaturized photoionization detector

Miniaturized photoionization detectors (PIDs) are used in conjunction with gas chromatography systems to detect volatile compounds in gases by collecting the current from the photoionized gas analytes. PIDs should be inexpensive and compatible with a wide range of analyte species. One such PID is based on the formation of a He plasma in a dielectric barrier discharge (DBD), which generates vacuum UV (VUV) photons from excited states of He to photoionize gas analytes. There are several design parameters that can be leveraged to increase the ionizing photon flux to gas analytes to increase the sensitivity of the PID. To that end, the methods to maximize the photon flux from a pulsed He plasma in a DBD-PID were investigated using a two-dimensional plasma hydrodynamics model. The ionizing photon flux originated from the resonance states of helium, He(3P) and He(2 1 P), and from the dimer excimer He 2 *. While the photon flux from the resonant states was modulated over the voltage pulse, the photon flux from He 2 * persisted long after the voltage pulse passed. Several geometrical optimizations were investigated, such as using an array of pointed electrodes. However, increasing the capacitance of the dielectric enclosing the plasma chamber had the largest effect on increasing the VUV photon fluence to gas analytes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Energy Flow in Thin Shell Implosions and Explosions

Energy flow and balance in convergent systems beyond petapascal energy densities controls the fate of late-stage stars and the potential for controlling thermonuclear inertial fusion ignition. Time-resolved x-ray self-emission imaging combined with a Bayesian inference analysis is used to describe the energy flow and the potential information stored in the rebounding spherical shock at 0.22 PPa (2.2 Gbar or billions of atmospheres pressure). This analysis, together with a simple mechanical model, describes the trajectory of the shell and the time history of the pressure at the fuel-shell interface, ablation pressure, and energy partitioning including kinetic energy of the shell and internal energy of the fuel. Here, the techniques used here provide a fully self-consistent uncertainty analysis of integrated implosion data, a thermodynamic-path independent measurement of pressure in the petapascal range, and can be used to deduce the energy ow in a wide variety of implosion systems to petapascal energy densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗
NASA NTRS2005

Medical and Scientific Evaluations aboard the KC-135. Microgravity-Compatible Flow Cytometer

A spaceflight-compatible flow cytometer would be useful for the diagnosis of astronaut illness during long duration spaceflight and for conducting in-flight research to evaluate the effects of microgravity on human physiology. Until recently, the primary limitations preventing the development of a spaceflight compatible flow cytometer have been largely mechanical. Standard commercially available flow cytometers are large, complex instruments that use high-energy lasers and require significant training to operate. Standard flow cytometers function by suspending the particles to be analyzed inside a sheath fluid for analysis. This requires the presence of several liters of sheath fluid for operation, and generates a corresponding amount of liquid hazardous waste. The particles are then passed through a flow cell which uses the fluid mechanical property of hydrodynamic focusing to place the cells in single-file (laminar flow) as they pass through a laser beam for scanning and evaluation. Many spaceflight experiments have demonstrated that fluid physics is dramatically altered in microgravity (MSF [Manned Space Flight] Fluid Physics Data Sheet-August 1997) and previous studies have shown that sheath-fluid based hydrodynamic focusing may also be altered during microgravity (Crucian et al, 2000). For these reasons it is likely that any spaceflight compatible design for a flow cytometer would abandon the sheath fluid requirement. The elimination of sheath fluid would remove both the problems of weight associated with large volumes of liquids as well as the large volume of liquid waste generated. It would also create the need for a method to create laminar particle flow distinct from the standard sheath-fluid based method. The spaceflight prototype instrument is based on a recently developed commercial flow cytometer possessing a novel flow cell design that creates single-particle laser scanning and evaluation without the need for sheath-fluid based hydrodynamic focusing. This instrument also possesses a number of design advances that make it conditionally microgravity compatible: it is highly miniaturized and lightweight, uses a low energy diode laser, has a small number of moving parts, does not use sheath fluid and does not generate significant liquid waste. Although possessing certain limitations, the commercial cytometer functions operationally like a standard bench top laboratory flow cytometer, aspirating liquid particle samples and generating histogram or dot-plot data in standard FCS file format. In its current configuration however, the cytometer is limited to three parameter/two-color capability (two color PMTs + forward scatter), does not allow compensation between colors, does not allow linear analysis and is operated by rather inflexible software with limited capabilities. This is due to the fact that the cytometer has been designed and marketed as an instrument specific to a few particular assays, not as a multipurpose cytometer.

Crucian, Brian↗
DOE OSTI2015

Brownian dynamics of confined rigid bodies

We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

A new multi-energy neutrino radiation-hydrodynamics code in full general relativity and its application to the gravitational collapse of massive stars

We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inelastic neutrino–electron scattering, thermal neutrino production via pair annihilation, and nucleon–nucleon bremsstrahlung are included. While the Einstein field equations and the spatial advection terms in the radiation-hydrodynamics equations are evolved explicitly, the source terms due to neutrino–matter interactions and energy shift in the radiation moment equations are integrated implicitly by an iteration method. To verify our code, we first perform a series of standard radiation tests with analytical solutions that include the check of gravitational redshift and Doppler shift. A good agreement in these tests supports the reliability of the GR multi-energy neutrino transport scheme. We then conduct several test simulations of core-collapse, bounce, and shock stall of a 15M{sub ⊙} star in the Cartesian coordinates and make a detailed comparison with published results. Our code performs quite well to reproduce the results of full Boltzmann neutrino transport especially before bounce. In the postbounce phase, our code basically performs well, however, there are several differences that are most likely to come from the insufficient spatial resolution in our current 3D-GR models. For clarifying the resolution dependence and extending the code comparison in the late postbounce phase, we discuss that next-generation Exaflops class supercomputers are needed at least.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Validation of radiative transfer computation with Monte Carlo method for ultra-relativistic background flow

We developed a three-dimensional radiative transfer code for an ultra-relativistic background flow-field by using the Monte Carlo (MC) method in the context of gamma-ray burst (GRB) emission. For obtaining reliable simulation results in the coupled computation of MC radiation transport with relativistic hydrodynamics which can reproduce GRB emission, we validated radiative transfer computation in the ultra-relativistic regime and assessed the appropriate simulation conditions. The radiative transfer code was validated through two test calculations: (1) computing in different inertial frames and (2) computing in flow-fields with discontinuous and smeared shock fronts. The simulation results of the angular distribution and spectrum were compared among three different inertial frames and in good agreement with each other. If the time duration for updating the flow-field was sufficiently small to resolve a mean free path of a photon into ten steps, the results were thoroughly converged. The spectrum computed in the flow-field with a discontinuous shock front obeyed a power-law in frequency whose index was positive in the range from 1 to 10 MeV. The number of photons in the high-energy side decreased with the smeared shock front because the photons were less scattered immediately behind the shock wave due to the small electron number density. The large optical depth near the shock front was needed for obtaining high-energy photons through bulk Compton scattering. Even one-dimensional structure of the shock wave could affect the results of radiation transport computation. Although we examined the effect of the shock structure on the emitted spectrum with a large number of cells, it is hard to employ so many computational cells per dimension in multi-dimensional simulations. Therefore, a further investigation with a smaller number of cells is required for obtaining realistic high-energy photons with multi-dimensional computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Differentiable lagrangian shock hydrodynamics with application to stable shock acceleration of density interfaces

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. Here, 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.

Hydrophysics↗