Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation decomposition”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Determining the Cause of Reduced Concurrent Flame Spread over Thin Solid Fuels in Low Pressure and Low Gravity

The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.

microgravity↗
DOE OSTI2022

Investigations into Plasma-Mediated Decomposition of Organoiodide Species as a Pretreatment for Mitigation of Radioiodine Emissions

A plasma-mediated pretreatment step is proposed to decompose methyl iodide (CH 3 I) for facilitating the capture of the iodine radioisotopes present in off-gas streams from used nuclear fuel reprocessing operations. Simulated CH 3 I gas streams were exposed to dielectric barrier discharge (DBD) plasma in a quartz glass continuous-flow reactor using copper electrodes. The kinetics of decomposition of CH 3 I in nitrogen was investigated by varying the applied voltage, inlet concentration, and residence time. The rate of decomposition was proportional to CH 3 I concentration and increased as the applied potential was increased. This effect was incorporated in the kinetics through the dependence of the rate constant on the applied electric field. Decomposition of CH 3 I was not affected by moisture, and the DBD pretreatment is likely to be similarly effective in an air environment as well; however, the results were confounded by a plasma-mediated reaction between nitrogen and oxygen.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Non-Stationary Power System Forced Oscillation Analysis using Synchrosqueezing Transform

Non-stationary forced oscillations (FOs) have been observed in power system operations. However, most detection methods assume that the frequency of FOs is stationary. In this paper, we present a methodology for the analysis of nonstationary FOs. Firstly, Fourier synchrosqueezing transform (FSST) is used to provide a concentrated time-frequency representation of the signals that allows identification and retrieval of non-stationary signal components. To continue, the Dissipating Energy Flow (DEF) method is applied to the extracted components to locate the source of forced oscillations. The methodology is tested using simulated as well as real PMU data. In conclusion, the results show that the proposed FSST-based signal decomposition provides a systematic framework for the application of DEF Method to non-stationary FOs.

42 ENGINEERING↗
DOE OSTI2022

Effective radiative forcing of anthropogenic aerosols in E3SMv1: historical changes, causality, decomposition, and parameterization sensitivities

The effective radiative forcing of anthropogenic aerosols (ERF aer ) is an important measure of the anthropogenic aerosol effects simulated by a global climate model. Here we analyze ERF aer simulated by the E3SMv1 atmosphere model using both century-long free-running atmosphere-land simulations and short nudged simulations. We relate the simulated ERF aer to characteristics of the aerosol composition and optical properties, and evaluate the relationships between key aerosol and cloud properties. In terms of historical changes from the year 1870 to 2014, our results show that the global mean anthropogenic aerosol burden and optical depth increase during the simulation period as expected, but the regional averages show large differences in the temporal evolution. The largest regional differences are found in the emission-induced evolution of the burden and optical depth of the sulfate aerosol: a strong decreasing trend is seen in the Northern Hemisphere high-latitude region after around 1970, while a continued increase is simulated in the tropics. Consequently, although the global mean anthropogenic aerosol burden and optical depth increase from 1870 to 2014, the ERF aer magnitude does not increase after around year 1970. The relationships between key aerosol and cloud properties (relative changes between preindustrial and present-day conditions) also show evident changes after 1970, diverging from the linear relationships exhibited for the period from 1870 to 2014. The ERF aer in E3SMv1 is relatively large compared to the recently published multi-model estimates; the primary reason is the large indirect aerosol effect (i.e., through aerosol-cloud interactions). Compared to other models, E3SMv1 features a stronger sensitivity of the cloud droplet effective radius to changes in the cloud droplet number concentration. Large sensitivity is also seen in the liquid cloud optical depth, which is determined by changes in both the effective radius and liquid water path. Aerosol-induced changes in liquid and ice cloud properties in E3SMv1 are found to have a strong correlation, as the evolution of anthropogenic sulfate aerosols affects both the liquid cloud formation and the homogeneous ice nucleation in cirrus clouds. The ERF aer estimates in E3SMv1 for the shortwave and longwave components are sensitive to the parameterization changes in both liquid and ice cloud processes. When the parameterization of ice cloud processes is modified, the top-of-atmosphere forcing changes in the shortwave and longwave components largely offset each other, so the net effect is negligible. This suggests that, to reduce the magnitude of the net ERF aer , it would be more effective to reduce the anthropogenic aerosol effect through liquid or mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effective radiative forcing of anthropogenic aerosols in E3SM version 1: historical changes, causality, decomposition, and parameterization sensitivities

