Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterized dynamical systems”

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 55 records · Page 3

A model-independent data assimilation (MIDA) module and its applications in ecology

Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Earth system dynamics↗

A Multiscale Modeling System: Developments, Applications, and Critical Issues

The foremost challenge in parameterizing convective clouds and cloud systems in large-scale models are the many coupled dynamical and physical processes that interact over a wide range of scales, from microphysical scales to the synoptic and planetary scales. This makes the comprehension and representation of convective clouds and cloud systems one of the most complex scientific problems in Earth science. During the past decade, the Global Energy and Water Cycle Experiment (GEWEX) Cloud System Study (GCSS) has pioneered the use of single-column models (SCMs) and cloud-resolving models (CRMs) for the evaluation of the cloud and radiation parameterizations in general circulation models (GCMs; e.g., GEWEX Cloud System Science Team 1993). These activities have uncovered many systematic biases in the radiation, cloud and convection parameterizations of GCMs and have led to the development of new schemes (e.g., Zhang 2002; Pincus et al, 2003; Zhang and Wu 2003; Wu et al. 2003; Liang and Wu 2005; Wu and Liang 2005, and others). Comparisons between SCMs and CRMs using the same large-scale forcing derived from field campaigns have demonstrated that CRMs are superior to SCMs in the prediction of temperature and moisture tendencies (e.g., Das et al. 1999; Randall et al 2003b; Xie et al. 2005).

Tao, Wei-Kuo↗

Exploring Climatology and Long-Term Variations of Aerosols from NASA Reanalysis MERRA-2 with Giovanni

Dust plays important roles in energy cycle and climate variations. The dust deposition is the major source of iron in the open ocean, which is an essential micronutrient for phytoplankton growth and therefore may influence the ocean uptake of atmospheric CO2. Mineral dust can also act as fertilizer for forests over long time periods. Over 35 years of simulated global aerosol products from NASA atmospheric reanalysis, second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) are available from NASA Goddard Earth Science Data and Information Services Center (GES DISC). The MERRA-2 covers the period 1980-present, continuing as an ongoing climate analysis. Aerosol assimilation is included throughout the period, using MODIS, MISR, AERONET, and AVHRR (in the pre-EOS period). The aerosols are assimilated by using MERRA-2 aerosol model, which interact directly with the radiation parameterization, and radiatively coupled with atmospheric model dynamics in the Goddard Earth Observing System Model, Version 5 (GEOS-5). Dust deposition data along with other major aerosol compositions (e.g. black carbon, sea salt, and sulfate, etc.) are simulated as dry and wet deposition, respectively. The hourly and monthly data are available at spatial resolution of 0.5ox0.625o (latitude x longitude). Quick data exploration of climatology and interannual variations of MERRA-2 aerosol can be done through the online visualization and analysis tool, Giovanni. This presentation, using dust deposition as an example, demonstrates a number of MERRA-2 data services at GES DISC. Global distributions of dust depositions, and their seasonal and inter-annual variations are investigated from MERRA-2 monthly aerosol products.

data services↗

DAO's Finite-volume/Physical-space Data Assimilation System: Stratospheric Applications

In this talk we describe the next-generation data assimilation system being developed at NASA's Data Assimilation Office (DAO), with emphasis on the applications to stratospheric forecasts and stratospheric constituent transport. This data assimilation system includes the General Circulation Model jointly developed by the DAO and the Climate and Global Dynamics Division (CGDD) at NCAR. This model is based on the finite-volume dynamical core) developed at DAO with physical parameterizations from the NCAR Climate Community Model. The Physical-space Statistical Analysis System (PSAS) is used to combine a first guess from the NASA-NCAR GCM with observational data to provide an updated estimate of the state of the atmosphere. Case studies for the (northern) Winter of 2000 will be discussed.

daSilva, Arlindo↗

Cascading toward a Kilometer‐scale GCM: Impacts of a Scale‐aware Convection Parameterization in the Goddard Earth Observing System GCM.

