Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “adaptive meshing”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Characterizing laser preheat for laser-driven magnetized liner inertial fusion using soft x-ray emission

Laser heating of a gas in a cylindrical liner held in by a thin foil window is a critical process in the MagLIF (magnetized liner inertial fusion) concept. Window burn-through and gas heating for OMEGA scale MagLIF cylinders as a function of time have been determined using spectrally integrated soft x-ray diagnostics. Window laser absorption is classified in terms of the emitted x-rays from the window plasma as a function of laser energy, and shows that laser energy absorbed is weakly dependent on incident intensity. Radiation-hydrodynamic simulations overestimate the amount of laser energy absorbed by the window as evidenced by the increase in x-ray radiation across several photon energy bands compared to experiments. Gas temperatures inferred from soft x-ray emission from the front 1 mm of the liner are shown to evolve in time in a similar manner to simulation predictions. In conclusion, soft x-ray emission from the gas within the region of the liner that is normally imploded is shown to meet the 100 eV requirements set by the initial point design for laser-driven MagLIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental quantification of the impact of heterogeneous mix on thermonuclear burn

In inertial confinement fusion, deuterium–tritium (DT) fuel is brought to densities and temperatures where fusion ignition occurs. However, mixing of the ablator material into the fuel may prevent ignition by diluting and cooling the fuel. MARBLE experiments at the National Ignition Facility provide new insight into how mixing affects thermonuclear burn. These experiments use laser-driven capsules containing deuterated plastic foam and tritium gas. Embedded within the foam are voids of known sizes and locations, which control the degree of heterogeneity of the fuel. Initially, the reactants are separated, with tritium concentrated in the voids and deuterium in the foam. During the implosion, mixing occurs between the foam and gas materials, leading to DT fusion reactions in the mixed region. Here, it is shown that by measuring the ratios of DT and deuterium–deuterium neutron yields for different macropore sizes and gas compositions, the effects of mix heterogeneity on thermonuclear burn may be quantified, supporting an improved understanding of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A mechanism for reduced compression in indirectly driven layered capsule implosions

High-yield implosions on the National Ignition Facility rely on maintaining low entropy in the deuterium–tritium fuel, quantified by its adiabat, in order to efficiently couple energy to the hot spot through high compression of the fuel layer. We present very-high-resolution xRAGE simulation results that study the impacts of interfacial mixing and the jetting of materials due to surface defects, defects on internal interfaces, voids, and engineering features on fuel layer compression. Defects and voids are typically neglected in implosion simulations due to their small size and three-dimensional geometry. Our results showed that supersonic jets of material arise through weak spots in the shell at peak implosion velocity that prevent uniform compression of the fuel layer even when they do not introduce contaminant into the hot spot. This occurs despite maintaining low fuel entropy, since the formation of the weak spots involves nonradial displacement of fuel mass. In contrast, simulations show that fuel–ablator mixing due to interfacial instabilities has a much smaller impact on compression. We show that defects on interior interfaces of plastic capsules decrease compression by 15% to 25% and interfacial mixing between the ablator and fuel decreases compression by less than 1% for implosions with plastic or high-density carbon (HDC) ablators. For low adiabat implosions, the impact of jetting seeded by the support tent can also decrease the compression by 25%. We demonstrate that the inclusion of interior defects in simulations can explain the inferred compression in two fielded plastic capsule implosions and that the inclusion of voids, for which available characterization has large uncertainties, in simulations of HDC capsule implosions has a qualitatively consistent impact. This mechanism offers a potential explanation for persistently overestimated fuel compression in design simulations of layered implosions on the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Morphology of shocked lateral outflows in colliding hydrodynamic flows

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. Here, we introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules

Bayesian optimization has shown promise for the design optimization of inertial confinement fusion targets. Specifically, in Vazirani et al. [Phys. Plasmas 28 , 122709 (2021)], optimal designs for double shell capsules with graded inner shells were identified using one-dimensional xRAGE simulation yield calculations. While the machine learning models were able to accurately learn and predict one-dimensional simulation target performance, using simulations with higher fidelity would improve design optimization and better match with the expected experimental performance. However, higher fidelity physics modeling, i.e., two-dimensional xRAGE simulations, requires significantly larger computational time/cost, usually at least an order of magnitude, in comparison with one-dimensional simulations. This study presents a multi-fidelity Bayesian optimization, in which the machine learning model leverages low-fidelity (one-dimensional xRAGE) and high-fidelity (two-dimensional xRAGE) simulations to more accurately predict “pre-shot” target performance with respect to the expected experimental performance. By building a multi-fidelity Bayesian optimization framework coupled with xRAGE, the low-fidelity and high-fidelity simulations are able to inform one another, such that we have: (1) improved physics modeling in comparison with using low-fidelity simulations alone, (2) reduced computational time/cost in comparison with using high-fidelity simulations alone, and (3) more confidence in the expected performance of optimized targets during real-world experiments. In the future, we plan to use this robust multi-fidelity Bayesian optimization methodology to expedite the design of graded inner shells further and eventually full capsules as a part of the current double shell campaign at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the importance of three-dimensional modeling for high-energy-density physics experiments