Abstract. The effective radiative forcing of anthropogenic aerosols (ERFaer) is an important measure of the anthropogenic aerosol effects simulated by a global climate model. Here we analyze ERFaer simulated by the E3SM version 1 (E3SMv1) atmospheric model using both century-long free-running atmosphere–land simulations and short nudged simulations. We relate the simulated ERFaer to characteristics of the aerosol composition and optical properties, and we evaluate the relationships between key aerosol and cloud properties. In terms of historical changes from the year 1870 to 2014, our results show that the global mean anthropogenic aerosol burden and optical depth increase during the simulation period as expected, but the regional averages show large differences in the temporal evolution. The largest regional differences are found in the emission-induced evolution of the burden and optical depth of the sulfate aerosol: a strong decreasing trend is seen in the Northern Hemisphere high-latitude region after around 1970, while a continued increase is simulated in the tropics. The relationships between key aerosol and cloud properties (relative changes between pre-industrial and present-day conditions) also show evident changes after 1970, diverging from the linear relationships exhibited for the period of 1870–1969. In addition to the regional differences in the simulated relationships, a reduced sensitivity in cloud droplet number and other cloud properties to aerosol perturbations is seen when the aerosol perturbation is large. Consequently, the global annual mean ERFaer magnitude does not increase after around 1970. The ERFaer in E3SMv1 is relatively large compared to the recently published multi-model estimates; the primary reason is the large indirect aerosol effect (i.e., through aerosol–cloud interactions). Compared to other models, E3SMv1 features large relative changes in the cloud droplet effective radius in response to aerosol perturbations. Large sensitivity is also seen in the liquid cloud optical depth, which is determined by changes in both the effective radius and liquid water path. Aerosol-induced changes in liquid and ice cloud properties in E3SMv1 are found to have a strong correlation, as the evolution of anthropogenic sulfate aerosols affects both the liquid cloud formation and the homogeneous ice nucleation in cirrus clouds (that causes a large effect on longwave ERFaer). As suggested by a previous study, the large ERFaer appears to be one of the reasons why the model cannot reproduce the observed global mean temperature evolution in the second half of the 20th century. Sensitivity simulations are performed to understand which parameterization and/or parameter changes have a large impact on the simulated ERFaer. The ERFaer estimates in E3SMv1 for the shortwave and longwave components are sensitive to the parameterization changes in both liquid and ice cloud processes. When the parameterization of ice cloud processes is modified, the top-of-model forcing changes in the shortwave and longwave components largely offset each other, so the net effect is negligible. This suggests that, to reduce the magnitude of the net ERFaer, it would be more effective to reduce the anthropogenic aerosol effect through liquid or mixed-phase clouds.

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

Understanding Solid-Electrolyte Interphase Formation at the Lithium Metal Anode of Lithium-Sulfur Batteries

Lithium-sulfur (Li-S) batteries are one of the most promising chemistries for the next generation of energy storage devices. In a standard Li-S cell, the metallic lithium anode is an essential component of the battery due to its low density, extremely high theoretical specific capacity, and very low negative electrochemical potential. However, several challenges related to the use of Li metal have prevented this battery technology from becoming commercially available. For instance, Li-anodes are highly reactive, which results in the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. This, combined with the non-uniform deposition of Liions during plating at the anode surface, can seriously affect the performance, cycling, and safety of the battery. In addition, it has been suggested that a controlled SEI formation at the metallic Li anode can yield enhanced battery performance. Hence, a comprehensive understanding of how the SEI layer is formed can help elucidate improvements in this battery technology. In this work, the formation of the SEI layer in Li-S batteries is investigated using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. The stability of liquid organic-based electrolyte components, salt concentration, electron-rich environments, and the use of inorganic solid-state electrolyte (SSE) materials are explored in order to provide molecular-level fundamental insights into how the nature and composition of the electrolyte can alter the initial stages of the SEI formation. The products of electrolyte decomposition and electro-/chemical reaction pathways of liquid electrolyte component such as solvent and salts are predicted from AIMD simulations of electrolyte-anode interfaces. Energetics of reactions from DFT ground-state optimizations are also presented to confirm the decomposition mechanisms. Some significant differences are then drawn regarding the use of low and high salt concentrations. Simulations under electron-rich environments show additional multielectron electrochemical reactions of solvent and salt decomposition taking place due to the excess of electrons and the presence of radical anions in the solution. Finally, the initial stages of formation, stability, and main constituents of interfaces between Li-metal and Sbased electrodes with sulfide-based SSEs are characterized in detail.

