Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical models”

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

Techno-economic feasibility of borehole thermal energy storage system connected to geothermal heat pumps for seasonal heating load of two buildings in Fairbanks, Alaska

Borehole thermal energy storage (BTES) is an effective solution for managing imbalanced heating and cooling loads in cold regions. This study evaluated the long-term feasibility of a BTES system in the Fairbanks area, Alaska, through building energy modeling, resource characterization, and numerical modeling. The system was designed to store waste heat from a nearby coal power plant during summer and provide thermal energy during winter to geothermal heat pumps supplying heating loads in two buildings. Heating load profiles were modeled for the buildings using EnergyPlus, and the results indicated the annual heating load was 5.6 times greater than the cooling load. 40 borehole heat exchangers were pre-designed approximately 100 m away from the two buildings in terms of land availability and regulatorily optimized depth. The 20-year performance of the designed BTES system under two operational scenarios—one with a 5-year preheating period and one without—was numerically modeled using subsurface temperature and properties characterized through the literature review and thermal response tests. Both scenarios demonstrated that the BTES has the capacity to fully cover the heating loads in the two buildings throughout the 20-year lifetime. Production temperatures at central wells were 33 % higher on average than at outer wells in both scenarios. The 5-year preheating period increased subsurface and extraction temperatures, and correspondingly annual average and total thermal energy production was higher for 8 years than in the scenario without the preheating period. These results highlight the long-term reliability and sustainability of the BTES system in meeting heating demands over its lifetime, with the preheating period offering potential performance improvements. Implementing the BTES system in cold regions with high heating demand, such as Fairbanks, Alaska, could provide a long-term, sustainable energy solution for managing imbalanced heating and cooling loads.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Scalable foundation models for numerical simulations on HPC platforms

In recent years, foundation models (FMs) have begun to reshape numerical simulations on high-performance computing (HPC) platforms. These large, pre-trained AI models enable rapid predictions across a broad range of physical domains, including Earth system modeling, fluid dynamics, materials science, as well as complex multi-modal simulations in aerospace engineering and fusion research. By training on diverse datasets, FMs learn intricate relationships and underlying physical behavior while also enabling the quantification of uncertainty in their predictions. This capability allows simulations that once required days of numerical calculation to be completed in minutes (FM inference), supporting real-time design optimization, uncertainty-aware decision making, and more comprehensive exploration of complex scenarios.

AI↗
DOE OSTI2025

A two-fluid solar-wind model with intermittent Alfvénic turbulence

In one of the leading theories for the origin of the solar wind, photospheric motions launch Alfvén waves (AWs) that propagate along open magnetic-field lines through the solar atmosphere and into the solar wind. The radial variation in the Alfvén speed causes some of the AWs to reflect, and counter-propagating AWs subsequently interact to produce Alfveńic turbulence, in which AW energy cascades from long wavelengths to short wavelengths and dissipates, heating the plasma. In this paper we develop a one-dimensional two-fluid solar-wind model that includes Alfvénic turbulence, proton temperature anisotropy and a novel method for apportioning the turbulent heating rate between parallel proton heating, perpendicular proton heating and electron heating. We employ a turbulence model that accounts for recent observations from NASA’s Parker Solar Probe, which find that AW fluctuations in the near-Sun solar wind are intermittent and less anisotropic than in previous models of anisotropic magnetohydrodynamic turbulence. Our solar-wind model reproduces a wide range of remote observations of the corona and in-situ measurements of the solar wind, and our turbulent heating model consists of analytic equations that could be usefully incorporated into other solar-wind models and numerical models of more distant astrophysical plasmas.

astrophysical plasmas↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2024

Recent Progress on Surface Water Quality Models Utilizing Machine Learning Techniques

Surface waterbodies are heavily exposed to pollutants caused by natural disasters and human activities. Empowering sensor technologies in water quality monitoring, sufficient measurements have become available to develop machine learning (ML) models. Numerous ML models have quickly been adopted to predict water quality indicators in various surface waterbodies. This paper reviews 78 recent articles from 2022 to October 2024, categorizing water quality models utilizing ML into three groups: Point-to-Point (P2P), which estimates the current target value based on other measurements at the same time point; Sequence-to-Point (S2P), which utilizes previous time series data to predict the target value at one time point ahead; and Sequence-to-Sequence (S2S), which uses previous time series data to forecast sequential target values in the future. The ML models used in each group are classified and compared according to water quality indicators, data availability, and model performance. Widely used strategies for improving performance, including feature engineering, hyperparameter tuning, and transfer learning, are recognized and described to enhance model effectiveness. The interpretability limitations of ML applications are discussed. This review provides a perspective on emerging ML for surface water quality models.

