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 271 records · Page 15

Multi-Scale Simulation of Wind Farm Performance during a Frontal Passage

Predicting the response of wind farms to changing flow conditions is necessary for optimal design and operation. Here in this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine–flow interaction, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with 10 m grid spacing. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed background velocity, turbine wake structure, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD model framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗

Historical and Future Changes in Air Pollutants from CMIP6 Models

Poor air quality is currently responsible for large impacts on human health across the world. In addition, the air pollutants ozone (O3) and particulate matter less than 2.5 µm in diameter (PM2.5) are also radiatively active in the atmosphere and can influence Earth's climate. It is important to understand the effect of air quality and climate mitigation measures over the historical period and in different future scenarios to ascertain any impacts from air pollutants on both climate and human health. The Coupled Model Intercomparison Project Phase 6 (CMIP6) presents an opportunity to analyse the change in air pollutants simulated by the current generation of climate and Earth system models that include a representation of chemistry and aerosols (particulate matter). The shared socio-economic pathways (SSPs) used within CMIP6 encompass a wide range of trajectories in precursor emissions and climate change, allowing for an improved analysis of future changes to air pollutants. Firstly, we conduct an evaluation of the available CMIP6 models against surface observations of O3 and PM2.5. CMIP6 models consistently overestimate observed surface O3 concentrations across most regions and in most seasons by up to 16 ppb, with a large diversity in simulated values over Northern Hemisphere continental regions. Conversely, observed surface PM2.5 concentrations are consistently underestimated in CMIP6 models by up to 10 µg m−3, particularly for the Northern Hemisphere winter months, with the largest model diversity near natural emission source regions. The biases in CMIP6 models when compared to observations of O3 and PM2.5 are similar to those found in previous studies. Over the historical period (1850–2014) large increases in both surface O3 and PM2.5 are simulated by the CMIP6 models across all regions, particularly over the mid to late 20th century, when anthropogenic emissions increase markedly. Large regional historical changes are simulated for both pollutants across East and South Asia with an annual mean increase of up to 40 ppb for O3 and 12 µg m−3 for PM2.5. In future scenarios containing strong air quality and climate mitigation measures (ssp126), annual mean concentrations of air pollutants are substantially reduced across all regions by up to 15 ppb for O3 and 12 µg m−3 for PM2.5. However, for scenarios that encompass weak action on mitigating climate and reducing air pollutant emissions (ssp370), annual mean increases in both surface O3 (up 10 ppb) and PM2.5 (up to 8 µg m−3) are simulated across most regions, although, for regions like North America and Europe small reductions in PM2.5 are simulated due to the regional reduction in precursor emissions in this scenario. A comparison of simulated regional changes in both surface O3 and PM2.5 from individual CMIP6 models highlights important regional differences due to the simulated interaction of aerosols, chemistry, climate and natural emission sources within models. The projection of regional air pollutant concentrations from the latest climate and Earth system models used within CMIP6 shows that the particular future trajectory of climate and air quality mitigation measures could have important consequences for regional air quality, human health and near-term climate. Differences between individual models emphasise the importance of understanding how future Earth system feedbacks influence natural emission sources, e.g. response of biogenic emissions under climate change.

Air pollutants↗

Experience in the control of a continuous flow cyrogenic tunnel

The concept of using liquid nitrogen to cool the test gas of a wind tunnel to cryogenic temperatures lead to the achievement of very high Reynolds number flows in relatively small transonic tunnels. The economical operation of liquid nitrogen cooled cryogenic tunnels is critically dependent on fast and accurate control of the tunnel variables. The control problem of a continuous flow fan driven cryogenic tunnel was addressed, first by developing a lumped multivariable mathematical model of a tunnel and validating the model by reconciling the responses of the Langley 0.3 m transonic cryogenic tunnel to the responses of the mathematical model on a simulator. Finally, the development of laws for the closed loop control of the tunnel pressure and temperature and the successful implementation of a control system for the 0.3 m transonic cryogenic tunnel based on these laws are presented.

