Search NASA⌕ Search

SEARCH · Search NASA

Results for “comparative modeling”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Constraints Pb-210 and Be-7 on Wet Deposition and Transport in a Global Three-Dimensional Chemical Tracer Model Driven by Assimilated Meteorological Fields

The atmospheric distributions of the aerosol tracers Pb-210 and Be-7 are simulated with a global three-dimensional model driven by assimilated meteorological observations for 1991-1996 from the NASA Goddard Earth Observing System (GEOSl). The combination of terrigenic Pb-210 and cosmogenic Be-7 provides a sensitive test of wet deposition and vertical transport in the model. Our simulation of moist transport and removal includes scavenging in wet convective updrafts (40% scavenging efficiency per kilometer of updraft), midlevel entrainment and detrainment, first-order rainout and washout from both convective anvils and large-scale precipitation, and cirrus precipitation. Observations from surface sites in specific years are compared to model results for the corresponding meteorological years, and observations from aircraft missions over the Pacific are compared to model results for the days of the flights. Initial simulation of Be-7 showed that cross-tropopause transport in the GEOSl meteorological fields is too fast by a factor of 3-4. We adjusted the stratospheric Be-7 source to correct the tropospheric simulation. Including this correction, we find that the model gives a good simulation of observed Pb-210 and Be-7 concentrations and deposition fluxes at surface sites worldwide, with no significant global bias and with significant success in reproducing the observed latitudinal and seasonal distributions. We achieve several improvements over previous models; in particular, we reproduce the observed Be-7 minimum in the tropics and show that its simulation is sensitive to rainout from convective anvils. Comparisons with aircraft observations up to 12-km altitude suggest that cirrus precipitation could be important for explaining the low concentrations in the middle and upper troposphere.

Liu, Hong-Yu↗
NASA NTRS1985

Features of a reattaching turbulent shear layer in divergent channel flow

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step in a diverging channel flow. Mean velocities, Reynolds stresses, and triple products that were measured by a laser Doppler velocimeter are presented for two cases of tunnel wall divergence. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) terms are extracted from the data. These data are compared with various eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model; however, like the others, it dramatically overpredicts the experimentally determined dissipation rate.

Driver, D. M.↗
NASA NTRS2014

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗
NASA NTRS1992

Comment on: The present state and future direction of eddy viscosity models

Two turbulence models are compared under various flow conditions. Presented first is a k-omega model which is best suited for wall bounded flows. Presented second is a k-epsilon model which is best suited for free shear flow. Generalizations of the k-omega model are discussed, and various aerodynamic parameters such as skin friction, velocity distribution, and hypersonic cylinder-flare, are compared for the different models. All results are presented in viewgraph format.

Coakley, T. J.↗
DOE OSTI2020

Sensitivity of Thermal Predictions to Uncertain Surface Tension Data in Laser Additive Manufacturing

To understand the process-microstructure relationships in additive manufacturing (AM), it is necessary to predict the solidification characteristics in the melt pool. This study investigates the influence of Marangoni driven fluid flow on the predicted melt pool geometry and solidification conditions using a continuum finite volume model. A calibrated laser absorptivity was determined by comparing the model predictions (neglecting fluid flow) against melt pool dimensions obtained from single laser melt experiments on a nickel super alloy 625 (IN625) plate. Using this calibrated efficiency, predicted melt pool geometries agree well with experiments across a range of process conditions. When fluid mechanics is considered, a surface tension gradient recommended for IN625 tends to overpredict the influence of convective heat transfer, but the use of an intermediate value reported from experimental measurements of a similar nickel super alloy produces excellent experimental agreement. Despite its significant effect on the melt pool geometry predictions, fluid flow was found to have a small effect on the predicted solidification conditions compared to processing conditions. We find that this result suggests that under certain circumstances, a model only considering conductive heat transfer is sufficient for approximating process-microstructure relationships in laser AM. Extending the model to multiple laser passes further showed that fluid flow also has a small effect on the solidification conditions compared to the transient variations in the process. Furthermore, limitations of the current model and areas of improvement, including uncertainties associated with the phenomenological model inputs are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Modeling CMB lensing cross correlations with CLEFT

A new generation of surveys will soon map large fractions of sky to ever greater depths and their science goals can be enhanced by exploiting cross correlations between them. In this paper we study cross correlations between the lensing of the CMB and biased tracers of large-scale structure at high z. We motivate the need for more sophisticated bias models for modeling increasingly biased tracers at these redshifts and propose the use of perturbation theories, specifically Convolution Lagrangian Effective Field Theory (CLEFT). Since such signals reside at large scales and redshifts, they can be well described by perturbative approaches. We compare our model with the current approach of using scale independent bias coupled with fitting functions for non-linear matter power spectra, showing that the latter will not be sufficient for upcoming surveys. Furthermore, we illustrate our ideas by estimating σ 8 from the auto- and cross-spectra of mock surveys, finding that CLEFT returns accurate and unbiased results at high z. We report on uncertainties due to the redshift distribution of the tracers, and several avenues for future development.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Multi-Decadal Change of Atmospheric Aerosols and their Effect on Surface Radiation

