Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameterization”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
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 OSTI2024

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

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

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗
DOE OSTI2025

Regional-Scale Modeling Parameterizations for Secondary Organic Aerosol Formation from Isoprene Epoxydiols: Experimentally Based Evaluation and Optimization

Isoprene is an abundant volatile organic compound emitted from broadleaf forests. Under low nitric oxide concentrations, isoprene is photochemically oxidized to form gas-phase isoprene epoxydiols (IEPOX). In the presence of acidified sulfate aerosols, IEPOX enhances the secondary organic aerosol (SOA) formation. Predictions of IEPOX-SOA in regional-scale models, e.g., the Community Multiscale Air Quality Model (CMAQ), are uncertain due to homogeneous aerosol assumptions, underpredictions of water uptake (hygroscopicity), and aerosol surface area. Here, we used experimental measurements of IEPOX-SOA tracers, 2-methyltetrols (2-MT) and 2-methyltetrol sulfates (2-MTS), formed at initial IEPOX-to-inorganic sulfate ratios ranging from 1–10.5, at ∼50% relative humidity to constrain key IEPOX-SOA parameters: phase separation, organic shell diffusivity (D org ), acidity, hygroscopic growth, mass accommodation, and kinetics. The base CMAQ parametrization overpredicted experimental IEPOX-SOA with an average normalized mean bias (NMB average ) of 1.63. CMAQ with phase separation underpredicted IEPOX-SOA (NMB average = −0.71). Using the phase-separated model, CMAQ model performance was optimized (NMB average = 0.077) with an increased D org = 2 × 10 –16 m 2 s –1 and increased rate constants (k 2-MT = 1 × 10 –3 M 2 s –1 , k 2-MTS = 8.83 × 10 –3 M 2 s –1 ). The optimized model explicitly accounted for hygroscopic growth by utilizing experimentally derived growth rates, improving aerosol surface area predictions. Our model highlights the importance of the aerosol mixing state (homogeneous versus phase-separated), aerosol size dynamics, and hygroscopic growth in modeling heterogeneous reactive uptake of IEPOX.

aerosols↗
DOE OSTI2023

Sensitivity of deep convection and cross-tropopause water transport to microphysical parameterizations in WRF

Deep convective storms can overshoot the tropopause, thus altering the composition of the stratosphere by vertically transporting tropospheric air. The transport of water vapor and ice particles into a sub-saturated environment can hydrate the stratosphere, with implications for radiative forcing and ozone chemistry. Cloud-resolved models, if employed at high spatial resolutions, are used to probe process-level questions about cross-tropopause deep convective hydration and its controls. There is considerable diversity in model representations of processes associated with water transport and transformation, and the choice of a microphysics scheme affects model skill in simulating deep convective events. This motivates our evaluation of state-of-the-art, as well as widely used standard schemes, in a high spatial- and temporal-resolution framework. Six bulk microphysics schemes were employed in a WRF-LES setup, initialized with a sounding profile representative of a tropopause-overshooting storm. Here, we used an idealized framework to isolate the effect of microphysics on the dominant processes that control the reach of deep convection and stratospheric hydration. All schemes produced the highest reaching updrafts 8-12 hours into the simulation but the strength and persistence of updrafts varied across the schemes; maximum storm heights ranged 9.1-12.6 km across the schemes. Varying microphysics produced large differences in the vertical extent and horizontal aggregation of convection, and an order of magnitude spread in above-tropopause water vapor concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Response of the Current Climate to Land‐Ocean Contrasts in Parameterized Cumulus Entrainment

