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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Crystal and Electronic Structures of A 2 NaIO 6 Periodate Double Perovskites (A = Sr, Ca, Ba): Candidate Wasteforms for I-129 Immobilization

The synthesis, structure, and thermal stability of the periodate double perovskites A 2 NaIO 6 (A= Ba, Sr, Ca) were investigated in the context of potential application for the immobilization of radioiodine. A combination of X-ray diffraction and neutron diffraction, Raman spectroscopy, and DFT simulations were applied to determine accurate crystal structures of these compounds and understand their relative stability. The compounds were found to exhibit rock-salt ordering of Na and I on the perovskite B-site; Ba 2 NaIO 6 was found to adopt the Fm- 3 m aristotype structure, whereas Sr 2 NaIO 6 and Ca 2 NaIO 6 adopt the P 2 1 / n hettotype structure, characterized by cooperative octahedral tilting. DFT simulations determined the Fm- 3 m and P 2 1 / n structures of Ba 2 NaIO 6 to be energetically degenerate at room temperature, whereas diffraction and spectroscopy data evidence only the presence of the Fm- 3 m phase at room temperature, which may imply an incipient phase transition for this compound. The periodate double perovskites were found to exhibit remarkable thermal stability, with Ba 2 NaIO 6 only decomposing above 1050 °C in air, which is apparently the highest recorded decomposition temperature so far recorded for any iodine bearing compound. As such, these compounds offer some potential for application in the immobilization of iodine-129, from nuclear fuel reprocessing, with an iodine incorporation rate of 25–40 wt%. The synthesis of these compounds, elaborated here, is also compatible with both current conventional and future advanced processes for iodine recovery from the dissolver off-gas.

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

Personalized Tucker Decomposition: Modeling Commonality and Peculiarity on Tensor Data

In this paper, we propose a personalized Tucker decomposition (perTucker) to address the limitations of traditional tensor decomposition methods in capturing heterogeneity across different datasets. perTucker decomposes tensor data into shared global components and personalized local components. We introduce an order orthogonality assumption and develop a proximal gradient regularized block coordinate descent algorithm guaranteed to converge to a stationary point. The unique and common representations learned by perTucker reveal intrinsic statistical patterns in data and provide valuable information for a wide range of downstream analytics, including anomaly detection, source classification, and clustering. We demonstrate perTucker’s effectiveness through a simulation study and two case studies on solar flare detection and tonnage signal classification.

14 SOLAR ENERGY↗
DOE OSTI2024

Cold Plasma Measurements

We have continued the simulation campaign in support of our ongoing magnetospheric cold plasma research project. This project aims to develop the next-generation particle instruments to measure the properties of the cold particle populations in the Earth’s magnetosphere. For this purpose, simulations have been performed with a Particle-In-Cell (PIC) code called the Curvilinear PIC (CPIC). The code is formulated in curvilinear geometry and couples the standard PIC algorithm with algorithms for the generation and adaptation of the underlaying computational mesh. It conforms to complex objects like spacecraft and it can place more grid points in regions where higher resolution is needed. The code also features a scalable solver based on the multigrid algorithm and it is fully parallelized via domain decomposition and MPI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Aging studies of Dual functional materials for CO 2 direct air capture with in situ methanation under simulated ambient conditions: Ru thrifting for cost reduction

Dual function materials (DFMs) comprised of 0.25 %Ru, 6.1 %Na 2 O/γ-Al 2 O 3 //monolith were evaluated for about 250 hours time-on-stream (TOS) at various simulated ambient climate capture conditions followed by temperature swing methanation to 280 °C. Herein this paper focuses on the impact of thrifting Ru to low levels and the impact on performance. Results showed both stable CO 2 capture capacity and CH 4 production with 0.25 % Ru DFM deposited on a ceramic monolith. The CO 2 conversion to CH 4 production was decreased slightly by the Ru decrease from 1 % to 0.25 %. The capture capacity decreased since a lower Ru content reduces the complete decomposition of the Na 2 CO 3 precursor producing fewer active adsorption sites (“Na 2 O”). The deployment of a low Ru DFM monolith would significantly decrease the overall capital cost and give more potential for a large-scale direct air capture and methanation (DACM) application.