We present an investigation on multi-decadal changes of atmospheric aerosols and their effects on surface radiation using a global chemistry transport model, GOCART, along with the near-term to long-term data records. We focus on a 28-year time period of satellite era from 1980 to 2007 during which a suite of aerosol data from satellite observations, ground-based measurements, and intensive field experiments have become available. Particularly: (1) We compare the model calculated clear sky downward radiation at the surface with surface network data from BSRN and CMA (2) We compare the model and surface data with satellite derived downward radiation products from ISCCP and SRS (3) We analyze the long-term global and regional aerosol trends in major anthropogenic source regions (North America, Europe, Asia) that have been experiencing considerable changes of emissions during the three decades, dust and biomass burning regions that have large interannual variability, downwind regions that are directly affected by the changes in the source area, and remote regions that are considered to representing "background" conditions. The comparisons and methods from this study can be applied to multiple model analysis in the AeroCom framework.

Chin, Mian↗
DOE OSTI2015

Lower bound on the electroweak wall velocity from hydrodynamic instability

The subsonic expansion of bubbles in a strongly first-order electroweak phase transition is a convenient scenario for electroweak baryogenesis. For most extensions of the Standard Model, stationary subsonic solutions (i.e., deflagrations) exist for the propagation of phase transition fronts. However, deflagrations are known to be hydrodynamically unstable for wall velocities below a certain critical value. We calculate this critical velocity for several extensions of the Standard Model and compare with an estimation of the wall velocity. In general, we find a region in parameter space which gives stable deflagrations as well as favorable conditions for electroweak baryogenesis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1979

A comparison of two accretional heating models

Models of accretional heating that include large-body impact effects are reviewed. Various details of the models are compared, as well as the results of their calculations. Within the uncertainties in the models, the results compare reasonably with each other, and a post-accretion moon with some melting in its outer layers, but not much towards the center, seems to be described by these analyses.

Ransford, G. A.↗
NASA NTRS2012

Multi-decadal Change of Atmospheric Aerosols and their Effect on Surface Radiation

We present an investigation on multi-decadal changes of atmospheric aerosols and their effects on surface radiation using a global chemistry transport model GOCART along with the near-term to long-term data records. We focus on a 28-year time period of satellite era from 1980 to 2007 during which a suite of aerosol data from satellite observations, ground-based measurements, and intensive field experiments have become available. Particularly: (1) We compare the model calculated clear sky downward radiation at the surface with surface network data from Baseline Surface Radiation Network (BSRN) and CMA (2) We compare the model and surface data with satellite derived downward radiation products from ISCCP and SRB (3) We analyze the long-term global and regional aerosol trends in major anthropogenic source regions (North America, Europe, Asia) that have been experiencing considerable changes of emissions during the three decades, dust and biomass burning regions that have large interannual variability, downwind regions that are directly affected by the changes in the source area, and remote regions that are considered to representing "background" conditions.

Chin, Mian↗
DOE OSTI2018

A coupled rate theory-Monte Carlo model of helium bubble evolution in plasma-facing micro-engineered tungsten

A multiscale model of helium bubble evolution in plasma-facing materials is developed. The model links different stages of helium bubble evolution: deposition, nucleation, growth, motion, and coalescence. Helium deposition is simulated with the SRIM Monte Carlo program to give spatial information on helium and displacement damage distributions near the surface. This deposition profile is then introduced into a space-dependent rate theory of bubble nucleation and growth to describe the early stages of the distribution and size of helium bubbles. The coarsening stage of bubble evolution as a result of whole bubble motion, interaction, and coalescence is modeled by a new Object Kinetic Monte Carlo (OKMC) model, for which initial conditions are taken from the mean-field rate theory calculations. The model is compared to experimental data on low-energy helium plasma interaction with micro-engineered tungsten (W), and on high-energy helium ion deposition in flat W samples. Novel features of the multiscale model include: (1) space-dependent rate theory; (2) OKMC model of bubble motion in stress and temperature fields; and (3) application of the model to micro-engineered materials, and comparison with experiments on the same time-scale. At low helium ion energy, it is found that the mechanism of trap mutation is essential in achieving good agreement with experimental measurements. On the other hand, good agreement with experiments at high incident ion energy and temperature showed the importance of bubble coalescence and coarsening as main mechanisms. The results of the model are compared with experiments on flat W surfaces irradiated at high ion energy (30 keV), and with micro-engineered W, where the surface is coated with high-density micro-pillars at low ion energy (around 100 eV). The predicted average bubble radius and density are in qualitative agreement with experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Extending the applicability of Redfield theories into highly non-Markovian regimes

