Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

On the Dependence of Cloud Feedbacks on Physical Parameterizations in WRF Aquaplanet Simulations

We investigate the effects of physical parameterizations on cloud feedback uncertainty in response to climate change. For this purpose, we construct an ensemble of eight aquaplanet simulations using the Weather Research and Forecasting (WRF) model. In each WRF‐derived simulation, we replace only one parameterization at a time while all other parameters remain identical. By doing so, we aim to (i) reproduce cloud feedback uncertainty from state‐of‐the‐art climate models and (ii) understand how parametrizations impact cloud feedbacks. Our results demonstrate that this ensemble of WRF simulations, which differ only in physical parameterizations, replicates the range of cloud feedback uncertainty found in state‐of‐the‐art climate models. We show that microphysics and convective parameterizations govern the magnitude and sign of cloud feedbacks, mostly due to tropical low‐level clouds in subsidence regimes. Finally, this study highlights the advantages of using WRF to analyze cloud feedback mechanisms owing to its plug‐and‐play parameterization capability.

Cesana, Gregory↗
DOE OSTI2026

Fantômas unconfined: global QCD fits with Bézier parameterizations

Fantômas is a C++ toolkit for exploring the parametrization dependence of parton distribution functions (PDFs) and other correlator functions in quantum chromodynamics (QCD). Fantômas facilitates the generation of adaptable polynomial parametrizations for PDFs, called metamorphs, to find best-fit PDF solutions and quantify the epistemic uncertainty associated with the parametrizations during their fitting. The method employs Bézier curves as universal approximators for a variety of PDF shapes. Integrated into the xFitter framework for the global QCD analysis, Fantômas provides a foundation for general models of PDFs, while reducing the computational time compared to the approaches utilizing traditional polynomial parametrizations as well as providing an interpretable alternative to neural-network-based models. This paper outlines the structure and practical usage of the Fantômas toolkit, including its inputs, outputs, and implementation within xFitter. It also provides a practical example of using Fantômas for uncertainty quantification as well as the combination of PDF fits into a single ensemble.

Bézier curves↗
NASA NTRS2003

Jell-Molds and Cookie-cutters: Shrinkwrap Isn't Just for Leftovers Anymore

So what is Shrinkwrap all about? For those of you who may not know about it, Shrinkwrap is a type of data structure that can manifest itself as a feature or model. It is cleverly covered up, almost hidden, and doesn't get the press or widespread use of a solid or surface. The shrinkwrap feature is located under the data sharing submenu of the feature menu. The shrinkwrap feature, as described by PTC, is a collection of surfaces and datum features of a model that represents the exterior of the model . The advantages and applications of the shrinkwrap feature are in the creation of minimal memory guzzling representations of assemblies. These can be used to represent subassemblies in parent assemblies, and can handle control of dependency issues, geometry represented, and additional references through the use of the shrinkwrap feature options. The shrinkwrap model is an option available under the save as umbrella. Its function, as described by PTC, is to share data with internal and external design groups and improve performance in large assembly design . Some of the benefits of the shrinkwrap model include being able to represent complex assemblies with a single, lightweight part that protects design intent and parametric data, and the ability to improve performance of large assembly modeling in the area of less load time. The proper-scale models can be saved as IGES, STEP, and VRML (for fly-throughs).

Randazzo, John↗
DOE OSTI2025

Optimization problems governed by systems of PDEs with uncertainties

This paper reviews current theoretical and numerical approaches to optimization problems governed by partial differential equations (PDEs) that depend on random variables or random fields. Such problems arise in many engineering, science, economics and societal decision-making tasks. This paper focuses on problems in which the governing PDEs are parametrized by the random variables/fields, and the decisions are made at the beginning and are not revised once uncertainty is revealed. Examples of such problems are presented to motivate the topic of this paper, and to illustrate the impact of different ways to model uncertainty in the formulations of the optimization problem and their impact on the solution. A linear–quadratic elliptic optimal control problem is used to provide a detailed discussion of the set-up for the risk-neutral optimization problem formulation, study the existence and characterization of its solution, and survey numerical methods for computing it. Different ways to model uncertainty in the PDE-constrained optimization problem are surveyed in an abstract setting, including risk measures, distributionally robust optimization formulations, probabilistic functions and chance constraints, and stochastic orders. Furthermore, approximation-based optimization approaches and stochastic methods for the solution of the large-scale PDE-constrained optimization problems under uncertainty are described. Some possible future research directions are outlined.