03 NATURAL GAS↗
DOE OSTI2023

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hyperstoichiometric Uranium Dioxides: Rapid Synthesis and Irradiation-Induced Structural Changes

Uranium dioxide (UO 2 ), the primary fuel for commercial nuclear reactors, incorporates excess oxygen forming a series of hyperstoichiometric oxides. Thin layers of these oxides, such as UO 2.12 , form readily on the fuel surface and influence its properties, performance, and potentially geologic disposal. This work reports a rapid and straightforward combustion process in uranyl nitrate–glycine–water solutions to prepare UO 2.12 nanomaterials and thin films. We also report on the investigation of the structural changes induced in the material by irradiation. Despite the simple processing aspects, the combustion synthesis of UO 2.12 has a sophisticated chemical mechanism involving several exothermic steps. Raman spectroscopy and single-crystal X-ray diffraction (XRD) measurements reveal the formation of a complex compound containing the uranyl moiety, glycine, H 2 O, and NO 3 – groups in reactive solutions and dried combustion precursors. Combustion diagnostic methods, gas-phase mass spectroscopy, differential scanning calorimetry (DSC), and extracted activation energies from DSC measurements show that the rate-limiting step of the process is the reaction of ammonia with nitrogen oxides formed from the decomposition of glycine and uranyl nitrate, respectively. However, the exothermic decomposition of the complex compound determines the maximum temperature of the process. In situ transmission electron microscopy (TEM) imaging and electron diffraction measurements show that the decomposition of the complex compound directly produces UO 2 . The incorporation of oxygen at the cooling stage of the combustion process is responsible for the formation of UO 2.12 . Spin coating of the solutions and brief annealing at 670 K allow the deposition of uniform films of UO 2.12 with thicknesses up to 300 nm on an aluminum substrate. Irradiation of films with Ar 2+ ions (1.7 MeV energy, a fluence of up to 1 × 10 17 ions/cm 2 ) shows unusual defect-simulated grain growth and enhanced chemical mixing of UO 2.12 with the substrate due to the high uranium ion diffusion in films. As a result, the method described in this work allows the preparation of actinide oxide targets for fundamental nuclear science research and studies associated with stockpile stewardship.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Breaking the curse of dimensionality: Solving configurational integrals for crystalline solids by tensor networks

Accurately evaluating configurational integrals for dense solids remains a central and difficult challenge in the statistical mechanics of condensed systems. Here, we present a tensor network approach that reformulates the high-dimensional configurational integral for identical-particle crystals into a sequence of computationally efficient summations. We represent the integrand as a high-dimensional tensor and apply tensor-train (TT) decomposition together with a custom TT-cross interpolation. This approach circumvents the need to explicitly construct the full tensor. We introduce tailored rank-1 and rank-2 schemes optimized for sharply peaked Boltzmann probability densities, typical for identical-particle crystals. When applied to the calculation of internal energy and pressure-temperature curves for crystalline Cu and Ar at high (GPa) pressures, as well as the alpha-to-beta phase transition diagram of Sn, our method accurately reproduces molecular dynamics simulation results using tight-binding, machine learning, hierarchical interacting particle–neural network, and modified embedded atom method potentials,all within seconds of computation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Surprising similarities in model and observational aerosol radiative forcing estimates

