Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer simulation”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Calorimeter with Bayesian unfolding of spectra of high-flux broadband x rays

We report the development of a multipurpose differential x-ray calorimeter with a broad energy bandwidth. The absorber architecture is combined with a Bayesian unfolding algorithm to unfold high energy x-ray spectra generated in high-intensity laser–matter interactions. Particularly, we show how to extract absolute energy spectra and how our unfolding algorithm can reconstruct features not included in the initial guess. The performance of the calorimeter is evaluated via Monte Carlo generated data. The method accuracy to reconstruct electron temperatures from bremsstrahlung is shown to be 5% for electron temperatures from 1 to 50 MeV. We study bremsstrahlung generated in solid target interaction showing an electron temperature of 0.56 ± 0.04 MeV for a 700 μm Ti titanium target and 0.53 ± 0.03 MeV for a 50 μm target. We investigate bremsstrahlung from a target irradiated by laser-wakefield accelerated electrons showing an endpoint energy of 551 ± 5 MeV, inverse Compton generated x rays with a peak energy of 1.1 MeV, and calibrated radioactive sources. The total energy range covered by all these sources ranges from 10 keV to 551 MeV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

A scintillator attenuation spectrometer for intense gamma-rays

A new type of compact high-resolution high-sensitivity gamma-ray spectrometer for short-pulse intense gamma-rays (250 keV to 50 MeV) has been developed by combining the principles of scintillators and attenuation spectrometers. The first prototype of this scintillator attenuation spectrometer (SAS) was tested successfully in Trident laser experiments at LANL. Later versions have been used extensively in the Texas Petawatt laser experiments in Austin, TX, and more recently in OMEGA-EP laser experiments at LLE, Rochester, NY. The SAS is particularly useful for high-repetition-rate laser applications. Furthermore, we give a concise description of the design principles, capabilities, and sample preliminary results of the SAS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

A new class of variable-radii diffraction optics for high-resolution x-ray spectroscopy at the National Ignition Facility (invited)

A new class of crystal shapes has been developed for x-ray spectroscopy of point-like or small (a few mm) emission sources. These optics allow for dramatic improvement in both achievable energy resolution and total throughput of the spectrometer as compared with traditional designs. This class of crystal shapes, collectively referred to as the Variable-Radii Spiral (VR-Spiral), utilize crystal shapes in which both the major and minor radii are variable. A crystal using this novel VR-Spiral shape has now been fabricated for high-resolution Extended X-ray Absorption Fine Structure (EXAFS) experiments targeting the Pb-L 3 (13.0 keV) absorption edge at the National Ignition Facility. The performance of this crystal has been characterized in the laboratory using a microfocus x-ray source, showing that high-resolution high-throughput EXAFS spectra can be acquired using this geometry. Importantly, these successful tests show that the complex three-dimensional crystal shape is manufacturable with the required precision needed to realize the expected performance of better than 5 eV energy resolution while using a 30 mm high crystal. An improved generalized mathematical form for VR-Spiral shapes is also presented allowing improved optimization as compared to the first sinusoidal-spiral based design. Finally, this new formulation allows VR-Spiral spectrometers to be designed at any magnification with optimized energy resolution at all energies within the spectrometer bandwidth.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Erratum: “DFTB+, a software package for efficient approximate density functional theory based atomistic simulations” [J. Chem. Phys. 152, 124101 (2020)]

The implementation of the GPU support in DFTB+, as described in Sec. III C of the original publication,1 was developed based on a previous unpublished implementation by Jacek Jakowski. In order to acknowledge his work on this first implementation, the authors of the original publication wish to include J. Jakowski as co-author. The scientific content of the original publication is not affected.

42 ENGINEERING↗
DOE OSTI2022

Preface for frontiers of magnetic reconnection research in heliophysical, astrophysical, and laboratory plasmas

Magnetic reconnection—the topological rearrangement of magnetic field—underlies many explosive phenomena across a wide range of natural and laboratory plasmas.3 It plays a pivotal role in electron and ion heating, particle acceleration to high energies, energy transport, and self-organization. Reconnection can have a complex relationship with turbulence at both large and small scales, leading to various effects that are only beginning to be understood. In heliophysics, magnetic reconnection plays a key role in solar flares, coronal mass ejections, coronal heating, solar wind dissipation, the interaction of interplanetary plasma with magnetospheres, dynamics of planetary magnetospheres such as magnetic substorms, and the heliospheric boundary with the interstellar medium. Here, the magnetic reconnection is integral to the solar and planetary dynamo processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Monte Carlo N-Particle forward modeling for density reconstruction of double shell capsule radiographs

