SEARCH · Engineering Papers

Results for “consistent scheme”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Photochemically induced acousto-optics in gases

Acousto-optics consists of launching acoustic waves in a medium (usually a crystal) in order to modulate its refractive index and create a tunable optical grating. Here, in this article, we present the theoretical basis of an alternative scheme to generate acousto-optics in a gas, where the acoustic waves are initiated by the localized absorption (and thus gas heating) of spatially modulated UV light, as was demonstrated by Michine and Yoneda [Commun. Phys. 3, 24 (2020)]. We identify the chemical reactions initiated by the absorption of UV light via the photodissociation of ozone molecules present in the gas, and calculate the resulting temperature increase in the gas as a function of space and time. Solving the Euler fluid equations shows that the modulated, isochoric heating initiates a mixed acoustic-entropy wave in the gas, whose high-amplitude density (and thus refractive index) modulation can be used to manipulate a high-power laser. We calculate that diffraction efficiencies near 100% can be obtained using only a few millimeters of gas containing a few percent ozone fraction at room temperature, with UV fluences of less than 100 mJ/cm 2 —consistent with the experimental measurements. Our analysis suggests possible ways to optimize the diffraction efficiency by changing the buffer gas composition. Gases have optics damage thresholds 2–3 orders of magnitude beyond those of solids; these optical elements should therefore be able to manipulate kilojoule-class lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Self-interaction-corrected Kohn–Sham effective potentials using the density-consistent effective potential method

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn–Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. Here, we adapt and implement the density-consistent effective potential method of Kohut, Ryabinkin, and Staroverov (2014) to obtain SI-corrected local effective potentials from the SI-corrected Fermi–Löwdin orbitals and density in the Fermi–Löwdin orbital self-interaction correction scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO–LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Dark matter in a mirror solution to the strong C P problem

We study thermal production of dark matter (DM) in a realization of the minimal models of Bonnefoy [], where parity is used to solve the strong C P problem by transforming the entire Standard Model (SM) into a mirror copy. Although the mirror electron e ′ is a good DM candidate, its viability is mired by the presence of the mirror up quark u ′ , whose abundance is intimately related to the e ′ abundance and must be suppressed. This can be achieved through a sequential freeze-in mechanism, where mirror photons are first produced from SM gluons, and then the mirror photons produce e ′ . After computing the details of this double freeze-in, we discuss the allowed parameter space of the model, which lies at the threshold of experimental observations. We find that this origin of e ′ DM requires a low reheating temperature after inflation and is consistent with the baryon asymmetry arising from leptogenesis, providing mirror neutrinos have a significant degeneracy. Finally, we show that this e ′ DM is not compatible with Higgs parity, the simplest scheme with exact parity, unless SM parameters deviate significantly from their central values or the minimal model is extended. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Bayesian learning of orthogonal embeddings for multi-fidelity Gaussian Processes

Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. Here in this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The \projection" mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction

42 ENGINEERING↗

DOE OSTI2024

Gauge preheating with full general relativity

Abstract We study gauge preheating following pseudoscalar-driven inflation in full general relativity. We implement the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) scheme to solve the full nonlinear evolution of the metric alongside the dynamics of the pseudoscalar and gauge fields. The dynamics of the background and emission of gravitational waves are broadly consistent with simulations in a Friedmann-Lemaître-Robertson-Walker (FLRW) spacetime. We find large, localized overdensities in the BSSN simulations of orderδ=δρ/ρ∼ 30, and the dimensionless power spectrum ofδpeaks above unity. These overdense regions are seeded on length scales only slightly smaller than the horizon, and have a compactnessC∼ 0.1. The scale of peak compactness is shorter than the Jeans length, which implies that pressure of the matter fields plays an important role in the evolution of these objects.

Astronomy & Astrophysics↗

DOE OSTI2023

Nonperturbatively renormalized nucleon gluon momentum fraction in the continuum limit of N f = 2 + 1 + 1 lattice QCD

