Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimization for fusion”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2024

On the effect of flux-surface shaping on trapped-electron modes in quasi-helically symmetric stellarators

Using a novel optimization procedure, it has been shown that the Helically Symmetric eXperiment stellarator can be optimized for reduced trapped-electron-mode (TEM) instability [Gerard et al., Nucl. Fusion 63, (2023) 056004]. Presently, with a set of 563 experimental candidate configurations, gyrokinetic simulations are performed to investigate the efficacy of available energy EA, quasi-helical symmetry, and flux-surface shaping parameters as metrics for TEM stabilization. It is found that lower values of EA correlate with reduced growth rates, but only when separate flux-surface shaping regimes are considered. Moreover, configurations with improved quasi-helical symmetry demonstrate a similar reduction in growth rates and less scatter compared to EA. Regarding flux-surface shaping, a set of helical shaping parameters is introduced that show increased elongation is strongly correlated with reduced TEM growth rates, however, only when the quasi-helical symmetry is preserved. Using a newly derived velocity-space-averaged TEM resonance operator, these trends are analyzed to provide insights into the physical mechanism of the observed stabilization. For elongation, stabilization is attributed to geometric effects that reduce the destabilizing particle drifts across the magnetic field. Regarding quasi-helical symmetry, the TEM resonance in the maximally resonant trapping well is shown to increase as the quasi-helical symmetry is broken, and breaking quasi-helical symmetry increases the prevalence of highly resonant trapping wells. While these results demonstrate the limitations of using any single metric as a linear TEM proxy, it is shown that quasi-helical symmetry and plasma elongation are highly effective metrics for reducing TEM growth rates in helical equilibria.

Physics↗
DOE OSTI2025

High-Power Impulse Magnetron Sputter Deposition of Boron Carbide with Full-Face Erosion Magnetron and Mixed Ar-Ne Plasma

Boron carbide (B 4 C) is an attractive inertial confinement fusion ablator material. The fabrication of B4C ablators by magnetron sputtering requires process optimization. To increase process flexibility, here we explore high-power impulse magnetron sputter (HiPIMS) deposition of B 4 C in pure Ar and mixed Ar-Ne plasmas. Here, the results show that higher plasma discharge currents can be reached with a mixed Ar-Ne plasma in the entire working pressure range studied (5 𝑡𝑜 50 mTorr). At 45 mTorr with 10% of Ne in the Ar-Ne mix, high peak target current densities of ~1 A cm −2 were demonstrated. Films deposited with such a mixed Ar-Ne plasma with a full-face erosion magnetron source on substrates biased at −25 V exhibited higher density and improved mechanical properties, albeit with higher compressive residual stresses compared to the case of HiPIMS deposition in a pure Ar plasma. This work demonstrates additional process flexibility of the HiPIMS discharge mode for the deposition of B 4 C coatings.

Ablator capsule↗
DOE OSTI2025

Feasibility of fusion plasma burn control via real-time, sub-divertor neutral gas isotopic and compositional analysis

The ability to provide fusion burn control without requiring physical access through the first wall and fuel breeding blankets, would be vital for any future, magnetically confined fusion power reactor. A multi-sensor, fusion fuel cycle exhaust, neutral gas analysis system on JET, capable of delivering real time data, and accessing only the sub-divertor region, provides an excellent example of such capability. Optimized for and operated during the deuterium–tritium experimental campaigns 2 and 3 (DTE2, DTE3), it is proving valuable for planning to explore fusion reactor burn control in ITER with a comparable diagnostic system called the Diagnostic Residual Gas Analyzer (DRGA). This paper aims to show feasibility of developing model-based controllers for ITER and next generation, reactor-relevant devices, by building both on the empirical experience in JET-DTE2, and on the already emerging experience on developing such models specifically for ITER. The paper begins with a specific use-case from JET-DTE2, pertaining to the observed sensitivity of the fusion neutron yield on the concentration of isotopic helium-3 ( 3 He), with data from one of the high-performance DT shots exhibited with emphasis on the 3 He measurement via the sub-divertor. Then, a first model is developed and then explored with simulations that aim to discover how well the controllers in the model react to either insufficient levels of 3 He or excessive levels of 3 He. The simulations then explore potential impact from a delay in the measurement (or the response) that would be comparable to the ∼1 s, conductance limited response for the ITER DRGA system, currently in its final design. The simulations show that control is feasible, and that its effectiveness is not significantly impacted by such delay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Angular-spectral filtering of recoil protons for optimization of fast neutron imaging employing proton converters