Heinkenschloss, Matthias [Rice Univ., Houston, TX ↗
DOE OSTI2025

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗
NASA NTRS1973

An analysis of atmospheric entry trajectories for manned and unmanned missions to the planet Venus

A unified atmospheric density model is presented for the planet Venus, and this model is compatible with Russian and American data obtained during the 1967 Venus launch opportunity. Trajectory characteristics involved in atmospheric entry are discussed with reference to roll-control modulation and entry corridor. Subsequent to a trajectory-oriented discussion of potential Venus mission characteristics, a parametric analysis of manned and unmanned vehicle entry trajectories into the Venusian atmosphere is presented. A sensitivity analysis with reference to atmospheric density deviations is included to show the dependence of corridor depth on the atmospheric density profile.

Craig, M. K.↗
NASA NTRS2014

Ice Roughness in Short Duration SLD Icing Events

Ice accretion codes depend on models of roughness parameters to account for the enhanced heat transfer during the ice accretion process. While mitigating supercooled large droplet (SLD or Appendix O) icing is a significant concern for manufacturers seeking future vehicle certification due to the pending regulation, historical ice roughness studies have been performed using Appendix C icing clouds which exhibit mean volumetric diameters (MVD) much smaller than SLD clouds. Further, the historical studies of roughness focused on extracting parametric representations of ice roughness using multiple images of roughness elements. In this study, the ice roughness developed on a 21-in. NACA 0012 at 0deg angle of attack exposed to short duration SLD icing events was measured in the Icing Research Tunnel at the NASA Glenn Research Center. The MVD's used in the study ranged from 100 micrometer to 200 micrometers, in a 67 m/s flow, with liquid water contents of either 0.6 gm/cubic meters or 0.75 gm/cubic meters. The ice surfaces were measured using a Romer Absolute Arm laser scanning system. The roughness associated with each surface point cloud was measured using the two-dimensional self-organizing map approach developed by McClain and Kreeger (2013) resulting in statistical descriptions of the ice roughness.

Icing↗
DOE OSTI2022

Impact of Irradiation-induced Lattice Defects and Dopants on Thermal Transport in Oxide Nuclear Fuels

Advanced nuclear reactor concepts aim to use fuels that must withstand unprecedented temperature and radiation extremes. In these fuels, thermal energy transport under irradiation is directly related to fuel longevity, reactor safety, and is arguably one of the most important performance metrics. Here we provide a comprehensive, first-principles-informed treatment of phonon mediated thermal transport in a defect-bearing actinide oxide with direct comparison to experimental measurements. Pristine and proton irradiated thorium dioxide was chosen as a model system to treat the complexity of thermal transport in the presence of lattice defects. A thermal transport model is implemented using the linearized Boltzmann transport equation (LBTE) with input from first principles calculations and defect evolution models. The output of the LBTE is compared directly to mesoscopic measurements of thermal conductivity on length scales commensurate with defect accumulation. Parametric measurements of conductivity with irradiation dose and temperature reveal a monotonic decrease in conductivity with irradiation dose, in the range of 0.001 dpa to 0.1 dpa. Besides radiation-induced defects, doping the thorium dioxide crystal with small quantities of uranium atoms also results in a drastic reduction in thermal conductivity. Temperature-dependent measurements of thermal conductivity in uranium-doped thorium dioxide single crystals suggests a resonant scattering mechanism is responsible for the observed reduction. This comprehensive, atomistic- to meso-scale treatment provides the necessary basis to investigate thermal transport under irradiation in more complex systems that exhibit strong electron correlation. This is a poster presentation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Investigation of methane and ethane pyrolysis with highly active and durable iron-alumina catalyst to produce hydrogen and valuable nano carbons: Continuous fluidized bed tests and reaction rate analysis

The goal of this research work was to evaluate a novel iron-alumina catalyst for methane and ethane pyrolysis with high conversion rates to H 2 with a long life time to address two major issues reported with previous catalysts. Fluidized bed tests demonstrated continuous hydrogen production while collecting carbon containing fines. During the 230 h fluid bed test the methane conversion rate to hydrogen at 700 °C remained above 80% for 160 h while creating valuable nano carbons. Similarly high catalyst lifetimes with high methane conversions have not been reported before. Ethane was more active than methane for pyrolysis with 100% conversion to H 2 . Transmission electron microscopy and Raman analysis of the carbon showed single and multiwall nanotubes and conglomerate nano-fibers. Ethane did not impact the allotrope of carbon formed. TGA was used to create a data array that was used for kinetic analysis and reaction model parametric regressions. Ethane had a higher rate of pyrolysis that was 2.3–2.8 times that of methane. Ethane pyrolysis was found to be more sensitive to temperature with a higher activation energy of 62 kJ/mol as compared to methane at 43 kJ/mol. The frequency factor, A i , was significantly higher for ethane at 0.73 [[kmol/m 3 -sec]/[N/m2] n ] compared to 0.03 for methane on the same basis. Both gases had a functional partial pressure order dependence of ~0.6 based on the regression analysis.

08 HYDROGEN↗
DOE OSTI2021

Influence of gas expansion on the propagation of a premixed flame in a spatially periodic shear flow

It has been previously demonstrated that thermal gas expansion might have a role in boundary layer flashback of premixed turbulent flames, inducing local flow-reversal in the boundary layer's low-velocity streaks on the reactants’ side of the flame and facilitating its upstream propagation. Here, we perform a two-dimensional numerical investigation of the interaction between a periodic shear flow and a laminar premixed flame. The periodic shear is a simplified model for the oncoming prolonged streamwise velocity streaks with alternating regions of high and low velocities found in turbulent boundary layers in the vicinity of the walls. The parametric study focuses on the amplitude and wavelength of the periodic shear flow and on the gas expansion ratio (unburnt-to-burnt density ratio). With the increase of the amplitudes of the periodic shear flow and of the gas expansion, the curved flame velocity increases monotonically. The flame velocity dependence on the periodic shear wavelength is non-monotonic, which is consistent with previous theoretical studies of curved premixed flame velocity. The flame shape that is initially formed by the oncoming periodic shear appears to be metastable. At a later stage of the flame propagation, the flame shape transforms into the stationary one dominated by the Darrieus-Landau instability.

42 ENGINEERING↗
DOE OSTI2024

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Force Balance Model Assessment for Mechanistic Prediction of Sliding Bubble Velocity in Vertical Subcooled Boiling Flow

The bubble sliding after departing from nucleation site is frequently observed in flow boiling systems and, the crucial impact on wall heat transfer has been evidenced through many experiments. As a result, the heat transfer modeling associated with sliding bubble has become one of the subjects of great attention in CFD boiling heat transfer community. The modeling efforts are primarily aimed at improving the existing Heat Flux Partitioning (HFP) model via the implementation of sliding bubble-induced heat transfer. The performance of HFP model depends inherently on the fidelity of sub-models used to predict the fundamental bubble parameters (e.g., bubble departure/lift-off diameter). In the same context, the accurate prediction of sliding bubble parameters (e.g., sliding bubble growth, sliding bubble velocity) is essential to achieving the successful heat transfer modeling associated with sliding bubble. Of many sliding bubble parameters, this paper deals with the sliding bubble velocity. Specifically, the force balance model was assessed in view of the sliding bubble velocity prediction. The parametric effect of key sub-models (e.g., drag force, bubble growth models) used in the force balance equation was investigated. The experimental data from Maity (2000) and Yoo et al. (2016) were used for demonstrating the model predictive performance. It was found that the force balance model proposed in this study was able to predict well the bubble sliding velocity based on the accurate prediction of bubble growth during sliding. The predictive performance was proven with the experimental data measured under various subcooled boiling conditions of both water and refrigerant (NOVEC-7000) flowing upward in vertical channels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1974

A parameterization for the absorption of solar radiation in the earth's atmosphere

A method is described for rapidly computing the amount of solar energy absorbed at the earth's surface and in the atmosphere as a function of altitude. The method is a parametric treatment, but the form of the solution and the coefficients involved are based on accurate multiple-scattering computations. In this treatment the absorption varies with the amount and type of clouds, the humidity, the zenith angle of the sun, and the albedo of the earth's surface. Within the stratosphere the absorption also depends on the vertical distribution of ozone. This parameterization for solar radiation is being used in current versions of the global atmospheric circulation model developed at the Goddard Institute for Space Studies.

Lacis, A. A.↗
NASA NTRS1976

Development of relaxation turbulence models

Relaxation turbulence models have been intensively studied. The complete time dependent mass averaged Navier-Stokes equations have been solved for flow into a two dimensional compression corner. A new numerical scheme has been incorporated into the developed computed code with an attendant order of magnitude reduction in computation time. Computed solutions are compared with experimental measurements of Law for supersonic flow. Details of the relaxation process have been studied; several different relaxation models, including different relaxation processes and varying relaxation length, are tested and compared. Then a parametric study has been conducted in which both Reynolds number and wedge angle are varied. To assess effects of Reynolds number and wedge angle, the parametric study includes the comparison of computed separation location and upstream extent of pressure rise; numerical results are also compared with the measurements of surface pressure, skin friction and mean velocity field.

Hung, C. M.↗
NASA NTRS2018

Improved Hydrology over Peatlands in a Global Land Modeling System

Peatlands of the Northern Hemisphere represent an important carbon pool that mainly accumulated since the last ice age under permanently wet conditions in specific geological and climatic settings. The carbon balance of peatlands is closely coupled to water table dynamics. Consequently, the future carbon balance over peatlands is strongly dependent on how hydrology in peatlands will react to changing boundary conditions, e.g. due to climate change or regional water level drawdown of connected aquifers or streams. Global land surface modeling over organic-rich regions can provide valuable global-scale insights on where and how peatlands are in transition due to changing boundary conditions. However, the current global land surface models are not able to reproduce typical hydrological dynamics in peatlands well. We implemented specific structural and parametric changes to account for key hydrological characteristics of peatlands into NASA's GEOS-5 Catchment Land Surface Model (CLSM, Koster et al. 2000). The main modifications pertain to the modeling of partial inundation, and the definition of peatland-specific runoff and evapotranspiration schemes. We ran a set of simulations on a high performance cluster using different CLSM configurations and validated the results with a newly compiled global in-situ dataset of water table depths in peatlands. The results demonstrate that an update of soil hydraulic properties for peat soils alone does not improve the performance of CLSM over peatlands. However, structural model changes for peatlands are able to improve the skill metrics for water table depth. The validation results for the water table depth indicate a reduction of the bias from 2.5 to 0.2 m, and an improvement of the temporal correlation coefficient from 0.5 to 0.65, and from 0.4 to 0.55 for the anomalies. Our validation data set includes both bogs (rain-fed) and fens (ground and/or surface water influence) and reveals that the metrics improved less for fens. In addition, a comparison of evapotranspiration and soil moisture estimates over peatlands will be presented, albeit only with limited ground-based validation data. We will discuss strengths and weaknesses of the new model by focusing on time series of specific validation sites.

Bechtold, M.↗
DOE OSTI2022

Evidence of Quadrupole and Octupole Deformations in Zr 96 + Zr 96 and Ru 96 + Ru 96 Collisions at Ultrarelativistic Energies

In the hydrodynamic model description of heavy-ion collisions, the elliptic flow v 2 and triangular flow v 3 are sensitive to the quadrupole deformation β 2 and octupole deformation β 3 of the colliding nuclei. The relations between v n and β n have recently been clarified and were found to follow a simple parametric form. The STAR Collaboration has just published precision v n data from isobaric Ru 96 + Ru 96 and Zr 96 + Zr 96 collisions, where they observe large differences in central collisions v 2 , Ru > v 2 , Zr and v 3 , Ru < v 3 , Zr . Using a transport model simulation, we show that these orderings are a natural consequence of $\beta_{2,Ru} \gg \beta_{2,Zr}$ and $\beta_{3,Ru} \ll \beta_{3,Zr}$ We reproduce the centrality dependence of the v 2 ratio qualitatively and v 3 ratio quantitatively and extract values of β 2 and β 3 that are consistent with those measured at low-energy nuclear structure experiments. STAR data provide the first direct evidence of strong octupole correlations in the ground state of Zr 96 in heavy-ion collisions. Our analysis demonstrates that flow measurements in high-energy, heavy-ion collisions, especially using isobaric systems, are a new precision tool to study nuclear structure physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2023

Multiphysics Modeling of Electrode Dissolution under Hydrodynamic Conditions in Low Conductivity Water

Efficient plug-flow electrochemical reactors are characterized by their ability to effectively convert species with significant ionic strength and/or with the assistance of a well-behaved electrolyte. Because of the classic electroconductivity requirement, electrochemical reactors are not particularly attractive options for applications with low-conductivity conditions, especially for low-pressure feeds. However, these reactors can still find application in dedicated processes such as ultra-clean water conditioning and biological-coupled separation systems. Consequently, a multi-physics model was developed using COMSOL Multiphysics to characterize different reactor-design conditions based on the Tertiary Current Distribution which accounts for the effects of expected variations in electrolyte composition and ionic strength on the electrochemical process, as well as solution resistance and electrode kinetics. The modelling framework presented here employs single-phase laminar fluid flow, the Nernst-Planck equation, the water-based electroneutrality condition, and concentration-dependent overpotentials. This investigation particularly explores the electrolysis process of “anodic dissolution,” in which elemental species dissociate from the solid matrix of the electrode and enter the liquid phase of the electrolyte as soluble ions. The reactor configuration of interest has a rectangular parallel-electrode design and is operated galvanostatically. A series of parametric studies are carried out to inspect the response of the system when hydrodynamic, kinetic, and geometric conditions are changed. The parametric sweep resulted in variable effluent concentrations and system voltage, which reveal underlaying optimization patterns for electrode dissolution in low conductivity water.

Multiphysics↗