Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation 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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Global climate mode resonance due to rapidly intensifying El Niño-Southern Oscillation

The El Niño-Southern Oscillation (ENSO) influences climate variability globally, encompassing various other modes of variability, and thus represents a key predictable climate signal on seasonal timescales. Yet, its response to greenhouse warming remains uncertain, with models projecting a range of outcomes. Here, we demonstrate that in response to warming, a state-of-the-art high-resolution climate model simulates a rapid transition from a moderate-amplitude irregular regime, as observed in the current climate, to a highly regular oscillation with intensifying amplitude. This behaviour can be attributed to increasing air-sea feedbacks, which approach criticality in the second half of this century, and growing atmospheric noise. As ENSO intensifies in this model, it synchronizes with other prominent climate modes, such as the North Atlantic Oscillation and the Indian Ocean Dipole, thereby imprinting its regular, predictable variability on them. If realized, this global climate mode resonance would have wide-ranging whiplash impacts on regional hydroclimates.

Atmospheric dynamics↗
DOE OSTI2026

Flash flourishing of Northern Hemisphere vegetation and its drivers

Rapid surges in vegetation growth—defined by thresholds in growth rate and duration—are critical yet understudied indicators of ecosystem responses to environmental change. Here, we investigate spatiotemporal patterns of such abrupt, short-lived flash flourishing events across the northern extratropical latitudes (NEL) from 2003 to 2022. We find more frequent occurrence of flash flourishing events at high latitudes (≥45° N), where their incidence is 1.6 times higher than at mid-latitudes. Moreover, there is an increasing tendency in frequency, duration, and intensity of flash flourishing events over the past two decades, alongside consistent rises in vegetation indices across onset, post-onset, and entire phases. Model simulations attribute these multiyear increases primarily to elevated atmospheric CO 2 , while temperature and radiation predominantly control phase-specific variability, with onset traits strongly predicting subsequent phenological responses. Together, these findings identify the patterns and drivers of NEL flash flourishing and highlight their large-scale impacts on ecosystem dynamics, offering critical insights for model improvement and the assessment of ecological shifts.