The NASA Goddard Earth Observing System global model is evaluated through a cascade of simulations with increasing horizontal resolution. This model employs a non-hydrostatic dynamical core and includes a scale-aware, deep convection parameterization (DPCP). 40-daysimulations at six resolutions (100km to 3km) with unvarying model formulation were produced. At the highest resolution, extreme experiments were carried out; one with no DPCP, and one with its scale-awareness eliminated. Simulated precipitation, radiative balance, and atmospheric thermodynamic and dynamical variables are well reproduced with respect to both observational and reanalysis data. As model resolution increases, the convective precipitation smoothly transitions from being mostly produced by the convection parameterization to the cloud microphysics parameterization. However, contrary to current thought, these extreme cases argue for maintaining, to some extent, the scale-aware deep convection parameterization even at 3kmscale, as the run relying solely on explicit grid-scale production of rainfall performs more poorly at this resolution.

Saulo Freitas↗

Stochastic control system parameter identifiability

The parameter identification problem of general discrete time, nonlinear, multiple input/multiple output dynamic systems with Gaussian white distributed measurement errors is considered. The knowledge of the system parameterization was assumed to be known. Concepts of local parameter identifiability and local constrained maximum likelihood parameter identifiability were established. A set of sufficient conditions for the existence of a region of parameter identifiability was derived. A computation procedure employing interval arithmetic was provided for finding the regions of parameter identifiability. If the vector of the true parameters is locally constrained maximum likelihood (CML) identifiable, then with probability one, the vector of true parameters is a unique maximal point of the maximum likelihood function in the region of parameter identifiability and the constrained maximum likelihood estimation sequence will converge to the vector of true parameters.

Lee, C. H.↗

Studying Diurnal Variations of Aerosols with NASA MERRA-2 Reanalysis Data

Aerosols play an important role in atmospheric dynamics, climate variations, and Earth's energy cycle by altering the radiation balance in the atmosphere through interaction with clouds, providing fertilizer for forests and canopy, and as a supply of iron to the ocean over long time periods. Studies suggest that much of the feedback between dust aerosols and dynamics is associated with diurnal and synoptic scale variability. However, the lack of sub-daily resolution of aerosols from satellite observations makes it difficult to study the diurnal characteristics, especially over tropical and subtropical regions. Investigation of this topic utilizes over 37 years of simulated global aerosol products from NASA atmospheric reanalysis, in the second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) data set, available from NASA Goddard Earth Science Data and Information Services Center (GES DISC). MERRA-2 covers the period 1980-present, and is continuing as an ongoing climate analysis. Aerosol assimilation is included throughout the period, using data from MODIS, MISR, AERONET, and AVHRR (in the pre-EOS period). The aerosols are assimilated using the MERRA-2 aerosol model, which interacts directly with radiation parameterization, and is radiatively coupled with atmospheric model dynamics in the Goddard Earth Observing System Model, Version 5 (GEOS-5). Hourly, monthly, and monthly diurnal data are available at spatial resolution of 0.5o x 0.625o (latitude x longitude). By using MERRA-2 hourly and monthly diurnal products, different aerosol diurnal variabilities are observed over North America, Africa, Asia, and Australia, that may be due to different meteorological conditions and aerosol sources. The presentation will also provide an overview of MERRA-2 data services at GES DISC, such as how to find and download data, and how to quickly visualize and analyze data online with Giovanni.

aerosol↗

Alpha-canonical form representation of the open loop dynamics of the Space Shuttle main engine

A parameter and structure estimation technique for multivariable systems is used to obtain a state space representation of open loop dynamics of the space shuttle main engine in alpha-canonical form. The parameterization being used is both minimal and unique. The simplified linear model may be used for fault detection studies and control system design and development.

Duyar, Almet↗

Variational Autoencoders for Learning Nonlinear Dynamics of Physical Systems

We develop data-driven methods for incorporating physical information for priors to learn parsimonious representations of nonlinear systems arising from parameterized PDEs and mechanics. Our approach is based on Variational Autoencoders (VAEs) for learning nonlinear state space models from observations. We develop ways to incorporate geometric and topological priors through general manifold latent space representations. We investigate the performance of our methods for learning low dimensional representations for the nonlinear Burgers equation and constrained mechanical systems.

97 MATHEMATICS AND COMPUTING↗

Improving Time Step Convergence in an Atmosphere Model With Simplified Physics: The Impacts of Closure Assumption and Process Coupling