We present the nonperturbatively renormalized nucleon gluon momentum fraction using ensembles with 2+1+1 flavors of highly improved staggered quarks (HISQ), generated by the MILC Collaboration. The calculation is done using clover fermions for the valence action with three pion masses, 220, 310, and 690 MeV, and three lattice spacings, 0.09, 0.12, and 0.15 fm. The renormalization is done using RI/MOM nonperturbative renormalization and using cluster-decomposition error reduction (CDER) to enhance the signal-to-noise ratio of the renormalization constant. We find the CDER technique is particularly important to improve the signal at the finer lattice ensembles where the lattice volume is larger. We extrapolate the gluon momentum fraction to the continuum-physical limit and obtain $\langle{x}\rangle$ g = 0.502⁢(53) stat+NPR ⁢(50) mixing in the $\overline{MS}$ scheme at 2 GeV, where first error includes the statistical error and uncertainties in nonperturbative renormalization, while the latter systematic error accounts for ignoring quark mixing. Our gluon momentum fraction is consistent with other recent lattice-QCD results at physical pion mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Soils' dirty little secret: Depth–based comparisons can be inadequate for quantifying changes in soil organic carbon and other mineral soil properties

Quantifying changes in soil organic carbon (SOC) stocks and other soil properties is essential for understanding how soils will respond to land management practices and global change. Although they are widely used, comparisons of SOC stocks at fixed depth (FD) intervals are subject to errors when changes in bulk density or soil organic matter occur. The equivalent soil mass (ESM) method has been recommended in lieu of FD for assessing changes in SOC stocks in mineral soils, but ESM remains underutilized for SOC stocks and has rarely been used for other soil properties. In this paper, we draw attention to the limitations of the FD method and demonstrate the advantages of the ESM approach. We provide illustrations to show that the FD approach is susceptible to errors not only for quantifying SOC stocks but also for soil mass-based properties such as SOC mass percent, C:N mass ratio, and δ 13 C. Here we describe the ESM approach and show how it mitigates the FD method limitations. Using bulk density change simulations applied to an empirical dataset from bioenergy cropping systems, we show that the ESM method provides consistently lower errors than FD when quantifying changes in SOC stocks and other soil properties. To simplify the use of ESM, we detail how the method can be integrated into sampling schemes, and we provide an example R computer script that can perform ESM calculations on large datasets. We encourage future studies, whether temporal or comparative, to utilize sampling methods that are amenable to the ESM approach. Overall, we agree with previous recommendations that ESM should be the standard method for evaluating SOC stock changes in mineral soils, but we further suggest that ESM may also be preferred for comparisons of other soil properties including mass percentages, elemental mass ratios, and stable isotope composition.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

A description of multi-material hydrodynamics without reliance upon Pressure-Temperature-Equilibrium and its relevance to Eulerian codes [Slides]

This report contains "a description of the full multi-material Euler system without reliance on the assumption of the PTE closure; a discussion of the details that a numerical scheme discretizing the above system might need to have; a description of how that might play out in a PTE code after all; and, an attempt to lay out the maximally consistent set of equations."

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Cyclometallated Co(III) Complexes with Lowest-Energy Charge Transfer Excited States Accessible with Visible Light

The Co(III) complexes, cis-[Co(ppy) 2 (L)]PF 6 , where ppy = 2-phenylpyridine and L = bpy (2,2′-bipyridine; 1), phen (1,10-phenanthroline; 2), and DAP (1,12-diazaperylene; 3), are reported and their photophysical properties were investigated to evaluate their potential as sensitizers for applications that include solar energy conversion schemes and photoredox catalysis. Calculations show that cyclometallation in the cis-[Co(ppy) 2 (L)]PF 6 series affords strong Co(dπ)/ppy(π) orbital interactions that result in a Co/ppy(π*) highest occupied molecular orbital (HOMO) and a lowest unoccupied molecular orbital (LUMO) localized on the diimine ligand, L(π*). Complexes 1–3 exhibit relatively invariant oxidation potentials, whereas the reduction event is dependent on the identity of the diimine ligand, L, consistent with the theoretical predictions. For 3 a broad Co/ppy(π*) → L(π*) metal/ligand-to-ligand charge transfer (ML-LCT) absorption band is observed in CH 3 CN with a maxima at 507 nm, extending beyond 600 nm. Upon excitation of the 1ML-LCT transition, transient absorption features consistent with the population of a 3 ML-LCT excited state with lifetimes, τ, of 3.0 ps, 4.6 and 42 ps for 1, 2 and 3 in CH 3 CN respectively are observed. Upon irradiation with 505 nm, 3 is able to reduce methyl viologen (MV 2+ ), an electron acceptor commonly in photocatalytic schemes. Here, to our knowledge, 3 represents the first heteroleptic molecular Co(III) complex that combines cyclometallation with a diimine ligand with lowest-lying metal-to-ligand charge transfer excited states able to undergo photoinduced charge transfer with low-energy green light. As such, the structural design of 3 represents an important step toward d6 photosensitizers based on earth abundant metals.

