Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPUTER PROGRAM”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code x-point gyrokinetic code (XGC), in order to study density pump-out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of heat transport modeling using directly driven beryllium spheres

Recent experiments involving directly driven beryllium spheres are reported. Plasma conditions are measured using Thomson scattering with the probe beam pointed 200, 300, and 400 μ m from the surface of the sphere. Laser coupling is assessed using calorimeters that collect scattered light placed at various locations within the target chamber. Laser intensities of 1 0 14 W / c m 2 and 2.5 × 1 0 14 W / c m 2 are chosen to minimize unmodeled laser-plasma interactions (LPIs) that lead to laser-target decoupling. Two-dimensional simulations are compared to the interpreted data using the radiation-hydrodynamics code Lasnex. Heat transport is simulated using flux-limited Spitzer–Harm with both high ( f = 0.15) and low ( f = 0.03) flux limiters and the Schurtz–Nicolai-Busquet (SNB) model. At 1 0 14 W / c m 2 , all three heat transport models agree well with the measurement, demonstrating that the heat flux is local at low intensities near the measurement locations. At 2.5 × 1 0 14 W / c m 2 , the SNB and high flux model roughly match the plasma conditions but predict 2% uncoupled light compared to 10% measured. Additionally, the use of drive multipliers to match the measured coupled light does not alter the agreement between measured and simulated plasma conditions, suggesting that decoupling due to LPI is unlikely to alter this agreement. The low flux model cannot match the plasma conditions and results in 19% scattered light. The use of a resonant absorption model can be used to bring the simulated scattered light into agreement, but the simulated plasma conditions are still in disagreement with the measurement. For this reason, the low flux model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Deficiencies in compression and yield in x-ray-driven implosions

This paper analyzes x-ray–driven implosions that are designed to be less sensitive to 2-D and 3-D effects in hohlraum and capsule physics. Key performance metrics including the burn-averaged ion temperature, hot-spot areal density, and fusion yield are found to agree with simulations where the design adiabat (internal pressure) is multiplied by a factor of 1.4. Furthermore, these results motivate the development of a simple model for interpreting experimental data, which is then used to quantify how improvements in compression could help achieve ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

ITER cold VDEs in the limit of perfectly conducting walls

Recently, it has been shown that a vertical displacement event (VDE) can occur in ITER even when the walls are perfect conductors, as a consequence of the current quench (CQ) [A. H. Boozer, Phys. Plasmas 26, 114501 (2019)]. We used the extended-MHD code M3D-C1 with an ITER-like equilibrium and induced a CQ to explore cold VDEs in the limit of perfectly conducting walls, using different wall geometries. In the case of a rectangular first wall with the side walls far away from the plasma, we obtained very good agreement with the analytical model developed by Boozer that considers a top/bottom flat-plates wall. We show that the solution in which the plasma remains at the initial equilibrium position is improved when bringing the side walls closer to the plasma. When approximating the ITER first wall as a perfect conductor, the plasma remains stable at the initial equilibrium position far beyond the value predicted by the flat-plates wall limit. When considering an opposite limit in which only the inner shell of the ITER vacuum vessel acts as a perfect conductor, the plasma is displaced during the CQ, but the edge safety factor remains above 2 longer in the current decay compared to the flat-plates wall limit. In all the simulated cases, the VDE is found to be strongly dependent on the plasma current, in agreement with a similar finding in the flat-plates wall limit, showing an important difference with hot VDEs in which the CQ is not a necessary condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Sodium tracer measurements of an expanded dense aluminum plasma from e-beam isochoric heating

Spatially and temporally resolved visible absorption spectroscopy is performed on sodium D-lines present as surface contaminants on an expanded dense aluminum plasma plume. An 80-ns FWHM, intense, relativistic electron beam deposits 5.4 J into a 100- μ m-thick Al foil, which isochorically heats and subsequently hydrodynamically expands the material through the warm dense matter state and into a classical-like plasma state, with a coupling parameter of approximately 0.2 and a degeneracy parameter of approximately 270. The Na contamination, carried along with the expanding plume, shows saturated absorption features in the dense Al continuum for λ > 450 nm. X-ray photoelectron spectroscopy and laser-induced breakdown spectroscopy confirm Na is a surface contaminant with an atomic concentration of ~ 0.1 % when interrogating identical foil samples. Additioanlly, a spectroscopic-quality radiation transport model is used to post-process 2D hydrodynamic simulations to interpret the plasma conditions based on the measured Na 3p-3s doublet line profiles. A sodium number density of 3 × 10 15 cm -3 best matches the experimental spectra, which originate from a dense surface plasma with n e = 3.0 ± 0.8 × 10 18 cm -3 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The effects of multispecies Hohlraum walls on stimulated Brillouin scattering, Hohlraum dynamics, and beam propagation

