Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density perturbation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Control of plasmons in topological insulators via local perturbations

Here, we use a fully quantum mechanical approach to demonstrate control of plasmonic excitations in prototype models of topological insulators by molecule-scale perturbations. Strongly localized surface plasmons are present in the host systems, arising from the topologically nontrivial single-particle edge states. A numerical evaluation of the random phase approximation equations for the perturbed systems reveals how the positions and the internal electronic structure of the added molecules affect the degeneracy of the locally confined collective excitations, i.e., shifting the plasmonic energies of the host system and changing their spatial charge density profile. In particular, we identify conditions under which significant charge transfer from the host system to the added molecules occurs. Furthermore, the induced field energy density in the perturbed topological systems due to external electric fields is determined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Monte Carlo event generator for initial conditions of conserved charges in nuclear geometry

At top collider energies where baryon stopping is negligible, the initial state of heavy-ion collisions is overall charge neutral and predominantly composed of gluons. Nevertheless, there can also be significant local fluctuations of the baryon number, strangeness, and electric charge densities about zero, perturbatively corresponding to the production of quark/antiquark pairs. These previously ignored local charge fluctuations can permit the study of charge diffusion in the quark-gluon plasma (QGP), even at top collider energies. In this paper we present a new model denoted ICCING (initial conserved charges in nuclear geometry) which can reconstruct the initial conditions of conserved charges in the QGP by sampling a (g→q¯q) splitting probability over the initial energy density. We find that the new charge distributions generally differ from the bulk energy density; in particular, the strangeness distribution is significantly more eccentric than standard bulk observables and appears to be associated with the geometry of hot spots in the initial state. The new information provided by these conserved charges opens the door to studying a wealth of new charge- and flavor-dependent correlations in the initial state and ultimately the charge transport parameters of the QGP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Experimental and numerical evaluation of laser shinethrough in polystyrene at 3 ω and 4 ω laser frequencies

In laser-direct-drive inertial confinement fusion (ICF), laser shinethrough is the process associated with laser propagation through the target plastic ablator, polystyrene is often used in ICF experiments, and leads to energy deposition in the target until the critical density is established at the ablator front surface. This process causes material modifications that may change the shock propagation and lead to density and pressure perturbations at the interface between the polystyrene and cryogenic deuterium-tritium (DT) layer. This work reports on experimental observations and numerical investigations of laser shinethrough at the third and the fourth harmonic of ICF class laser systems (1 ω corresponds to a wavelength of ∼1.06 μm ). The experiments evidence that, at 4 ω , the laser is absorbed over a distance of about a few microns and laser shinethrough is completely eliminated in contrast with the behavior at 3 ω . These experimental observations are reproduced by numerical simulations. The laser fluence required to reach the critical density at the ablator front surface is also found to be more than one order of magnitude higher at the third harmonic. Finally, modeling suggests that the back surface of the ablator exhibits ion temperatures greater than the melting temperature of cryogenic DT and electron pressure larger than the bulk and shear moduli of cryogenic DT for the case of using a 3 ω laser.

