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 325 records · Page 18

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Long‐Term Large‐Scale Atmospheric Forcing Data From Three‐Dimensional Constrained Variational Analysis for the ARM SGP Site

Here, this study presents a long‐term three‐dimensional large‐scale forcing data set (VARANAL3D) derived from the three‐dimensional constrained variational analysis (3DCVA) method at the Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) site from 2004 to 2018. Building on the same input data sets as the conventional continuous forcing data set (VARANAL), VARANAL3D maintains overall consistency in domain‐averaged fields while introducing spatial variability, offering critical insights into the influence of mesoscale synoptic systems on cloud‐related processes. Evaluations are conducted across four cloud and precipitation regimes: Clear‐sky, Shallow‐clouds, Afternoon‐precipitation, and Nocturnal‐precipitation, presenting high consistency of the domain‐mean forcing data sets while emphasizing the role of subdomain forcing variability particularly in precipitating regimes. Single column model (SCM) simulations demonstrate that subdomain VARANAL3D forcing improves cloud and precipitation representation, with the ensemble outperforming domain‐mean forcing in three cloudy and precipitating regimes. Overall, these results highlight VARANAL3D's value for investigating the impacts of spatial variability of large‐scale forcing on atmospheric processes. The VARANAL3D data set provides new opportunities for evaluating model physics, advancing the development of scale‐aware parameterizations and deepening our understanding of cloud and precipitation dynamics.

Environmental sciences↗

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗

ERF: Energy Research and Forecasting

The Energy Research and Forecasting (ERF) code is a new model that simulates the mesoscale and microscale dynamics of the atmosphere using the latest high-performance computing architectures. It employs hierarchical parallelism using an MPI+X model, where X may be OpenMP on multicore CPU-only systems, or CUDA, HIP, or SYCL on GPU-accelerated systems. ERF is built on AMReX (Zhang et al., 2019, 2021), a block-structured adaptive mesh refinement (AMR) software framework that provides the underlying performance-portable software infrastructure for block-structured mesh operations. The "energy" aspect of ERF indicates that the software has been developed with renewable energy applications in mind. In addition to being a numerical weather prediction model, ERF is designed to provide a flexible computational framework for the exploration and investigation of different physics parameterizations and numerical strategies, and to characterize the flow field that impacts the ability of wind turbines to extract wind energy. The ERF development is part of a broader effort led by the US Department of Energy's Wind Energy Technologies Office.

17 WIND ENERGY↗

Numerical modeling of the hybrid plume plasma rocket

This paper reports on ongoing research in the area of numerical simulations of hybrid plume dynamics, and describes the basic mechanisms for energy and particle transport. The solution approach assumes cylindrical geometry and includes a multiplicity of atomic reactions, as well as the presence of a strong magnetic field. The principal reactions are electron impact ionization, and resonant charge exchange between hot ions and cold neutrals. Radial particle and energy transport is mainly by diffusion. The present description includes a modified Bohm diffusion model for plasma in the core of the plume, and classical neutral particle diffusion in the cooler regions of the flow. Neutrals are allowed to 'free stream' in the low density regions, where the collision mean-free-path becomes comparable, or even larger than the characteristic dimensions of the system. A parameteric scan of the various regimes of interest has been completed and the most significant results are presented.

Chang, F. R.↗

Offline GCSS Intercomparison of Cloud-Radiation Interaction and Surface Fluxes

Simulations of deep tropical clouds by both cloud-resolving models (CRMs) and single-column models (SCMs) in the GEWEX Cloud System Study (GCSS) Working Group 4 (WG4; Precipitating Convective Cloud Systems), Case 2 (19-27 December 1992, TOGA-COARE IFA) have produced large differences in the mean heating and moistening rates (-1 to -5 K and -2 to 2 grams per kilogram respectively). Since the large-scale advective temperature and moisture "forcing" are prescribed for this case, a closer examination of two of the remaining external types of "forcing", namely radiative heating and air/sea hear and moisture transfer, are warranted. This paper examines the current radiation and surface flux of parameterizations used in the cloud models participating in the GCSS WG4, be executing the models "offline" for one time step (12 s) for a prescribed atmospheric state, then examining the surface and radiation fluxes from each model. The dynamic, thermodynamic, and microphysical fluids are provided by the GCE-derived model output for Case 2 during a period of very active deep convection (westerly wind burst). The surface and radiation fluxes produced from the models are then divided into prescribed convective, stratiform, and clear regions in order to examine the role that clouds play in the flux parameterizations. The results suggest that the differences between the models are attributed more to the surface flux parameterizations than the radiation schemes.