machine learning (ML)↗
NASA NTRS1985

Diagnostics of severe convection and subsynoptic scale ageostrophic circulations

Diagnostics of severe convection and subsynoptic scale ageostrophic circulations are reported. Mesoscale circulations through forcing of ageostrophic motion by adiabatic, diabatic and frictional processes were studied. The development and application of a hybrid isentropic sigma coordinate numerical model was examined. The numerical model simulates mesoscale ageostrophic circulations associated with propagating jet streaks and severe convection. A complete list of publications and these completed through support of the NASA severe storms research project is included.

Source record↗
NASA NTRS1989

Low level remote sensing: Orographic winds

The orographic flow data set was obtained from a flight program to measure the influence of orographic features on turbulence momentum, heat, and moisture fluxes. The NASA B-57 aircraft instrumented with probes for measuring the three fluctuating wind speed components, temperature, and humidity was the primary measuring vehicle. Ancillary measurements were made with several ground-based sensors. These include the NOAA radar wind profilers, the Boulder wind network, the PROFS mesoscale surface network, the Boulder Atmospheric Observatory 300 m tower, special rawinsonde observations, and the NOAA/WPL Doppler lidar. The major objective of the flight program was to provide planetary boundary layer parameter information for new and current general circulation computer models. A numerical code, WINDER, based on a discrete element technique was run to numerically model the water tunnel simulated flow. Comparison of the analytical model with the experimental results is very good. Physical fluid dynamic principles embedded in the computational model and visual and hot wire anemometer measurements from the simulation are being rationalized to develop a physical understanding of the vortex flow. The results will be interpreted as they pertain to full scale atmospheric flows.

Frost, Walter↗
NASA NTRS1993

The 24 GHz measurements of 2.2 lambda conical horn antennas illuminating a conducting sheet

Monostatic reflection-coefficient magnitude, absolute value of Gamma, measurements occurring between a radiating horn and a metal reflecting plate are presented for a family of three 2.2 lambda diameter conical horn antennas. The three horns have different aperture phase deviations: 6 deg, 22.5 deg, and 125 deg. Measurements of the magnitude of absolute value of Gamma as a function of horn-plate separation (d) extend from an effective antenna aperture short (d = O) to beyond the far-field boundary (d = 2D(sup 2)/lambda, where D is the antenna diameter). Measurement data are presented with various physical environments for each of the horns. Measured scalar data are compared with theoretical data from two models, a numerical model for a circular waveguide aperture in a ground plane and a scalar diffraction theory model. This work was conducted in support of the development effort for a spaceborne multifrequency microwave reflectometer designed to accurately determine the distance from a space vehicle's surface to a reflecting plasma boundary. The metal reflecting plate was used to simulate the RF reflectivity of a critically dense plasma. The resulting configuration, a ground plane mounted aperture facing a reflecting plane in close proximity, produces a strong interaction between the ground plane and the reflecting plate, especially at integral half-wavelength separations. The transition coefficient is characterized by large amplitude variations.

Cross, A. E.↗
NASA NTRS1993

Research on Earth's rotation and the effect of atmospheric pressure on vertical deformation and sea level variability