Laser-driven cylindrical implosion experiments enable direct measurements of hydrodynamic instability growth in convergent geometries, providing a wealth of validation data in the high-energy-density regime. These experiments are designed to be nearly axially invariant, allowing for modeling with complementary two-dimensional slices of the cylinder. Two distinct hydrodynamics codes are employed to model a subset of these experiments, and the results are shown to be in very good agreement with each other and the available experimental data. While this 2D modeling approach adequately captures most of the physics of the implosion and ensuing instability growth, there are crucial aspects from the three-dimensional nature of the experiments that are missed in 2D. The first fully 3D simulations of these experiments are presented, and small but significant differences are found to arise from both the axial and azimuthal non-uniformity in the laser drive. Recent experimental results confirming the drive asymmetry are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulated impact of fill tube geometry on recent high-yield implosions at the National Ignition Facility

Inertial confinement fusion capsules fielded at the National Ignition Facility are filled with deuterium and tritium fuel by means of a fill tube. The fill tube introduces a low-density pathway into the fuel region of the capsule that allows high Z contaminant to invade the hot spot during the course of the implosion. A recent series of nominally identical high-yield implosions on the NIF has exhibited significant variability in performance. We evaluate the impact of the fill tube in these implosions computationally to determine whether variations in fill tube geometry could have contributed to this variability. The main contrast between the fill tube geometry in the six shots was the outer diameter of the capsule bore hole, a conical hole into which the fill tube is inserted. In our simulations, the geometry of the bore hole can play a significant role in the development of nonlinear flows seeded by the fill tube. Additionally, we find that the amount of space between the bore hole and the fill tube is the primary factor that determines the amount of contaminant jetted into the hot spot by the fill tube and, in turn, the level of yield reduction due to the fill tube in our simulations. As a consequence, some capsules with 5 μm fill tubes are predicted to outperform capsules with 2 μm fill tubes. We also find that micrometer-scale changes to bore hole size can impact fusion yields by up to four times near the ignition threshold. Nevertheless, simulation trends do not reproduce experimental yield trends, suggesting that the fill tube geometry was not the primary factor contributing to the observed variability in performance and that the fill tube could be masking sensitivity to other asymmetries such as other micrometer-scale capsule defects like voids that were not included in our simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Laser-generated plasmas in length scales relevant for thin film growth and processing: simulation and experiment

In pulsed laser deposition, thin film growth is mediated by a laser-generated plasma, whose properties are critical for controlling the film microstructure. The advent of 2D materials has renewed the interest in how this ablation plasma can be used to manipulate the growth and processing of atomically thin systems. For such purpose, a quantitative understanding of the density, charge state, and kinetic energy of plasma constituents is needed at the location where they contribute to materials processes. Here, we study laser-induced plasmas over expansion distances of several centimeters from the ablation target, which is the relevant length scale for materials growth and modification. The study is enabled by a fast implementation of a laser ablation/plasma expansion model using an adaptive Cartesian mesh solver. Simulation outcomes for KrF excimer laser ablation of Cu are compared with Langmuir probe and optical emission spectroscopy measurements. Furthermore, simulation predictions for the plasma-shielding threshold, the ionization state of species in the plasma, and the kinetic energy of ions, are in good correspondence with experimental data. For laser fluences of 1–4 J cm -2 , the plume is dominated by Cu 0 , with small concentrations of Cu+ and electrons at the expansion front. Higher laser fluences (e.g. 7 J cm -2 ) lead to a Cu + -rich plasma, with a fully ionized leading edge where Cu 2+ is the dominant species. In both regimes, simulations indicate the presence of a low-density, high-temperature plasma expansion front with a high degree of ionization that may play a significant role in doping, annealing, and kinetically-driven phase transformations in 2D materials.

2D materials↗
DOE OSTI2014

EVOLUTION OF COLD STREAMS AND THE EMERGENCE OF THE HUBBLE SEQUENCE