Tao, W.-K.↗

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗

Specifications of a Plasmasphere Modeling Code for GGCM

The Dynamic Global Core Plasma Model (DGCPM) is a parameterized model for core or thermal plasma in the magnetosphere. The model accounts for dayside ionospheric outflow and nightside inflow. It accounts for the global pattern of convection and corotation. The model is capable of being coupled to ring current and superthermal electron models for the purpose of providing thermal plasma spatial distributions and for the purpose of accepting the dynamic influences of these plasma populations back upon the thermal plasma. The DGCPM is designed to operate alone or to operate as part of a larger integrated package. The convection electric field and magnetic field used within the DGCPM can be shared with models of other plasma populations, in addition to the exchange of parameters important to the collective modeling of whole plasma systems in the inner magnetosphere. This talk will present the features of the DGCPM model code and the various forms of information that can be exchanged with other cooperating codes.

Gallagher, D. L.↗

The effect of latent heat release on synoptic-to-planetary wave interactions and its implication for satellite observations: Theoretical modeling

The project objectives are to develop process models to investigate the interaction of planetary and synoptic-scale waves including the effects of latent heat release (precipitation), nonlinear dynamics, physical and boundary-layer processes, and large-scale topography; to determine the importance of latent heat release for temporal variability and time-mean behavior of planetary and synoptic-scale waves; to compare the model results with available observations of planetary and synoptic wave variability; and to assess the implications of the results for monitoring precipitation in oceanic-storm tracks by satellite observing systems. Researchers have utilized two different models for this project: a two-level quasi-geostrophic model to study intraseasonal variability, anomalous circulations and the seasonal cycle, and a 10-level, multi-wave primitive equation model to validate the two-level Q-G model and examine effects of convection, surface processes, and spherical geometry. It explicitly resolves several planetary and synoptic waves and includes specific humidity (as a predicted variable), moist convection, and large-scale precipitation. In the past year researchers have concentrated on experiments with the multi-level primitive equation model. The dynamical part of that model is similar to the spectral model used by the National Meteorological Center for medium-range forecasts. The model includes parameterizations of large-scale condensation and moist convection. To test the validity of results regarding the influence of convective precipitation, researchers can use either one of two different convective schemes in the model, a Kuo convective scheme or a modified Arakawa-Schubert scheme which includes downdrafts. By choosing one or the other scheme, they can evaluate the impact of the convective parameterization on the circulation. In the past year researchers performed a variety of initial-value experiments with the primitive-equation model. Using initial conditions typical of climatological winter conditions, they examined the behavior of synoptic and planetary waves growing in moist and dry environments. Surface conditions were representative of a zonally averaged ocean. They found that moist convection associated with baroclinic wave development was confined to the subtropics.

Branscome, Lee E.↗

DEM Solutions Develops Answers to Modeling Lunar Dust and Regolith