Fast neutron imaging is an important capability for diverse applications such as inertial confinement fusion diagnostics, cargo security, nuclear nonproliferation and arms control, and industrial inspection. Traditional phosphor image plates can be enhanced for fast neutron imaging using hydrogenous plastic converters which allow fast neutrons to scatter off hydrogen nuclei to produce energetic protons that can be recorded by the image plate. However, protons emitted by image plates are not constrained in their emission angle, which contributes to the blur of the resulting image. Here, we investigate two methods that can alter the spatial extent of converted protons that deposit energy in the image plate: reducing the converter thickness, and introducing a proton filter between the plastic converter and image plate to reduce the contribution of lower-energy, off-axis protons to the image. Here we determine the optimal plastic converter thickness for maximizing the signal intensity to be 2–3 mm through Monte Carlo simulations, and we benchmark this result against experimental measurements with a deuterium-tritium (DT) neutron generator. Next, we evaluate the image smearing and signal loss for various converters to show that solely reducing the converter thickness has the expected effect of reducing the blur from proton image smearing of the sharpness of an edge recorded on the image plate at the cost of reducing the signal intensity. The use of a proton filter is shown to achieve a similar improvement of edge sharpness as reducing the converter thickness while also sacrificing the signal intensity. We conclude that the use of proton energy filtering can improve the sharpness of fast neutron images in situations where the converter thickness cannot be reduced below some practical minimum. For more intense neutron sources, the signal intensity is of less concern, and optimizing the resolution of the image plate and therefore of the imaging system could have greater value. In these applications, proton filters may allow for improved fast neutron imaging measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2017

The Pulsed Fission-Fusion (PUFF) Concept for Deep Space Exploration and Terrestrial Power Generation

This team is exploring a modified Z-pinch geometry as a propulsion system, imploding a liner of liquid lithium onto a pellet containing both fission and fusion fuel. The plasma resulting from the fission and fusion burn expands against a magnetic nozzle, for propulsion, or a magnetic confinement system, for terrestrial power generation. There is considerable synergy in the concept; the lithium acts as a temporary virtual cathode, and adds reaction mass for propulsion. Further, the lithium acts as a radiation shield against generated neutrons and gamma rays. Finally, the density profile of the column can be tailored using the lithium sheath. Recent theoretical and experimental developments (e.g. tailored density profile in the fuel injection, shear stabilization, and magnetic shear stabilization) have had great success in mitigating instabilities that have plagued previous fusion efforts. This paper will review the work in evaluating the pellet sizes and z-pinch conditions for optimal PuFF propulsion. Trades of pellet size and composition with z-pinch power levels and conditions for the tamper and lithium implosion are evaluated. Current models, both theoretical and computational, show that a z-pinch can ignite a small (~1 cm radius) fission-fusion target with significant yield. Comparison is made between pure fission and boosted fission targets. Performance is shown for crewed spacecraft for high speed Mars round trip missions and near interstellar robotic missions. The PuFF concept also offers a solution for terrestrial power production. PuFF can, with recycling of the effluent, achieve near 100% burnup of fission fuel, providing a very attractive power source with minimal waste. The small size of PuFF relative to today's plants enables a more distributed power network and less exposure to natural or man-made disruptions.

Adams, Robert↗
DOE OSTI2024

A Gaussian Process-Based extended Goldak heat source model for finite element simulation of laser powder bed fusion additive manufacturing process

