Engineering Papers⌕ Search

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.

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

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Testing of Cloud Microphysics Scheme with Snow Events

One of the grand challenges of the Global Precipitation Measurement (GPM) mission is to improve precipitation measurements in mid- and high-latitudes during cold seasons through the use of high-frequency passive microwave radiometry. For this, the Weather Research Forecast (WRF) model with the Goddard microphysics scheme is coupled with the Satellite Data Simulation Unit (WRF-SDSU) that has been developed to facilitate the over-snowfall retrieval algorithm by providing virtual cloud library and microwave brightness temperature (To) measurements consistent to the GPM Microwave Imager (GMI). This study tested the Goddard cloud microphysics scheme in WRF in snowstorm events (January 20-22, 2007) over the Canadian CloudSat/CALIPSO Validation Project (C3VP) site up in Ontario, Canada. In this meeting, we will present the performance of the Goddard cloud microphysics scheme both with 2ice (ice and snow) and 3ice (ice, snow and graupel) as well as other WRF microphysics schemes. Results will be compared with the King Radar data. We will also use the WRF model outputs to drive the Goddard SDSU to calculate radiances and backscattering signals consistent to satellite direct observations. These simulated radiance are evaluated against the measurement from A-Train satellites.

Tao, Wei-Kuo↗
NASA NTRS2016

The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure-temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

planets and satellites↗
DOE OSTI2025

The Principle of Maximum Conformality Correctly Resolves the Renormalization-Scheme-Dependence Problem

In this paper, we clarify a serious misinterpretation and consequent misuse of the Principle of Maximum Conformality (PMC), which also can serve as a mini-review of PMC. In a recently published article, P. M. Stevenson has claimed that “the PMC is ineffective and does nothing to resolve the renormalization-scheme-dependence problem”, concluding incorrectly that the success of PMC predictions is due to the PMC being a “laborious, ad hoc, and back-door” version of the Principle of Minimal Sensitivity (PMS). We show that such conclusions are incorrect, deriving from a misinterpretation of the PMC and an overestimation of the applicability of the PMS. The purpose of the PMC is to achieve precise fixed-order pQCD predictions, free from conventional renormalization schemes and scale ambiguities. We demonstrate that the PMC predictions satisfy all the self-consistency conditions of the renormalization group and standard renormalization-group invariance; the PMC predictions are thus independent of any initial choice of renormalization scheme and scale. The scheme independence of the PMC is also ensured by commensurate scale relations, which relate different observables to each other. Moreover, in the Abelian limit, the PMC dovetails into the well-known Gell-Mann–Low framework, a method universally revered for its precision in QED calculations. Due to the elimination of factorially divergent renormalon terms, the PMC series not only attains a convergence behavior far superior to that of its conventional counterparts but also deftly curtails any residual scale dependence caused by the unknown higher-order terms. This refined convergence, coupled with its robust suppression of residual uncertainties, furnishes a sound and reliable foundation for estimating the contributions from unknown higher-order terms. Anchored in the bedrock of standard renormalization-group invariance, the PMC simultaneously eradicates the factorial divergences and eliminates superfluous systematic errors, which inversely provides a good foundation for achieving high-precision pQCD predictions. Consequently, owing to its rigorous theoretical underpinnings, the PMC is eminently applicable to virtually all high-energy hadronic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2008

The a(3) Scheme--A Fourth-Order Space-Time Flux-Conserving and Neutrally Stable CESE Solver

The CESE development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To initiate a systematic CESE development of high order schemes, in this paper we provide a thorough discussion on the structure, consistency, stability, phase error, and accuracy of a new 4th-order space-time flux-conserving and neutrally stable CESE solver of an 1D scalar advection equation. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables (the numerical analogues of the dependent variable and its 1st-order and 2ndorder spatial derivatives, respectively) and three equations per mesh point, the new scheme is referred to as the a(3) scheme. Through the von Neumann analysis, it is shown that the a(3) scheme is stable if and only if the Courant number is less than 0.5. Moreover, it is established numerically that the a(3) scheme is 4th-order accurate.

Chang, Sin-Chung↗
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
NASA NTRS1997