With the proposed return to the Moon, scientists like NASA-KSC's Dr. Calle are concerned for a number of reasons. We will be staying longer on the planet's surface, future missions may include dust-raising activities, such as excavation and handling of lunar soil and rock, and we will be sending robotic instruments to do much of the work for us. Understanding more about the chemical and physical properties of lunar dust, how dust particles interact with each other and with equipment surfaces and the role of static electricity build-up on dust particles in the low-humidity lunar environment is imperative to the development of technologies for removing and preventing dust accumulation, and successfully handling lunar regolith. Dr. Calle is currently working on the problems of the electrostatic phenomena of granular and bulk materials as they apply to planetary surfaces, particularly to those of Mars and the Moon, and is heavily involved in developing instrumentation for future planetary missions. With this end in view, the NASA Kennedy Space Center's Innovative Partnerships Program Office partnered with OEM Solutions, Inc. OEM Solutions is a global leader in particle dynamics simulation software, providing custom solutions for use in tackling tough design and process problems related to bulk solids handling. Customers in industries such as pharmaceutical, chemical, mineral, and materials processing as well as oil and gas production, agricultural and construction, and geo-technical engineering use OEM Solutions' EDEM(TradeMark) software to improve the design and operation of their equipment while reducing development costs, time-to-market and operational risk. EDEM is the world's first general-purpose computer-aided engineering (CAE) tool to use state-of-the-art discrete element modeling technology for the simulation and analysis of particle handling and manufacturing operations. With EDEM you'can quickly and easily create a parameterized model of your granular solids system. Computer-aided design (CAD) models of real particles can be imported to obtain an accurate representation of their shape. EDEM(TradeMark) uses particle-scale behavior models to simulate bulk solids behavior. In addition to particle size and shape, the models can account for physical properties of particles along with interaction between particles and with equipment surfaces and surrounding media, as needed to define the physics of a particular process.

Dunn, Carol Anne↗

A long short-term memory embedding for hybrid uplifted reduced order models

In this paper, we introduce an uplifted reduced order modeling (UROM) approach through the integration of standard projection based methods with long short-term memory (LSTM) embedding. Our approach has three modeling layers or components. In the first layer, we utilize an intrusive projection approach to model dynamics represented by the largest modes. The second layer consists of an LSTM model to account for residuals beyond this truncation. This closure layer refers to the process of including the residual effect of the discarded modes into the dynamics of the largest scales. However, the feasibility of generating a low rank approximation tails off for higher Kolmogorov n -width systems due to the underlying nonlinear processes. The third uplifting layer, called super-resolution, addresses this limited representation issue by expanding the span into a larger number of modes utilizing the versatility of LSTM. Therefore, our model integrates a physics-based projection model with a memory embedded LSTM closure and an LSTM based super-resolution model. In several applications, we exploit the use of Grassmann manifold to construct UROM for unseen conditions. We performed numerical experiments by using the Burgers and Navier-Stokes equations with quadratic nonlinearity. Finally, our results show robustness of the proposed approach in building reduced order models for parameterized systems and confirm the improved trade-off between accuracy and efficiency.

42 ENGINEERING↗

Carbon–Nutrient Economy of the Rhizosphere: Improving Biogeochemical Prediction and Scaling Feedbacks from Ecosystem to Regional Scales

Our project has advanced the science of plant-soil-microbial dynamics across these areas: i) nutrient cycling and plant uptake; ii) root exudation and priming; and, iii) mycorrhizal dynamics. Our project has accomplished 5 main developments: 1) Incorporation of phosphorus cycling into the Fixation & Uptake of Nutrients (FUN 3.0) model. 2) Coupling of FUN 3.0 into the E3SM Land Model (ELM). 3) Data collection across a large mycorrhizal gradient in the US, as well data in the tropics, to parameterize, test, and validate the model. 4)Scaling up mycorrhizal association measurements across landscapes using airborne hyperspectral remote sensing data. 5) Evaluation of global carbon and nutrient cycle impacts in the Community Land Model (CLM5.0) from a suite of new global mycorrhizal association maps. Over 25 publications resulted from this project, with more continuing past the project funded lifetime. Paper highlights from most of these publications have already been submitted to the DOE paper submission online system. These publications include journals such as Science and PNAS, as well as top disciplinary journals from the Nature journals, Global Change Biology, New Phytologist, and Ecology Letters, for example. Our project also contributed to improving the process representation, capabilities, and accuracy of the DOE ELM. Overall, this project significantly advanced the science of belowground plant-soil-microbial interactions as well as technical capabilities from remote sensing to modeling.

59 BASIC BIOLOGICAL SCIENCES↗