In this study, laser powder bed fusion (L-PBF) additive manufacturing (AM) is a key enabling technology to manufacture highly complex and integrated metallic structures. In L-PBF AM process, the melting of the metal powders and the layers underneath can be governed by either “conduction mode” or “keyhole mode”, with the keyhole mode reportedly leading to porosity and decreased strength and ductility by many studies. In part scale simulations, finite element (FE) model is often used to study the temperature distribution during printing and to predict the residual stress, where a volumetric heat flux with a Gaussian or a double ellipsoidal (Goldak) distribution is often applied as the laser heat source. However, the above heat source models can only capture the melt pool shape in the conduction mode, and fail to capture the transition to keyhole melting mode when the process parameters change. To overcome this inaccuracy, an extended Goldak heat source model is proposed by introducing a laser penetration term as a function of laser parameters obtained from a Gaussian-Process (GP) model. The model is validated by “2D pad” AlSi10Mg L-PBF experiments under a wide range of laser power, scan speed, and laser focus offset, and the results show the model successfully captures the measured melt pool shape in all conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fatigue behavior of low-cost, non-spherical Ti-6Al-4V powder processed via laser powder bed fusion

Hydride-dehydride Ti-6Al-4V powder with particle size distribution of 75–175 μm were used to manufacture fatigue samples. Optimized processing parameters were used to maintain a relative density of >99.5% in the samples. Fatigue tests were carried out in R = σ min /σ max = -1 condition under stresses ranging from 150 to 500 MPa and the results were compared to samples that were additively manufactured using spherical powders. The fatigue lives were similar between the two types of powder; i.e., powder morphology has no effect on fatigue performance. However, the high surface roughness and martensitic microstructure resulting from the laser powder bed fusion process result in poor performance compared to samples machined from conventionally produced materials. In conclusion, crack initiation occurred consistently from surface defects because of high surface roughness of the additively manufactured parts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Test of CP-invariance of the Higgs boson in vector-boson fusion production and in its decay into four leptons

A search for CP violation in the decay kinematics and vector-boson fusion production of the Higgs boson is performed in the H → ZZ * → 4 ℓ ( ℓ = e, μ ) decay channel. The results are based on proton-proton collision data produced at the LHC at a centre-of-mass energy of 13 TeV and recorded by the ATLAS detector from 2015 to 2018, corresponding to an integrated luminosity of 139 fb –1 . Matrix element-based optimal observables are used to constrain CP-odd couplings beyond the Standard Model in the framework of Standard Model effective field theory expressed in the Warsaw and Higgs bases. Differential fiducial cross-section measurements of the optimal observables are also performed, and a new fiducial cross-section measurement for vector-boson-fusion production is provided. All measurements are in agreement with the Standard Model prediction of a CP-even Higgs boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Precipitation hardening of laser powder bed fusion Ti-6Al-4V

Here, the laser powder bed fusion (PBF-L) additive manufacturing (AM) community has dedicated significant efforts into process optimization and control for defect-free Ti-6Al-4V. As defects become less of an issue for PBF-L Ti-6Al-4V, the processing-structure-properties (PSP) relationships between AM microstructures can now be explored to optimize mechanical properties. Lower temperature aging treatments around 550 °C in wrought Ti-6A-4V have been historically understood to precipitation harden the α phase with an ordered, hexagonal close packed (HCP) α 2 phase. The α 2 phase existed as nanoscale Ti 3 Al precipitates coherent with the parent α phase. The goal of the present investigation was to implement a vacuum heat treatment of 545 °C for 100 hours on PBF-L Ti-6Al-4V. This vacuum heat treatment took place after an initial hot isostatic pressure (HIP) treatment that decomposed the as-built, martensitic microstructure. The vacuum aging successfully produced nanoscale precipitates of α 2 phase within α-laths, confirmed via atom probe tomography (APT). Microstructural-length scale and quasi-static mechanical properties were investigated by nanoindentation and uniaxial tensile tests of miniaturized test specimens. Given the sensitivity of Ti-6Al-4V mechanical properties to small changes in chemistry, all test specimens originated from the same build. The α 2 precipitation resulted in significantly harder α-laths (≈ 2 GPa) as measured via nanoindentation, as well as a relative yield and ultimate tensile strength increase of 60 MPa and 38 MPa, respectively. Analysis of Variance (ANOVA) of the datasets revealed no statistical differences in total elongation between the HIPed and HIPed + aged specimens, indicative of the aging treatment producing a net benefit of strength with no loss in ductility.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Sensor fusion display evaluation using information integration models in enhanced/synthetic vision applications