25 ENERGY STORAGE↗
DOE OSTI2021

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A phase-field model for void and gas bubble superlattice formation in irradiated solids

A phase-field model to simulate the formation of both void and gas bubble superlattices is derived from a grand potential functional, assuming 1D diffusion of self-interstitial atoms. The model is capable of accounting for superlattice formation by either a nucleation and growth or spinodal decomposition mechanism; in this work, we focus on the nucleation and growth mechanism, using a discrete nucleation approach in the phase-field model. In simulations of void formation, short aligned rows of voids were initially formed, followed by growth in the size of the aligned rows, and finally leading to superlattice formation, consistent with experimental observations. Void superlattice spacing as a function of model parameters was quantified. Increased nucleation rates led to smaller superlattice spacing, while increased vacancy diffusivity led to larger superlattice spacing. In simulations of gas bubble superlattice formation, superlattice spacing decreased with increasing gas atom flux, consistent with experimental observations. Here, simulations with increasing ratio of gas atom flux to vacancy-interstitial pair production led to increased superlattice spacing, which is inconsistent with experimental observations at high gas atom flux rates; this discrepancy is attributed to the fact that the model does not account for bound gas atom-vacancy pairs that may dominate at high gas flux rates.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Computer Simulation of the far Infrared Collision Induced Absorption Spectra of Gaseous CO2

Far infrared collision induced absorption spectra of gaseous CO2 were computed using molecular dynamics simulations. The quadrupole and hexadecapole multipolar induction, through the trace, and the anisotropy of the molecular polarizability were found to be insufficient to represent properly the dipole induction mechanism. For a detailed analysis of the induction process the spectra obtained were decomposed into components resulting from different terms of the induced dipole. Based on this decomposition, all additional overlap contribution for each term was proposed. When spectra were recomputed including such overlap, good agreement between experiment and simulation was achieved over the temperature range at which measurements exist (233-400 K). The use of an anisotropic intermolecular potential was found to be of critical importance for obtaining the right shape of the far wings of the spectra.

Gruszka, Marcin↗
NASA NTRS1976

A model for the simulation of atmospheric turbulence

A method that produces realistic simulations of atmospheric turbulence is developed and analyzed. The procedure makes use of a generalized spectral analysis, often called a proper orthogonal decomposition or the Karhunen-Loeve expansion. A set of criteria, emphasizing a realistic appearance, a correct spectral shape, and non-Gaussian statistics, is selected in order to evaluate the model turbulence. An actual turbulence record is analyzed in detail, providing both a standard for comparison and input statistics for the generalized spectral analysis. The simulation utilizes the eigenfunction-expansion procedure to produce preliminary time histories of the three velocity components simultaneously. It is found that important cross-statistical features are reasonably well-behaved.

Petersen, E. L.↗
DOE OSTI2021

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Polar decomposition for attitude determination from vector observations