Kong, Xiangxu [Nanjing University of Information S↗
DOE OSTI2021

Delocalization Transition of a Disordered Axion Insulator

The axion insulator is a higher-order topological insulator protected by inversion symmetry. We show that, under quenched disorder respecting inversion symmetry on average, the topology of the axion insulator stays robust, and an intermediate metallic phase in which states are delocalized is unavoidable at the transition from an axion insulator to a trivial insulator. We derive this conclusion from general arguments, from classical percolation theory, and from the numerical study of a 3D quantum network model simulating a disordered axion insulator through a layer construction. We find the localization length critical exponent near the delocalization transition to be ν = 1.42 ± 0.12. We further show that this delocalization transition is stable even to weak breaking of the average inversion symmetry, up to a critical strength. Finally, we also quantitatively map our quantum network model to an effective Hamiltonian and we find its low-energy k ∙ p expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Theoretical Control-Centric Modeling for Precision Model-Based Sliding Mode Control of a Hydraulic Artificial Muscle Actuator

Artificial muscles (AMs) typically rely on pneumatic sources of fluid power. The use of hydraulics can increase the power and force to weight and volume ratios of AM actuators. This paper develops a control-centric third-order single-input single-output (SISO) lumped-parameter dynamic model and sliding mode position controller based on Filippov's principle of equivalent dynamics for a braided hydraulic artificial muscle (HAM) actuator. The model predicts the nonlinear behavior of the HAM free contraction and captures the fluid and actuator nonlinear dynamic interactions in addition to the braid deformation. Model simulations are compared to experimental results for quasi-static pressurization, isometric pressurization, and open-loop square wave commands at 0.25, 0.5, and 1 Hz. Experiments of sine wave tracking at 0.25, 0.5, and 1 Hz and continuous square wave tracking at 0.067 Hz are conducted using a sliding mode controller (SMC) derived from the model. The SMC achieves a steady-state error of 6 μm at multiple setpoints within the actuator's 17 mm stroke. Compared to a proportional-integral-derivative (PID) controller, the SMC root-mean-square (RMS) error, mean error, and absolute maximum error are reduced on average by 53%, 61%, and 44%, respectively, demonstrating the benefit of model-based approaches for controlling HAMs.

42 ENGINEERING↗
DOE OSTI2023

Scalable self attraction and loading calculations for unstructured ocean tide models

Self attraction and earth-loading effects are important for accurately modeling global tides. A common approach of handling this forcing is to expand mass anomalies into spherical harmonics, which are scaled by load Love numbers to account for elastic earth deformation. We investigate two different approaches to perform these calculations for ocean models that employ unstructured meshes and distributed memory parallelization. The first approach leverages a highly efficient spherical harmonics library, but requires all-to-one and one-to-all communications and interpolation operations between the unstructured and a structured mesh. This approach is compared to a parallel algorithm that computes the spherical harmonic transformations directly on the unstructured mesh with an all-reduce communication. Here, our results show that although the unstructured mesh calculations are more expensive, the scalability of the unstructured mesh approach allows for more efficient spherical harmonics transforms for high-resolution meshes and large processor counts. This methodology enables the efficient inclusion of tidal dynamics large-scale Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Observed humidity trends in dry regions contradict climate models

Arid and semi-arid regions of the world are particularly vulnerable to greenhouse gas–driven hydroclimate change. Climate models are our primary tool for projecting the future hydroclimate that society in these regions must adapt to, but here, we present a concerning discrepancy between observed and model-based historical hydroclimate trends. Over the arid/semi-arid regions of the world, the predominant signal in all model simulations is an increase in atmospheric water vapor, on average, over the last four decades, in association with the increased water vapor–holding capacity of a warmer atmosphere. In observations, this increase in atmospheric water vapor has not happened, suggesting that the availability of moisture to satisfy the increased atmospheric demand is lower in reality than in models in arid/semi-arid regions. This discrepancy is most clear in locations that are arid/semi-arid year round, but it is also apparent in more humid regions during the most arid months of the year. Here, it indicates a major gap in our understanding and modeling capabilities which could have severe implications for hydroclimate projections, including fire hazard, moving forward.

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

A Bottom-up Approach to Improve the Representation of Deep Convective Clouds in Weather and Climate Models (Final Report)

A fundamental yet outstanding question about deep cumulus convection regards how convective updrafts, downdrafts, and surface-based cold pools are inter-related within organized storms, and how these three convective components are modulated by external and internal factors. We considered this question through numerical model simulations as well as data collected during the MC3E field campaign. One of our basic hypotheses was that deep cold pools are forced by wide downdrafts which arise from wide updrafts. This “bottom-up” perspective required consideration of the external controls on updraft width, and then the internal controls on cold-pool characteristics including depth. Model experiments showed that environmental vertical wind shear (an external factor) exerts a particularly strong control on updraft width; analyses of linear and nonlinear dynamic pressure forcing readily explain this control. The experiments also showed that increases in updraft width are accompanied by corresponding increases in downdraft width and cold-pool depth. The updraft–downdraft connection is through precipitation: the wide area of precipitation that will be initiated in a wide updraft will in turn lead to a wide downdraft. The downdraft–cold pool connection is through mass continuity: at the downdraft terminus, the large downdraft mass flux through a wide downdraft must lead to a deep cold pool. Model experiments also showed that the cold-pool depth and other characteristics can be significantly modulated by microphysical processes (internal factors). Specifically, the cold-pool expansion rate, depth, and intensity was found to depend most strongly on the amount of rain evaporation near the ground. However, the sustenance of the cold pool itself depended upon the downdrafts strengthened by sublimation of graupel higher in the storm. Finally, field campaign data revealed compelling inter-relationships between the cold-pool characteristics, and downdraft transience in terms of the microphysics. For example, pulses of hail and graupel (also internal factors) within downdrafts resulted in short-lived increases in the cold-pool depth. However, the cold-pool speed appeared to be relatively independent of the cold-pool depth and intensity, suggesting limitations of density current theory to convective cold pools in nature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Comparison of Deterministic and Statistical Models for Water Quality Compliance Forecasting in the San Joaquin River Basin, California

Model selection for water quality forecasting depends on many factors including analyst expertise and cost, stakeholder involvement and expected performance. Water quality forecasting in arid river basins is especially challenging given the importance of protecting beneficial uses in these environments and the livelihood of agricultural communities. In the agriculture-dominated San Joaquin River Basin of California, real-time salinity management (RTSM) is a state-sanctioned program that helps to maximize allowable salt export while protecting existing basin beneficial uses of water supply. The RTSM strategy supplants the federal total maximum daily load (TMDL) approach that could impose fines associated with exceedances of monthly and annual salt load allocations of up to $1 million per year based on average year hydrology and salt load export limits. The essential components of the current program include the establishment of telemetered sensor networks, a web-based information system for sharing data, a basin-scale salt load assimilative capacity forecasting model and institutional entities tasked with performing weekly forecasts of river salt assimilative capacity and scheduling west-side drainage export of salt loads. Web-based information portals have been developed to share model input data and salt assimilative capacity forecasts together with increasing stakeholder awareness and involvement in water quality resource management activities in the river basin. Two modeling approaches have been developed simultaneously. The first relies on a statistical analysis of the relationship between flow and salt concentration at three compliance monitoring sites and the use of these regression relationships for forecasting. The second salt load forecasting approach is a customized application of the Watershed Analysis Risk Management Framework (WARMF), a watershed water quality simulation model that has been configured to estimate daily river salt assimilative capacity and to provide decision support for real-time salinity management at the watershed level. Analysis of the results from both model-based forecasting approaches over a period of five years shows that the regression-based forecasting model, run daily Monday to Friday each week, provided marginally better performance. However, the regression-based forecasting model assumes the same general relationship between flow and salinity which breaks down during extreme weather events such as droughts when water allocation cutbacks among stakeholders are not evenly distributed across the basin. A recent test case shows the utility of both models in dealing with an exceedance event at one compliance monitoring site recently introduced in 2020.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

An evolving Coupled Model Intercomparison Project phase 7 (CMIP7) and Fast Track in support of future climate assessment

The Coupled Model Intercomparison Project (CMIP) coordinates community-based efforts to answer key and timely climate science questions, facilitate delivery of relevant multi-model simulations through shared infrastructure, and support national and international climate assessments. Generations of CMIP have evolved through extensive community engagement from punctuated phasing into more continuous support for the design of experimental protocols, infrastructure for data publication and access, and public delivery of climate information. We identify four fundamental research questions motivating a seventh phase of coupled model intercomparison relating to patterns of sea surface temperature change, changing weather, the water–carbon–climate nexus, and tipping points. Key CMIP7 advances include an expansion of baseline experiments, a focus on CO 2 -emissions-driven experiments, sustained support for community MIPs, periodic updating of historical forcings and diagnostics requests, and a collection of prioritized experiments, or the “Assessment Fast Track”, drawn from community MIPs to support climate research, assessment, and service goals across prediction and projection, characterization, attribution, and process understanding.

Environmental sciences↗
DOE OSTI2021

Mathematical analysis of laboratory microbial experiments demonstrating deterministic chaotic dynamics

Presented is a system of four ordinary differential equations and a mathematical analysis of microbiological experiments in a four-component chemostat—nutrient n, rods r, cocci c, and predators p. The analysis is consistent with the conclusion that previous experiments produced features of deterministic chaotic and classical dynamics depending on dilution rate. The surrogate model incorporates as much experimental detail as possible, but necessarily contains unmeasured parameters. The objective is to understand better the differences between model simulations and experimental results in complex microbial populations. The key methodology for simulation of chaotic dynamics, consistent with the measured dilution rate and microbial volume averages, was to cause the preference of p for r vs. c to vary with the r and c concentrations, to make r more competitive for nutrient than c, and to recycle some dying p biomass, leading to a modified version of the Monod kinetics model. Our mathematical model demonstrated that the occurrence of chaotic dynamics requires a predator, p, preference for r versus c to increase significantly with increases in r and c populations. Also included is a discussion of several generalizations of the existing model and a possible involvement of the minimum energy dissipation principle. This principle appears fundamental to thermodynamic systems including living systems. Several new experiments are suggested.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Quantifying UAS Observation Error Variance Used in Data Assimilation Systems and Its Impact on Predictive Skill

Observation error determines the weights of the observations and background state used in data assimilation to generate analyses. Quantifying observation error is critical for the optimal assimilation of observational data sets. Uncrewed Aircraft System (UAS) observations have shown potential benefits in filling observational gaps in the lower atmosphere; however, characterization of their error characteristics has been limited. To optimize the use of UAS observations in numerical weather prediction, UAS observation error is estimated based on the 3‐cornered hat diagnostic approach which uses three independent estimates of the atmospheric state. This approach is applied to data from the 2018 Lower Atmospheric Profiling Studies at Elevation‐a Remotely‐piloted Aircraft Team Experiment field campaign using collocated UAS and rawinsonde observations along with output from a set of convection‐permitting model simulations. The estimated observation error values for UAS temperature, wind, and relative humidity measurements were found to be only weakly dependent on height AGL with mean values equal to 0.5°C, 0.8 m s −1 , and 3%, respectively. Only the newly estimated observation error for temperature differed from that previously used to assimilate commercial aircraft observations into global models (1.0°C). However, using this reduced temperature observation error produced more accurate mesoscale analyses and forecasts of both terrain‐driven flows and convection initiation generated by colliding outflow boundaries within the San Luis Valley of Colorado.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sparse Linear Solvers for Large-scale Electromagnetic Transient Simulations

Linear solvers form the basis for electromagnetic transient (EMT) simulations. There is a need to speed up EMT simulations as larger regions are analyzed using EMT simulations. For the same, the performance of linear solvers plays an important role. Exploiting the sparsity of the matrices generated in EMT simulations could assist with speed-up. Scalability is also crucial as power grids expand, demanding solutions capable of accommodating the increasing system size. Recent studies from the North American Electric Reliability Corporation (NERC) increasingly emphasize that EMT simulation models of the power grid will grow larger with the inclusion of power electronics components. Parallelisms in sparsity patterns exploit modern central processing units (CPUs), multi-core CPUs, and graphics processing units (GPUs) architectures in sparse solver designs. Therefore, this paper explores publicly available existing linear solvers and investigates their efficiency in large-scale power grid simulations. A large-scale power grid is developed by increasing the size of the IEEE 39 bus test system to up to 39000 bus systems.

Hsu, Kuan-Chieh↗
DOE OSTI2021

Improving Fission Products at CARIBU: Near Field Detection (Q4/FY21 Quarterly Progress Report)

Since the last quarterly, we have been working on improving the GEANT4 simulation model, which we currently use for the 4-pi gas counter efficiency calculations, to streamline the work when calculating multi-level decay properties. We have also been preparing for the next branching-ratio measurement, which will be for 161 Tb. Because travel is still limited and restricted, we have worked with a summer student on performing simulations to estimate our detector’s efficiency for different thresholds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Analysis of the HAE activity in the TJ-II stellarator using a Landau closure model

The aim of this study is to analyze the stability of helical Alfvén eigenmodes (HAEs) in TJ-II discharges and the stabilizing effect of the energetic particles generated by the neutral beam injector (NBI) on pressure gradient-driven modes (PGDMs). HAE and PGDM stability is studied using the linear version of the gyro-fluid code FAR3d and the continuous structure by the STELLGAP code. First, Alfvén eigenmode (AE) and PGDM activity observed in the experiments is reproduced by the simulations, identifying unstable m/n = 4/7 − 2/3 and 7/12 − 5/8 HAEs triggered around ρ = 0.66 showing a frequency of 209 and 204 kHz, respectively, as well as 5/3 PGDM. Next, a parametric study is performed with respect to the thermal ion density and iota profile in the middle-outer plasma region to verify the robustness of the simulation results with respect to the uncertainty of experimental profiles. The analysis confirms that experimental uncertainty does not cause large deviations in the simulation results, showing the destabilization of the same HAEs for all the configurations tested. The simulations also indicate the decay of the 5/3 PGDM growth rate as the energetic particle (EP) population in the plasma increases, consistent with the experiment. Stability analysis of the n = 3, 7, 11, n = 5, 9, 13, n = 6, 10, 14, and n = 8, 12 helical families is performed with respect to the NBI operational regime for different EP energies, β as well as deposition profiles. The most unstable configuration is the radially localized on-axis NBI operation (stiff EP density profile gradients nearby the magnetic axis). Using the simulation model that reproduces the observed Alfvén activity, we extend the study to analyze NBI performance within a theoretical framework. It shows that increasing NBI voltage (which raises EP energy) leads to a degradation in NBI performance for a given power (related to EP β and their density). To achieve better NBI operation, higher voltage must be balanced with lower injection power, ensuring stable AEs while keeping the same EP β.

Plasma confinement↗
DOE OSTI2020

Identification of the molecular-weight growth reaction network in counterflow flames of the C 3 H 4 isomers allene and propyne

The reaction networks responsible for aromatics formation in counterflow flames of the C 3 H 4 isomers allene and propyne are identified through a combined experimental and modeling study. Mole fraction profiles of near-atmospheric pressure (933 mbar) diffusion flames fueled by the C 3 H 4 isomers are analyzed by means of a newly assembled, chemically detailed kinetic mechanism. The experiment consists of a counterflow burner system that is coupled to a high-resolution time-of-flight molecular-beam mass spectrometer with single-photon ionization via synchrotron-generated vacuum-ultraviolet photons. Flame-sampled, mass-specific photoionization efficiency curves are used to identify the presence of aliphatically substituted aromatic species in addition to the commonly considered pericondensed ring structures. The new mechanism describes the formation and growth of aromatics through repetitive sequences of radical–radical and radical–molecule reactions that include C 1 6 intermediates. Higher concentrations of aromatic species are observed in the allene flame and the new mechanism captures the observed experimental trends very accurately. The results indicate the importance of the aliphatically substituted aromatics and of ring-enlargement reactions for the growth reactions. According to the model simulations, radical+radical recombination and PAH-radical+molecule reactions play an important role in PAH growth.

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

A new metrics framework for quantifying and intercomparing atmospheric rivers in observations, reanalyses, and climate models

We present a new atmospheric river (AR) analysis and benchmarking tool, namely Atmospheric River Metrics Package (ARMP). It includes a suite of new AR metrics that are designed for quick analysis of AR characteristics via statistics in gridded climate datasets such as model output and reanalysis. This package can be used for climate model evaluation in comparison with reanalysis and observational products. Integrated metrics such as mean bias and spatial pattern correlation are efficient for diagnosing systematic AR biases in climate models. For example, the package identifies the fact that, in CMIP5 and CMIP6 (Coupled Model Intercomparison Project Phases 5 and 6) models, AR tracks in the South Atlantic are positioned farther poleward compared to ERA5 reanalysis, while in the South Pacific, tracks are generally biased towards the Equator. For the landfalling AR peak season, we find that most climate models simulate a completely opposite seasonal cycle over western Africa. This tool can also be used for identifying and characterizing structural differences among different AR detectors (ARDTs). For example, ARs detected with the Mundhenk algorithm exhibit systematically larger size, width, and length compared to the TempestExtremes (TE) method. The AR metrics developed from this work can be routinely applied for model benchmarking and during the development cycle to trace performance evolution across model versions or generations and set objective targets for the improvement of models. They can also be used by operational centers to perform near-real-time climate and extreme event impact assessments as part of their forecast cycle.

58 GEOSCIENCES↗
DOE OSTI2025

Emulation With Uncertainty Quantification of Regional Sea‐Level Change Caused by the Antarctic Ice Sheet

Abstract Projecting regional sea‐level change under various climate‐change scenarios typically involves running forward simulations of the Earth's gravitational, rotational and deformational (GRD) response to ice‐mass change, which requires substantial computational cost if applied to probabilistic frameworks requiring thousands to millions of samples. Here we build emulators of regional sea‐level change at 27 coastal locations, due to the GRD effects associated with future Antarctic Ice Sheet mass change over the 21st century. The emulators are evaluated against a numerical sea‐level model applied to an ensemble of ice‐sheet model simulations of the Antarctic Ice Sheet through 2100. We build a physics‐based emulator using a recent sensitivity kernel approach and compare it to machine learning based emulators (neural network and conditional variational autoencoder methods). In order to quantify uncertainty, we derive well‐calibrated prediction intervals for regional sea‐level change via split‐conformal inference and linear regression, and show that Monte Carlo dropout does not yield well‐calibrated uncertainties in this instance. We also demonstrate substantial gains in computational efficiency using both the physics‐based emulator and neural networks in comparison to the numerical model for the complete regional sea‐level solution. Overall, we find the physics‐based emulator modestly outperforms the machine learning emulators for this problem.

58 GEOSCIENCES↗