Kilgore, R. A.↗

A continuum damage mechanics (CDM) based Wilshire model for creep deformation, damage, and rupture prediction

In this report, a new constitutive model is derived that combines the Wilshire equations with continuum damage mechanics (CDM) to enable the long-term prediction of creep deformation, damage, and rupture. The Wilshire equations have been demonstrated to accurately extrapolate the stress-rupture, minimum-creep-strain-rate, and time-to-strain of various alloys across decade of life. Recently, the time-to-creep-strain equation has been rearranged to predict creep deformation curves; however, the resulting equation is difficult to calibrate and does not track damage. The CDM-based Sinh constitutive model accurately predicts creep deformation, damage, and stress-rupture for a variety of alloys; however, it lacks an explicit description of stress and temperature co-dependence which limits interpolation and extrapolation ability. In this study, the Wilshire and the CDM-based Sinh equations are combined to create the “WCS” model. The WCS model combines the best features of each model while eliminating their deficiencies. Experimental data for alloy P91 is gathered and the WCS model is calibrated to the data. The WCS model accurately predicts the stress rupture, minimum-creep-strain-rate, creep deformation, and damage of alloy P91 across the experimental stress and temperature range as well as consistency in parametric simulations.

36 MATERIALS SCIENCE↗

A time dependent anatomically detailed model of cardiac conduction

In order to understand the determinants of transitions in cardiac electrical activity from normal patterns to dysrhythmias such as ventricular fibrillation, we are constructing an anatomically and physiologically detailed finite element simulation of myocardial electrical propagation. A healthy human heart embedded in paraffin was sectioned to provide a detailed anatomical substrate for model calculations. The simulation of propagation includes anisotropy in conduction velocity due to fiber orientation as well as gradients in conduction velocities, absolute and relative refractory periods, action potential duration and electrotonic influence of nearest neighbors. The model also includes changes in the behaviour of myocardial tissue as a function of the past local activity. With this model, we can examine the significance of fiber orientation and time dependence of local propagation parameters on dysrhythmogenesis.

Computer Simulation↗

Meso- to microscale modeling of atmospheric stability effects on wind turbine wake behavior in complex terrain

Abstract. Terrain-induced flow phenomena modulate wind turbine performance and wake behavior in ways that are not adequately accounted for in typical wind turbine wake and wind plant design models. In this work, we simulate flow over two parallel ridges with a wind turbine on one of the ridges, focusing on conditions observed during the Perdigão field campaign in 2017. Two case studies are selected to be representative of typical flow conditions at the site, including the effects of atmospheric stability: a stable case where a mountain wave occurs (as in ∼ 50 % of the nights observed) and a convective case where a recirculation zone forms in the lee of the ridge with the turbine (as occurred over 50 % of the time with upstream winds normal to the ridgeline). We use the Weather Research and Forecasting Model (WRF), dynamically downscaled from the mesoscale (6.75 km resolution) to microscale large-eddy simulation (LES) at 10 m resolution, where a generalized actuator disk (GAD) wind turbine parameterization is used to simulate turbine wakes. We compare the WRF–LES–GAD model results to data from meteorological towers, lidars, and a tethered lifting system, showing good qualitative and quantitative agreement for both case studies. Significantly, the wind turbine wake shows different amounts of vertical deflection from the terrain and persistence downstream in the two stability regimes. In the stable case, the wake follows the terrain along with the mountain wave and deflects downwards by nearly 100 m below hub height at four rotor diameters downstream. In the convective case, the wake deflects above the recirculation zone over 40 m above hub height at the same downstream distance. Overall, the WRF–LES–GAD model is able to capture the observed behavior of the wind turbine wakes, demonstrating the model's ability to represent wakes over complex terrain for two distinct and representative atmospheric stability classes, and, potentially, to improve wind turbine siting and operation in hilly landscapes.

17 WIND ENERGY↗

The Impact of the Atlantic Cold Tongue on West African Monsoon Onset in Regional Model Simulations for 1998-2002

