Search NASA⌕ Search

SEARCH · Search NASA

Results for “ENERGY EQUATION”

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

Practical Effects of Integrating Temperature with Strang Split Reactions

For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A whole-core steady-state thermal-hydraulic model for annular fuel type fluoride-salt-cooled reactors

A whole-core, steady-state thermal-hydraulic model is developed for the fluoride-salt-cooled small modular advanced high-temperature reactor (SmAHTR) that employs an annular fuel configuration. This pre-conceptual reactor design by Oak Ridge National Laboratory (ORNL) has the annular fuel and moderator pins arranged in a hexagonal layout. The FLiBe coolant flows from the bottom to the top of the core, parallel to the hexagonal bundle. The fuel and moderator pins in the core are discretized into finite volumes and the 3-D heat conduction equation is solved to obtain the temperature profile. Inter-fuel assembly conduction is also addressed. For this fuel assembly configuration, the coolant flows through two distinct regions – the hexagonal pin bundle and the annulus between the fuel pin and the tie rod. The fluid flow through the hexagonal bundles is modeled using the subchannel approach, in which the coolant region is discretized into corner, edge and interior subchannels and the resulting conservation equations are systematically solved. The 1-D mass, momentum and energy equations are solved for the annulus channels between the fuel pin and the tie rod. Pertinent closure models from the literature are employed to close the system of equations. We also performed a preliminary code-to-code comparison between the present model and a CFD model.. The resulting thermal-hydraulic model can provide temperature, flow rate and pressure drop profiles for the different solid and fluid regions throughout the entire core. Whole-core thermal-hydraulic results for a representative power profile are presented and discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

ANTS

The ANTS code (Alternate Non-Linear Two-phase Solver) is based on a novel non-linear solution algorithm for the solution of the two-phase, subchannel fluid equations. It achieves its performance through decoupling of the two-phase momentum equations (axial and transverse) from the axial phasic mass and energy equations which allows for a nested non-linear iteration scheme. This enables a plane-by-plane solution which the inner iteration focuses on a reduced non-linear equation set for the primitives in phasic mass flow rate, enthalpy and void for each node edge. Single node edges are coupled as part of the outer iteration via surface mass fluxes which appear as source terms in the inner iteration scheme. The outer iteration readily accommodates two-phase flow phenomena closure relationships for subchannel mixing and void drift. A primary feature is the use of a non-staggered mesh computational mesh and steady-state iterative solver in contrast to all existing subchannel codes.

Kropaczek, David J↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Quantifying scalar field dynamics with DESI 2024 Y1 BAO measurements

Quintessence scalar fields are a natural candidate for evolving dark energy. Unlike the phenomenological w 0 ⁢w a parameterization of the dark energy equation of state, they cannot accommodate the phantom regime of dark energy w⁡(z) <–1, or crossings into the phantom regime. Recent baryon acoustic oscillation (BAO) measurements by the Dark Energy Spectroscopic Instrument (DESI) indicate a preference for evolving dark energy over a cosmological constant, ranging from 2.6⁢σ –3.9⁢σ when fitting to w 0⁢ w a , and combining the DESI BAO measurements with other cosmological probes. In this work, we directly fit three simple scalar field models to the DESI BAO data, combined with cosmic microwave background anisotropy measurements and supernova datasets. We find the best fit model to include a 2–4% kinetic scalar field energy Ω scf,k , for a canonical scalar field with a quadratic or linear potential. However, only the DESY-Y5 supernova dataset combination shows a preference for quintessence over Λ cold dark matter (CDM) at the 95% confidence level. Fitting to the supernova datasets Pantheon, Pantheon+, DES-Y5, and Union3, we show that the mild tension (n σ < 3.4) under Λ ⁢CDM emerges from a BAO preference for smaller values of fractional mass-energy density Ω m < 0.29, while all supernova datasets, except for Pantheon, prefer larger values, Ω m > 0.3. The tension under Λ ⁢CDM remains noticeable (n σ < 2.8), when replacing two of the DESI BAO redshift bins with effective redshifts z eff = 0.51, and z eff = 0.706 with comparable BOSS DR 12 BAO measurements at z eff = 0.51, and z eff = 0.61. Furthermore, canonical scalar fields as dark energy are successful in mitigating that tension.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2020