Pineau, A. [Univ. of Rochester, NY (United States)↗
DOE OSTI2023

Mallat Scattering Transformation based surrogate for Magnetohydrodynamics

Abstract A Machine and Deep Learning (MLDL) methodology is developed and applied to give a high fidelity, fast surrogate for 2D resistive MagnetoHydroDynamic (MHD) simulations of Magnetic Liner Inertial Fusion (MagLIF) implosions. The resistive MHD code is used to generate an ensemble of implosions with different liner aspect ratios, initial gas preheat temperatures (that is, different adiabats), and different liner perturbations. The liner density and magnetic field as functions of x , y , and z were generated. The Mallat Scattering Transformation (MST) is taken of the logarithm of both fields and a Principal Components Analysis (PCA) is done on the logarithm of the MST of both fields. The fields are projected onto the PCA vectors and a small number of these PCA vector components are kept. Singular Value Decompositions of the cross correlation of the input parameters to the output logarithm of the MST of the fields, and of the cross correlation of the SVD vector components to the PCA vector components are done. This allows the identification of the PCA vectors vis-a-vis the input parameters. Finally, a Multi Layer Perceptron (MLP) neural network with ReLU activation and a simple three layer encoder/decoder architecture is trained on this dataset to predict the PCA vector components of the fields as a function of time. Details of the implosion, stagnation, and the disassembly are well captured. Examination of the PCA vectors and a permutation importance analysis of the MLP show definitive evidence of an inverse turbulent cascade into a dipole emergent behavior. The orientation of the dipole is set by the initial liner perturbation. The analysis is repeated with a version of the MST which includes phase, called Wavelet Phase Harmonics (WPH). While WPH do not give the physical insight of the MST, they can and are inverted to give field configurations as a function of time, including field-to-field correlations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Spinor $GW$ Bethe-Salpeter calculations in BerkeleyGW: Implementation, symmetries, benchmarking, and performance

Computing the GW quasiparticle band structure and Bethe-Salpeter equation (BSE) absorption spectra for materials with spin-orbit coupling have commonly been done by treating GW corrections and spin-orbit coupling (SOC) as separate perturbations to density-functional theory. However, accurate treatment of materials with strong spin-orbit coupling (such as many topological materials of recent interest, and thermoelectrics) often requires a nonperturbative approach using spinor wave functions in the Kohn-Sham equation and GW/BSE. Such calculations have only recently become available, in particular for the BSE. Here, we have implemented this approach in the plane-wave pseudopotential GW/BSE code BerkeleyGW, which is highly parallelized and widely used in the electronic-structure community. We present reference results for quasiparticle band structures and optical absorption spectra of solids with different strengths of spin-orbit coupling, including Si, Ge, GaAs, GaSb, CdSe, Au, and Bi 2 Se 3 . The calculated quasiparticle band gaps of these systems are found to agree with experiment to within a few tens of meV. SOC splittings are found to be generally in better agreement with experiment, including quasiparticle corrections to band energies. The absorption spectrum of GaAs is not significantly impacted by the inclusion of spin-orbit coupling due to its relatively small value (0.2 eV) in the Λ direction, while the absorption spectrum of GaSb calculated with the spinor GW/BSE captures the large spin-orbit splitting of peaks in the spectrum. For the prototypical topological insulator Bi 2 Se 3 , we find a drastic change in the low-energy band structure compared to that of DFT, with the spinorial treatment of the GW approximation correctly capturing the parabolic nature of the valence and conduction bands after including off-diagonal self-energy matrix elements. We present the detailed methodology, approach to spatial symmetries for spinors, comparison against other codes, and performance compared to spinless GW/BSE calculations and perturbative approaches to SOC. This work aims to spur further development of spinor GW/BSE methodology in excited-state research software and enables a more accurate and detailed exploration of electronic and optical properties of materials containing elements with large atomic numbers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

First principles thermal transport modeling in GaN and related materials

Gallium nitride is a wide bandgap material utilized in a variety of technologies, including high-power electronics and light-emitting diodes, partly due to its favorable thermal properties. This chapter describes modern first-principles-based modeling of phonons and lattice thermal conductivity (k) of GaN, III-nitrides and related materials. In particular, we describe the theoretical underpinnings of calculating phonon dispersions, intrinsic phonon interactions, and other lattice dynamical properties from quantum perturbation theory and density functional theory (DFT) methods. Description of how these methods are then coupled with the Peierls-Boltzmann transport (PBT) equation to determine phonon distributions and lifetimes relevant for thermal transport is given. These theoretical and numerical methods have demonstrated quantitative accuracy and predictive power for calculating pristine k and defect-limited k from first principles for a variety of materials. We present a review of the literature utilizing DFT-PBT methods to understand novel k behaviors in III-nitrides and related materials.

Lindsay, Lucas↗
DOE OSTI2022

Ion dynamics driven by a strongly nonlinear plasma wake

Abstract In plasma wakefield accelerators, the wave excited in the plasma eventually breaks and leaves behind slowly changing fields and currents that perturb the ion density background. We study this process numerically using the example of a Facility for Advanced aCcelerator Experimental Tests (FACET) experiment where the wave is excited by an electron bunch in the bubble regime in a radially bounded plasma. Four physical effects underlie the dynamics of ions: (1) attraction of ions toward the axis by the fields of the driver and the wave, resulting in formation of a density peak, (2) generation of ion-acoustic solitons following the decay of the density peak, (3) positive plasma charging after wave breaking, leading to acceleration of some ions in the radial direction, and (4) plasma pinching by the current generated during the wave-breaking. The interplay of these effects results in the formation of various radial density profiles, which are difficult to produce in any other way.

Physics↗
DOE OSTI2025

Azimuthal and axial structures in 3D particle-in-cell simulation of Penning discharge

We report the results of a 3D particle-in-cell simulation of cylindrical Penning discharge with the axial magnetic field in the so-called reflex configuration, where the cathode and anti-cathode are biased to the same negative potential. The discharge is supported by thermal electron emission from the cathode. Electron and ion collisions, including ionization, are fully accounted for. The emphasis is on a specific regime in which the plasma potential at the center of the discharge is positive with respect to the chamber walls, serving as an anode. Spatial and temporal scales of the observed azimuthal and axial fluctuations and structures are characterized. It is suggested that azimuthal structures are caused by the dissipative gradient-drift instability. We find that the axial fluctuations related to the plasma-beam instabilities are weakly correlated with the azimuthal perturbations of the density, so that the azimuthal modes rotate as a whole and do not show any axial shear. The behavior of the electric potential is more involved, demonstrating intermittent standing wave and propagating structures in the axial direction that modulate the electron transport, producing the standing wave pattern (along the z-direction) in the radial electron flux.

3D PIC simulation of discharge↗
DOE OSTI2023

Preliminary Monte Carlo and Thermal Hydraulic Analysis using a Hybrid ETF-Corrected-Diffusion Prediction Block

This paper builds upon previous work to accelerate the Picard iteration (PI) method typically applied for coupled Monte Carlo-Thermal hydraulic (MC-TH) solutions. Previously, the use of the generalized transfer functions (GTFs) to predict variation in macroscopic cross sections following a perturbation in TH properties was demonstrated for a subset of simple 3D problems. In addition, the reduced-order transport prediction block relied on the first order perturbation (FOP) method, which was shown to have computational overheads. Recent work replaced the FOP block with a 1-group nodal diffusion solver to eliminate these overheads. While the use of diffusion is desirable for large-scale problems, the new solver introduces significant homogenization error. This work aims to address this issue by using the Jacobian-Free Newton Krylov (JFNK) method to generate a set of super homogenization (SPH) factors to improve the accuracy of the diffusion solution. The SPH factors will be used in conjunction with an improved cross section prediction method – the expanded transfer function (ETF) method – to produce a highly accurate flux prediction for an axial 1D boiling water reactor (BWR) pincell following a large perturbation in moderator density. The ETF-corrected diffusion (ETF-CD) block is shown to be highly accurate for the 1D test case. Future work will investigate the accuracy of the method for a realistic 3D pressurized water reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Electron Energization by Inertial Alfvén Waves in Density Depleted Flux Tubes at Jupiter

Juno satellite observations have illustrated substantial examples of broadband electron energization up to 10 5 –10 6 eV levels. In order to explain these observations, we use a hybrid gyrofluid kinetic-electron model in an untilted dipolar topology to illustrate energization to high levels in weak current conditions by inertial Alfvén waves close to the Jupiter ionosphere for ambient plasma densities and magnetic field perturbations inferred from Juno satellite observations. The key to the high energization is the extremely low densities evident in the observations which necessitates the acceleration of electrons to very high velocities in order to carry the field-aligned current.

58 GEOSCIENCES↗
DOE OSTI2021

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
DOE OSTI2022

Investigations of plasma response associated with resonant magnetic perturbation fields using perturbation method in KSTAR H-mode plasmas

Abstract The plasma response associated with the resonant magnetic perturbation (RMP) field was investigated using the small edge perturbations induced by a modulated supersonic molecular beam injection (SMBI) in KSTAR. The modulated SMBI provides a time-varying perturbation of the plasma density source in the region just inside the last closed flux surface and a modulated flow damping rate. Radial propagation of the toroidal rotation perturbation induced by SMBI from the q = 3 surface to the q = 2 surface was observed. Theoretical analysis using the general perturbed equilibrium code of the RMP intensity profiles of the RMP field is consistent with the phase profile of the toroidal rotation perturbation.

Physics↗
DOE OSTI2021

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Seeding the explosion of a high-current-density conductor in a controlled manner through the addition of micron-scale surface defects

Inhomogeneities in a current-carrying conductor promote non-uniform heating and expansion through the complex feedback between current density, electrical resistivity, Ohmic heating, temperature, and hydrodynamics. Three-dimensional-magnetohydrodynamic (3D-MHD) simulations suggest that μm-scale resistive inclusions or voids seed local overheating and through hydrodynamic explosion generate continuously growing craters which become several times larger than the initial perturbation. The ejected mass is the genesis of an electrothermally driven plasma filament which develops at lower current than plasmas on uniform surfaces adjacent to the defect. This result suggests that 1D or even 2D treatments are largely inadequate for detailed prediction of plasma formation. To test computational predictions, z-pinch experiments driven to 1 MA studied ultra-high-purity aluminum rods which were then machined to include pairs of quasi-hemispherical voids or “engineered defects (ED)” on the current-carrying surface. ED are the dominant current-density perturbation and reproducibly drive local overheating which can be compared with 3D-MHD simulation. Data from high-resolution-gated imagers of visible surface emissions confirm many simulation predictions, including the surface topography of local overheating, and the propensity for neighboring ED to prematurely source plasmas which then connect to form a plasma filament. Results also provide conditional support of theory which suggests heating similarity; that is, heating is independent of ED size for geometrically scaled ED.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wiggler-based electron beam microbunching amplifier study

The wiggler-based amplifier of the microbunching instability in the electron beam is considered and compared to a longitudinal space charge amplifier consisting of the drifts and chicanes. The gain in the magnitude of a seed perturbation of the electron beam density is evaluated analytically as a function of the period of the perturbation, the electron beam parameters, the wiggler parameters, and the amplifier lattice parameters. The theory is compared with the results of the numerical calculations of the gain based on simulations of the electron beam transport through the amplifier lattice. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Microscopic constraints for the equation of state and structure of neutron stars: A Bayesian model mixing framework

Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒⁢EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20⁢𝑛 0 . The uncertainties of the 𝜒⁢EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2⁢𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.

Bayesian methods↗