Aero-Propulsive Modeling for eVTOL Aircraft Using Wind Tunnel Testing with Multisine Inputs

A novel approach for modeling the aero-propulsive characteristics of an electric vertical takeoff and landing (eVTOL) aircraft was developed and demonstrated in wind tunnel testing. The approach was applied to the NASA LA-8 tandem tilt-wing eVTOL aircraft, using an efficient, hybrid experiment design composed of a static I-optimal response surface design for slowly-varying test variables, and dynamic orthogonal phase-optimized multisine excitations for the control surfaces and electric propulsors. Both the static and dynamic experiment designs were executed simultaneously to collect informative data for model identification. Statistically-weighted multivariate orthogonal function modeling was used to combine local modeling results computed in the frequency domain using data collected with dynamic excitations operating on the control effectors to form an aggregate aero-propulsive model. The final identified model exhibited good predictive capability when compared to validation data acquired separately from the data used to develop the model. The required test time using these new techniques was reduced by at least a factor of five compared to previous static wind tunnel testing for the LA-8 aircraft, while providing more informative data, greater parameterization flexibility, and high-quality models.

system identification↗

Building Intelligent Cyberinfrastructure to Learn Iteratively from both Observations and Models for Understanding Watershed Dynamics

Focal Area(s): Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising of a hierarchy of models (e.g., AI-driven model/component/parameterization selection). Science Challenge: Watershed processes, such as the fate and transport of sediment, carbon and nutrients across landscapes and their fluxes to water bodies (e.g., streams, rivers and lakes), have important implications for global and regional carbon and nutrient dynamics, biogeochemical functioning of terrestrial ecosystems, and soil functions. The magnitude of lateral surface/subsurface transport and fluxes of sediment, carbon and nutrients are key factors controlling the vulnerability of watersheds to climate extremes such as droughts, wildfires, and floods. Recent field observations and other scientific evidence suggest that the magnitudes of lateral transport and fluxes of sediment, carbon and nutrients are governed primarily by the spatial and vertical heterogeneity of landscape and soil properties and by pedogenic processes. However, the current generation of land surface and watershed models do not mechanistically couple the terrestrial and hydrologic systems, nor do they represent sufficiently the spatial and vertical heterogeneity of land surface and subsurface properties. On the other hand, increasing complexity of coupled watershed and land surface models requires more data to parameterize, calibrate and validate. Remote sensing (RS) provides a means to acquire spatial data and characterize their heterogeneity at the watershed scale, overcoming a major limitation associated with conventional point measurements. To improve the representation of land-surface and surface/subsurface process coupling and sub-grid heterogeneity in watershed models, it is essential to build our predictive understanding by learning from both the multi-scale multi-process modeling and diverse multi-scale data while leveraging powerful artificial intelligence (AI) techniques.

54 ENVIRONMENTAL SCIENCES↗

Stratiform clouds and their interaction with atmospheric motion

During 1989 and 1990, the researchers saw the publication of two papers and the submission of a third for review on work supported primarily by the previous contract, NAS8-36150; the delivery of an invited talk at the SIAM Conference on Dynamical Systems in Orlando, Florida; and the start of two new projects on the radiative effects of stratocumulus on the large-scale flow. The published papers discuss aspects of stratocumulus circulations (Laufersweiler and Shirer, 1989) and the Hadley to Rossby regime transition in rotating spherical systems (Higgins and Shirer, 1990). The submitted paper (Haack and Shirer, 1990) discusses a new nonlinear model of roll circulations that are forced both dynamically and thermally. The invited paper by H. N. Shirer and R. Wells presented an objective means for determining appropriate truncation levels for low-order models of flows involving two incommensurate periods; this work has application to the Hadley to Rossby transition problem in quasi-geostrophic flows (Moroz and Holmes, 1984). The new projects involve the development of a multi-layered quasi-geostrophic channel model for study of the modulation of the large-scale flow by stratocumulus clouds that typically develop off the coasts of continents. In this model the diabatic forcing in the lowest layer will change in response to the (parameterized) development of extensive fields of stratocumulus clouds. To guide creation of this parameterization scheme, researchers are producing climatologies of stratocumulus frequency and the authors correlate these frequencies with the phasing and amplitude of the large-scale flow pattern. Researchers discuss the above topics in greater detail.