Based on existing integration models in the psychological literature, an evaluation framework is developed to assess sensor fusion displays as might be implemented in an enhanced/synthetic vision system. The proposed evaluation framework for evaluating the operator's ability to use such systems is a normative approach: The pilot's performance with the sensor fusion image is compared to models' predictions based on the pilot's performance when viewing the original component sensor images prior to fusion. This allows for the determination as to when a sensor fusion system leads to: poorer performance than one of the original sensor displays, clearly an undesirable system in which the fused sensor system causes some distortion or interference; better performance than with either single sensor system alone, but at a sub-optimal level compared to model predictions; optimal performance compared to model predictions; or, super-optimal performance, which may occur if the operator were able to use some highly diagnostic 'emergent features' in the sensor fusion display, which were unavailable in the original sensor displays.

Foyle, David C.↗
DOE OSTI2025

Enhancing Unknown Waveform Detection by Learning Intra and Inter-domain Dependencies with Advanced Attention Fusion Mechanisms

Detection of unknown waveforms in mission-critical communications is a crucial area of interest for the Department of Energy (DoE). Traditional methods and recent deep learning-based approaches often assume that the training set includes all possible classes, which is impractical for detecting new waveforms. This limitation gives rise to the problem of open-set recognition (OSR), which involves correctly identifying known classes while detecting and rejecting unknown or unseen classes. To address this limitation, we propose a novel dual-domain complex-valued neural architecture that jointly processes time-domain and frequency-domain signal representations using transformer mechanisms. A transformer model is a deep learning architecture that uses self-attention mechanisms to process and learn relationships in sequential data. Our model employs a cosine similarity loss to extract domain-specific features and incorporates a transformer architecture in the latent space to weigh the importance of different features from the time and frequency domains. The transformer layer includes stacked self-attention and cross-attention modules to learn intra-domain and inter-domain dependencies, creating a more holistic signal representation. An attention-based fusion module intelligently combines the time and frequency-domain features using multi-head attention, enabling the network to learn the optimal feature for each domain in each input signal. Quantitative results demonstrate the impact of these architectural choices on overall performance, showing significant improvement after incorporating self and cross-attention modules and using complex attention fusion over simple weighted fusion. Our ongoing work will focus on addressing the limitations of threshold-based OSR methods by developing a novel generative framework that integrates a conditional diffusion probabilistic model (DPM). DPM is a generative framework that learns to synthesize complex data by reversing a gradual noising process using a neural network trained to denoise step-by-step. Our goal is to leverage the inherent strengths of DPMs for identifying unknown signals more robustly. One primary advantage of using a DPM is its ability to provide a more reliable anomaly score based on the model's reconstruction error, rather than relying solely on classifier confidence. Additionally, the iterative denoising process of DPMs makes this approach naturally resilient to low Signal-to-Noise Ratio (SNR) conditions, where traditional methods often fail. By implementing this generative framework, we aim to enhance the model's capability to accurately detect unknown waveforms and maintain performance in challenging environments.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

ARC physics basis–magnetohydrodynamics

ARC is designed to produce 400 ⁢MW of net electricity and prove the commercial feasibility of a fusion power plant. In order to achieve this goal ARC has to operate with optimal core performance in a stationary scenario that minimises wear on the first wall and divertor. This requires avoiding or mitigating magnetohydrodynamic (MHD) instabilities which have the potential to not only degrade the plasma core but also lead to deleterious transient heat loads on plasma facing components. Therefore, this work aims at characterising the MHD stability of the high performance ARC scenario and inform the design of error field correction coils. Firstly, simulations of vertical displacement events show that an in-vessel coil is not needed and instead the poloidal shaping coils can be used to control vertical stability. These simulations also inform the demands on the corresponding coil power supplies. Stability analysis of the ideal kink mode with or without a conducting wall and kinetic effects suggests that the ARC baseline scenario operates deeply in the stable region. Using RDCON, tearing modes at the 𝑚/𝑛 =2/1 and 3/2 surfaces (with poloidal mode number 𝑚, and toroidal mode number 𝑛) are shown to be linearly stable, and including thermal transport effects in the rational surfaces lead to further stabilisation. However, other transient plasma instabilities can seed neoclassical tearing modes (NTMs). The marginally stable width of NTMs in ARC strongly depends on the internal inductance and can fall below 0.1% of the normalised poloidal flux. Furthermore, an empirical cross-machine model of the 𝑛 =1 error field leading to a disruption predicts a critical error field larger than SPARC but smaller than ITER. Three-dimensional coils can be designed with the Generalised Purturbed Equilbium Code based on a simple model that calculates the maximum correctable error field that is limited by the neoclassical toroidal viscosity torque. Broad scans of different coil geometries identify a set of 2 rows of off-midplane coils to be a suitable solution. It is also determined that such a set of three-dimensional coils is capable of correcting 𝑛 =2 error fields to some degree and creating strong enough 𝑛 =2 or 𝑛 =3 edge resonant perturbation fields for the suppression of edge-localised modes at reasonable coil currents. The final design of the first ARC will be further informed by results from SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Very high cycle fatigue characteristics of laser beam powder bed fused AlSi10Mg: A systematic evaluation of part geometry