We present a new, computationally inexpensive method for the calculation of reduced density matrix dynamics for systems with a potentially large number of subsystem degrees of freedom coupled to a generic bath. The approach consists of propagation of weak-coupling Redfield-like equations for the high-frequency bath degrees of freedom only, while the low-frequency bath modes are dynamically arrested but statistically sampled. We examine the improvements afforded by this approximation by comparing with exact results for the spin-boson model over a wide range of parameter space. We further generalize the method to multi-site models and compare with exact results for a model of the Fenna–Matthews–Olson complex. The results from the method are found to dramatically improve Redfield dynamics in highly non-Markovian regimes, at a similar computational cost. Relaxation of the mode-freezing approximation via classical (Ehrenfest) evolution of the low-frequency modes results in a dynamical hybrid method. We find that this Redfield-based dynamical hybrid approach, which is computationally more expensive than bare Redfield dynamics, yields only a marginal improvement over the simpler approximation of complete mode arrest.

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

Comparison of structurally diverse simulation models for prediction of epidemic outcomes caused by a long-distance dispersed pathogen

Long-distance dispersal (LDD) pathogens pose substantial challenges for epidemic control due to their ability to generate new infection foci at great distances. While various modeling approaches have been developed to understand and manage such outbreaks, little work has compared how models of different structures behave under shared conditions. Here, in this study, we compare four structurally distinct epidemiological models — EPIMUL, GEMF, PoPS, and Warwick — each adapted to simulate the spread of wheat stripe rust (WSR), a wind-dispersed LDD pathogen, under identical epidemiological parameters and dispersal kernel. Using data from a controlled field experiment, we evaluate the ability of each model to replicate disease prevalence under nine intervention scenarios that vary in timing and culling area. While the models differ substantially in design — ranging from spatial grid-based to network-based and raster-based frameworks — the shared dispersal kernel allowed for close alignment in their predictions. All models accurately captured general epidemic trends, particularly the strong effect of early intervention on disease suppression. We qualitatively compared their behavioral responses across scenarios and also evaluated an ensemble prediction by averaging across model outputs. Our findings highlight how integrating shared epidemiological components into distinct modeling frameworks can improve consistency and accuracy, while reinforcing the importance of early culling in managing LDD pathogen outbreaks.

Dispersal kernel↗
DOE OSTI2022

Machine Learning with Gradient-Based Optimization of Nuclear Waste Vitrification with Uncertainties and Constraints

Gekko is an optimization suite in Python that solves optimization problems involving mixed-integer, nonlinear, and differential equations. The purpose of this study is to integrate common Machine Learning (ML) algorithms such as Gaussian Process Regression (GPR), support vector regression (SVR), and artificial neural network (ANN) models into Gekko to solve data based optimization problems. Uncertainty quantification (UQ) is used alongside ML for better decision making. These methods include ensemble methods, model-specific methods, conformal predictions, and the delta method. An optimization problem involving nuclear waste vitrification is presented to demonstrate the benefit of ML in this field. ML models are compared against the current partial quadratic mixture (PQM) model in an optimization problem in Gekko. GPR with conformal uncertainty was chosen as the best substitute model as it had a lower mean squared error of 0.0025 compared to 0.018 and more confidently predicted a higher waste loading of 37.5 wt% compared to 34 wt%. The example problem shows that these tools can be used in similar industry settings where easier use and better performance is needed over classical approaches. Future works with these tools include expanding them with other regression models and UQ methods, and exploration into other optimization problems or dynamic control.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2009

Validation of PICA Ablation and Thermal-Response Model at Low Heat Flux

Phenolic Impregnated Carbon Ablator (PICA) was the forebody heatshield material on the Stardust sample-return capsule and is also a primary candidate material for the Mars Science Lander (MSL), the Orion Crew Module, and the SpaceX Dragon vehicle. As part of the heatshield qualification for Orion, physical and thermal properties of virgin and charred PICA were measured, and an ablation and thermal response model was developed. We validated the model by comparing it with recession and temperature data from stagnation arcjet tests conducted over a wide range of stagnation heat flux of 107 to 1102 W/sq cm. The effect of orthotropic thermal conductivity was evident in the thermal response of the arcjet models. In general, model predictions compared well with the data; however, the uncertainty of the recession prediction was greatest for heat fluxes below 200 W/sq cm. More recent MSL testing focused on the low heat flux regime of 45 to 250 W/sq cm. The new results confirm the recession uncertainty, especially for pressures below 6 kPa. In this work we focus on improving the model predictions for MSL and Orion tests below 250 W/sq cm.

Milos, Frank S.↗
DOE OSTI2017

