Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling & Simulation (M&S)”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Simulated interannual variability of the Greenland Sea deep water formation and its connection to surface forcing

A fully prognostic Arctic ice-ocean model is used to study the interannual variability of deepwater formation in the Greenland Sea Gyre based on the simulations for the Arctic ice-ocean system for the period 1955 and 1960 - 1985. The model uses monthly climatology for thermodynamic forcing components (such as air temperature and cloudiness), together with constant annual net precipitation and river runoff. The daily wind forcing is derived from analyzed sea level air pressures from the National Center for Atmospheric Research (NCAR). In summary, the model shows that the occurence of deep convection in the Greenland Sea Gyre is controlled by the extensive Fram Strait ice export and/or local wind conditions in the Greenland Sea. In the latter case the weakening of the local wind curl allows the Polar Front to move eastward. The movement of the Polar Front causes adverse ice conditions, often together with much larger than normal ice export from the Arctic, such as in 1968, which can block convection in the gyre. The density difference between upper and lower layers is investigated as an indication of water mass formation through convection, occurring as strong diffusion in the model. The model-simulated density difference between the average top 100 m and deep levels reveals that the period 1960 - 1985 had only a few distinct years with weak stratification, and, especially, the model predicts no deep convection since the nid-1970s. The common factor for the years of the weakest decrease of the model-predicted heat content of the upper 2000 m which can, to a high degree, be explained by local heat loss.

Haekkinen, Sirpa↗
DOE OSTI2022

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H 2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories. This is particularly the case for bentonite, an often-preferred barrier material. Gas flow in preferential pathways is characterized by localization and spontaneous behavior, which is challenging to simulate in numerical models due to strong hydro-mechanical coupling. To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project, this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations. The proposed model utilized hydro-mechanical spatial distributions, namely Young's modulus and gas entry pressure, and elastoplasticity combined with a linear swelling model. A conceptual strain-dependent permeability approach simulated dilatancy-controlled gas flow based on hydro-mechanical coupling. To test the effectiveness of the presented approach, a gas injection test in a compacted, saturated bentonite sample was simulated using the open-source code OpenGeoSys 5.8 and compared with experimental observations. The presented methodology is capable of simulating localized gas flow in preferential pathways. The spatial distributions of Young's modulus and gas entry pressure affect the swelling pressure, relative permeability and, in combination with the strain-dependent permeability model, also the intrinsic permeability.

42 ENGINEERING↗
DOE OSTI2022

Barometric Pumping Through Fractured Rock: A Mechanism for Venting Deep Methane to Mars' Atmosphere

Abstract Both the source of methane on Mars and the mechanism for transmission from the subsurface to the atmosphere are not fully understood. Previous seepage simulations have invoked relatively shallow subsurface sources to explain observed methane signatures on Mars. We propose that barometric‐pressure pumping through fracture networks could be an effective mechanism for methane transport from the deep subsurface on Mars. Using atmospheric pressure data gathered by Curiosity as input, we simulate methane gas transport from depths of 200 m to the surface. Even with such a deep source, our model reproduces the observed seasonality of methane, and the simulated surface methane fluxes fall within the range of previous estimates derived from atmospheric observations. Because 200 m is the likely minimum hospitable depth for living methanogenic microbes, our fracture network model indirectly reinvigorates the possibility of a microbial source of methane on Mars.

58 GEOSCIENCES↗
NASA NTRS1994

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.↗
DOE OSTI2021

Complementary Study of Radiative Heat Transfer and Flow Physics from Moderate-scale Hydrocarbon Pool Fire Simulations

As part of the Advanced Simulation and Computing Verification and Validation (ASCVV) program, a 0.3-m diameter hydrocarbon pool fire with multiple fuels was modeled and simulated. In the study described in this report, systematic examination was performed on the radiation model used in a series of coupled Fuego/Nalu simulations. A calibration study was done with a medium-scale methanol pool fire and the effect of calibration traced throughout the radiation model. This analysis provided a more detailed understanding of the effect of radiation model parameters on each other and on other quantities in the simulations. Heptane simulation results were also examined using this approach and possible areas for further improvement of the models were identified. The effect of soot on radiative losses was examined by comparing heptane and methanol results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