Wake Statistics of Different-Scale Wind Turbines under Turbulent Boundary Layer Inflow

Subscale wind turbines can be installed in the field for the development of wind technologies, for which the blade aerodynamics can be designed in a way similar to that of a full-scale wind turbine. However, it is not clear whether the wake of a subscale turbine, which is located closer to the ground and faces different incoming turbulence, is also similar to that of a full-scale wind turbine. In this work we investigate the wakes from a full-scale wind turbine of rotor diameter 80 m and a subscale wind turbine of rotor diameter of 27 m using large-eddy simulation with the turbine blades and nacelle modeled using actuator surface models. The blade aerodynamics of the two turbines are the same. In the simulations, the two turbines also face the same turbulent boundary inflows. The computed results show differences between the two turbines for both velocity deficits and turbine-added turbulence kinetic energy. Such differences are further analyzed by examining the mean kinetic energy equation.

17 WIND ENERGY↗
DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Sensitivity of Arctic Surface Temperature to Including a Comprehensive Ocean Interior Reflectance to the Ocean Surface Albedo Within the Fully Coupled CESM2

Abstract Almost all current climate models simplify the ocean surface albedo (OSA) by assuming the reflected solar energy without the ocean interior contribution. In this study, an improved ocean surface albedo scheme is incorporated into the Community Earth System Model version 2 (CESM2) to assess the sensitivity of Arctic surface temperature to including ocean interior reflectance to the OSA. Fully coupled CESM2 simulations with and without ocean interior reflectance are subsequently performed, we focus on the analysis of Arctic surface temperature responses. Incorporating ocean interior reflectance increases absorbed solar radiation and warms the ocean, enhancing seasonal heat storage and release across the Arctic Ocean, and increasing sea ice reduction and positive climate feedbacks that elevates Arctic surface temperature. Seasonal variations in air‐surface temperature differences induce changes in turbulent heat flux patterns, concurrently modifying dynamic advection and moisture processes that affect boundary layer humidity and low clouds, especially in winter. Based on partitioning physical processes in the thermodynamic energy equation, surface air warming is induced primarily through positive heating anomalies of vertical advection, latent heat release, and longwave radiative forcing. Through an examination of the surface energy budget, skin temperature warming is driven predominantly by increased downward longwave radiation, positive surface albedo feedback in summer, and increased conductive heat transport from the ocean particularly in winter. Significant effects of ocean interior reflectance on the Arctic Ocean, including sea surface warming and sea ice reduction, justify the importance of ocean interior reflectance in climate models for better understanding of ongoing Arctic climate changes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Reconciling and Improving Formulations for Thermodynamics and Conservation Principles in Earth System Models (ESMs)

Abstract This paper provides a comprehensive derivation of the total energy equations for the atmospheric components of Earth System Models (ESMs). The assumptions and approximations made in this derivation are motivated and discussed. In particular, it is emphasized that closing the energy budget is conceptually challenging and hard to achieve in practice without resorting to ad hoc fixers. As a concrete example, the energy budget terms are diagnosed in a realistic climate simulation using a global atmosphere model. The largest total energy errors in this example are spurious dynamical core energy dissipation, thermodynamic inconsistencies (e.g., coupling parameterizations with the host model) and missing processes/terms associated with falling precipitation and evaporation (e.g., enthalpy flux between components). The latter two errors are not, in general, reduced by increasing horizontal resolution. They are due to incomplete thermodynamic and dynamic formulations. Future research directions are proposed to reconcile and improve thermodynamics formulations and conservation principles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Semi-classical approximation description of static properties of nuclei

In this report, we present recent development in semi-classical description of static properties of nuclei. By employing the Wigner transform, we derive simple semiclassical approximations for evaluating properties of nuclear system. We present results of calculations of static properties, using the energy density functional associated with Skyrme effective nucleon–nucleon interaction. In particular, we consider properties of the Wigner distribution function, line of beta stability, Coulomb and symmetry energies, deformation energies, equation of state, nuclear radii, and single-particle level density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The physical origin of dark energy constraints from rubin observatory and CMB-S4 lensing tomography