The work done under NASA grant NAG5-485 included modelling the deformation of the earth caused by variations in atmospheric pressure. The amount of deformation near coasts is sensitive to the nature of the oceanic response to the pressure. The PSMSL (Permanent Service for Mean Sea Level) data suggest the response is inverted barometer at periods greater than a couple months. Green's functions were constructed to describe the perturbation of the geoid caused by atmospheric and oceanic loading and by the accompanying load-induced deformation. It was found that perturbation of up to 2 cm are possible. Ice mass balance data was used for continental glaciers to look at the glacial contributions to time-dependent changes in polar motion, the lod, the earth's gravitational field, the position of the earth's center-of-mass, and global sea level. It was found that there can be lateral, non-hydrostatic structure inside the fluid core caused by gravitational forcing from the mantle, from the inner core, or from topography at the core/mantle or inner core/outer core boundaries. The nutational and tidal response of a non-hydrostatic earth with a solid inner core was modeled. Monthly, global tide gauge data from PSMSL was used to look at the 18.6-year ocean tide, the 14-month pole tide, the oceanic response to pressure, the linear trend and inter-annual variability in the earth's gravity field, the global sea level rise, and the effects of post glacial rebound. The effects of mantle anelasticity on nutations, earth tides, and tidal variation in the lod was modeled. Results of this model can be used with Crustal Dynamics observations to look at the anelastic dissipation and dispersion at tidal periods. The effects of surface topography on various components of crustal deformation was also modeled, and numerical models were developed of post glacial rebound.

Wahr, John↗
NASA NTRS1999

Solidification of II-VI Compounds in a Rotating Magnetic Field

This project is aimed at using a rotating magnetic field (RMF) to control fluid flow and transport during directional solidification of elemental and compound melts. Microgravity experiments have demonstrated that small amounts of residual acceleration of less than a micro-g can initiate and prolong fluid flow, particularly when there is a static component of the field perpendicular to the liquid solid interface. Thus a true diffusion boundary layer is not formed, and it becomes difficult to verify theories of solidification or to achieve diffusion controlled solidification. The RMF superimposes a stirring effect on an electrically conducting liquid, and with appropriate field strengths and frequencies, controlled transport of material through a liquid column can be obtained. As diffusion conditions are precluded and complete mixing conditions prevail, the technique is appropriate for traveling solvent zone or float zone growth methods in which the overall composition of the liquid can be maintained throughout the growth experiment. Crystals grown by RMF techniques in microgravity in previous, unrelated missions have shown exceptional properties. The objective of the project is two-fold, namely (1) using numerical modeling to simulate the behavior of a solvent zone with applied thermal boundary conditions and demonstrate the effects of decreasing gravity levels, or an increasing applied RMF, or both, and (2) to grow elements and II-VI compounds from traveling solvent zones both with and without applied RMFs, and to determine objectively how well the modeling predicts solidification parameters. Numerical modeling has demonstrated that, in the growth of CdTe from a tellurium solution, a rotating magnetic field can advantageously modify the shape of the liquid solid interface such that the interface is convex as seen from the liquid. Under such circumstances, the defect structure is reduced as any defects which are formed tend to grow out and not propagate. The flow of liquid, however, is complex due to the competing flow induced by the rotating magnetic field and the buoyancy driven convection. When the acceleration forces are reduced to one thousandth of gravity, the flow pattern is much simplified and well controlled material transport through the solvent zone can be readily achieved. Triple axis diffractometry and x-ray synchrotron topography have demonstrated that there is no significant improvement in crystal quality for HgCdTe grown on earth from a tellurium solution when a rotating magnetic field is applied. However, modeling shows that the flow in microgravity with a rotating magnetic field would produce a superior product.

Gillies, D. C.↗
NASA NTRS1998

Volume Diffusion Growth Kinetics and Step Geometry in Crystal Growth

The role of step geometry in two-dimensional stationary volume diff4sion process used in crystal growth kinetics models is investigated. Three different interface shapes: a) a planar interface, b) an equidistant hemispherical bumps train tAx interface, and c) a train of right angled steps, are used in this comparative study. The ratio of the super-saturation to the diffusive flux at the step position is used as a control parameter. The value of this parameter can vary as much as 50% for different geometries. An approximate analytical formula is derived for the right angled steps geometry. In addition to the kinetic models, this formula can be utilized in macrostep growth models. Finally, numerical modeling of the diffusive and convective transport for equidistant steps is conducted. In particular, the role of fluid flow resulting from the advancement of steps and its contribution to the transport of species to the steps is investigated.

Mazuruk, Konstantin↗
NASA NTRS2000

Analysis of Post-Support and Wind-Tunnel Wall Interference on Flow Field About Subsonic High-Lift High-Speed Research Configuration