Cloud Radiation Forcings and Feedbacks: General Circulation Model Tests and Observational Validation

Using an atmospheric general circulation model (the National Center for Atmospheric Research Community Climate Model: CCM2), the effects on climate sensitivity of several different cloud radiation parameterizations have been investigated. In addition to the original cloud radiation scheme of CCM2, four parameterizations incorporating prognostic cloud water were tested: one version with prescribed cloud radiative properties and three other versions with interactive cloud radiative properties. The authors' numerical experiments employ perpetual July integrations driven by globally constant sea surface temperature forcings of two degrees, both positive and negative. A diagnostic radiation calculation has been applied to investigate the partial contributions of high, middle, and low cloud to the total cloud radiative forcing, as well as the contributions of water vapor, temperature, and cloud to the net climate feedback. The high cloud net radiative forcing is positive, and the middle and low cloud net radiative forcings are negative. The total net cloud forcing is negative in all of the model versions. The effect of interactive cloud radiative properties on global climate sensitivity is significant. The net cloud radiative feedbacks consist of quite different shortwave and longwave components between the schemes with interactive cloud radiative properties and the schemes with specified properties. The increase in cloud water content in the warmer climate leads to optically thicker middle- and low-level clouds and in turn to negative shortwave feedbacks for the interactive radiative schemes, while the decrease in cloud amount simply produces a positive shortwave feedback for the schemes with a specified cloud water path. For the longwave feedbacks, the decrease in high effective cloudiness for the schemes without interactive radiative properties leads to a negative feedback, while for the other cases, the longwave feedback is positive. These cloud radiation parameterizations are empirically validated by using a single-column diagnostic model. together with measurements from the Atmospheric Radiation Measurement program and from the Tropical Ocean Global Atmosphere Combined Ocean-Atmosphere Response Experiment. The inclusion of prognostic cloud water produces a notable improvement in the realism of the parameterizations, as judged by these observations. Furthermore, the observational evidence suggests that deriving cloud radiative properties from cloud water content and microphysical characteristics is a promising route to further improvement.

Lee,Wan-Ho↗
NASA NTRS2009

The 20-22 January 2007 Snow Events over Canada: Microphysical Properties

One of the grand challenges of the Global Precipitation Measurement (GPM) mission is to improve precipitation measurements in mid- and high-latitudes during cold seasons through the use of high-frequency passive microwave radiometry. Toward this end, the Weather Research and Forecasting (WRF) model with the Goddard microphysics scheme is coupled with a Satellite Data Simulation Unit (WRF-SDSU) that has been developed to facilitate over-land snowfall retrieval algorithms by providing a virtual cloud library and microwave brightness temperature (Tb) measurements consistent with the GPM Microwave Imager (GMI). This study tested the Goddard cloud microphysics scheme in WRF for snowstorm events (January 20-22, 2007) that took place over the Canadian CloudSAT/CALIPSO Validation Project (C3VP) ground site (Centre for Atmospheric Research Experiments - CARE) in Ontario, Canada. In this paper, the performance of the Goddard cloud microphysics scheme both with 2ice (ice and snow) and 3ice (ice, snow and graupel) as well as other WRF microphysics schemes will be presented. The results are compared with data from the Environment Canada (EC) King Radar, an operational C-band radar located near the CARE site. In addition, the WRF model output is used to drive the Goddard SDSU to calculate radiances and backscattering signals consistent with direct satellite observations for evaluating the model results.

Tao. W.K.↗
DOE OSTI2021

Using spin alignment of inelastically excited nuclei in fast beams to assign spins: The spectroscopy of 13 O as a test case

Excited states in O 13 were investigated using inelastic scattering of an E / A = 69.5 MeV O 13 beam off of a Be 9 target. The excited states were identified in the invariant-mass spectra of the decay products. Both single-proton and sequential two-proton decays of the excited states were examined. Furthermore, for a number of the excited states, the protons were emitted with strong anisotropy where emissions transverse to the beam axis are favored. The measured proton-decay angular distributions were compared to predictions from distorted-wave Born-approximation calculations of the spin alignment which was shown to be largely independent of the excitation mechanism. The deduced O 13 level scheme is compared to ab initio no-core shell model with continuum predictions. The lowest-energy excited states decay isotropically consistent with predictions of strong proton 1 s 1 / 2 structure. Above these states in the level scheme, we observed a number of higher-spin states not predicted within the model. Possibly these are associated with rotational bands built on deformed cluster configurations predicted by antisymmetrized molecular dynamics calculations. The spin alignment mechanism is shown to be useful for making spin assignments and may have widespread use.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2007