Clark, John H. E.↗

Mass Gains of the Antarctic Ice Sheet Exceed Losses

During 2003 to 2008, the mass gain of the Antarctic ice sheet from snow accumulation exceeded the mass loss from ice discharge by 49 Gt/yr (2.5% of input), as derived from ICESat laser measurements of elevation change. The net gain (86 Gt/yr) over the West Antarctic (WA) and East Antarctic ice sheets (WA and EA) is essentially unchanged from revised results for 1992 to 2001 from ERS radar altimetry. Imbalances in individual drainage systems (DS) are large (-68% to +103% of input), as are temporal changes (-39% to +44%). The recent 90 Gt/yr loss from three DS (Pine Island, Thwaites-Smith, and Marie-Bryd Coast) of WA exceeds the earlier 61 Gt/yr loss, consistent with reports of accelerating ice flow and dynamic thinning. Similarly, the recent 24 Gt/yr loss from three DS in the Antarctic Peninsula (AP) is consistent with glacier accelerations following breakup of the Larsen B and other ice shelves. In contrast, net increases in the five other DS of WA and AP and three of the 16 DS in East Antarctica (EA) exceed the increased losses. Alternate interpretations of the mass changes driven by accumulation variations are given using results from atmospheric-model re-analysis and a parameterization based on 5% change in accumulation per degree of observed surface temperature change. A slow increase in snowfall with climate waRMing, consistent with model predictions, may be offsetting increased dynamic losses.

Zwally, H. Jay↗

Optimization of an Impedance-Matched Test Fixture with the Modal Projection Error

Across many industries and engineering disciplines, components and systems are designed and deployed into their operational environment of intended use. It is the desire of the design agency to be able to predict whether their component or system will function in its shock and vibration environments or if it will fail due to mechanical stresses. One method to determine if the component will survive the shock and vibration environments is to expose the component to the operational environment in a laboratory. One difficulty in executing a representative laboratory test is that the component may not have the same boundary condition in the laboratory as in the operational configuration. This paper examines the use of parameterized optimization to design the test fixture in order to better match the operational configuration. Several frequency and modal-based objective functions are examined for the optimization. The study shows that the Modal Projection Error objective function performs the best of the functions studied. The efficacy of the Modal Projection Error is demonstrated with respect to dynamic test fixture design on several analytical and experimental exemplars.

dynamic↗

Adaptive pruning-based optimization of parameterized quantum circuits

Abstract Variational hybrid quantum–classical algorithms are powerful tools to maximize the use of noisy intermediate-scale quantum devices. While past studies have developed powerful and expressive ansatze, their near-term applications have been limited by the difficulty of optimizing in the vast parameter space. In this work, we propose a heuristic optimization strategy for such ansatze used in variational quantum algorithms, which we call ‘parameter-efficient circuit training (PECT)’. Instead of optimizing all of the ansatz parameters at once, PECT launches a sequence of variational algorithms, in which each iteration of the algorithm activates and optimizes a subset of the total parameter set. To update the parameter subset between iterations, we adapt the Dynamic Sparse Reparameterization scheme which was originally proposed for training deep convolutional neural networks. We demonstrate PECT for the Variational Quantum Eigensolver, in which we benchmark unitary coupled-cluster ansatze including UCCSD and k -UpCCGSD, as well as the Low-Depth Circuit Ansatz (LDCA), to estimate ground state energies of molecular systems. We additionally use a layerwise variant of PECT to optimize a hardware-efficient circuit for the Sycamore processor to estimate the ground state energy densities of the one-dimensional Fermi-Hubbard model. From our numerical data, we find that PECT can enable optimizations of certain ansatze that were previously difficult to converge and more generally can improve the performance of variational algorithms by reducing the optimization runtime and/or the depth of circuits that encode the solution candidate(s).

Physics↗