We report experiments and simulations have been conducted to investigate the efficacy of Ta 2 O 5 -lined Hohlraum walls at reducing stimulated Brillouin backscattering (SBS) as well as any subsequent effects on the Hohlraum dynamics and capsule implosions in indirect drive experiments at the National Ignition Facility. Using a 1.1 MJ 400 TW, 351 nm, shaped laser pulse, we measure a 5× reduction in SBS power in the peak of the pulse from the wall on the outer 50° cone beams. The SBS spectrum indicates a reduction in the high-Z spectral signature when using multispecies wall materials. Detailed hydrodynamic simulations were performed using different heat conduction models with flux limiters. Additional simulations were performed on the plasma maps using the 3D parallel paraxial code pF3D to compare backscatter powers between the pure Au and Ta 2 O 5 -lined Hohlraums. Further analysis, using hydrodynamically equivalent plasmas, shows that the SBS reduction is clearly a result of the added ion Landau damping caused by the oxygen ions and not from differences in plasma conditions. The experimental and simulation results also show an increase in the wall plasma expansion when using the Ta 2 O 5 liner leading to a 70% more oblate implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fusion gain from cylindrical liner-driven implosions of field reversed configurations

MagLIF experiments [M.R. Gomez et al., Phys. Plasmas 22, 056306 (2015)] on Z have demonstrated the basic principles of Magneto-Inertial Fusion (MIF) for wall confined plasmas. Other MIF schemes have been proposed based on the liner implosion of closed field magnetically confined plasmas such as Field Reversed Configurations (FRCs) [T. P. Intrator et al., Phys. Plasmas 15, 042505 (2008)]. We present a semi-analytical model of liner driven FRC implosions that predicts the fusion gain of such systems. The model predicts a fusion gain near unity for an FRC imploded by a liner driven by the Z Machine. We show that FRCs could be formed and imploded at the Z facility using the AutoMag liner concept [S. A. Slutz et al., Phys. Plasmas 24, 012704 (2017)]. Furthermore, an initial bias magnetic field can be supplied by the external magnets used in MagLIF experiments. The reverse field is then supplied by an AutoMag liner, which has helical conducting paths imbedded in an insulating substance. Experiments [Shipley et al., Phys. Plasmas 26, 052705 (2019)] have demonstrated that AutoMag can generate magnetic fields greater than 30 Tesla inside of the liner. We have performed 2D Radiation MHD simulations of the formation and implosion of an FRC, which are in good agreement with the analytical model. The FRC formation process could be studied on small pulsed power machines delivering about 1 MA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dense hydrogen layers for high performance MagLIF

We report that Magnetized Liner Inertial Fusion (MagLIF) experiments driven by the Z machine produce >10 13 deuterium-deuterium fusion reactions. Simulations indicate high yields and gains (1000) with increased current and deuterium-tritium layers for burn propagation. Such a coating also isolates the metal liner from the gaseous fuel, which should reduce mixing of liner material into the fuel. However, the vapor density at the triple point is only 0.3 kg/m3, which is not high enough for MagLIF operation. We present two solutions to this problem. First, a fuel wetted low-density plastic foam can be used to form a layer on the inside of the liner. The desired vapor density can be obtained by controlling the temperature. This does however introduce carbon into the layer which will enhance radiation losses. Simulations indicate that this wetted foam layer can significantly contribute to the fusion yield when the foam density is less than 35 kg/m 3 . Second, we show that a pure frozen fuel layer can first be formed on the inside of the liner and then low temperature gaseous fuel can be introduced just before the implosion without melting a significant amount of the ice layer. This approach is the most promising for MagLIF to produce high yield and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Understanding and controlling capsule symmetry in near vacuum hohlraums at the National Ignition Facility