The Atlantic cold tongue (ACT) develops during spring and early summer near the Equator in the Eastern Atlantic Ocean and Gulf of Guinea. The hypothesis that the ACT accelerates the timing of West African monsoon (WAM) onset is tested by comparing two regional climate model (RM3) simulation ensembles. Observed sea surface temperatures (SST) that include the ACT are used to force a control ensemble. An idealized, warm SST perturbation is designed to represent lower boundary forcing without the ACT for the experiment ensemble. Summer simulations forced by observed SST and reanalysis boundary conditions for each of five consecutive years are compared to five parallel runs forced by SST with the warm perturbation. The article summarizes the sequence of events leading to the onset of the WAM in the Sahel region. The representation of WAM onset in RM3 simulations is examined and compared to Tropical Rainfall Measuring Mission (TRMM), Global Precipitation Climatology Project (GPCP) and reanalysis data. The study evaluates the sensitivity of WAM onset indicators to the presence of the ACT by analysing the differences between the two simulation ensembles. Results show that the timing of major rainfall events and therefore theWAM onset in the Sahel are not sensitive to the presence of the ACT. However, the warm SST perturbation does increase downstream rainfall rates over West Africa as a consequence of enhanced specific humidity and enhanced northward moisture flux in the lower troposphere.

Regional model simulation↗

Comparing multi-source urban flood indicators: satellite, simulation, and citizen-reported data

Urban flooding arises from complex mechanisms, making it challenging to capture accurately with a single detection method. This study evaluates three complementary approaches to detect flooding across three Chicago neighborhoods: (i) Sentinel-1 synthetic aperture radar (SAR), offering weather-independent, high-resolution (10 m) imagery of surface inundation; (ii) the storm water management model (SWMM), simulating combined sewer overflow and drainage performance; and (iii) citizen-generated 311 service requests, capturing observed flooding impacts. By analyzing six storms ranging from severe to mild, we examine how each source uniquely contributes to identifying urban flood events. SAR imagery effectively identifies standing water but can miss brief flooding due to satellite revisit constraints. SWMM provides detailed insights into system-wide drainage behavior yet may underestimate localized street-level flooding. Meanwhile, 311 calls reflect real-world flooding impacts but are vulnerable to underreporting. Statistical overlap analysis highlights chronic flood hotspots repeatedly identified across multiple detection methods, indicating persistent infrastructure and topographic vulnerabilities. Temporal analysis further reveals that while SWMM flooding aligns closely with rainfall peaks, 311 calls typically precede or persist beyond these peaks. Our findings emphasize the value of using satellite observations, hydrological modeling, and resident-reported data in a complementary manner to better interpret patterns in flood timing, severity, and spatial distribution—providing insights that can inform targeted infrastructure improvements and contribute to urban flood resilience planning.

311↗

ExaWind: A multifidelity modeling and simulation environment for wind energy

We introduce the open-source ExaWind modeling and simulation environment for wind energy. The primary physics codes of ExaWind are Nalu-Wind and OpenFAST. Nalu-Wind is a wind-focused computational fluid dynamics (CFD) code that is coupled to the whole-turbine simulation code OpenFAST. The ExaWind environment was created under U.S. Department of Energy funding to achieve the highest-fidelity simulations of wind turbines and wind farms to date, with the goal of enabling disruptive changes to turbine and plant design and operation. Innovation will be gleaned through better understanding of the complex flow dynamics in wind farms, including wake evolution and the impact of wakes on downstream turbines and turbulent flow from complex terrain. High-fidelity predictive simulations employ hybrid turbulence models, geometry/boundary-layer-resolving CFD meshes, atmospheric turbulence, nonlinear structural dynamics, and fluid-structure interaction. While there is an emphasis on very high-fidelity simulations (e.g., blade resolved with full fluid-structure coupling), the ExaWind environment supports lower-fidelity modeling capabilities including actuator-line and -disk methods. Important in the development of ExaWind codes is that the codes scale well on today's largest petascale supercomputers and on the next-generation platforms that will enable exascale computing.

