Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling and design of compact, permanent-magnet transport systems for highly divergent, broad energy spread laser-driven proton beams

Laser-driven (LD) ion acceleration has been explored in a newly constructed short focal length laser beamline at the BELLA petawatt facility (interaction point 2, iP2). For applications utilizing such LD ion beams, a beam transport system is required, which for reasons of compactness be ideally contained within 3 m. While they are generated from a micron-scale source, large divergence and energy spread of LD ion beams present a unique challenge to transporting them compared to beams from conventional accelerators. This study gives an overview of proposed compact transport designs using permanent magnets satisfying different requirements depending on the application for the iP2 laser beamline such as radiation biology, material science, and high-energy density science. These designs are optimized for different parameters such as energy spread and peak proton density according to the application’s need. The various designs consist solely of permanent magnet elements, which can provide high magnetic field gradients on a small footprint. While the field strengths are fixed, we have shown that the beam size is able to be tuned effectively by varying the placement of the magnets. The performance of each design was evaluated based on high-order particle tracking simulations of typical LD proton beams. We also examine the ability of certain configurations to tune and select beam energies, critical for specific applications. A more detailed investigation was carried out for a design to deliver 10 MeV LD accelerated ions for radiation biology applications. With these transport system designs, the iP2 laser beamline is ready to house various application experiments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

First-principles electronic structure in second-moment calculation of mode frequencies: Failure of quasiharmonic approximation in silicon

Using first-principles electronic structure calculations (local density approximation), we present second-moment calculations of the temperature dependence of the vibrational normal mode frequencies in Si up to 1500 K . The method is based on simple ensemble averages of displacements and forces and is easy to implement. These efficient calculations incorporate all interactions nonperturbatively, implicitly including all high order force constants, in contrast to typical perturbative approaches. Also in this work, we propose and apply to our classical calculations an estimated quantum correction that could be used for classical approaches other than the moments method. We compare our results to available experiments, including recent inelastic neutron scattering data and conclude that the second-moment gives the frequency shift to within a factor of 2 in contrast to quasiharmonic approaches that fail by an order of magnitude to predict the anharmonic shift. Noting the relative simplicity and speed of the calculations, we conclude that the second-moment approximation reasonably accounts for anharmonicity in a system that has been very closely studied experimentally and may reliably be applied in an exploratory way to more complex systems. Finally, we suggest that further improvement on the results could be obtained by incorporating the fourth moment, which can be implemented with little additional cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First-principles predictions of temperature-dependent infrared dielectric function of polar materials by including four-phonon scattering and phonon frequency shift

Recently, first-principles calculations based on density functional theory have been widely used to predict the temperature-dependent infrared spectrum of polar materials, but the calculations are usually limited to the harmonic frequency (0 K) and three-phonon scattering damping for the zone-center infrared-active optical phonon modes, and fail to predict the high-temperature infrared optical properties of materials such as sapphire (α–Al 2 O 3 ), GaAs, TiO 2 , etc., due to the neglect of high-order phonon scattering damping and phonon frequency shift. In this work, we implemented first-principles calculations to predict the temperature-dependent infrared dielectric function of polar materials by including four-phonon scattering and phonon frequency shift. The temperature-dependent phonon damping by including three- and four-phonon scattering as well as the phonon frequency shift by including cubic and quartic anharmonicity and the thermal expansion effect are calculated based on anharmonic lattice dynamics method. The infrared dielectric function of α–Al 2 O 3 is parameterized, and then the temperature-dependent infrared optical reflectance is determined. We find that our predictions agree better with the experimental data than the previous density functional theory-based methods. Lastly, this work will help to effectively predict the thermal radiative properties of polar materials at elevated temperature, which is generally difficult to measure, and will enable predictive design of new materials for radiative applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Large spin splittings due to the orbital degree of freedom and spin textures in a ferroelectric nitride perovskite