The near vacuum hohlraum platform is an inertial confinement fusion design at the National Ignition Facility (NIF) that uses the lowest practical density of helium gas of 30 μg/cc to fill the hohlraum, which is ten times lower than now used routinely. This has several advantages, such as high laser coupling; however, the inability to understand and simulate the symmetry of the imploded capsule has limited the use of this platform. This work presents the first simulations that are able to accurately capture the highly prolate implosion seen experimentally without unphysical, ad hoc model changes. While previous investigations attributed this asymmetry to multi-species interpenetration in the hohlraum, we find that this alone has little effect on symmetry. Instead, it is the presence of crossed-beam energy transfer (CBET), occurring with no applied wavelength shift between the laser beams, that increases the laser power to the inner cones and causes a more prolate implosion. The effect of CBET is increased in the simulation model when the hohlraum laser entrance hole hardware is included. Using this understanding, CBET is exploited by shifting the inner-beam wavelength by –0.75 Å (at 1ω) with respect to the outer-beams. This transfers laser power to the outer-beams in contrast to positive wavelength shifts as done routinely on NIF and produces a round capsule implosion in our simulations. Furthermore, this work shows the possibility of the near vacuum hohlraum as a viable experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Density measurements for the National Ignition Facility (NIF) opacity platform

The Opacity Platform on the National Ignition Facility (NIF) has been developed to measure opacities at varying densities and temperatures relevant to the solar interior and thermal cooling rates in white dwarf stars. The typical temperatures reached at NIF range between 150 and 210 eV, which allow these measurements to be performed experimentally. The captured opacities are crucial to validating radiation-hydrodynamic models that are used in astrophysics. The NIF opacity platform has a unique new capability that allows in situ measurement of the sample expansion. The sample expansion data are used to better understand the plasma conditions in our experiments by inferring the sample density throughout the duration of the laser drive. Here, we present the details of the density measurement technique, data analysis, and recent results for Fe and MgO.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

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

Measurement of early time outer laser beam reflection inside a cylindrical hohlraum at the National Ignition Facility

In a National Ignition Facility (NIF) cylindrical hohlraum, any specular reflection (“glint”) off the wall from the outer cone (incident angle of 50° and 44° relative to the surface of the wall) laser beam will irradiate the capsule poles. If the glint power is sufficiently large during the picket (early time) of the laser pulse, it may seed high-mode perturbations on the capsule surface that can grow during the implosion. To quantify the glint power on the capsule during the picket by the outer beams, we performed dedicated experiments on NIF using a flat witness foil as a surrogate for the capsule in a half-hohlraum target. We found that the measured glint power is approximately comparable to simulations using a nominal electron conduction flux limiter of f = 0.15, and over an order of magnitude lower than that predicted using f = 0.03 in the wall. Based on our current understanding, we conclude that the glint from the outer beams plays an insignificant role in capsule drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Design and analysis of dudded fuel experiments at the National Ignition Facility

Recent experiments conducted at the National Ignition Facility (NIF) within the past 2 years have achieved the burning plasma state and exceeded the Lawson criterion for the first time in the laboratory. Here, we report on a set of experiments where the deuterium and tritium (DT) ice layers were replaced with dudded tritium, hydrogen, and deuterium (THD) fuel mixtures to remove the influence of alpha-heating on hot spot dynamics. The hot spot compression and yield in the absence of alpha particle self-heating were measured to assess the proximity of NIF implosions toward the ignition cliff. We find that the “burn-off” Lawson parameters χnoα inferred from the THD experiments are in good agreement with the inferences from postshot simulations of the DT-layered implosions. The THD for burning plasma shot N210307 yielded χnoα≈0.88±0.03 while the THD for ignition shot N210808 yielded χnoα≈1.04±0.04. These results also provide important context for the observed variability in the repeat attempts of ignition shot N210808 since implosions on the ignition cliff are expected to exhibit very large variations in the fusion yield from small changes in the initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluating nonlocal heat transport in directly driven chromium spheres using x-ray spectroscopy

We report on experiments investigating heat transport in laser-generated plasmas using directly driven chromium spheres. The spheres are fielded at the OMEGA laser facility and are driven with laser intensities of 5×10 14 Wcm −2 . Plasma conditions in the corona and scattered light are measured experimentally and compared against predictions from two-dimensional (2D) radiation-hydrodynamic simulations using different heat transport models. Spectroscopic analysis of x-ray self-emission is used as an additional diagnostic. X-ray emission is integrated over a large region of the plasma, probing regions that are not observed by localized optical Thomson scattering. In particular, x-ray emission peaks near the plasma critical density, so emission from optically thin lines provides information on plasma conditions where nonlocal transport is most likely to be significant. Three common heat transport models are considered: local transport with flux limiters f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. Consistent with previous work, both the high-flux (f = 0.15) and SNB models show good agreement with experimentally measured plasma conditions in the corona despite overpredicting laser absorption, whereas the low-flux (f = 0.03) model fails to match any experimental data. Conditions inferred from x-ray self-emission line ratios support this conclusion during the period of laser peak power, although synthetic spectra for all models fail to match the experiment during the transient portions of the pulse. For these reasons, the low-flux model is again rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