17 WIND ENERGY↗

Magnetic Field Orientation Effects on the Standoff Distance of Earth's Bow Shock

Three-dimensional, global MHD simulations of solar wind flow onto a prescribed magnetopause obstacle are used to show that a bow shock's nose location a(sub s), and the relative subsolar magnetosheath thickness Delta(sub ms)/a(sub mp) are strong functions of the IMF cone angle theta (between v(sub sw) and B(sub sw)) and the Alfven Mach number M(sub A). For a given M(sub A) the shock is more distant for higher theta (restricted to the interval 0-90deg by symmetries), while a(sub s)/a(sub mp) and Delta(sub ms/a(sub mp) increase with decreasing M(sub A) for theta greater than or approximately 20deg but decrease with decreasing M(sub A) for theta approximately Odeg. Large differences in Delta(sub ms/a(sub mp) are predicted between theta = Odeg and 90deg at low M(sub A), with smaller differences remaining even at M(sub A) approximately 10. The theta = Odeg results confirm and extend the previous work of Spreiter and Rizzi [1974]. The simulations show that successful models for the subsolar shock location cannot subsume the dependences on M(sub A) and theta into a sole dependence on M(ms). Instead, they confirm a recent prediction [Cairns and Grabbe, 1994] that a(sub s)/a(sub mp) and Delta(sub ms)/a(sub mp) should depend strongly on theta and M(sub A) for M(sub A) less than or approximately 10 (as well as on other MHD variables). Detailed comparisons between theory and data remain to be done. However, preliminary comparisons show good agreement, with distant shock locations found for low M(sub A) and large theta greater than or approximately 45deg and closer locations found for theta less than or approximately 20deg even at M/A approximately 8.

Cairns, Iver H.↗

Aero-Propulsive Damping Characterization for eVTOL Aircraft Using Free Motion Wind-Tunnel Testing

This paper describes an electric vertical takeoff and landing (eVTOL) aircraft system identification method applied using three degree-of-freedom (3DOF) free motion wind-tunnel testing. The approach, similar to flight-test system identification, allows for efficient mathematical model development of the aero-propulsive moments applied on an eVTOL vehicle, including aerodynamic damping effects. The approach is demonstrated using a subscale tiltrotor eVTOL aircraft mounted on a new 3DOF wind-tunnel apparatus. To execute the test, a model-based 3DOF control system is designed to track attitude commands and transition the aircraft based on the freestream dynamic pressure. While the flight controller is active, orthogonal phase-optimized multisine inputs are injected into the attitude command and control effector command signals to enable collection of informative data for model identification. Aero-propulsive models are then identified at several reference conditions in the transition flight envelope using the equation-error method in the frequency domain. The identified models are shown to have a good fit to the modeling data and good prediction capability of data not used for model identification. The method yields aerodynamic damping estimates using less wind-tunnel test time compared to traditional forced oscillation experiments and supplements static wind-tunnel testing to produce a comprehensive transition aero-propulsive model suitable for use in flight dynamics simulations.

RAVEN↗

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape

Arctic warming is altering vegetation and carbon dynamics with global implications, yet Earth System Model (ESM) predictions in the Arctic remain highly uncertain, in part due to historically limited data for model parameterization and validation. As such, ESMs typically represent Arctic ecosystems in an oversimplified manner. Recently, nine plant functional types (PFTs) designed to realistically represent tundra vegetation were integrated into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) and parameterized using plot-scale observations from a single site. Additional evaluation was needed to determine their transferability across the Arctic. Here, in this study, we evaluated whether refined representation of tundra vegetation improved model accuracy by conducting spatially explicit 100 × 100 m resolution ELM simulations on Alaska's Seward Peninsula. Simulations with the default two-PFT configuration and with the nine Arctic-specific PFTs were benchmarked against observations of net ecosystem exchange, gross primary production, and aboveground biomass from multiple data streams including an eddy covariance flux tower, flux chambers, and aircraft and unoccupied aerial system hyperspectral remote sensing. Evaluation revealed that Arctic-specific PFT simulations produced more realistic landscape-level carbon exchanges, and better captured observed heterogeneity in biomass and productivity, explaining 60%–70% of spatial variance (R 2 = 0.6–0.7) compared to just 12%–18% (R 2 = 0.12–0.18) with the default configuration. However, the refined model failed to reproduce observed aboveground biomass for highly productive alder-willow communities, requiring further evaluation of carbon allocation parameterizations for tall shrubs that are increasingly expanding across tundra landscapes. Our results demonstrate that enhanced representation of vegetation heterogeneity boosts predictive understanding of tundra carbon dynamics, facilitating regional to pan-Arctic model and remote-sensing scaling.