Energy↗

DOE OSTI2020

Focusing a round coherent beam by spatial filtering the horizontal source

This paper illustrates the use of spatial filtering with a horizontal slit near the source to enlarge the horizontal coherence in an experimental station and produce a diffraction-limited round focus at an insertion device beamline for X-ray photon correlation spectroscopy experiments. Simple expressions are provided to guide the optical layout, and wave propagation simulations confirm their applicability. The two-dimensional focusing performance of Be compound refractive lenses to produce a round diffraction-limited focus at 11 keV capable of generating a high-contrast speckle pattern of an aerogel sample is demonstrated. The coherent scattering patterns have comparable speckle sizes in both horizontal and vertical directions. The focal spot sizes are consistent with hybrid ray-tracing calculations. Producing a two-dimensional focus on the sample can be helpful to resolve speckle patterns with modern pixel array detectors with high visibility. This scheme has now been in use since 2019 for the 8-ID beamline at the Advanced Photon Source, sharing the undulator beam with two separate beamlines, 8-ID-E and 8-ID-I at 7.35 keV, with increased partially coherent flux, reduced horizontal spot sizes on samples, and good speckle contrast.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

A Computationally Efficient Algorithm for Computing Convex Hull Prices

Electricity markets worldwide allow participants to bid non-convex production offers. While non-convex offers can more accurately reflect a resource's capabilities, they create challenges for market clearing processes. For example, system operators may execute side payments when a participant’s cost is not covered through energy sale settlements from locational marginal pricing schemes, or when a participant incurs lost opportunity costs to follow the dispatch signal. Convex hull pricing minimizes these and other types of side payments while providing uniform (i.e., locationally and temporally consistent) prices. However, computing convex hull prices involves solving either a large-scale linear program - which in turn requires explicit descriptions of market participants’ convex hulls -or the Lagrangian dual of the corresponding non-convex scheduling problem. Here, we propose a computationally feasible and industrially scalable Benders decomposition approach to computing convex hull prices at least an order of magnitude faster than the current state-of-the-art while leveraging recent advances in convex hull formulations for thermal generating units.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2024

Conceptual Design of a Thermal–to-14-MeV Neutron Conversion Device for Use in the Advanced Test Reactor

Here, this paper details the conceptual design of a thermal-to-14-MeV neutron converter consisting of a mixture of lithium and deuterium in a blanket material. Such a device operates a two-step reaction, first generating tritons via thermal neutron absorption in the tritium breeding material, and in the second step, high-energy neutrons are produced either via deuterium-tritium fusion reaction or with tritium reacting with lithium. A thermal-to-14-MeV neutron converter significantly hardens the neutron spectrum by virtually removing thermal neutrons and adding a high-energy 14 MeV component to the neutron spectrum. While similar concepts have been previously proposed and tested in other reactors, the unique characteristics of the Advanced Test Reactor (ATR), namely its important thermal flux (up to $10^{15}$ n.cm$^{-2}$.s$^{-1}$) makes it markedly attractive for obtaining a very large fast neutron flux, usable for irradiation studies under neutron flux conditions prototypical of fusion reactors. The paper provides a description of a new computational scheme developed for handling the coupled neutron-triton transport mechanism using the Geant4 toolkit. Resulting neutron spectra and high-energy neutron yields are summarized for different irradiation positions and potential neutron breeder materials. Maximum predicted thermal-to-14-MeV neutron yields are on the order of $2.10^{-4}$, which is consistent with previous studies found in the literature. Thus, when placed inside the ATR, such a neutron converter will be providing the largest high-energy neutron source available for activation and irradiation studies of materials foreseen for use in fusion reactors. Future steps will involve qualifying the computational scheme using the ATR critical facility using activation foil measurements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method

Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.

42 ENGINEERING↗

DOE OSTI2023

Generation of entangled-photons by a quantum dot cascade source in polarized cavities: Using cavity resonances to boost signals and preserve the entanglements