High sensitivity of simulated fog properties to parameterized aerosol activation in case studies from ParisFog

Aerosols influence fog properties such as visibility and lifetime by affecting fog droplet number concentrations (N d ). Numerical weather prediction (NWP) models often represent aerosol–fog interactions using highly simplified approaches. Incorporating prognostic size-resolved aerosol microphysics from climate models could allow them to simulate N d and aerosol–fog interactions without incurring excessive computational expense. However, microphysics code designed for coarse spatial resolution may struggle with sub-kilometer-scale grid spacings. Here, we test the ability of the UK Met Office Unified Model to simulate aerosol and fog properties during case studies from the ParisFog field campaign in 2011. We examine the sensitivity of fog properties to variations in N d caused by modifications to simulated aerosol activation. Our model, with a 500 m horizontal resolution and interactive aerosol and cloud microphysics, significantly underpredicts N d , although it only slightly underestimates the cloud condensation nuclei concentration. With an updated version of the Abdul-Razzak and Ghan (2000) activation scheme, we produce N d that are more consistent with those predicted by a cloud parcel model under fog-like conditions. We activate droplets only by adiabatic cooling. We incorporate more realistic hygroscopicities for sulfate and organic aerosols and explore the sensitivity of simulated N d to unresolved updrafts. We find that both N d and simulated fog liquid water content are very sensitive to the updated activation scheme but remain less affected by the update to hygroscopicities. Our improvements offer insights into the physical processes regulating N d in stable conditions, potentially laying foundations for improved operational fog forecasts that incorporate interactive aerosol simulations or aerosol climatologies.

Ghosh, Pratapaditya [Carnegie Mellon University, P↗
DOE OSTI2026

Simulation budgeting for hybrid effective field theories

In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model. Our accuracy goals, determined with careful consideration of statistical and systematic uncertainties, are 1% accurate in the high-likelihood range of cosmological parameters, and 2% accurate over a broader parameter space volume for k < 1 h Mpc -1 and z < 3. Focusing in part on the 8-parameter w 0 w a CDM+ m ν cosmological model, we find that < 225 simulations are required to meet our error goals over our wide parameter space, including models with rapidly evolving dark energy, given our simulation and emulator recommendations. For a more restricted parameter space volume, as few as 80 simulations are sufficient. We additionally present simulation forecasts for example use cases, and make the code used in our analyses publicly available. These results offer practical guidance for efficient emulator design and simulation budgeting in future cosmological analyses.

cosmological parameters from LSS↗
NASA NTRS2024

Improved Simulations of Biomass Burning Aerosol Optical Properties and Lifetimes in the NASA GEOS Model during the ORACLES-I Campaign

In order to improve aerosol representation in the NASA Goddard Earth Observing System (GEOS) model, we evaluated simulations of the aerosol properties and transport over the southern African biomass burning source and outflow region using observations made during the first deployment of the ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) field campaign in September 2016. An example case study of September 24 was analyzed in detail, during which aircraft-based in-situ and remote sensing observations showed the presence of a multi-layered smoke plume structure with significant vertical variation in single scattering albedo (SSA). Our baseline GEOS simulations were not able to represent the observed SSA variation, nor the observed organic aerosol-to-black carbon ratio (OA : BC). Analyzing the simulated smoke age suggests that the higher altitude, less absorbing smoke plume was younger (~4 days), while the lower altitude and more absorbing smoke plume was older (~7 days). We hypothesize a chemical or microphysical loss process exists to explain the change in aerosol absorption as the smoke plume ages, and we apply a simple 6-day e-folding loss rate to the model hydrophilic biomass burning OA to mimic this process. Adding this loss process required some adjustment to the model assumed scaling factors of aerosol emissions to conserve the regional aerosol loading and further improve the simulated OA : BC ratio. Accordingly, we have increased our biomass burning emissions of OA by 60 % and biomass burning BC by 15 %. We also utilized the ORACLES airborne observations to better constrain the simulation of aerosol optical properties, adjusting the assumed particle size, hygroscopic growth, and absorption. Our final GEOS model simulation with additional OA loss and updated optics showed a better performance in simulating AOD and SSA compared to independent ground and space-based retrievals for the entire month of September 2016, including the OMI Aerosol Index. In terms of radiative implications of our model adjustments, the final GEOS simulation suggested a decreased atmospheric warming of about 10 % (~2 W m -2 ) over the south-east Atlantic region and above the stratocumulus cloud decks compared to the model baseline simulations. These results improve the representation of the smoke age, transport, and optical properties in Earth system models.