Abstract Cumulus entrainment substantially regulates the earth's climate but remains poorly constrained in global climate models. Recent studies have shown that cumulus bulk entrainment (or dilution) is particularly sensitive to continentality, with the entrainment rate in simulated maritime cumuli nearly double that of continental cumuli. This study examines the impacts of such land–ocean entrainment contrasts on the current climate using 21‐year simulations with the Geophysical Fluid Dynamics Laboratory's High‐Resolution Atmospheric Model (HIRAM). In response to a 25% reduction in the HIRAM entrainment parameter c 0 over land, precipitation over tropical land regions increases by up to 40%. Along with directly facilitating enhanced convective precipitation, this c 0 reduction induces an increase in soil moisture, which may contribute to a further enhancement of convective precipitation over land. A 25% c 0 reduction over the oceans leads to more widespread modifications of convection patterns, with the strongest signal in the tropical Pacific. Deep convection shifts upstream (eastward) there, inducing enhanced large‐scale ascent over the central Pacific with compensating subsidence and reduced humidity and precipitation over the western Pacific (WP). Land–ocean variations in c 0 project onto the Pacific Walker circulation, with the 25% land reduction strengthening it by 4% and the 25% ocean reduction weakening it by 14%. These changes are driven by variations in convective and large‐scale stratiform heating over the Pacific. While reduced c 0 over land enhances diabatic heating in the Maritime Continent to strengthen the Walker circulation, reduced c 0 over the oceans decreases diabatic heating in the WP to weaken the Walker circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

Marine cold-air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed-phase cloud properties during CAO conditions, the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (~1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (~3–5 km) and cold (–30 to –50 °C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ~1,000 g m –2 ) associated with strong updrafts and enhanced turbulence. We use a high-resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ~10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well-correlated with a cloud transition. Finally, the eddy-diffusivity/mass-flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. Furthermore, we recommend the use of a non-local turbulence closure approach to better capture turbulent processes in intense CAOs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Sensitivity of Anthropogenic Cloud Droplet Number Change to Preindustrial Emission Inventories and Physics Parameterizations

Understanding the preindustrial (PI) to present-day (PD) change in cloud droplet number concentration (ΔNd) is important to constrain the anthropogenic influence on clouds and aid future model projections. Perturbed parameter ensemble simulations of PI and PD conditions in the Energy Exascale Earth System Model (E3SMv3) reveal two locations with consistent negative ΔNd signal across parameter space–the Southeastern US and the United Kingdom–despite a trend of positive ΔNd globally from increased fossil fuel consumption. These negative ΔNd signals are driven by higher anthropogenic biofuel emissions in the PI. By varying the PI aerosol mass and/or number emissions across these data sets' uncertainty ranges, we can change the magnitude and sign of ΔNd and the effective radiative forcing from aerosol-cloud interactions (ERFaci) regionally and globally, with near-zero ΔNd leading to near-zero ERFaci. These results highlight the need to constrain PI emission uncertainty to better understand Earth system responses to aerosol-cloud interactions.

Brown, Hunter Y. [University of Wyoming, Laramie, ↗
DOE OSTI2024

A Lie algebraic theory of barren plateaus for deep parameterized quantum circuits

Variational quantum computing schemes train a loss function by sending an initial state through a parametrized quantum circuit, and measuring the expectation value of some operator. Despite their promise, the trainability of these algorithms is hindered by barren plateaus (BPs) induced by the expressiveness of the circuit, the entanglement of the input data, the locality of the observable, or the presence of noise. Up to this point, these sources of BPs have been regarded as independent. In this work, we present a general Lie algebraic theory that provides an exact expression for the variance of the loss function of sufficiently deep parametrized quantum circuits, even in the presence of certain noise models. Our results allow us to understand under one framework all aforementioned sources of BPs. This theoretical leap resolves a standing conjecture about a connection between loss concentration and the dimension of the Lie algebra of the circuit’s generators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The ellipticity parameterization for an NFW profile: An overlooked angular structure in strong lens modeling

Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is a Navarro Frenk White (NFW) profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because an unmodeled azimuthal structure has been shown to be able to bias lens model results, we explored the degree to which this azimuthal structure that was introduced can affect the model accuracy. We constructed input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculated the value of the Hubble parameter H 0 one would determine from the lensing fit. We found that the fits to the NFWp input return H 0 values that are systematically biased by about 3% lower than the NFWm counterparts. We explored whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, ξ 2 . We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias was found beyond the MST due to the azimuthal structure exterior to the Einstein radius. We recommend modelers use an NFWm prescription in the future, such that the azimuthal structure can be introduced explicitly rather than implicitly.

79 ASTRONOMY AND ASTROPHYSICS↗