The present study was performed to determine how significant the interference effects of the wind-tunnel model support system and tunnel walls can be for a high-speed configuration during takeoff and landing conditions. A 5-percent scale model of the Technology Concept Airplane was recently tested in the Langley 14- by 22-Foot Sub-sonic Tunnel. The model was numerically modeled with and without the support and tunnel walls and compared with experimental data. Detailed analysis of the flow provided additional insight concerning what effects the post support and tunnel walls had on the flow field. This study revealed that although the overall forces and moments could be experimentally accounted for, the detailed flow features, such as the surface pressure distributions, could not be accurately simulated without including the post support in the computations.

Lessard, Wendy B.↗
NASA NTRS2009

Numerical Analysis of Film Cooling at High Blowing Ratio

Computational Fluid Dynamics is used in the analysis of a film cooling jet in crossflow. Predictions of film effectiveness are compared with experimental results for a circular jet at blowing ratios ranging from 0.5 to 2.0. Film effectiveness is a surface quantity which alone is insufficient in understanding the source and finding a remedy for shortcomings of the numerical model. Therefore, in addition, comparisons are made to flow field measurements of temperature along the jet centerline. These comparisons show that the CFD model is accurately predicting the extent and trajectory of the film cooling jet; however, there is a lack of agreement in the near-wall region downstream of the film hole. The effects of main stream turbulence conditions, boundary layer thickness, turbulence modeling, and numerical artificial dissipation are evaluated and found to have an insufficient impact in the wake region of separated films (i.e. cannot account for the discrepancy between measured and predicted centerline fluid temperatures). Analyses of low and moderate blowing ratio cases are carried out and results are in good agreement with data.

El-Gabry, Lamyaa↗
NASA NTRS2013

Large/Complex Antenna Performance Validation for Spaceborne Radar/Radiometeric Instruments

Over the past decade, Earth observing missions which employ spaceborne combined radar & radiometric instruments have been developed and implemented. These instruments include the use of large and complex deployable antennas whose radiation characteristics need to be accurately determined over 4 pisteradians. Given the size and complexity of these antennas, the performance of the flight units cannot be readily measured. In addition, the radiation performance is impacted by the presence of the instrument's service platform which cannot easily be included in any measurement campaign. In order to meet the system performance knowledge requirements, a two pronged approach has been employed. The first is to use modeling tools to characterize the system and the second is to build a scale model of the system and use RF measurements to validate the results of the modeling tools. This paper demonstrates the resulting level of agreement between scale model and numerical modeling for two recent missions: (1) the earlier Aquarius instrument currently in Earth orbit and (2) the upcoming Soil Moisture Active Passive (SMAP) mission. The results from two modeling approaches, Ansoft's High Frequency Structure Simulator (HFSS) and TICRA's General RF Applications Software Package (GRASP), were compared with measurements of approximately 1/10th scale models of the Aquarius and SMAP systems. Generally good agreement was found between the three methods but each approach had its shortcomings as will be detailed in this paper.

Antenna Modeling↗
NASA NTRSDate not supplied

From Exploration Flight Test-1 to Artemis II--A NASA Langley's Orion Aerosciences Overview

The Orion Aerosciences program at NASA Langley has played a central role in advancing the aerodynamic and aeroheating prediction capabilities required for the Orion crew vehicle’s return from deep space. This presentation provides a technical overview of aerosciences contributions spanning Exploration Flight Test-1 (EFT-1), Artemis I, and the ongoing post-flight analysis of Artemis II. EFT-1 provided the first high-energy entry dataset for Orion, enabling critical validation of aerodynamic force and moment predictions, static and dynamic stability characteristics, and aeroheating environments at relevant flight Mach and Reynolds numbers. Flight-derived pressure data were used to refine the Flush Air Data System (FADS) methodology for atmospheric density reconstruction and to improve Best Estimated Trajectory (BET) solutions. The EFT-1 data also offered key insights into heat shield performance, including char layer recession, in-depth thermal response, and material retention behavior under flight conditions, informing updates to both thermal response models and uncertainty quantification practices. Building on EFT-1, Artemis I extended the database to true lunar-return conditions. Observations of heat shield performance, including localized char loss, bondline response, and recession variability, provided an unprecedented opportunity to reassess Thermal Protection System (TPS) and aeroheating modeling assumptions. Aerodynamic reconstruction efforts incorporated improved FADS calibration, enhanced atmospheric modeling, and refined force and moment databases to reduce trajectory and load uncertainties. Aeroheating comparisons between pre-flight predictions and flight data enabled targeted model updates, particularly in transitional flow environments and wake heating regions. For Artemis II, these lessons were systematically incorporated into the pre-flight prediction process. Updates included refined aerodynamic databases anchored to flight-validated corrections, improved density estimation and BET methodologies using enhanced database interpolation algorithm and FADS modeling, and revised aeroheating design environments informed by Artemis I material response observations. By the time of the workshop, Artemis II post-flight analysis will be underway, and preliminary findings will be presented where available, including early comparisons of aerodynamic reconstruction, atmospheric density estimation, and thermal protection system performance relative to updated predictions. Collectively, this body of work is a testament to the dedicated and multidisciplinary team whose sustained efforts have contributed to the program’s success and to the progressive maturation of Orion aerosciences modeling through numerical modeling, ground and flight data assimilation. The integrated advancement of aerodynamics, trajectory reconstruction, FADS-based density estimation, and aeroheating analysis has reduced predictive uncertainty and strengthened confidence for future crewed lunar and deep-space missions.