A comprehensive iso-octane combustion model with improved thermochemistry and chemical kinetics

Iso-Octane (2,2,4-trimethylpentane) is a primary reference fuel and an important component of gasoline fuels. Furthermore, it is a key component used in surrogates to study the ignition and burning characteristics of gasoline fuels. This paper presents an updated chemical kinetic model for iso-octane combustion. Specifically, the thermodynamic data and reaction kinetics of iso-octane have been re-assessed based on new thermodynamic group values and recently evaluated rate coefficients from the literature. The adopted rate coefficients were either experimentally measured or determined by analogy to theoretically calculated values. New alternative isomerization pathways for peroxy-alkyl hydroperoxide ($\dot{O}$OQOOH) radicals were added to the reaction mechanism. The updated kinetic model was compared against new ignition delay data measured in rapid compression machines (RCM) and a high-pressure shock tube. Our experiments were conducted at pressures of 20 and 40 atm, at equivalence ratios of 0.4 and 1.0, and at temperatures in the range of 632–1060 K. The updated model was further compared against shock tube ignition delay times, jet-stirred reactor oxidation speciation data, premixed laminar flame speeds, counterflow diffusion flame ignition, and shock tube pyrolysis speciation data available in the literature. Finally, the updated model was used to investigate the importance of alternative isomerization pathways in the low temperature oxidation of highly branched alkanes. When compared to available models in the literature, the present model represents the current state-of-the-art in fundamental thermochemistry and reaction kinetics of iso-octane; and thus provides the best prediction of wide ranging experimental data and fundamental insights into iso-octane combustion chemistry.

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

Effective hopping in holographic Bose and Fermi-Hubbard models

In this paper, we analyze a proposed gravity dual to a SU(N) Bose-Hubbard model, as well as construct a holographic dual of a SU(N) Fermi-Hubbard model from D-branes in string theory. In both cases, the SU(N) is dynamical, i.e. the hopping degrees of freedom are strongly coupled to SU(N) gauge bosons which themselves are strongly interacting. The vacuum expectation value (VEV) of the hopping term (i.e. the hopping energy) is analyzed in the gravity dual as a function of the bulk mass of the field dual to the hopping term, as well as of the coupling constants of the model. The bulk mass controls the anomalous dimension (i.e. the critical exponent) of the hopping term in the SU(N) Bose-Hubbard model. We compare the hopping energy to the corresponding result in a numerical simulation of the ungauged SU(N) Bose-Hubbard model. We find agreement when the hopping parameter is smaller than the other couplings. Our analysis shows that the kinetic energy increases as the bulk mass increases, due to increased contributions from the IR. The holographic Bose-Hubbard model is then compared with the string theory construction of a SU(N) Fermi-Hubbard model. The string theory construction makes it possible to describe fluctuations around a half-filled state in the supergravity limit, which map to O(1) occupation number fluctuations in the Fermi-Hubbard model at half filling. Finally, the VEV of the Bose-Hubbard model is shown to agree with the one of the fermionic Hubbard model with the help of a two-site version of the Jordan-Wigner transformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

The Separate Physics and Dynamics Experiment (SPADE) framework for determining resolution awareness: A case study of microphysics

The ability to use multi-resolution dynamical cores for weather and climate modeling is pushing the atmospheric community towards developing scale aware or, more specifically, resolution aware parameterizations that will function properly across a range of grid spacings. Determining the resolution dependence of specific model parameterizations is difficult due to strong resolution dependencies in many pieces of the model. This study presents the Separate Physics and Dynamics Experiment (SPADE) framework that can be used to isolate the resolution dependent behavior of specific parameterizations without conflating resolution dependencies from other portions of the model. To demonstrate the SPADE framework, the resolution dependence of the Morrison microphysics from the Weather Research and Forecasting model and the Morrison-Gettelman microphysics from the Community Atmosphere Model are compared for grid spacings spanning the cloud modeling gray zone. It is shown that the Morrison scheme has stronger resolution dependence than Morrison-Gettelman, and that the ability of Morrison-Gettelman to use partial cloud fractions is not the primary reason for this difference. This study also discusses how to frame the issue of resolution dependence, the meaning of which has often been assumed, but not clearly expressed in the atmospheric modeling community. It is proposed that parameterization resolution dependence can be expressed in terms of "resolution dependence of the first type," RA1, which implies that the parameterization behavior converges towards observations with increasing resolution, or as "resolution dependence of the second type," RA2, which requires that the parameterization reproduces the same behavior across a range of grid spacings when compared at a given coarser resolution. RA2 behavior is considered the ideal, but brings with it serious implications due to limitations of parameterizations to accurately estimate reality with coarse grid spacing. The type of resolution awareness developers should target in their development depends upon the particular modeler’s application.

Gustafson, William I.↗