In the Double Shell Inertial Confinement Fusion concept, characterizing the shape asymmetry of imploding metal shells is vital for understanding energy-efficient compression and radiative losses of the thermonuclear fuel. The Monte Carlo N-Particle MCNP® code forward models radiography of Double Shell capsule implosions using the Advanced Radiographic Capability at the National Ignition Facility. A procedure is developed for using MCNP to reconstruct density profiles from the radiograph image intensity. For a given Double Shell imploding target geometry, MCNP radiographs predict image contrast, which can help guide experimental design. Finally, in future work, the calculated MCNP synthetic radiographs will be compared with experimental radiographs to determine the radial and azimuthal density profiles of the Double Shell capsules.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Scaling patch analysis of planar turbulent mixing layers

Proper scales for the mean flow and Reynolds shear stress in planar turbulent mixing layers are determined from a scaling patch analysis of the mean continuity and momentum equations. By seeking an admissible scaling of the mean continuity equation, a proper scale for the mean transverse flow is determined as V ref =(dδ/dx)U ref , where dδ/dx is the growth rate of the mixing layer width and U ref =U h –U l is the difference between the velocity of the high speed stream U h and the velocity of the low speed stream U l . By seeking an admissible scaling for the mean momentum equation, a proper scale for the kinematic Reynolds shear stress is determined as R uv,ref =U avg V ref =[$\frac{1}{2A_u}$$\frac{dδ}{dx}$]$U^{2}_{ref}$ where A u $_{=}^{def}$(U h –U l )/(U h +U l ) is the normalized velocity difference that emerges naturally in the admissible scaling of the mean momentum equation. Self-similar equations for the scaled mean transverse flow V* and Reynolds shear stress $R^{*}_{uv}$=R uv /R uv,ref are derived from the mean continuity and mean momentum equations. Finally, approximate equations for V* and $R^{*}_{uv}$ are developed and found to agree well with experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Metal–ceramic composite structures for fabrication of high power density plasmonic devices

The recent decade brought many advances to plasmonics, but high power density plasmonic antennas designed to behave as heaters or operate in high temperature environments are still facing material stability challenges preventing their ultimate use. Gold has been the optimal choice among plasmonic materials but experiences morphology changes at temperature that result in device efficiency reduction and failure. Bulk titanium nitride has been explored as a solution but has deal-breaking tradeoffs in device quality factor. In this paper, we explore via proof-of-concept the use of a metal–ceramic composite structure to determine whether a bulk Au nanorod can provide strong plasmonic resonances while coated with an ultrathin conformal layer of titanium nitride or silica to provide morphological stability and sufficient plasmonic activity without excessive resonance quality degradation. We show SEM-level morphological stability for temperatures up to 500 °C with coatings below 4 nm. Computer modeling suggests the ultrathin titanium nitride has detrimental effects on the strong plasmonic resonances of a Au nanorod. Here, we then looked into other possible coatings for solutions to stabilize high power density plasmonic antennas including plasmonic oxides, metal adhesion layers, and silica, the latter appearing to be the best option while lowering the overall peak electric field intensity, the silica increases the electric field intensity at its boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗
DOE OSTI2023

Temperature effects in plasma-based positron acceleration schemes using electron filaments

Preserving the quality of a positron beam in a plasma-based accelerator, where a wakefield suitable for positron transport and acceleration is generated by means of an electron filament, is challenging. This is due to the nature of the wakefields, characterized by focusing fields that vary nonlinearly in the transverse direction, and by accelerating fields that are non-uniform. These fields also change slice-by-slice along the beam. Maintaining a high beam quality is pivotal for application of positron beams in a plasma-based collider. In this paper, we show that an initial background plasma temperature can help mitigate the positron beam quality degradation in plasma-based accelerators that rely on electron filaments. We show that temperature effects broaden the electron filament and smooth radially both the non-linear transverse and the non-uniform longitudinal wakefields. Using warm plasmas opens up new possibilities to improve beam quality in several plasma-based positron acceleration concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Investigation of boosted proton energies through proton radiography of target normal sheath acceleration fields in the multi-ps regime

