Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average model”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

On the Quality of the Nimbus 7 LIMS Version 6 Ozone for Studies of the Middle Atmosphere

The Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) radiance profile dataset of 1978/79 was reconditioned and reprocessed to Version 6 (V6) profiles of temperature and species that are improved significantly over those from Version 5 (V5). The LIMS V6 dataset was archived for public use in 2002. Improvements for its ozone include: (1) a more accurate accounting for instrument and spacecraft motion effects in the radiances, (2) the use of better spectroscopic line parameters for its ozone forward model, (3) retrievals of all its scans, (4) more accurate and compatible temperature versus pressure profiles (or T(p)), which are needed for the registration of the ozone radiances and for the removal of temperature effects from them, and (5) a better accounting for interfering species in the lower stratosphere. The retrieved V6 ozone profiles extend from near cloud top altitudes to about 80 km and from 64S to 84N latitude with better sampling along the orbit than for the V5 dataset. Calculated estimates of the single-profile precision and accuracy are provided for the V6 ozone from this study. Precision estimates based on the data themselves are of order 3% or better from 1 to 30 hPa. Estimates of total systematic error for a single profile are hard to generalize because the separate sources of error may not all be of the same sign and they depend somewhat on the atmospheric state. It is estimated that the V6 zonal mean ozone distributions are accurate to within 9% to 7% from 50 hPa to 3 hPa, respectively. Effects of a temperature bias can be significant and may be present at 1 to 2 hPa though. There may be ozone biases of order 10% at those levels due to possible biases of up to +2 K, but there is no indication of a similar problem elsewhere in the stratosphere. Simulation studies show that the LIMS retrievals are also underestimating slightly the small amplitudes of the atmospheric temperature tides, which affect its retrieved day/night ozone differences. There are small biases in the middle to lower stratosphere for the ascending versus descending node LIMS ozone, due principally to not accounting for the asymmetric weighting of its radiances across the tangent layer. The estimates of total accuracy were assessed by comparing the daily zonal mean LIMS ozone distributions against those from the Nimbus 7 SBUV Version 8 (V8) dataset for the same period. Generally, the LIMS V6 ozone agrees well with SBUV, except perhaps in the tropical lower stratosphere where the LIMS ozone is less. Still, the accuracy for LIMS V6 ozone in the lower stratosphere is improved over that found for LIMS V5, as indicated by several LIMS comparisons with ECC ozonesonde profiles. The LIMS V6 ozone is considered especially suitable for detailed studies of large-scale stratospheric processes above the 100-hPa level. Comparison of diurnal, photochemical model calculations with the monthly-averaged, upper stratospheric ozone observed from LIMS V6 indicates only a slight ozone deficit for the model at about 2 hPa. However, that deficit exhibits little to no seasonal variation and is in good agreement with similar model comparisons for a seasonal time series of ozone obtained with ground-based microwave instruments. Because the LIMS V6 ozone in the lower stratosphere has improved accuracy and sampling versus that of V5, it should now be possible to conduct quantitative studies of ozone transport and chemistry for the northern hemisphere, polar stratospheric winter of 1978/79 a time period when the catalytic loss of ozone due to reactive chlorine should not have been a major factor for the Arctic region.

Remsberg, Ellis↗
NASA NTRS2019

The Badhwar-O'Neill 2020 Model

The Badhwar-O'Neill (BON) model has been used for some time to describe the galactic cosmic ray (GCR) environment encountered in deep space by astronauts and sensitive electronics. The most recent version of the model, BON2014, was calibrated to available measurements to reduce model errors for particles and energies of significance to astronaut exposure. Although subsequent studies showed the model to be reasonably accurate for such applications, modifications to the sunspot number (SSN) classification system and a large number of new high precision measurements suggested the need to develop an improved and more capable model. In this work, the BON2020 model is described. The new model relies on daily integral flux from the Advanced Composition Explorer Cosmic Ray Isotope Spectrometer (ACE/CRIS) to describe solar activity. For time periods not covered by ACE/CRIS, the updated international SSN database is used. Parameters in the new model are calibrated to available data, which includes the new Alpha Magnetic Spectrometer (AMS-02) and Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) high-precision measurements. It is found that the BON2020 model is a significant improvement over BON2014. Systematic errors associated with BON2014 have been removed. The average relative error of the BON2020 model compared to all available measurements is found to be <1%, and BON2020 is found to be within ±15% of a large fraction of the available measurements (26,269 of 27,646 → 95%).