Orion↗
NASA NTRS1989

Development of high resolution simulations of the atmospheric environment using the MASS model

Numerical simulations were performed with a very high resolution (7.25 km) version of the MASS model (Version 4.0) in an effort to diagnose the vertical wind shear and static stability structure during the Shuttle Challenger disaster which occurred on 28 January 1986. These meso-beta scale simulations reveal that the strongest vertical wind shears were concentrated in the 200 to 150 mb layer at 1630 GMT, i.e., at about the time of the disaster. These simulated vertical shears were the result of two primary dynamical processes. The juxtaposition of both of these processes produced a shallow (30 mb deep) region of strong vertical wind shear, and hence, low Richardson number values during the launch time period. Comparisons with the Cape Canaveral (XMR) rawinsonde indicates that the high resolution MASS 4.0 simulation more closely emulated nature than did previous simulations of the same event with the GMASS model.

Kaplan, Michael L.↗
NASA NTRS1993

Computation of flows in a turn-around duct and a turbine cascade using advanced turbulence models

Numerical investigation has been carried out to evaluate the capability of the Algebraic Reynolds Stress Model (ARSM) and the Nonlinear Stress Model (NLSM) to predict strongly curved turbulent flow in a turn-around duct (TAD). The ARSM includes the near-wall damping term of pressure-strain correlation phi(sub ij,w), which enables accurate prediction of individual Reynolds stress components in wall flows. The TAD mean flow quantities are reasonably well predicted by various turbulence models. The ARSM yields better predictions for both the mean flow and the turbulence quantities than the NLSM and the k-epsilon (k = turbulent kinetic energy, epsilon = dissipation rate of k) model. The NLSM also shows slight improvement over the k-epsilon model. However, all the models fail to capture the recovery of the flow from strong curvature effects. The formulation for phi(sub ij,w) appears to be incorrect near the concave surface. The hybrid k-epsilon/ARSM, Chien's k-epsilon model, and Coakley's q-omega (q = the square root of k, omega = epsilon/k) model have also been employed to compute the aerodynamics and heat transfer of a transonic turbine cascade. The surface pressure distributions and the wake profiles are predicted well by all the models. The k-epsilon model and the k-epsilon/ARSM model provide better predictions of heat transfer than the q-omega model. The k-epsilon/ARSM solutions show significant differences in the predicted skin friction coefficients, heat transfer rates and the cascade performance parameters, as compared to the k-epsilon model. The k-epsilon/ARSM model appears to capture, qualitatively, the anisotropy associated with by-pass transition.

Lakshminarayana, B.↗
NASA NTRS1990

Three real-time architectures - A study using reward models

Numerous applications in the area of computer system analysis can be effectively studied with Markov reward models. These models describe the evolutionary behavior of the computer system by a continuous-time Markov chain, and a reward rate is associated with each state. In reliability/availability models, upstates have reward rate 1, and down states have reward rate zero associated with them. In a combined model of performance and reliability, the reward rate of a state may be the computational capacity, or a related performance measure. Steady-state expected reward rate and expected instantaneous reward rate are clearly useful measures which can be extracted from the Markov reward model. The diversity of areas where Markov reward models may be used is illustrated with a comparative study of three examples of interest to the fault tolerant computing community.

Sjogren, J. A.↗