Convergence testing is a common practice in the development of dynamical cores of atmospheric models but is not as often exercised for the parameterization of sub-grid physics. An earlier study revealed that the stratiform cloud parameterizations in several predecessors of the Energy Exascale Earth System Model (E3SM) showed strong time-step sensitivity and slower-than-expected convergence when the model's time step was systematically refined. In this work, a simplified atmosphere model is configured that consists of the spectral-element dynamical core of the E3SM atmosphere model coupled with a large-scale condensation parameterization based on commonly used assumptions. This simplified model also resembles E3SM and its predecessors in the numerical implementation of process coupling and shows poor time-step convergence in short ensemble tests. We present a formal error analysis to reveal the expected time-step convergence rate and the conditions for obtaining such convergence. Numerical experiments are conducted to investigate the root causes of convergence problems. We show that revisions in the process coupling and closure assumption help to improve convergence in short simulations using the simplified model; the same revisions applied to a full atmosphere model lead to significant changes in the simulated long-term climate. This work demonstrates that causes of convergence issues in atmospheric simulations can be understood by combining analyses from physical and mathematical perspectives. Addressing convergence issues can help to obtain a discrete model that is more consistent with the intended representation of the physical phenomena.

54 ENVIRONMENTAL SCIENCES↗

Structural inference of networked dynamical systems with universal differential equations

Networked dynamical systems are common throughout science in engineering; e.g., biological networks, reaction networks, power systems, and the like. For many such systems, nonlinearity drives populations of identical (or near-identical) units to exhibit a wide range of nontrivial behaviors, such as the emergence of coherent structures (e.g., waves and patterns) or otherwise notable dynamics (e.g., synchrony and chaos). Here, we seek to infer (i) the intrinsic physics of a base unit of a population, (ii) the underlying graphical structure shared between units, and (iii) the coupling physics of a given networked dynamical system given observations of nodal states. These tasks are formulated around the notion of the Universal Differential Equation, whereby unknown dynamical systems can be approximated with neural networks, mathematical terms known a priori (albeit with unknown parameterizations), or combinations of the two. We demonstrate the value of these inference tasks by investigating not only future state predictions but also the inference of system behavior on varied network topologies. The effectiveness and utility of these methods are shown with their application to canonical networked nonlinear coupled oscillators.

97 MATHEMATICS AND COMPUTING↗

NN-OpInf

SAND2026-18878O The NN-OpInf tool is a PyTorch-based approach to operator inference that uses composable, structure-preserving neural networks to represent nonlinear operators. Operator inference is a machine learning method for inferring low-dimensional systems from data and polynomial models for system dynamics. However, many systems do not conform to polynomial structures, which NN-OpInf addresses by parameterizing operators with neural networks. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

3D multi-system Bayesian calibration with energy conservation to study rapidity-dependent dynamics of nuclear collisions

Considerable information about the early-stage dynamics of heavy ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems — large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late-time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Precise Orbit and the Challenge of Long Term Stability

The computation of a precise orbit reference is a fundamental component of the altimetric measurement. Since the dawn of the modern altimeter age, orbit accuracy has been determined by the quality of the GPS, SLR, and DORIS tracking systems, the fidelity of the measurement and force models, and the choice of parameterization for the orbit solutions, and whether a dynamic or a reduced-dynamic strategy is used to calculate the orbits. At the start of the TOPEX mission, the inaccuracies in the modeling of static gravity, dynamic ocean tides, and the nonconservative forces dominated the orbit error budget. Much of the error due to dynamic mismodeling can be compensated by reduced-dynamic tracking techniques depending on the measurement system strength. In the last decade, the launch of the GRACE mission has eliminated the static gravity field as a concern, and the background force models and the terrestrial reference frame have been systematically refined. GPS systems have realized many improvements, including better modeling of the forces on the GPS spacecraft, large increases in the ground tracking network, and improved modeling of the GPS measurements. DORIS systems have achieved improvements through the use of new antennae, more stable monumentation, and of satellite receivers that can track multiple beacons, and as well as through improved modeling of the nonconservative forces. Many of these improvements have been applied in the new reprocessed time series of orbits produced for the ERS satellites, Envisat, TOPEX/Poseidon and the Jason satellites, and as well as for the most recent Cryosat-2 and HY2A. We now face the challenge of maintaining a stable orbit reference for these altimetric satellites. Changes in the time-variable gravity field of the Earth and how these are modelled have been shown to affect the orbit evolution, and the calibration of the altimetric data with tide gauges. The accuracy of the reference frame realizations, and their projection into the future remains a source of error. Other sources of omission error include the geocenter for which no consensus model is as of yet applied. Although progress has been made in nonconservative force modeling through the use of detailed satellite-specific models, radiation pressure modeling, and atmospheric density modeling remain a potential source of orbit error. The longer term influence of variations in the solar and terrestrial radiation fields over annual and solar cycles remains principally untested. Also the long term variation in optical and thermal properties of the space vehicle surfaces would contribute to biases in the orbital frame if ignored. We review the status of altimetric precision orbit determination as exemplified by the recent computations undertaken by the different analysis centers for ERS, Envisat, TOPEX/Poseidon, Jason, Cryosat2 and HY2A, and we provide a perspective on the challenges for future missions such as the Jason-3, SENTINEL-3 and SWOT.