Slaba, Tony C.↗
NASA NTRS2020

The Badhwar-O'Neill 2020 Model

The Badhwar-O’Neill (BON) model has been used for some time to describe the galactic cosmic ray (GCR) environment encountered in deep space by astronauts and sensitive electronics. The most recent version of the model, BON2014, was calibrated to available measurements to reduce model errors for particles and energies of significance to astronaut exposure. Although subsequent studies showed the model to be reasonably accurate for such applications, modifications to the sunspot number (SSN) classification system and a large number of new high-precision measurements suggested the need to develop an improved and more capable model. In this work, the BON2020 model is described. The new model relies on daily integral flux from the Advanced Composition Explorer Cosmic Ray Isotope Spectrometer (ACE/CRIS) to describe solar activity. For time periods not covered by ACE/CRIS, the updated international SSN database is used. Parameters in the new model are calibrated to available data, which includes the new Alpha Magnetic Spectrometer (AMS-02) and Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) high-precision measurements. It is found that the BON2020 model is a significant improvement over BON2014. Systematic errors associated with BON2014 have been removed. The average relative error of the BON2020 model compared to all available measurements is found to be <1%, and BON2020 is found to be within +15% of a large fraction of the available measurements (26,269 of 27,646  95%.

T C Slaba↗
DOE OSTI2026

TCAD-Machine Learning Enabled TID Compact Model Development for Commercial SiC MOSFET

We propose a TCAD (Technology Computer Aided Design)-machine learning coupled approach that combines a TCAD tool (Charon), optimization/uncertainty quantification tool (Dakota), surrogate models, and Bayesian learning capabilities. The coupling approach is used for accurate modeling and calibration of total ionizing dose (TID) induced threshold voltage (V th ) shifts in Commercial-Off-The-Shelf (COTS) semiconductor devices and to develop physics-informed TID compact models. This versatile approach is applied to model the TID effect in an exemplar COTS 3.3 kV SiC power MOSFET (Metal-Oxide-Semiconductor Field-Effect Transistor). With the Charon-Dakota coupling, we can determine key device geometry and doping values based on device physics, which are difficult to obtain or not available for COTS devices but important for TCAD simulation; additionally, we can efficiently generate thousands of simulation results in a large parameter space, which makes it possible to develop data-driven surrogate models and perform Bayesian calibration. Utilizing the full tool-coupling approach, we achieve calibrated TCAD simulation models that accurately capture the average TID-induced V th shifts behavior with total doses and V th shifts saturation at high doses as observed in experimental data. More importantly, the calibrated TCAD simulations are obtained with determined TID model parameters (e.g., hole trap density and capture cross section) values that contain well quantified uncertainties. Furthermore, we can isolate and quantify the noises that are not captured by the TCAD models but exist in the measured data due to measurements and devices variabilities. Lastly, the calibrated surrogate models are used to develop physics-informed TID compact models. The method is generalizable to other devices and/or radiation conditions with few modifications and can provide well-determined uncertainties.

COTS↗
DOE OSTI2024

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

A Correction for the IRI Topside Electron Density Model Based on Alouette/ISIS Topside Sounder Data

The topside segment of the International Reference Ionosphere (IRI) electron density model (and also of the Bent model) is based on the limited amount of topside data available at the time (40,OOO Alouette 1 profiles). Being established from such a small database it is therefore not surprising that the models have well-known shortcomings, for example, at high solar activities. Meanwhile a large data base of close to 200,000 topside profiles from Alouette 1,2, and ISIS I, 2 has become available online. A program of automated scaling and inversion of a large volume of digitized ionograms adds continuously to this data pool. We have used the currently available ISIs/Alouette topside profiles to evaluate the IRI topside model and to investigate ways of improving the model. The IRI model performs generally well at middle latitudes and shows discrepancies at low and high latitudes and these discrepancies are largest during high solar activity. In the upper topside IRI consistently overestimates the measurements. Based on averages of the data-model ratios we have established correction factors for the IRI model. These factors vary with altitude, modified dip latitude, and local time.

Bilitza, D.↗
NASA NTRSDate not supplied

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.

turbulent boundary layers↗
DOE OSTI2026

Two-population Rouse models for polymer segmental dynamics in nanocomposites

Segmental dynamics of polymer chains in a model nanocomposite of poly(ethylene oxide) and silica nanoparticles (NPs) was investigated using quasielastic neutron scattering. The dynamics can be accurately described with the Rouse model. The bulklike polymer far from the NP surface behaves as the neat polymer. However, the slower polymer in the interface close to the NP surface is described either with a second Rouse population with different relaxation time or using the suppressed Rouse model. These simple two-population models accurately reproduce the experimental data, with the suppressed Rouse model describing topological constraints, on average, every 12 beads with an interfacial thickness up to 13.5 beads, and the effect of the interface extending to a layer of thickness comparable to the end-to-end distance of the polymer. This modeling provides an explanation for the observed reinforcement in PNC even at low loadings, consistent with current understanding of the relevance of the interphase.

Composite materials↗
NASA NTRS1973

A radiation model for calculating atmospheric corrections to remotely sensed infrared measurements, version 2

A numerical model is developed which calculates the atmospheric corrections to infrared radiometric measurements due to absorption and emission by water vapor, carbon dioxide, and ozone. The corrections due to aerosols are not accounted for. The transmissions functions for water vapor, carbon dioxide, and water are given. The model requires as input the vertical distribution of temperature and water vapor as determined by a standard radiosonde. The vertical distribution of carbon dioxide is assumed to be constant. The vertical distribution of ozone is an average of observed values. The model also requires as input the spectral response function of the radiometer and the nadir angle at which the measurements were made. A listing of the FORTRAN program is given with details for its use and examples of input and output listings. Calculations for four model atmospheres are presented.

Boudreau, R. D.↗
NASA NTRS1980

Morphology of the Venus ultraviolet night airglow

Images of the nightside of Venus in the (0,1)delta band of nitric oxide have been obtained by the Pioneer Venus orbiter ultraviolet spectrometer (OUVS). The emission, which is produced by radiative association of N and O, shows a bright spot reaching 5 kR and located at 2 a.m. local solar time just south of the equator. The emitting layer is at 111 + or - 7-km altitude. A one dimensional vertical transport model shows that the hemispheric average brightness of 0.8 kR is consistent with the orbiter neutral mass spectrometer (ONMS) measurements of N and O near 167 km, and that the altitude of the emitting layer is consistent with the eddy mixing model proposed to explain the dayside helium profile measured by the bus neutral mass spectrometer. In the model, N reaches a peak of 7 x 10 to the 8th per cu cm at 114 km, and O reaches a peak of 2.6 x 10 to the 11th per cu cm at 106 km. There is a fair degree of consistency between the ONMS, OUVS, and other airglow measurements, except as regards the local time dependence.

Stewart, A. I. F.↗
DOE OSTI2022

FLOW Estimation and Rose Superposition (FLOWERS): an integral approach to engineering wake models

Abstract. Annual energy production (AEP) is often the objective function in wind plant layout optimization studies. The conventional method to compute AEP for a wind farm is to first evaluate power production for each discrete wind direction and speed using either computational fluid dynamics simulations or engineering wake models. The AEP is then calculated by weighted-averaging (based on the wind rose at the wind farm site) the power produced across all wind directions and speeds. We propose a novel formulation for time-averaged wake velocity that incorporates an analytical integral of a wake deficit model across every wind direction. This approach computes the average flow field more efficiently, and layout optimization is an obvious application to exploit this benefit. The clear advantage of this new approach is that the layout optimization produces solutions with comparable AEP performance yet is completed 2 orders of magnitude faster. The analytical integral and the use of a Fourier expansion to express the wind speed and wind direction frequency create a relatively smooth solution space for the gradient-based optimizer to excel in comparison to the existing weighted-averaging power calculation.

17 WIND ENERGY↗
DOE OSTI2022

Spatial patterns of snow distribution in the sub-Arctic

Abstract. The spatial distribution of snow plays a vital role in sub-Arctic and Arctic climate, hydrology, and ecology due to its fundamental influence on the water balance, thermal regimes, vegetation, and carbon flux. However, the spatial distribution of snow is not well understood, and therefore, it is not well modeled, which can lead to substantial uncertainties in snow cover representations. To capture key hydro-ecological controls on snow spatial distribution, we carried out intensive field studies over multiple years for two small (2017–2019; ∼ 2.5 km2) sub-Arctic study sites located on the Seward Peninsula of Alaska. Using an intensive suite of field observations (> 22 000 data points), we developed simple models of the spatial distribution of snow water equivalent (SWE) using factors such as topographic characteristics, vegetation characteristics based on greenness (normalized different vegetation index, NDVI), and a simple metric for approximating winds. The most successful model was random forest, using both study sites and all years, which was able to accurately capture the complexity and variability of snow characteristics across the sites. Approximately 86 % of the SWE distribution could be accounted for, on average, by the random forest model at the study sites. Factors that impacted year-to-year snow distribution included NDVI, elevation, and a metric to represent coarse microtopography (topographic position index, TPI), while slope, wind, and fine microtopography factors were less important. The characterization of the SWE spatial distribution patterns will be used to validate and improve snow distribution modeling in the Department of Energy's Earth system model and for improved understanding of hydrology, topography, and vegetation dynamics in the sub-Arctic and Arctic regions of the globe.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

An investigation of turbulence models

The accuracy to which a turbulent boundary layer or wake can be predicted numerically depends on the validity of the turbulence closure model used. The modeling of turbulence physics is one of the most difficult problems in computational fluid dynamics (CFD). In fact, it is one of the pacing factors in the development of CFD. In general, there are three main approaches to the description of trubulence physics. First is turbulence modeling in which the Reynolds averaged Navier-Stokes equations are used and some closure approximation is made for the Reynolds stresses. A second approach to turbulence is large eddy simulation (LES) in which the computational mesh is taken to be fine enough that the large scale structure of the turbulence can be calculated directly. An empirical assumption must be made for the small scale sub-grid turbulence. The third approach is direct simulation. In this technique the Navier-Strokes equations are solved directly on a mesh which if fine enough to resolve the smallest length scale of the turbulence. The Reynolds averaged equations are not used and no closure assumption is required. These last two approaches require extensive computer resources and as such are not engineering tools. The purpose of the work was to investigate the various engineering turbulence models for accuracy and ease of programming. This involved comparison of the models with each other and with experimental data.

Doria, Michael L.↗
DOE OSTI2022

Method to account for natural fracture induced elastic anisotropy in geomechanical characterization of shale gas reservoirs

Shale has been usually recognized as a transverse isotropic (TI) medium in conventional geomechanical log interpretation due to its laminated nature. However, when natural fractures exist in the shale rock, additional elastic anisotropy is introduced, converting laminated Shale to an orthorhombic (OB) medium. Previous studies illustrate that neglecting the natural fracture induced anisotropy in shale geomechanical log interpretation could lead to inaccurate evaluations of elastic moduli and in-situ stresses. In this paper, a new method is developed to account for the natural fracture induced anisotropy in geomechanical log interpretation based upon the TI acoustic model developed by the author and a characterization technique of elastic wave anisotropy (Sayers, 1991). The new OB model incorporates the four acoustic log data inputs and five modeling constraints in a nonlinear optimization algorithm to solve for the nine independent stiffness coefficients of an OB rock, and further to solve for the geomechanical properties and in-situ stress profiles in an OB formation. The new method was validated with a Marcellus Gas Shale field case. Both the new OB model and the conventional TI model were applied to interpret the minimum horizontal stress profile for the same formation. By comparing the results, the OB model is more robust than the TI from two aspects. First, the average stress magnitude predicted by the OB model is closer to the one measured by the Diagnostic Fracture Injection Test (DFIT). Second, the OB model predicts a more obvious stress barrier between the lower Marcellus and upper Onondaga Limestone than the TI model does. Finally, the predicted stress barrier is consistent with the observation of the microseismic events of a horizontal well drilled and completed nearby, which reveals that no hydraulic fracture propagates downward through the bottom boundary of Marcellus Shale into the underlying Onondaga Limestone.

03 NATURAL GAS↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

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

leray regularization↗
NASA NTRS1986

The dynamics of warm and cold climates

Changes occurring as a result of climatic warming or cooling are presently investigated by comparing the atmospheric dynamics of five different climate simulations encompassing two ice ages, current and doubled CO2 climates, and a simulation of the warm Cretaceous; a global average surface air temperature range of 13 C is covered by the models. It is noted that many zonally averaged processes do not change systematically as climate changes: the January Hadley cell, the jet stream, mean precipitation patterns, and total atmospheric transport, show very little variation among the different climate simulations. These features result from balances between competing factors; while individual processes differ in the cold and warm climates, there is considerable compensation.

Rind, D.↗