Abstract. The radiative forcing from aerosols (particularly through their interaction with clouds) remains one of the most uncertain components of the human forcing of the climate. Observation-based studies have typically found a smaller aerosol effective radiative forcing than in model simulations and were given preferential weighting in the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). With their own sources of uncertainty, it is not clear that observation-based estimates are more reliable. Understanding the source of the model and observational differences is thus vital to reduce uncertainty in the impact of aerosols on the climate. These reported discrepancies arise from the different methods of separating the components of aerosol forcing used in model and observational studies. Applying the observational decomposition to global climate model (GCM) output, the two different lines of evidence are surprisingly similar, with a much better agreement on the magnitude of aerosol impacts on cloud properties. Cloud adjustments remain a significant source of uncertainty, particularly for ice clouds. However, they are consistent with the uncertainty from observation-based methods, with the liquid water path adjustment usually enhancing the Twomey effect by less than 50 %. Depending on different sets of assumptions, this work suggests that model and observation-based estimates could be more equally weighted in future synthesis studies.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

An Efficient and Robust Singular Value Method for Star Pattern Recognition and Attitude Determination

A new star pattern recognition method is developed using singular value decomposition of a measured unit column vector matrix in a measurement frame and the corresponding cataloged vector matrix in a reference frame. It is shown that singular values and right singular vectors are invariant with respect to coordinate transformation and robust under uncertainty. One advantage of singular value comparison is that a pairing process for individual measured and cataloged stars is not necessary, and the attitude estimation and pattern recognition process are not separated. An associated method for mission catalog design is introduced and simulation results are presented.

Juang, Jer-Nan↗
NASA NTRS1986

Towards a general object-oriented software development methodology

Object diagrams were used to design a 5000 statement team training exercise and to design the entire dynamics simulator. The object diagrams are also being used to design another 50,000 statement Ada system and a personal computer based system that will be written in Modula II. The design methodology evolves out of these experiences as well as the limitations of other methods that were studied. Object diagrams, abstraction analysis, and associated principles provide a unified framework which encompasses concepts from Yourdin, Booch, and Cherry. This general object-oriented approach handles high level system design, possibly with concurrency, through object-oriented decomposition down to a completely functional level. How object-oriented concepts can be used in other phases of the software life-cycle, such as specification and testing is being studied concurrently.

Seidewitz, ED↗
NASA NTRS2007

Wavefront Control Toolbox for James Webb Space Telescope Testbed

We have developed a Matlab toolbox for wavefront control of optical systems. We have applied this toolbox to the optical models of James Webb Space Telescope (JWST) in general and to the JWST Testbed Telescope (TBT) in particular, implementing both unconstrained and constrained wavefront optimization to correct for possible misalignments present on the segmented primary mirror or the monolithic secondary mirror. The optical models implemented in Zemax optical design program and information is exchanged between Matlab and Zemax via the Dynamic Data Exchange (DDE) interface. The model configuration is managed using the XML protocol. The optimization algorithm uses influence functions for each adjustable degree of freedom of the optical mode. The iterative and non-iterative algorithms have been developed to converge to a local minimum of the root-mean-square (rms) of wavefront error using singular value decomposition technique of the control matrix of influence functions. The toolkit is highly modular and allows the user to choose control strategies for the degrees of freedom to be adjusted on a given iteration and wavefront convergence criterion. As the influence functions are nonlinear over the control parameter space, the toolkit also allows for trade-offs between frequency of updating the local influence functions and execution speed. The functionality of the toolbox and the validity of the underlying algorithms have been verified through extensive simulations.

Shiri, Ron↗
NASA NTRS2017

Order-Reduced Solution of the Nonlinear High-Fidelity Generalized Method of Cells Micromechanics Relations

The High-Fidelity Generalized Method of Cells (HFGMC) is one technique for accurately simulating nonlinear composite material behavior. The HFGMC uses a higher-order approximation for the subcell displacement field that allows for a more accurate determination of the subcell stressstrain fields at the cost of some computational efficiency. In order to reduce computational costs associated with the solution of the ensuing system of simultaneous equations, the HFGMC global system of equations for doubly-periodic repeating unit cells with nonlinear constituents was reduced in size through the use of a Petrov-Galerkin-based Proper Orthogonal Decomposition order-reduction scheme. A number of cases were presented that address the computational feasibility of using order-reduction techniques to solve solid mechanics problems involving complex microstructures.

order-reduction↗
NASA NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗
NASA NTRSDate not supplied