ABSTRACT We seek to clarify the origin of constraints on the dark energy equation of state parameter from CMB lensing tomography, that is the combination of galaxy clustering and the cross-correlation of galaxies with CMB lensing in a number of redshift bins. We focus on the analytic understanding of the origin of the constraints. Dark energy information in these data arises from the influence of three primary relationships: distance as a function of redshift (geometry), the amplitude of the power spectrum as a function of redshift (growth), and the power spectrum as a function of wavenumber (shape). We find that the effects from geometry and growth play a significant role and partially cancel each other out, while the shape effect is unimportant. We also show that Dark Energy Task Force figure of merit forecasts from the combination of LSST galaxies and CMB-S4 lensing are comparable to the forecasts from cosmic shear in the absence of the CMB lensing map, thus providing an important independent check. Compared to the forecasts with the LSST galaxies alone, combining CMB lensing and LSST clustering information increases the FoM by roughly a factor of 3–4 in the optimistic scenario where systematics are fully under control. We caution that achieving these forecasts will likely require a full analysis of higher-order biasing, photometric redshift uncertainties, and stringent control of other systematic limitations, which are outside the scope of this work, whose primary purpose is to elucidate the physical origin of the constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Testing cosmology with double source lensing

Double source lensing provides a dimensionless ratio of distance ratios, a "remote viewing" of cosmology through distances relative to the gravitational lens, beyond the observer. We use this to test the cosmological framework, particularly with respect to spatial curvature and the distance duality relation. We derive a consistency equation for constant spatial curvature, allowing not only the investigation of flat vs curved but of the Friedmann-Lemaître-Robertson-Walker framework itself. For distance duality, we demonstrate that the evolution of the lens mass profile slope must be controlled to ≳ 5 times tighter fractional precision than a claimed distance duality violation. Furthermore, using LensPop forecasts of double source lensing systems in Euclid and LSST surveys we also explore constraints on dark energy equation of state parameters and any evolution of the lens mass profile slope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Connecting primordial gravitational waves and dark energy

Cosmic acceleration manifested in the early universe as inflation, generating primordial gravitational waves detectable in the cosmic microwave background (CMB) radiation. Cosmic acceleration is occurring again at present as dark energy, detectable in cosmic distance and structure surveys. We explore the intriguing idea of connecting the two occurrences through quintessential inflation by an α-attractor potential without a cosmological constant. For this model we demonstrate robustness of the connection 1 + w 0 ≈ 4/(3N 2 r) between the present day dark energy equation of state parameter w 0 and the primordial tensor to scalar ratio r for a wide range of initial conditions. Analytic and numerical solutions produce current thawing behavior, resulting in a tight relation w a ≈ -1.53(1 + w 0 )≈ -0.2 (4 × 10 -3 /r). Upcoming CMB and galaxy redshift surveys can test this consistency condition. In conclusion, within this model, lack of detection of a dark energy deviation from Λ predicts a higher r, and lack of detection of r predicts greater dark energy dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Comparison of Artificial Compressibility Methods for Coupled Laminar Fluid Flow and Heat Transfer

Artificial compressibility (AC) methods for incompressible flows perturb the mass continuity equation by the time derivative of the pressure field to achieve direct coupling between pressure and velocity without the need for a projection method. Several variations of the original AC method have been proposed (e.g., the entropically damped artificial compressibility [EDAC] method) that can, in some scenarios, improve on the original model. In this paper, we methodically compare the results of the original AC model and the EDAC model for both steady-state and transient problems with coupling between the momentum and energy equations. We demonstrate that the EDAC model and the AC model yield similar numerical results, and the EDAC model is 10%-20% faster than the AC model for a given problem.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗
DOE OSTI2023

Dark energy with a triplet of classical U(1) fields

Here, we present a new mechanism for cosmic acceleration consisting of a scalar field coupled to a triplet of classical U(1) gauge fields. The gauge fields are arranged in a homogeneous, isotropic configuration, with both electric- and magneticlike vacuum expectation values. The gauge fields provide a masslike term via a Chern-Simons interaction that suspends the scalar away from its potential minimum, thereby enabling potential-dominated evolution. We show this mechanism can drive a brief period of acceleration, such as dark energy, without the need for fine-tunings. We obtain simple analytic results for the dark energy equation of state and dependence on model parameters. In this model, the presence of the gauge field generically leads to a suppression of long-wavelength gravitational waves, with implications for the experimental search for cosmic microwave background B modes and direct detection of a stochastic gravitational wave background.

79 ASTRONOMY AND ASTROPHYSICS↗