Lemoine, Frank G.↗

Stabilized Hyperfoam Modeling of the General Plastics EF4003 (3 PCF) Flexible Foam

Constitutive model parameterizations for the General Plastics EF4003 low density 3 pound per cubic foot are needed for design and qualification purposes in normal and abnormal mechanical simulations. The material is expected to be deformed in two ways: first during preloading, and second under impact conditions of the system (transient dynamic). All analyses are to be performed at room temperature. The goal is to provide the analysis community a robust constitutive model parameterization to represent the compression behavior of the EF4003 foam from small deformations up to massive compressive deformations when the foam is densifying. It is worth noting the EF4003 exhibits anisotropy in its stress-strain behavior between the rise and transverse directions (See figure 2.8c-d) as well as plateau behavior that is very likely to cause material stability issues, due to the buckling transition, (and has historically done so) when using Sandia’s current workhorse models for flexible foams, Hyperfoam and Flex Foam. A Stability-informed Hyperfoam parameterization procedure is developed and executed to calibrate a hyperfoam model for the EF4003 room temperature, transversely loaded data. A rise orientation parameterization was not attempted due to localization in the experiments.

36 MATERIALS SCIENCE↗

Development of a rotorcraft. Propulsion dynamics interface analysis, volume 2

A study was conducted to establish a coupled rotor/propulsion analysis that would be applicable to a wide range of rotorcraft systems. The effort included the following tasks: (1) development of a model structure suitable for simulating a wide range of rotorcraft configurations; (2) defined a methodology for parameterizing the model structure to represent a particular rotorcraft; (3) constructing a nonlinear coupled rotor/propulsion model as a test case to use in analyzing coupled system dynamics; and (4) an attempt to develop a mostly linear coupled model derived from the complete nonlinear simulations. Documentation of the computer models developed is presented.

Hull, R.↗

A model-independent data assimilation (MIDA) module and its applications in ecology

Abstract. Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

58 GEOSCIENCES↗

Implementation and Evaluation of Physics-Driven Dynamic Entrainment-Mixing Parameterization in a Climate Model and Its Impact on Low-Cloud Simulation

The turbulent entrainment-mixing process in the Community Earth System Model version 1.2 (CESM1.2) is assumed to follow the extremely inhomogeneous entrainment-mixing. However, different entrainment-mixing scenarios can occur in real clouds. To address this deficiency, a unifying parameterization that represents different entrainment-mixing processes is implemented and evaluated in CESM1.2. The results indicate that the homogeneous mixing degree values simulated by the new parameterization in CESM1.2 are predominantly greater than 50%, suggesting a tendency toward homogeneous mixing. Compared to the extremely inhomogeneous mixing mechanism, the new parameterization increases the cloud droplet number concentration (Nc). More importantly, the new parameterization improves low-cloud fraction (CLDLOW) simulation in Northwest Pacific (NWP) and Southeast Pacific (SEP) regions, with relative improvements of 2.95% and 4.17%, respectively. Furthermore, the improvements reach up to 44.6% and 16.2% in the NWP and SEP regions, respectively, when considering the relationship between N c and CLDLOW. Further analysis reveals that the new parameterization enhances cloud optical depth, longwave radiative cooling effect, net condensation rate, cloud water mixing ratio, lower-troposphere stability, and CLDLOW by increasing N c . Additionally, these results underscore the importance of improving entrainment-mixing parameterization in climate models.

54 ENVIRONMENTAL SCIENCES↗