First-principles simulations are conducted to predict that ferroelectric nitride perovskite LaWN 3 not only exhibits large spin splittings (2.7 eV Å) but also possesses unique spin textures for some of its conduction levels. Such spin splittings around the Γ and L points cannot be interpreted as a typical mixture of Rashba or Dresselhaus configurations but rather require the development of four-band $\textit{k ∙ p}$ models with high-order terms. We further identify that, for some bands, spin splittings can be greatly contributed by the pure orbital degree of freedom (PODF), a unique character of our four-band Hamiltonian compared to the traditional two-band version. The concept of PODF-enhanced spin splittings paves a way for designing materials with large spin splittings. Moreover, the energy levels possessing such large splittings and complex spin textures can be brought close to the conduction-band minimum by applying epitaxial strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Composite fermion mass: Experimental measurements in ultrahigh quality two-dimensional electron systems

Composite fermions (CFs), exotic quasiparticles formed by pairing an electron and an even number of magnetic flux quanta, emerge at high magnetic fields in an interacting electron system, and can explain phenomena such as the fractional quantum Hall state (FQHS) and other many-body phases. CFs possess an effective mass (m CF ) whose magnitude is inversely related to the most fundamental property of a FQHS, namely its energy gap. Here we present here experimental measurements of m CF in ultrahigh quality two-dimensional electron systems confined to GaAs quantum wells of varying thickness. An advantage of measuring m CF over gap measurements is that mass values are insensitive to disorder and are therefore ideal for comparison with theoretical calculations, especially for high-order FQHS. Our data reveal that m CF increases with increasing well width, reflecting a decrease in the energy gap as the electron layer becomes thicker and the in-plane Coulomb energy softens. Comparing our measured masses with available theoretical results, we find significant quantitative discrepancies, highlighting that more rigorous and accurate calculations are needed to explain the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Convergence of nuclear effective field theory with perturbative pions