Cartesian Mesh Simulations for the Third AIAA Sonic Boom Prediction Workshop

Simulation results are presented for all cases from the Third AIAA Sonic Boom Prediction Workshop. An inviscid, embedded-boundary Cartesian-mesh flow solver is used in conjunction with adjoint-based mesh adaptation to compute nearfield pressure signatures. Specialized techniques are applied to maximize accuracy and minimize cost on Cartesian meshes. The Richardson-based error estimate highlights regions of the signatures most sensitive to mesh refinement. Timing results and coarse, medium, and fine mesh sizes for nearfield cases demonstrate that the parallel decomposition approach is efficient in both computational time and wall-clock. Pressure signals are propagated to the ground using an augmented Burgers’ equation solver to predict boom carpets. Ground signatures and loudness metrics are presented for a standard atmosphere as well as more realistic atmospheric profiles, which affect overall noise levels and can significantly widen the boom carpet. Mesh convergence studies show that high sampling frequencies, around 500 kHz, are required for propagation, and the sampling frequency increases at large off-track angles with longer acoustic ray paths and propagation times. The numerical methods yield accurate results for predicting low sonic boom signatures while being among the least computationally expensive of the workshop.

ARMD↗
DOE OSTI2024

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗
DOE OSTI2022

Liquid and Glass Phases of an Alkylguanidinium Sulfonate Hydrogen-Bonded Organic Framework

Glassy phases of framework materials feature unique and tunable properties that are advantageous for gas separation membranes, solid electrolytes, and phase-change memory applications. However, the structural and chemical diversity of porous frameworks that can be liquified and quenched into a glass has been limited by thermal decomposition at-or below-the high temperatures required to induce a melting transition. Utilizing a desymmetrization strategy, in this work we report a new guanidinium organosulfonate hydrogen-bonded organic framework (HOF) that melts and vitrifies below 100 °C. In this low-temperature regime, non-covalent interactions between guest molecules and the porous framework become a dominant contributor to the overall stability of the structure, resulting in unusual phase behavior such as guest dependent melting, glass, and recrystallization transitions. Through molecular dynamics simulations and pair distribution function analysis, we show that the local structure of the amorphous liquid and glass phase resembles that of the parent crystalline framework. Access to molten phases of framework materials at moderate temperatures should permit the use of more thermally sensitive functional groups and enhance the structural control and tunability that can be realized in network-forming glasses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

DRDMannTurb: A Python package for scalable, data-driven synthetic turbulence

Synthetic turbulence models (STMs) are used in wind engineering to generate realistic flow fields and are employed as inputs to industrial wind simulations. Examples include prescribing inlet conditions in large eddy simulations that model loads on wind turbines and tall buildings. We are interested in STMs capable of generating fluctuations based on prescribed second-moment statistics since such models can simulate environmental conditions that closely resemble on-site observations. To this end, the widely used Mann model (see Mann, 1994, 1998) is the inspiration for DRDMannTurb. The Mann model is described by three physical parameters: a magnitude parameter influencing the global variance of the wind field and corresponding to the Kolmogorov constant multiplied by the rate of viscous dissipation of the turbulent kinetic energy to the two-thirds, αϵ 2/3 , a turbulence length scale parameter L, and a nondimensional parameter Γ related to the lifetime of the eddies. A number of studies, as well as international standards (e.g., those by the International Electrotechnical Commission (IEC)), include recommended values for these three parameters with the goal of standardizing wind simulations according to observed energy spectra. Yet, having only three parameters, the Mann model faces limitations in accurately representing the diversity of observable spectra. This Python package enables users to extend the Mann model and more accurately fit field measurements through flexible neural network models of the eddy lifetime function. Following Keith et al. (2021), we refer to this class of models as Deep Rapid Distortion (DRD) models. DRDMannTurb also includes a general module implementing an efficient method for synthetic turbulence generation based on a domain decomposition technique. This technique is also described in Keith et al. (2021).

17 WIND ENERGY↗