Multi-kilojoule, multi-picosecond short-pulse lasers, such as the National Ignition Facility-Advanced Radiographic Capability laser and the OMEGA-Extended Performance laser, which have been constructed over the last two decades, enable exciting opportunities to produce high-brightness, high-energy laser-driven proton sources for applications in high-energy-density science like proton fast ignition for inertial fusion energy, particle radiography, and materials science studies. Results on these platforms have demonstrated enhanced accelerated proton energies and electron temperatures when compared to established scaling laws. Recent work has developed a new scaling for proton TNSA in the multi-ps regime. However, this new physics in the multi-ps regime motivates the need to understand the origin of the enhancement in proton energies. Toward this goal, here this work presents the first measurements of the TNSA accelerating sheath field in the multi-ps regime for pulse durations of 0.6, 5, and 10 ps. This measurement was achieved by using a separate TNSA proton source to radiograph the spatiotemporal profile of the accelerating sheath that is responsible for proton acceleration. The use of stacked radiochromic film detectors allows for a discrete time profile of the radiographs, thus enabling the measurement of the temporal and spatial evolution of the accelerating field. In performing this measurement, we extract quantities such as the sheath strength as a function of time and pulse duration, which shows that longer pulse durations sustain a stronger electric field for a longer duration when compared to sub-ps laser pulses, which may enable the observed boosted proton energies and proton conversion efficiencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

X-ray continuum spectroscopy of inertial confinement fusion implosions at the National Ignition Facility

A methodology for measuring x-ray continuum spectra of inertial confinement fusion (ICF) implosions is described. The method relies on the use of ConSpec, a high-throughput spectrometer using a highly annealed pyrolytic graphite crystal , which measures the spectra in the ≃20–30 keV range. Due to its conical shape, the crystal is sagittally focusing a Bragg-reflected x-ray spectrum into a line, which enhances the recorded x-ray emission signal above the high neutron-induced background accompanying ICF implosions at the National Ignition Facility. To improve the overall measurement accuracy, the sensitivity of the spectrometer measured in an off-line x-ray laboratory setting was revised. The error analysis was expanded to include the accuracy of the off-line measurements, the effect of the neutron-induced background, as well as the influence of possible errors in alignment of the instrument to the ICF target. Here we demonstrate how the improved methodology is applied in the analysis of ConSpec data with examples of a relatively low-neutron-yield implosion using a tritium–hydrogen–deuterium mix as a fuel and a high-yield deuterium–tritium (DT) implosion producing high level of the background. In both cases, the shape of the measured spectrum agrees with the exponentially decaying spectral shape of bremsstrahlung emission to within ±10%. In the case of the high-yield DT experiment, non-monotonic deviations slightly exceeding the measurement uncertainties are observed and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A nonlinear journey from structural phase transitions to quantum annealing

Motivated by an exact mapping between equilibrium properties of a one-dimensional chain of quantum Ising spins in a transverse field (the transverse field Ising (TFI) model) and a two-dimensional classical array of particles in double-well potentials (the “$\phi$ 4 model”) with weak inter-chain coupling, we explore connections between the driven variants of the two systems. Here, we argue that coupling between the fundamental topological solitary waves in the form of kinks between neighboring chains in the classical $\phi$ 4 system is the analog of the competing effect of the transverse field on spin flips in the quantum TFI model. As an example application, we mimic simplified measurement protocols in a closed quantum model system by studying the classical $\phi$ 4 model subjected to periodic perturbations. This reveals memory/loss of memory and coherence/decoherence regimes, whose quantum analogs are essential in annealing phenomena. In particular, we examine regimes where the topological excitations control the thermal equilibration following perturbations. This paves the way for further explorations of the analogy between lower-dimensional linear quantum and higher-dimensional classical nonlinear systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Inference of the neutron down-scatter ratio using neutron time-of-flight data for DT-layered implosions on the OMEGA laser