This study explores the influence of geometry and part size on defect distribution, melt pool size, and mechanical characteristics in laser beam powder bed fused (LB-PBF) AlSi10Mg. Five distinct geometries—hourglass, small rod, small block, large block, and large rod—were fabricated under identical process parameters. Fully reversed ultrasonic fatigue testing, operating at a frequency of 20 kHz, was conducted to assess the very high cycle fatigue properties. The findings indicated that part geometry had a major impact on the fatigue properties of the material in the very high cycle fatigue regime. Specimens machined from large rods and large blocks had the lowest porosity and highest fatigue resistance. Microstructural analysis indicated that hourglass, small rod, and small block specimens had shallower melt pools and overlap depths compared to other geometries. This observation suggests a higher cooling rate in specimens with smaller cross-sectional areas, leading to the increased presence of entrapped gas pores and a lack of fusion defects. Understanding the relationship between part geometry and fatigue properties in LB-PBF components offers insights for optimizing design and manufacturing processes in additive manufacturing applications.

AlSi10Mg↗
NASA NTRS1992

Proposed evaluation framework for assessing operator performance with multisensor displays

Despite aggressive work on the development of sensor fusion algorithms and techniques, no formal evaluation procedures have been proposed. Based on existing integration models in the literature, an evaluation framework is developed to assess an operator's ability to use multisensor, or sensor fusion, displays. The proposed evaluation framework for evaluating the operator's ability to use such systems is a normative approach: The operator's performance with the sensor fusion display can be compared to the models' predictions based on the operator's performance when viewing the original sensor displays prior to fusion. This allows for the determination as to when a sensor fusion system leads to: 1) poorer performance than one of the original sensor displays (clearly an undesirable system in which the fused sensor system causes some distortion or interference); 2) better performance than with either single sensor system alone, but at a sub-optimal (compared to the model predictions) level; 3) optimal performance (compared to model predictions); or, 4) super-optimal performance, which may occur if the operator were able to use some highly diagnostic 'emergent features' in the sensor fusion display, which were unavailable in the original sensor displays. An experiment demonstrating the usefulness of the proposed evaluation framework is discussed.

Foyle, David C.↗
DOE OSTI2026

Primary Heat Transport System Design Considerations for Xcimer Energy’s Athena Fusion Pilot Plant

Fusion energy promises a reliable, carbon-free source of power; however, significant challenges remain before it can be deployed as an economical energy source. In addition to achieving fusion conditions, power plants must operate under extreme temperatures, radiation, and mechanical loads while maintaining high efficiency and availability. These requirements place strong demands on engineering design and plant operation. This work focuses on the engineering challenges associated with balance of plant analysis for inertial fusion energy systems. In particular, this paper examines the design considerations for primary heat transfer systems in fusion pilot plants employing molten fluoride salt coolants, with particular emphasis on system layout optimization and the balance between competing design objectives using the Xcimer Energy Athena inertial pilot plant design as a case study. Through systematic analysis of candidate system configurations and parametric sensitivity studies, we identify key engineering trade-offs governing salt inventory, pumping power requirements, and operational flexibility. The analysis employs system-level modeling tools to explore the design space and establish relationships between geometric parameters and system performance metrics.

Greenwood, Scott [ORNL] (ORCID:0000000333480736)↗
DOE OSTI2025

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗
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↗