Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First Order Second Moment”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Deep Direct-Use Feasibility Study Numerical Modeling and Uncertainty Analysis using iTOUGH2 for West Virginia University

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. As a result, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided to use single permeability model for further analysis; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the perdition, therefore, is considered in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis are used for economic analysis. The dataset submitted here support the described study. Manuscript is submitted to Geothermics, will be linked once paper is accepted.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2020

Modeling study of deep direct use geothermal on the West Virginia university campus-morgantown, WV

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use (DDU) system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. Consequently, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided only the single permeability model is needed; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the prediction, therefore, heterogeneity is included in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis will be used for further economic studies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

A surrogate-model-based approach for estimating the first and second-order moments of offshore wind power

Power curve, the functional relationship that governs the process of converting a set of weather variables experienced by a wind turbine into electric power, is widely used in the wind industry to estimate power output for planning and operational purposes. Existing methods for power curve estimation have three main limitations: (i) they mostly rely on wind speed as the sole input, thus ignoring the secondary, yet possibly significant effects of other environmental factors, (ii) they largely overlook the complex marine environment in which offshore turbines operate, potentially compromising their value in offshore wind energy applications, and (ii) they solely focus on the first-order properties of wind power, with little (or null) information about the variation around the mean behavior, which is important for ensuring reliable grid integration, asset health monitoring, and energy storage, among others. In light of that, this study investigates the impact of several wind-and wave-related factors on offshore wind power variability, with the ultimate goal of accurately predicting its first two moments. Further, our approach couples OpenFAST—a multi-physics wind turbine simulator—with Gaussian Process (GP) regression to reveal the underlying relationships governing offshore weather-to-power conversion. We first find that a multi-input power curve which captures the combined impact of wind speed, direction, and air density, can provide double-digit improvements, in terms of prediction accuracy, relative to univariate methods which rely on wind speed as the sole explanatory variable (e.g. the standard method of bins). Wave-related variables are found not important for predicting the average power output, but interestingly, appear to be extremely relevant in describing the fluctuation of the offshore power around its mean. Tested on real-world data collected at the New York/New Jersey bight, our proposed multi-input models demonstrate a high explanatory power in predicting the first two moments of offshore wind generation, testifying their potential value to the offshore wind industry.

17 WIND ENERGY↗
DOE OSTI2020

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Lagrangian analysis for turbulent transport in variable-density turbulence

In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Second Order Closures for the Radiative Transfer Equation: Some Are Unstable

The largest existing simulations of cosmic reionization model radiative transfer with moment methods that require a closure relation. The two most commonly used closure relations are M1 and OTVET; both close the moment hierarchy at the first moment. We explore the properties of a higher, second-order closure. We show that direct generalizations of M1 and OTVET to one higher order are physically unstable - i.e., the closure equations themselves result in unstable solutions, not just their numerical implementation. In fact, a generalization of OTVET to any order higher than the first one is unstable. We are also able to show that any local (i.e., depending only on the local moments of the radiation field, like M1) second-order closure that depends only on the radiation intensity and radiation flux, but does not explicitly depend on the radiation pressure, is physically unstable. This result restricts the choice of possible second-order closure relations.