aerosols↗
DOE OSTI2026

Performance of reanalysis and mesoscale models off the coast of Hawai'i

The eastern Hawai'i coast in the United States is characterized by considerable wind resource fuelled by persistent trade winds, making it an important area for energy research. The need is strong for reanalyses and higher-resolution regional simulations where observations have been historically limited, such as Hawai'i's offshore environments. However, studies using offshore observations in other parts of the world have shown that significant errors can occur in reanalyses and wind datasets, which can lead to inaccurate estimates of wind energy generation, payback periods, and extreme weather risks at project locations. The degree of such errors is influenced by a number of factors, including spatial resolution and the handling of processes within the planetary boundary layer (PBL). In this work, we provide a wind resource characterization from year-long lidar buoy measurements off the eastern coast of O'ahu, Hawai'i, an environment previously unobserved at the rotor level, and use the characterization to evaluate the performance of two simulation datasets. The O'ahu deployment location is meteorologically unique and less complex than land-based wind resource characterizations, being strongly characterized by trade winds with minimal land–atmosphere interaction influences. Despite the unique and fairly consistent meteorological conditions, we hypothesize that distinct simulation datasets will exhibit diverse ranges of errors similar to those that have been seen for other offshore locations. We find the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) to strongly underestimate observed wind speeds at the O'ahu location (bias = −1.54 m s −1 at a height of 140 m above sea level), while a regional Weather Research and Forecasting Model (WRF) simulation produced by the University of Hawai'i (UH-WRF) provides a significantly smaller wind speed bias (−0.25 m s −1 ), highlighting the value of running regional, higher-resolution simulations. The large bias noted for ERA5 is driven by significant underestimation of fast wind speeds (>9 m s −1 ), which the study site is largely characterized by, along with discontinuities in the ERA5 diurnal cycle. We also speculate that the relative sparsity of observations for data assimilation in this remote part of the world could influence the performance of ERA5 and that challenges with characterizing island effects could impact the performance of both datasets.

17 WIND ENERGY↗
DOE OSTI2024

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Cross-scale modeling of storm surge, tide, and inundation in Mid-Atlantic Bight and New York City during Hurricane Sandy, 2012

Driven by high-resolution NAM (North America Mesoscale) and ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric fields, a 3-D unstructured grid SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model) was applied to simulate total water level in the Mid-Atlantic Bight during Hurricane Sandy in 2012. The simulated storm surge, tide, significant wave height, and wave peak period results compared favorably with NOAA observations along US East Coast, Long Island Sound, and New York Harbor. The maximum total water level at The Battery in the New York Harbor was accurately simulated with an absolute error of less than 0.08 m (out of 2.9 m peak surge) and a timing difference less than 10 min, which includes the effects on the order of 0.1–0.3 m (5–10%) from the coastal setup induced by the surface waves. The scenario comparisons of (1) “NAM” versus “ECMWF”, (2) “2-D” versus “3-D”, and (3) “with” versus “without” wind wave model were examined. The 3-D barotropic model forced by ECMWF including the effects of wind waves performs the best, attributed to wave-induced radiation stress and reduction of bottom stress when the 3-D version is used. Simultaneously, the ELCIRC-Sub model, using a 5-m topography/bathymetry sub-grid and a regular 200-m resolution finite volume computational model grid, was developed to simulate the street-level inundation in New York City. The momentum and mass fluxes calculated by the coarser base grid model were effectively coupled with the sub-grid so that running a full-blown high-resolution base model is not required. The ELCIRC-Sub model uses MPI (Message Passing Interface) parallel computing to enlarge the coverage of the land surface and an efficient non-linear solver to improve the accuracy of the wetting-and-drying processes. The temporal comparisons of modeled water level with NOAA's tidal stations and USGS' rapid-deployed gauges showed overall performance with an average error of 0.1 m. Particularly, ELCIRC-Sub captured the profile of the highest peak surge (3.9 m), which rose 2 m within 3 hour at Kings Point, NY. The spatial comparisons of the modeled peak water level at 80 surveyed locations showed an average error less than 0.13 m. The modeled maximum inundation extent also matched well with an 80% hit rate against FEMA's Hurricane Sandy maximum flooding extent.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Artificial stereo presentation of meteorological data fields