The alpha(3) Scheme - A Fourth-Order Neutrally Stable CESE Solver

The conservation element and solution element (CESE) development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To provide a solid foundation for a systematic CESE development of high order schemes, in this paper we describe a new 4th-order neutrally stable CESE solver of the advection equation Theta u/Theta + alpha Theta u/Theta x = 0. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables u(sup n) (sub j), (u(sub x))(sup n) (sub j) , and (uxz)(sup n) (sub j) (the numerical analogues of u, Theta u/Theta x, and Theta(exp 2)u/Theta x(exp 2), respectively) and four equations per mesh point, the new scheme is referred to as the alpha(3) scheme. As in the case of other similar CESE neutrally stable solvers, the alpha(3) scheme enforces conservation laws in space-time locally and globally, and it has the basic, forward marching, and backward marching forms. These forms are equivalent and satisfy a space-time inversion (STI) invariant property which is shared by the advection equation. Based on the concept of STI invariance, a set of algebraic relations is developed and used to prove that the alpha(3) scheme must be neutrally stable when it is stable. Moreover it is proved rigorously that all three amplification factors of the alpha(3) scheme are of unit magnitude for all phase angles if |v| <= 1/2 (v = alpha delta t/delta x). This theoretical result is consistent with the numerical stability condition |v| <= 1/2. Through numerical experiments, it is established that the alpha(3) scheme generally is (i) 4th-order accurate for the mesh variables u(sup n) (sub j) and (ux)(sup n) (sub j); and 2nd-order accurate for (uxx)(sup n) (sub j). However, in some exceptional cases, the scheme can achieve perfect accuracy aside from round-off errors.

Chang, Sin-Chung↗
NASA NTRS2008

Using Bayesian Networks for Candidate Generation in Consistency-based Diagnosis

Consistency-based diagnosis relies heavily on the assumption that discrepancies between model predictions and sensor observations can be detected accurately. When sources of uncertainty like sensor noise and model abstraction exist robust schemes have to be designed to make a binary decision on whether predictions are consistent with observations. This risks the occurrence of false alarms and missed alarms when an erroneous decision is made. Moreover when multiple sensors (with differing sensing properties) are available the degree of match between predictions and observations can be used to guide the search for fault candidates. In this paper we propose a novel approach to handle this problem using Bayesian networks. In the consistency- based diagnosis formulation, automatically generated Bayesian networks are used to encode a probabilistic measure of fit between predictions and observations. A Bayesian network inference algorithm is used to compute most probable fault candidates.

Narasimhan, Sriram↗
NASA NTRS1993

A two-dimensional model with coupled dynamics, radiation, and photochemistry. 1: Simulation of the middle atmosphere

Spatial and seasonal distribution of ozone and other trace gases are simulated using a interactive two-dimensional model of the stratosphere updated to include full gas phase chemistry. The model consists of a primitive equation dynamics module, a full radiative transfer scheme, and a comprehensive gas phase chemistry module. The circulation is derived from heating rates in the stratosphere that are calculated using model-generated ozone. In the troposphere, parameterized heating rates are adopted. The eddy momentum flux divergence in the zonal mean momentum equation is given by the eddy fluxes of potential vorticity. Eddy fluxes of potential vorticity and tracers are parameterized using a set of predetermined diffusion coefficients. The adopted values for K(sub yy), show a hemispherical asymmetry in that the values in the lower stratosphere are consistently smaller in the southern hemisphere. The asymmetry in K(sub yy) and in the parameterization of the tropospheric heating rate results in an asymmetry in the circulation giving rise to unique signatures in the trace gas distributions. The model successfully simulates the observed asymmetry in the column abundance of the springtime ozone maxima between northern and southern hemisphere. Results for other trace gases are in agreement with the gross observed features although specific differences exist.