Motivated by recent advances in the development of single photon emitters for quantum information sciences, here we design and formulate a quantum cascade model that describes cascade emission by a quantum dot (QD) in a cavity structure while preserving entanglement that stores information needed for single photon emission. The theoretical approach is based on a photonic structure that consists of two orthogonal cavities in which resonance with either the first or second of the two emitted photons is possible, leading to amplification and rerouting of the entangled light. The cavity–QD scheme uses a four-level cascade emitter that involves three levels for each polarization, leading to two spatially entangled photons for each polarization. By solving the Schrodinger equation, we identify the characteristic properties of the system, which can be used in conjunction with optimization techniques to achieve the “best” design relative to a set of prioritized criteria or constraints in our optical system. The theoretical investigations include an analysis of emission spectra in addition to the joint spectral density profile, and the results demonstrate the ability of the cavities to act as frequency filters for the photons that make up the entanglements and to modify entanglement properties. The results provide new opportunities for the experimental design and engineering of on-demand single photon sources.

Chemistry↗

DOE OSTI2024

Perturbatively including inhomogeneities in axion inflation

Axion inflation, i.e. an axion-like inflaton coupled to an Abelian gauge field through a Chern-Simons interaction, comes with a rich and testable phenomenology. This is particularly true in the strong backreaction regime, where the gauge field production heavily impacts the axion dynamics. Lattice simulations have recently demonstrated the importance of accounting for inhomogeneities of the axion field in this regime. We propose a perturbative scheme to account for these inhomogeneities while maintaining high computational efficiency. Our goal is to accurately capture deviations from the homogeneous axion field approximation within the perturbative regime as well as self-consistently determine the onset of the nonperturbative regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: cosmological implications from multitracer BAO analysis with galaxies and voids

We construct cosmic void catalogues with the dive void finder upon SDSS BOSS DR12 and eBOSS DR16 galaxy samples with BAO reconstruction, and perform a joint BAO analysis using different types of galaxies and the corresponding voids. The BAO peak is evident for the galaxy–galaxy, galaxy–void, and void–void correlation functions of all data sets, including the ones cross-correlating LRG and ELG samples. Two multitracer BAO fitting schemes are tested, one combining the galaxy and void correlation functions with a weight applied to voids, and the other using a single BAO dilation parameter for all clustering measurements. Both methods produce consistent results with mock catalogues, and on average ~10 per cent improvements of the BAO statistical uncertainties are observed for all samples, compared to the results from galaxies alone. By combining the clustering of galaxies and voids, the uncertainties of BAO measurements from the SDSS data are reduced by 5–15 per cent, yielding 0.9 per cent, 0.8 per cent, 1.1 per cent, 2.3 per cent, and 2.9 per cent constraints on the distance $D_{_{\rm V}}(z)$, at effective redshifts 0.38, 0.51, 0.70, 0.77, and 0.85, respectively. When combined with BAO measurements from SDSS MGS, QSO, and Lyα samples, as well as the BBN results, we obtain $H_0 = 67.58 \pm 0.91\, {\rm km}\, {\rm s}^{-1}\, {\rm Mpc}^{-1}$, Ω m = 0.290 ± 0.015, and $\Omega _\Lambda h^2 = 0.3241 \pm 0.0079$ in the flat-ΛCDM framework, where the 1σ uncertainties are around 6 per cent, 6 per cent, and 17 per cent smaller respectively, compared to constraints from the corresponding anisotropic BAO measurements without voids and LRG–ELG cross-correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Generalized parton distributions from symbolic regression

AI/ML informed Symbolic Regression is the next stage of scientific modeling. We utilize a highly customizable symbolic regression package “PySR” to model the x and t dependence of the flavor isovector combination Hu-d(x, t, ) at = 0. These PySR models were trained on GPD results provided by both Lattice QCD and phenomenological sources GGL, GK, and VGG. We demonstrate, for the first time, the consistency and systematic convergence of Symbolic Regression by quantifying the disparate models through their Taylor expansion coefficients. In addition to PySR penalizing models with higher complexity and mean-squared error, we implement schemes that test specific physics hypotheses, including force-factorized x and t dependence and Regge behavior in PySR GPDs. We show that PySR can identify factorizing GPD sources based on their response to the Force-Factorized model. Here, knowing the precise behavior of the GPDs, and their uncertainties in a wide range in x and t, crucially impacts our ability to concretely and quantitatively predict hadronic spatial distributions and their derived quantities.

Dotson, Andrew [New Mexico State University, Las C↗