A new physical framework for the emergence of the Hubble sequence is outlined, based on novel analyses performed to quantify the evolution of cold streams of a large sample of galaxies from a state-of-the-art ultra-high resolution, large-scale adaptive mesh-refinement hydrodynamic simulation in a fully cosmological setting. It is found that the following three key physical variables of galactic cold inflows crossing the virial sphere substantially decrease with decreasing redshift: the number of streams N {sub 90} that make up 90% of concurrent inflow mass flux, average inflow rate per stream M-dot {sub 90} and mean (mass flux weighted) gas density in the streams n {sub gas}. Another key variable, the stream dimensionless angular momentum parameter λ, is found to instead increase with decreasing redshift. Assimilating these trends and others naturally leads to a physically coherent scenario for the emergence of the Hubble sequence, including the following expectations: (1) the predominance of a mixture of disproportionately small irregular and complex disk galaxies at z ≥ 2 when most galaxies have multiple concurrent streams, (2) the beginning of the appearance of flocculent spirals at z ∼ 1-2 when the number of concurrent streams are about 2-3, (3) the grand-design spiral galaxies appear at z ≤ 1 when galaxies with only one major cold stream significantly emerge. These expected general trends are in good accord with observations. Early-type galaxies are those that have entered a perennial state of zero cold gas stream, with their abundance increasing with decreasing redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Parallel DTFE Surface Density Field Reconstruction

We improve the interpolation accuracy and efficiency of the Delaunay tessellation field estimator (DTFE) for surface density field reconstruction by proposing an algorithm that takes advantage of the adaptive triangular mesh for line-of-sight integration. The costly computation of an intermediate 3D grid is completely avoided by our method and only optimally chosen interpolation points are computed, thus, the overall computational cost is significantly reduced. The algorithm is implemented as a parallel shared-memory kernel for large-scale grid rendered field reconstructions in our distributed-memory framework designed for N-body gravitational lensing simulations in large volumes. We also introduce a load balancing scheme to optimize the efficiency of processing a large number of field reconstructions. Our results show our kernel outperforms existing software packages for volume weighted density field reconstruction, achieving ~ 10x speedup, and our load balancing algorithm gains an additional ~ 3.6x speedup at scales with ~ 16k processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗
DOE OSTI2014

SU-E-J-87: Lung Deformable Image Registration Using Surface Mesh Deformation for Dose Distribution Combination

Purpose: To allow a reliable deformable image registration (DIR) method for dose calculation in radiation therapy. This work proposes a performance assessment of a morphological segmentation algorithm that generates a deformation field from lung surface displacements with 4DCT datasets. Methods: From the 4DCT scans of 15 selected patients, the deep exhale phase of the breathing cycle is identified as the reference scan. Varian TPS EclipseTM is used to draw lung contours, which are given as input to the morphological segmentation algorithm. Voxelized contours are smoothed by a Gaussian filter and then transformed into a surface mesh representation. Such mesh is adapted by rigid and elastic deformations to match each subsequent lung volumes. The segmentation efficiency is assessed by comparing the segmented lung contour and the TPS contour considering two volume metrics, defined as Volumetric Overlap Error (VOE) [%] and Relative Volume Difference (RVD) [%] and three surface metrics, defined as Average Symmetric Surface Distance (ASSD) [mm], Root Mean Square Symmetric Surface Distance (RMSSD) [mm] and Maximum Symmetric Surface Distance (MSSD) [mm]. Then, the surface deformation between two breathing phases is determined by the displacement of corresponding vertices in each deformed surface. The lung surface deformation is linearly propagated in the lung volume to generate 3D deformation fields for each breathing phase. Results: The metrics were averaged over the 15 patients and calculated with the same segmentation parameters. The volume metrics obtained are a VOE of 5.2% and a RVD of 2.6%. The surface metrics computed are an ASSD of 0.5 mm, a RMSSD of 0.8 mm and a MSSD of 6.9 mm. Conclusion: This study shows that the morphological segmentation algorithm can provide an automatic method to capture an organ motion from 4DCT scans and translate it into a volume deformation grid needed by DIR method for dose distribution combination.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Enhancing Least-Squares Finite Element Methods Through a Quantity-of-Interest

Here, we introduce an approach that augments least-squares finite element formulations with user-specified quantities-of-interest. The method incorporates the quantity-of-interest into the least-squares functional and inherits the global approximation properties of the standard formulation as well as increased resolution of the quantity-of-interest. We establish theoretical properties such as optimality and enhanced convergence under a set of general assumptions. Central to the approach is that it offers an element-level estimate of the error in the quantity-of-interest. As a result, we introduce an adaptive approach that yields efficient, adaptively refined approximations. Several numerical experiments for a range of situations are presented to support the theory and highlight the effectiveness of our methodology. Notably, the results show that the new approach is effective at improving the accuracy per total computational cost.

97 MATHEMATICS AND COMPUTING↗