Ko, Malcolm K. W.↗
DOE OSTI2021

Federated Quantum Machine Learning

Distributed training across several quantum computers could significantly improve the training time and if we could share the learned model, not the data, it could potentially improve the data privacy as the training would happen where the data is located. One of the potential schemes to achieve this property is the federated learning (FL), which consists of several clients or local nodes learning on their own data and a central node to aggregate the models collected from those local nodes. However, to the best of our knowledge, no work has been done in quantum machine learning (QML) in federation setting yet. In this work, we present the federated training on hybrid quantum-classical machine learning models although our framework could be generalized to pure quantum machine learning model. Specifically, we consider the quantum neural network (QNN) coupled with classical pre-trained convolutional model. Our distributed federated learning scheme demonstrated almost the same level of trained model accuracies and yet significantly faster distributed training. It demonstrates a promising future research direction for scaling and privacy aspects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Data for Impact of Vertical and Seasonal Variation in Leaf Traits on Simulating Soybean Canopy Photosynthesis via 1D and 3D Modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax25), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax25 profile achieved comparable performance to the scheme with the measured Vcmax25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

Photosynthesis↗
DOE OSTI2024

$t\overline{t}b\overline{b}$ at NLO precision in a variable flavor number scheme

Top-quark pair production in association with two $b$-jets is computed at next-to-leading order QCD precision, including effects of the $b$-quark mass, and matched to a $t\overline{t}$+jets simulation in a variable flavor number scheme. The Monte Carlo realization of this method, called fusing, consistently embeds the four-flavor calculation in a particle-level event generator. As a first phenomenological application, we present observables relevant to the data-driven estimation of irreducible backgrounds to $t\overline{t}H$ -production.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-Löwdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSICKLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Lastly the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Levels in 125 Cd populated by the β decay of 125 m Ag and 125 Ag

Here, the β decay of 125 m, 125 Ag into levels in 125 Cd was investigated at the Holifield Radioactive Ion Beam Facility (HRIBF). Uranium-238 targets were bombarded with 50-MeV protons with an intensity of 15 μ A, and the induced fission products were mass separated and deposited on a moving tape in the center of the VANDLE array consisting of γ detectors and plastic scintillators. A partial decay scheme has been assigned for both β decay of the (9/2 + ) ground state of 125 Ag and its low-lying (1/2 - ) isomer, with the energy of the low-lying (11/2 - ) isomeric state in 125 Cd assigned as 188.5 keV. In addition, β -delayed neutron emission probabilities were also determined to be 1.2(2)% for the (9/2 + ) 125 Ag ground state and 4.6(10)% for the (1/2 - ) isomer, which are substantially lower than the previously reported value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Demonstration of controlled spatial incoherence for beam smoothing

A concept for beam smoothing of high-energy laser systems is demonstrated. Multiple coherent beams are multiplexed so that their incoherent mix decreases the far-field nonuniformity after random spatial phase modulation introduced by a phase plate. Modeling shows that the incoherent mix effectively leads to the expected scaling in terms of the number of uncorrelated beams, and that spatial multiplexing can be combined with smoothing by spectral dispersion. The results of an experimental demonstration of this concept with up to six beams, leading to more than 40 modes with polarization multiplexing and angular dispersion, are consistent with expectations. The potential practical advantages of this smoothing scheme and future developments are discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Single-distance nano-holotomography with coded apertures

High-resolution phase-contrast 3D imaging using nano-holotomography typically requires collecting multiple tomograms at varying sample-to-detector distances, usually 3 to 4. This multi-distance approach limits temporal resolution, making it impractical for operando studies. Moreover, shifting the sample complicates reconstruction, requiring precise alignment, registration, and interpolation to correct for shift-dependent magnification on the detector. In response, we propose and validate through simulations a novel, to the best of our knowledge, single-distance approach that leverages coded apertures to structure beam illumination while the sample rotates. Finally, this approach relies on a joint reconstruction scheme, which integrates phase retrieval with 3D tomography, ensuring data consistency and achieving artifact-free reconstructions from a single distance, unlocking dynamic experiments.

Nikitin, Viktor [Argonne National Laboratory (ANL)↗