The classic paper by Fleming, Mehen, and Stewart (FMS) cast doubts on the convergence of spin-triplet nucleon-nucleon partial wave scattering amplitudes when following the proposal of Kaplan, Savage, and Wise to construct nuclear effective field theory (EFT) around the unitary fermion limit with perturbative pion exchange. FMS identified the subclass of iterated one-pion exchange potential graphs as the cause of this poor convergence, which they showed persisted in the chiral limit. Theoretical tools are developed in this work to compute these Feynman graphs analytically to high order in all angular momentum channels simultaneously, examining the amplitudes computed to seven loops in the $\textit{L = J}$ channels and three loops in the coupled $\textit{L = J}$ ± 1 channels. One finds that there is nothing pathological about the perturbative expansion of a 1/$r^3$ potential in general and that the expansion converges satisfactorily in all partial waves except those with the lowest angular momentum, particularly the $^3P_0$ and the coupled $^3S_1 - ^3D_1$ channels. The results corroborate work by Birse, which suggests possible avenues to explore for improving the range of validity of the EFT expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory (χ EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χ EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. Furthermore, we then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium χ EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

X 17 boson and the H 3 ( p , e + e – ) He 4 and He 3 ( n , e + e – ) He 4 processes: A theoretical analysis

The present work deals with e + – e – pair production in the four-nucleon system. We first analyze the process as a purely electromagnetic one in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory ( χ EFT ) . Next, we examine how the exchange of a hypothetical low-mass boson would impact the cross section for such a process. We consider several possibilities, that this boson is either a scalar, pseudoscalar, vector, or axial particle. The ab initio calculations use exact hyperspherical-harmonics methods to describe the bound state and low-energy spectrum of the A = 4 continuum, and they fully account for initial state interaction effects in the 3 + 1 clusters. While electromagnetic interactions are treated to high orders in the chiral expansion, the interactions of the hypothetical boson with nucleons are modeled in leading-order χ EFT (albeit, in some instances, selected subleading contributions are also accounted for). We also provide an overview of possible future experiments probing pair production in the A = 4 system at a number of candidate facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Magnetic structure of few-nucleon systems at high momentum transfers in a chiral effective field theory approach

The five low-energy constants (LECs) in the electromagnetic current derived in chiral effective field theory ( χ EFT ) up to one loop are determined by a simultaneous fit to the A = 2 – 3 nuclei magnetic moments and to the deuteron magnetic form factor and threshold electrodisintegration at backward angles over a wide range of momentum transfers. The resulting parametrization then yields predictions for the 3 He/ 3 H magnetic form factors in excellent accord with the experimental values for momentum transfers ranging up to ≈0.8 GeV/c, beyond the expected regime of validity of the χEFT approach. The calculations are based on last-generation two-nucleon interactions including high orders in the chiral expansion and derived by Entem, Machleidt, and Nosyk [Phys. Rev. C 96, 024004 (2017)] and by Piarulli et al. [Phys. Rev. C 94, 054007 (2016)], using different χEFT formulations. In the A = 3 calculations, (chiral) three-nucleon interactions are also accounted for. The model dependence resulting from these different formulations of the interactions is found to be mild for momentum transfer below ≈0.8 GeV/c . Further, an analysis of the convergence of the chiral expansion is also provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multiparticle integral and differential correlation functions

This paper formalizes the use of integral and differential cumulants for measurements of multiparticle event-by-event transverse momentum fluctuations, rapidity fluctuations, as well as net-charge fluctuations. This enables the introduction of multiparticle balance functions, defined based on differential correlation functions (factorial cumulants), that suppress two- and three-prong resonance decays effects and enable measurements of underlying long-range correlations obeying quantum number conservation constraints. These multiparticle balance functions satisfy simple sum rules determined by quantum number conservation. It is additionally shown that these multiparticle balance functions arise as an intrinsic component of high-order net-charge cumulants. This implies that the magnitude of these cumulants, measured in a specific experimental acceptance, is strictly constrained by charge conservation and primarily determined by the rapidity and momentum width of these balance functions. The paper also presents techniques to reduce the computation time of differential correlation functions up to order n = 10 based on the methods of moments. Published by the American Physical Society 2024

Physics↗
DOE OSTI2021

Deep learning simulations of the microwave sky

Here we present 500 high-resolution, full-sky millimeter-wave deep learning (DL) simulations that include lensed CMB maps and correlated foreground components. We find that these MillimeterDL simulations can reproduce a wide range of non-Gaussian summary statistics matching the input training simulations, while only being optimized to match the power spectra. The procedure we develop in this work enables the capability to mass produce independent full-sky realizations from a single expensive full-sky simulation, when ordinarily the latter would not provide enough training data. We also circumvent a common limitation of high-resolution DL simulations that they be confined to small sky areas, often due to memory or GPU issues; we do this by developing a “stitching” procedure that can recover the large-scale, high-order statistics and avoid discontinuities or repeated features in the maps. In addition, since our network takes as input a full-sky lensing convergence map, it can in principle take a full-sky lensing convergence map from any large-scale structure (LSS) simulation and generate the corresponding lensed CMB and correlated foreground components at millimeter wavelengths; this is especially useful in the current era of combining results from both CMB and LSS surveys, which require a common set of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Simulating magnetized neutron stars with discontinuous Galerkin methods

Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod /ΛΠ N limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinuous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized Tolman-Oppenheimer-Volkoff stars without any fine-tuning of parameters. Among the limiters explored in the paper, the only limiting strategy we can endorse is a discontinuous Galerkin-finite-difference hybrid method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Krylov spaces for truncated spectrum methodologies

We propose herein an extension of truncated spectrum methodologies, a nonperturbative numerical approach able to elucidate the low energy properties of quantum field theories. TSMs, in their various flavors, involve a division of a computational Hilbert space, H , into two parts, one part, H 1 that is “kept” for the numerical computations, and one part, H 2 , that is discarded or “truncated.” Even though H 2 is discarded, truncated spectrum methodologies will often try to incorporate the effects of H 2 in some effective way. In these terms, we propose to keep the dimension of H 1 small. We pair this choice of H 1 with a Krylov subspace iterative approach able to take into account the effects of H 2 . This iterative approach can be taken to arbitrarily high order and so offers the ability to compute quantities to arbitrary precision. In many cases it also offers the advantage of not needing an explicit UV cutoff. To compute the matrix elements that arise in the Krylov iterations, we employ a Feynman diagrammatic representation that is then evaluated with Monte Carlo techniques. Each order of the Krylov iteration is variational and is guaranteed to improve upon the previous iteration. The first Krylov iteration is akin to the next-to-leading order approach of Elias-Miró [NLO renormalization in the Hamiltonian truncation, ]. To demonstrate this approach, we focus on the ( 1 + 1 d )-dimensional ϕ 4 model and compute the bulk energy and mass gaps in both the Z 2 -broken and unbroken sectors. We estimate the critical ϕ 4 coupling in the broken phase to be g c = 0.2645 ± 0.002 . Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Renormalons and power corrections in pseudo- and quasi-GPDs

High-order behavior of the perturbative expansion for short-distance observables in QCD is intimately related to the contributions of small momenta in the corresponding Feynman diagrams and this correspondence provides one with a useful tool to investigate power-suppressed nonperturbative corrections. We use this technique to study the structure of power corrections to parton quasi- and pseudo-GPDs which are used in lattice calculations of generalized parton distributions. As the main result, we predict the functional dependence of the leading power corrections to quasi(pseudo)-GPDs on x variable for nonzero skewedness parameter ξ . The kinematic point x = ± ξ turns out to be special. We find that the nonperturbative corrections to quasi-GPDs at this point are suppressed by the first power of the hard scale only. These contributions come from soft momenta and have nothing to do with the known UV renormalon in the Wilson line. We also show that power corrections can be strongly suppressed by the normalization procedure. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Velocity circulation intermittency in finite-temperature turbulent superfluid helium

We study intermittency of circulation moments in turbulent superfluid helium by using experimental grid turbulence and numerical simulations of the Hall-Vinen-Bekarevich-Khalatnikov model. More precisely, we compute the velocity circulation Γ r in loops of size r laying in the inertial range. For both experimental and numerical data, the circulation variance shows a clear Kolmogorov scaling < $Γ$$^2_r$ > ~ r 8/3 in the inertial range, independently of the temperature. Scaling exponents of high-order moments are comparable, within error bars, to previously reported anomalous circulation exponents in classical turbulence and low-temperature quantum turbulence numerical simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Pressure fields produced by single-bubble collapse near a corner

Damage produced by repeated bubble collapse to neighboring rigid objects in hydraulic systems is an important consequence of cavitation. Although bubble collapse near a single wall has received significant attention in the past, few studies exist on the dynamics of bubbles collapsing near a corner, i.e., two flat rigid surfaces intersecting at a right angle. In this work we quantify the pressure fields produced by a single bubble collapsing near two perpendicular rigid walls. Using a high-order accurate shock- and interface-capturing method to solve the three-dimensional compressible Navier-Stokes equations for gas and liquid flows, we simulate the dynamics of a single bubble collapsing at different initial stand-off distances from the two walls. In contrast to a bubble collapsing near a single wall, the collapse of bubbles within a critical stand-off distance is not symmetric about the bisecting plane due to the interaction between the bubble and the second wall. The second wall affects the pressure produced during the collapse in the following ways: (i) For bubbles initially located sufficiently close to both walls, the reentrant jet produced during collapse no longer points in the direction normal to the closest wall but at an angle toward the corner, (ii) the part of the emitted shock with the highest amplitude propagates in line with the jet, and (iii) the bubble migrates in that same direction during its collapse with a dependence on the stand-off distance, consistent with predictions made using Kelvin impulse. The location of maximum pressure along the walls is measured for the different initial stand-off distances. Using acoustic arguments, we find a semiempirical relationship to predict the initial stand-off distances for which the maximum pressure occurs in the corner. Here we find that when the bubble is sufficiently close to equidistant from each boundary, the maximum pressure is observed in the corner due to the water-hammer and implosion shocks reflecting off the boundaries and intersecting in the corner. We also show that when bubbles are initially attached to either wall the wall pressure produced can be significantly increased compared to bubbles detached from either wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