The innate capability to perceive three-dimensional stereo imagery has been exploited to present multidimensional meteorological data fields. Variations on an artificial stereo technique first discussed by Pichel et al. (1973) are used to display single and multispectral images in a vivid and easily assimilated manner. Examples of visible/infrared artificial stereo are given for Hurricane Allen and for severe thunderstorms on 10 April 1979. Three-dimensional output from a mesoscale model also is presented. The images may be viewed through the glasses inserted in the February 1981 issue of the Bulletin of the American Meteorological Society, with the red lens over the right eye. The images have been produced on the interactive Atmospheric and Oceanographic Information Processing System (AOIPS) at Goddard Space Flight Center. Stereo presentation is an important aid in understanding meteorological phenomena for operational weather forecasting, research case studies, and model simulations.

Hasler, A. F.↗
NASA NTRS2010

Premonsoon Aerosol Characterization and Radiative Effects Over the Indo-Gangetic Plains: Implications for Regional Climate Warming

The Himalayas have a profound effect on the South Asian climate and the regional hydrological cycle, as it forms a barrier for the strong monsoon winds and serves as an elevated heat source, thus controlling the onset and distribution of precipitation during the Indian summer monsoon. Recent studies have suggested that radiative heating by absorbing aerosols, such as dust and black carbon over the Indo-Gangetic Plains (IGP) and slopes of the Himalayas, may significantly accelerate the seasonal warming of the Hindu Kush-Himalayas-Tibetan Plateau (HKHT) and influence the subsequent evolution of the summer monsoon. This paper presents a detailed characterization of aerosols over the IGP and their radiative effects during the premonsoon season (April-May-June) when dust transport constitutes the bulk of the regional aerosol loading, using ground radiometric and spaceborne observations. During the dust-laden period, there is a strong response of surface shortwave flux to aerosol absorption indicated by the diurnally averaged forcing efficiency of -70 W/sq m per unit optical depth. The simulated aerosol single-scattering albedo, constrained by surface flux and aerosol measurements, is estimated to be 0.89+/- 0.01 (at approx.550 nm) with diurnal mean surface and top-of-atmosphere forcing values ranging from -11 to -79.8 W/sq m and +1.4 to +12 W/sq m, respectively, for the premonsoon period. The model-simulated solar heating rate profile peaks in the lower troposphere with enhanced heating penetrating into the middle troposphere (5-6 km), caused by vertically extended aerosols over the IGP with peak altitude of approx.5 km as indicated by spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization observations. On a long-term climate scale, our analysis, on the basis of microwave satellite measurements of tropospheric temperatures from 1979 to 2007, indicates accelerated annual mean warming rates found over the Himalayan-Hindu Kush region (0.21 C/decade+/-0.08 C/decade) and underscores the potential role of enhanced aerosol solar absorption in the maximum warming localized over the western Himalayas (0.26 C/decade f 0.09 C/decade) that significantly exceed the entire HKHT and global warming rates. We believe the accelerated warming rates reported here are critical to both the South Asian summer monsoon and hydro-glaciological resource variability in the Himalayan-Hindu Kush snowpack and therefore to the densely populated downstream regions.

Gautam, Ritesh↗