Murphy, Bailey A. [Oak Ridge National Laboratory (↗

Quantifying vertical wind shear effects in shallow cumulus clouds over Amazonia

Abstract. This study analyses and quantifies the effects of vertical wind shear (VWS) on the properties of shallow cumulus cloud fields over Central Amazonia. We perform idealised simulations with high resolution (50 m horizontally and 20 m vertically) using the Dutch Atmospheric Large-Eddy Simulation (DALES) model, changing the initial conditions and large-scale forcing of VWS. The resulting cloud field is analysed by applying a cloud tracking algorithm to generate Lagrangian datasets of the life cycle of individual clouds as well as their time-varying core and margin dimensions. The reference run has no wind speed or directional shear and represents a typical day in the local dry season. Numerical experiments with moderate and high wind speed shear are simulated by adding linear increases in the wind speed of 1.2 and 2.4 m s−1 km−1, respectively. Three additional runs are made by adding 90∘ of wind rotation between the surface and the top of the domain (5 km) on top of the three wind speed shear conditions. We find that clouds developing in a sheared environment have horizontal equivalent diameter increased by up to 100 m on average, but the cloud depth is reduced. Our quantification shows that VWS tends to increase the size of the cloud cores but reduces their relative area, volume, and mass fractions compared to the overall cloud dimensions. The addition of 2.4 m s−1 km−1 of VWS decreases the relative core area by about 0.03 (about 10 % of the overall average) and its volume and mass ratios by about 0.05 (10 %–25 % in relative terms). Relevant for the cloud transport properties is that the updraught speed and the liquid water content are lower within the cores, and consequently so is the upward mass flux. All quantifications of mean cloud properties point to the inhibition of convective strength by VWS, therefore hampering the shallow-to-deep transition. However, open questions still remain given that the individually deepest clouds were simulated under high environmental shear, even though they occur in small numbers. This could indicate other indirect effects of VWS that have opposite effects on cloud development if found to be significant in the future.

54 ENVIRONMENTAL SCIENCES↗

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ICRF near-field effects sensitivity to magnetic tilt-angle

This paper complements the study reported in G. Urbanczyk et al 2025 Nucl. Fusion 65 046018 by investigating how near-field effects from Ion Cyclotron Range of Frequency (ICRF) antennas depend on the magnetic geometry, particularly focusing on the sensitivity of rectified sheath potentials and associated impurity production to variations in magnetic field tilt-angle. Simulations were conducted using Petra-M, employing a flat model of the ASDEX Upgrade (AUG) three-strap antenna, and COMSOL, utilizing a curved model of the WEST antenna. Results demonstrate indications of sensitivity of the rectified sheath potentials and impurity sputtering to the alignment of Faraday screen bars with the magnetic field lines, emphasizing the necessity for careful consideration of geometric alignment in antenna design for fusion devices. This trend is consistent with experimental data measured with various diagnostics.

Urbanczyk, Guillaume [Universite de Lorraine, Nanc↗

Simulations and models for the Richtmyer–Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P ∝ km, where k is the wavenumber. The hydrodynamic condition is the same as in the θ-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ∼1.86 shock and gamma-law = 5/3 fluids with Atwood number A = −0.5. The bubble amplitude hB is found to grow in two stages. Initially, hB undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, hB grows as a power law ∼tθ when k|hB| > O(1), where k is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