Gnedin, Nickolay Y. [Fermilab; Chicago U., Astron.↗
DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: High order kinematic and dynamic relations

Dark matter, if it exists, accounts for five times as much as ordinary baryonic matter. To better understand the self-gravitating collisionless dark matter flow on different scales, a statistical theory involving kinematic and dynamic relations must be developed for different types of flow, e.g., incompressible, constant divergence, and irrotational flow. This is mathematically challenging because of the intrinsic complexity of dark matter flow and the lack of a self-closed description of flow velocity. Here, this paper extends our previous work on second-order statistics Xu to kinematic relations of any order for any type of flow. Dynamic relations were also developed to relate statistical measures of different orders. The results were validated by N-body simulations. On large scales, we found that (i) third-order velocity correlations can be related to density correlation or pairwise velocity; (ii) the pth-order velocity correlations follow ∝ a (p+2)/2 for odd p and ∝ a p/2 for even p, where a is the scale factor; (iii) the overdensity δ is proportional to density correlation on the same scale, $\langle$δ$\rangle$∝$\langle$δδ'$\rangle$; (iv) velocity dispersion on a given scale r is proportional to the overdensity on the same scale. On small scales, (i) a self-closed velocity evolution is developed by decomposing the velocity into motion in haloes and motion of haloes; (ii) the evolution of vorticity and enstrophy are derived from the evolution of velocity; (iii) dynamic relations are derived to relate second- and third-order correlations; (iv) while the first moment of pairwise velocity follows $\langle$Δu L $\rangle$=-Har (H is the Hubble parameter), the third moment follows $\langle$(Δu L ) 3 $\rangle$ ∝ ε u ar that can be directly compared with simulations and observations, where ε u ≈ 10 -7 m 2 /s 3 is the constant rate for energy cascade; (v) the pth order velocity correlations follow ∝ a (3p-5)/4 for odd p and ∝ a 3p/4 for even p. Finally, the combined kinematic and dynamic relations lead to exponential and one-fourth power-law velocity correlations on large and small scales, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Field moment expansion method for interacting bosonic systems

In this work, we introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2021

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole–dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

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

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Anisotropic magnetism and electronic structure of trigonal EuAl 2 Ge 2 single crystals

Understanding the interplay between magnetic and electronic degrees of freedom is of profound recent interest in different Eu-based magnetic topological materials. Here, we studied the magnetic and electronic properties of the layered Zintl-phase compound EuAl 2 Ge 2 crystallizing in the trigonal CaAl 2 Si 2 -type structure. We report zero-field neutron diffraction, temperature T - and magnetic-field H-dependent magnetic susceptibility χ(T, H), isothermal magnetization M(T, H), heat capacity C p (T, H), and electrical resistivity ρ(T, H) measurements, together with T-dependent angle-resolved photoemission spectroscopy (ARPES) measurements complemented with first-principle calculations. EuAl 2 Ge 2 undergoes second-order A-type antiferromagnetic (AFM) ordering below T N = 27.5(5) K, with the Eu moments (Eu 2+ , S = 7/2) aligned ferromagnetically in the ab plane while these layers are stacked antiferromagnetically along the c axis. The critical fields at which all moments become parallel to the field are 37.5(5) and 52.5(5) kOe for H ∥ ab and H ∥ c, respectively. The H = 0 magnetic structure consists of trigonal AFM domains associated with ab-plane magnetic anisotropy and a field-induced reorientation of the Eu spins in the domains is also evident at T = 2 K below the critical field H c1 = 2.5(1) kOe. The ρ(T) measurements reveal metallic behavior transforming into a slight resistivity increase on cooling towards T N . A pronounced loss of spin-disorder scattering is observed below T N . The ARPES results show that EuAl 2 Ge 2 is metallic both above and below T N , and the Fermi surface is anisotropic with two hole pockets at the zone center and one small electron pocket at each M point. In the AFM phase, we directly observe folded bands in ARPES due to the doubling of the magnetic unit cell along the c axis with an enhancement of quasiparticle weight due to the complex change in the coupling between the magnetic moments and itinerant electrons on cooling below T N . The observed electronic structure is well reproduced by first-principle calculations, which also predict the presence of nontrivial electronic states near the Fermi level in the AFM phase with Z 2 topological numbers 1;(000)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-order antiferromagnetic transitions of SrMn 2 P 2 and CaMn 2 P 2 single crystals containing corrugated-honeycomb Mn sublattices

Significance With rare exceptions, an antiferromagnetic (AFM) transition in zero magnetic field is thermodynamically of second order where the thermal-average magnetic moments of the magnetic atoms (ordered moments) vary continuously on cooling through the AFM ordering temperature T N with no latent heat at the transition. Such materials include the AFM pnictides CaMn 2 As 2 , SrMn 2 As 2 , CaMn 2 Sb 2 , SrMn 2 Sb 2 , and CaMn 2 Bi 2 . Here we demonstrate that the closely related SrMn 2 P 2 and CaMn 2 P 2 insulators instead exhibit first-order AFM transitions at T N = 53 and 70 K, respectively, where the heat capacity exhibits a latent heat at T N . The mechanism causing these first-order transitions remains to be explained, but its understanding may lead to the development of novel magnetic materials of technological interest.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Vacancy-tuned magnetism in LaMn x Sb 2

The layered AMP 2 (A = alkali-earth or rare-earth atom, M = transition metal, P = Sb, Bi) compounds are widely studied for their rich magnetism and electronic structure topology. We provide a detailed characterization of the magnetic and transport properties of LaMn x Sb 2 , an understudied member of the AMP 2 family. LaMn x Sb 2 forms with intrinsic Mn vacancies, and we demonstrate that by varying the starting ratio of La, Mn, and Sb, we can synthetically control the Mn occupancy and produce single crystals with x = 0.74 – 0.97. Magnetization and transport measurements indicate LaMn x Sb 2 has a rich temperature-composition (T–x) magnetic phase diagram with physical properties strongly influenced by the Mn occupancy. LaMn x Sb 2 orders antiferromagnetically at T 1 = 130–180 K, where T 1 increases with x. Below T 1 , the T–x phase diagram is complicated. At high x, there is a second transition T 2 that decreases in temperature as x is lowered, vanishing below x ≤ 0.85. A third, first-order, transition T 3 is detected at x ≈ 0.92, and the transition temperature increases as x is lowered, crossing above T 2 near x ≈ 0.9. On moving below x < 0.79, here we find the crystal structure changes from the P 4/nmm arrangement to an I$\bar{4}$2m structure with partially ordered Mn vacancies. The change in crystal structure results in the sudden appearance of two new low-temperature phases and a crossover between regimes of negative and positive magnetoresistance when x ≤ 0.78. Finally, we provide powder neutron diffraction for x = 0.93, and find that the high-x compositions first adopt a G-type antiferromagnetic structure with the Mn moments aligned within the ab plane, which is followed upon further cooling by a second transition to a different, noncollinear structure where the moments are rotated within the basal plane. Our results demonstrate that LaMn x Sb 2 is a highly tunable material with six unique magnetically ordered phases, depending on T and x.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A dissipation-induced superradiant transition in a strontium cavity-QED system

Driven-dissipative many-body systems are ubiquitous in nature and a fundamental resource for quantum technologies. However, they are also complex and hard to model because they cannot be described by the standard tools in equilibrium statistical mechanics. Probing nonequilibrium critical phenomena in pristine setups can illuminate fresh perspectives on these systems. Here, we use an ensemble of cold 88 Sr atoms coupled to a driven high-finesse cavity to study the cooperative resonance fluorescence (CRF) model, a classic driven-dissipative model describing coherently driven dipoles superradiantly emitting light. We observe its nonequilibrium phase diagram characterized by a second-order phase transition. Below a critical drive strength, the atoms quickly reach the so-called superradiant steady state featuring a macroscopic dipole moment; above the critical point, the atoms undergo persistent Rabi-like oscillations. At longer times, spontaneous emission transforms the second-order transition into a discontinuous first-order transition. Our observations pave the way for harnessing robust entangled states and exploring boundary time crystals in driven-dissipative systems.

Science & Technology - Other Topics↗
DOE OSTI2021

Measurements of multiplicity fluctuations of identified hadrons in inelastic proton–proton interactions at the CERN Super Proton Synchrotron

Measurements of multiplicity fluctuations of identified hadrons produced in inelastic p+p interactions at 31, 40, 80, and 158 $\text {Ge}\text {V}/c$ beam momentum are presented. Three different measures of multiplicity fluctuations are used: the scaled variance $\omega $ and strongly intensive measures $\Sigma $ and $\Delta $. These fluctuation measures involve second and first moments of joint multiplicity distributions. Data analysis is preformed using the Identity method which corrects for incomplete particle identification. Strongly intensive quantities are calculated in order to allow for a direct comparison to corresponding results on nucleus–nucleus collisions. The results for different hadron types are shown as a function of collision energy. A comparison with predictions of string-resonance Monte-Carlo models: Epos, Smash and Venus, is also presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