This work treats the problem of weighted least squares fitting of a 3D Euclidean-coordinate transformation matrix to a set of unit vectors measured in the reference and transformed coordinates. A closed-form analytic solution to the problem is re-derived. The fact that the solution is the closest orthogonal matrix to some matrix defined on the measured vectors and their weights is clearly demonstrated. Several known algorithms for computing the analytic closed form solution are considered. An algorithm is discussed which is based on the polar decomposition of matrices into the closest unitary matrix to the decomposed matrix and a Hermitian matrix. A somewhat longer improved algorithm is suggested too. A comparison of several algorithms is carried out using simulated data as well as real data from the Upper Atmosphere Research Satellite. The comparison is based on accuracy and time consumption. It is concluded that the algorithms based on polar decomposition yield a simple although somewhat less accurate solution. The precision of the latter algorithms increase with the number of the measured vectors and with the accuracy of their measurement.

Bar-Itzhack, Itzhack Y.↗
DOE OSTI2021

Modeling interval trendlines: Symbolic singular spectrum analysis for interval time series

Abstract In this article we propose an extension of singular spectrum analysis for interval‐valued time series. The proposed methods can be used to decompose and forecast the dynamics governing a set‐valued stochastic process. The resulting components on which the interval time series is decomposed can be understood as interval trendlines, cycles, or noise. Forecasting can be conducted through a linear recurrent method, and we devised generalizations of the decomposition method for the multivariate setting. The performance of the proposed methods is showcased in a simulation study. We apply the proposed methods so to track the dynamics governing the Argentina Stock Market (MERVAL) in real time, in a case study over a period of turbulence that led to discussions of the government of Argentina with the International Monetary Fund.

de Carvalho, Miguel↗
DOE OSTI2021

Insight into SEI Growth in Li-Ion Batteries using Molecular Dynamics and Accelerated Chemical Reactions

The performance of lithium-ion batteries (LIB) using organic electrolytes strongly depends on the formation of a stable solid electrolyte interphase (SEI) film. Elucidating the dynamic evolution and spatial composition of the SEI can be very useful to study the stability of SEI components and help optimize the formation cycles of LIB. We propose a classical molecular dynamics simulation protocol for predicting the first stages of SEI formation using a reaction method involving the decomposition of EC and LiPF 6 molecules in the electrolyte. We accelerate the formation of SEI components near the anode surface by increasing the probability of reactions, implemented through a geometry matching scheme, followed by a force-field reconfiguration. We observe the formation of gases (C 2 H 4 ), inorganic (Li 2 CO 3 and LiF) and organic (LEDC) components. Finally, this protocol shows promise to be able to evaluate the effects of varying electrolyte compositions and additives on SEI layer structure and composition.

25 ENERGY STORAGE↗
DOE OSTI2020

Characterizing magnetized plasmas with dynamic mode decomposition

Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear and computationally intensive and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the helicity injected torus with steady inductive (HIT-SI) experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An Algorithmic and Software Pipeline for Very Large Scale Scientific Data Compression with Error Guarantees

Efficient data compression is becoming increasingly critical for storing scientific data because many scientific applications produce vast amounts of data. This paper presents an end-to-end algorithmic and software pipeline for data compression that guarantees both error bounds on primary data (PD) and derived data, known as Quantities of Interest (QoI).We demonstrate the effectiveness of the pipeline by compressing fusion data generated by a large-scale fusion code, XGC, which produces tens of petabytes of data in a single day. We demonstrate that the compression is conducted by setting aside computational resources known as staging nodes, and does not impact the simulation performance. For efficient parallel I/O, the pipeline uses ADIOS2, which many codes such as XGC already use for their parallel I/O. We show that our approach can compress the data by two orders of magnitude while guaranteeing high accuracy on both the PD and the QoIs. Further, the amount of resources required by compression is a few percent of the resources required by simulation while ensuring that the compression time for each stage is less than the corresponding simulation time.This pipeline consists of three main steps. The first step decomposes the data using domain decomposition into small subdomains. Each subdomain is then compressed independently to achieve a high level of parallelism. The second step uses existing techniques that guarantee error bounds on the primary data for each subdomain. The third step uses a post-processing optimization technique based on Lagrange multipliers to reduce the QoI errors for data corresponding to each subdomain. The Lagrange multipliers generated can be further quantized or truncated to increase the compression level. All of the above characteristics of our approach make it highly practical to apply on-the-fly compression while guaranteeing errors on QoIs that are critical to the scientists.

Banerjee, Tania↗