First analysis of the neutron time of flight (nTOF) data is presented to infer the areal density using the neutron down-scatter ratio (DSR) for OMEGA deuterium–tritium (DT)-layered (cryogenic) implosions. The required very-high-dynamic range (nTOF) signal is constructed using an nTOF detector [Forrest et al., Rev. Sci. Ins. 83, 10D919 (2012)] with multiple time-gated photomultiplier tubes. The nTOF data are analyzed using a forward-fit technique to account for the detector responses and infer the DSR and areal density for DT-layered implosions. The areal densities inferred using this analysis are found to be in agreement with the areal densities measured using the magnetic recoil spectrometer detector that is located in a similar line-of-sight. This work establishes the feasibility of the DSR measurement using nTOF data on OMEGA and provides an additional areal-density measurement that will benefit the assessment of the implosion performance and 3D reconstruction of the imploded core.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Many roads to the seam: How conformational flexibility drives nonadiabatic relaxation in a prototypical tetrapyrrolic chromophore

Large and structurally flexible chromophores pose challenges for in silico modeling of photodeactivation due to the many vibrational modes that can funnel the system toward energy degeneracy. In this work, we examine how the multiple degrees of freedom in biliverdin, a prototypical tetrapyrrolic chromophore, cooperate to drive access to the S 1 /S 0 intersection seam in vacuo. We begin by mapping the ground-state potential energy surface to identify representative biliverdin conformers relevant to photoexcitation. We then use a CASSCF-based framework to map the excited-state landscape and characterize the intersection seam, identifying distinct conical-intersection types. Finally, we employ ab initio multiple spawning to resolve the dynamical pathways by which the system accesses these regions. DFT potential-energy and free-energy mappings indicate that, although several conformers are relevant, the “locked-helix” ZsZsZs conformer predominates in the ground state. The intersection seam comprises numerous geometrically distinct regions characterized by varying degrees and combinations of dihedral torsion, pyramidalization, and bond-length alternation. Yet only select regions lie within energetic reach, and moderate barriers separate them from the S 1 minimum. Nonadiabatic dynamics combined with multivariate analyses show that, despite extensive mode coupling during deactivation that guides the system toward multiple regions of the seam, a single dihedral torsion, together with bond-length alternation, predominantly drives energy degeneracy. This work offers new insight into biliverdin’s intrinsic photochemical response and underscores a general feature of flexible chromophores: many modes may participate during photorelaxation, but only a limited subset ultimately dictates seam accessibility.

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

An idealized model for the crystal structure of intermetallic compounds isostructural with Mg 3 Cr 2 Al 18

Mg 3 Cr 2 Al 18 (abbreviated in this report as MCA) is the parent phase for a large class of intermetallic compounds that belong to the cubic crystal space group, $Fd\overline{3}m$. The purpose of this paper is to introduce an ideal, unrelaxed crystal structure for compounds isostructural with MCA. There are five distinct atomic sublattices in MCA compounds, which can be denoted, $A, B, C, D,$ and $E$. With this, a general description for MCA structures can be written as $A^{8a}_{1}B^{16c}_{2}C^{16d}_{2}D^{48f}_{6}E^{96g}_{12}$, where the superscripts represent the Wyckoff special equipoints associated with the various sublattices in MCA, and the subscripts indicate the contributions of each sublattice to the stoichiometry of one formula unit in an any given MCA structured compound. Sublattices D and E are where deviations from ideality occur in real, MCA-like compounds. This paper examines MCA bond lengths, nearest-neighbour polyhedral arrangements, 3-D sublattice crystal structures, 2-D atom tessellation patterns, and crystal chemical effects associated with atomic relaxations on the $D$ and $E$ sublattices. The ideal MCA crystal structure developed in this report provides an appropriate initial structure for use as input to crystal structure refinements of diffraction data for MCA-like phases being examined experimentally, or as input for computational, atomistic simulations of the structures of such compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Orbit-modulated transport and sources in time-dependent plasmas

The orbit-modulation framework for steady-state ion-orbit boundary fluxes is reformulated to apply to time-dependent plasmas. The new formulation can quantify the extent to which orbit fluxes are affected by the time dependence of the distribution function, and may be applied both to time-dependent physical scenarios and to computationally expensive simulations that are unable to reach a true steady state. As the system tends to a time-independent state, the orbit-loss fluxes approach their steady-state values, which are determined entirely by upstream transport (turbulent and collisional) and upstream sources and sinks. This result holds for the fluxes of any moments that are expressible in terms of the invariants that determine the orbits, including parallel toroidal angular momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