Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural Model 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

Subsystem Estimation: A Modal Decomposition Approach

In structural analysis it is often the requirement that the response of a critical subsystem must be extracted from the response of the overall structure even when not directly excited or sometimes even measured. Based on this requirement, a number of signal processing approaches evolve that can provide a meaningful solution to this troubling problem; however, there is a variety of prior information (parameters or data) that must be available. For instance, knowledge of subsystem bandwidth may be enough to design a simple filter if there is little or no overlap with other component systems. If this is not satisfactory and a specific spectral response is required, then a model-based approach can be applied using an available free-free response of the subsystem from calibration or acceptance tests. However, when assembled and coupled with the overall system, the subsystem boundary conditions have been altered and therefore its response is altered as well.

42 ENGINEERING↗

Progressive Hedging Decomposition for Solutions of Large-Scale Process Family Design Problems

In previous work, we have introduced a mathematical model for solving a discretized version of the process family design problem. This involves two sets of decision variables. One set selects which unit module designs are included in the process platform out of a candidate set of options; the other set determines which of these unit module designs are assigned to each variant. In this work, we exploit a parallelized Progressive Hedging (PH) algorithm to solve even larger scale design problems. PH is a well-known algorithm traditionally used to solve stochastic programming problems. While our problem is not a two-stage stochastic programming problem, the structure is similar, and it can be directly mapped to the PH approach, which we employ here to solve this deterministic optimization problem. We decompose our problem by process variant. We treat the platform unit module design variables as first-stage and the assignment of unit module designs to variants as second-stage, solving the problem using mpi-sppy. We demonstrate this approach on case studies of CC, water desalination, and refrigeration.

Stinchfield, Georgia↗

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗

Seismic fault weakening via CO 2 pressurization enhanced by mechanical deformation of dolomite fault gouges

Carbon dioxide emissions from dolomite decarbonation play an essential role in the weakening of carbonate faults by lowering the effective normal stress, which is thermally activated at temperatures above 600–700 °C. However, the mechanochemical effect of low-crystalline ultrafine fault gouge on the decarbonation and slip behavior of dolomite-bearing faults remains unclear. In this study, we obtained a series of artificial dolomite fault gouges with systematically varying particle sizes and dolomite crystallinities using a high-energy ball mill. The laboratory-scale pulverization of dolomite yielded MgO at temperatures below 50 °C, indicating that mechanical decarbonation without significant heating occurred due to the collapse of the crystalline structure, as revealed by X-ray diffraction and solid-state nuclear magnetic resonance results. Furthermore, the onset temperature of thermal decarbonation decreased to ~400 °C. Numerical modeling reproduced this two-stage decarbonation, where the pore pressure increased due to low-temperature thermal decarbonation, leading to slip weakening on the fault plane even at 400–500 °C; i.e., 200–300 °C lower than previously reported temperatures. Thus, the presence of small amounts of low-crystalline dolomite in a fault plane may lead to a severely reduced shear strength due to thermal decomposition at ~400 °C with a small slip weakening distance.

58 GEOSCIENCES↗

Convergence of Weak-SINDy Surrogate Models

In this paper, we give an in-depth error analysis for surrogate models generated by a variant of the Sparse Identification of Nonlinear Dynamics (SINDy) method. We start with an overview of a variety of nonlinear system identification techniques, namely SINDy, weak-SINDy, and the occupation kernel method. Under the assumption that the dynamics are a finite linear combination of a set of basis functions, these methods establish a linear system to recover coefficients. We illuminate the structural similarities between these techniques and establish a projection property for the weak-SINDy technique. Following the overview, we analyze the error of surrogate models generated by a simplified version of weak-SINDy. In particular, under the assumption of boundedness of a composition operator given by the solution, we show that (i) the surrogate dynamics converges towards the true dynamics and (ii) the solution of the surrogate model is reasonably close to the true solution. Finally, as an application, we discuss the use of a combination of weak-SINDy surrogate modeling and proper orthogonal decomposition (POD) to build a surrogate model for partial differential equations (PDEs).

97 MATHEMATICS AND COMPUTING↗

Tensor force role in β decays analyzed within the Gogny-interaction shell model

The half-life of the famous C 14 β decay is anomalously long, with different mechanisms: the tensor force, cross-shell mixing, and three-body forces, proposed to explain the cancellations that lead to a small transition matrix element. In this study, we revisit and analyze the role of the tensor force for the β decay of C 14 as well as of neighboring isotopes. We add a tensor force to the Gogny interaction, and derive an effective Hamiltonian for shell-model calculations. The calculations were carried out in a p – s d model space to investigate cross-shell effects. Furthermore, we decompose the wave functions according to the total orbital angular momentum L in order to analyze the effects of the tensor force and cross-shell mixing. The inclusion of the tensor force significantly improves the shell-model calculations of the β -decay properties of carbon isotopes. In particular, the anomalously slow β decay of C 14 can be explained by the isospin T = 0 part of the tensor force, which changes the components of N 14 with the orbital angular momentum L = 0 , 1 , and results in a dramatic suppression of the Gamow-Teller transition strength. At the same time, the description of other nearby β decays are improved. Decomposition of wave function into L components illuminates how the tensor force modifies nuclear wave functions, in particular suppression of β -decay matrix elements. Cross-shell mixing also has a visible impact on the β -decay strength. Inclusion of the tensor force does not seem to significantly change, however, binding energies of the nuclei within the phenomenological interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Metabolic tradeoffs and heterogeneity in microbial responses to temperature determine the fate of litter carbon in simulations of a warmer world

Abstract. Climate change has the potential to destabilize the Earth's massive terrestrial carbon (C) stocks, but the degree to which models project this destabilization to occur depends on the kinds and complexities of microbial processes they simulate. Of particular note is carbon use efficiency (CUE), which determines the fraction of C processed by microbes that is anabolized into microbial biomass rather than lost to the atmosphere and soil as carbon dioxide and extracellular products. The temperature sensitivity of CUE is often modeled as an intrinsically fixed (homogeneous) property of the community, which contrasts with empirical data and has unknown impacts on projected changes to the soil C cycle under global warming. We used the Decomposition Model of Enzymatic Traits (DEMENT) – which simulates taxon-level litter decomposition dynamics – to explore the effects of introducing organism-level heterogeneity into the CUE response to temperature for decomposition of leaf litter under 5 ∘C of warming. We found that allowing the CUE temperature response to differ between taxa facilitated increased loss of litter C, unless fungal taxa were specifically restricted to decreasing CUE with temperature. Litter C loss was exacerbated by variable and elevated CUE at higher temperature, which effectively lowered costs for extracellular enzyme production. Together these results implicate a role for diversity of taxon-level CUE responses in driving the fate of litter C in a warmer world within DEMENT, which should be explored within the framework of additional model structures and validated with empirical studies.

54 ENVIRONMENTAL SCIENCES↗

Electronic Structure Calculations of the Fundamental Interactions in Solvent Extraction Desalination

Solvent extraction desalination (SED) is a promising membrane-free desalination technology, but new solvents are still needed that can improve the molecular-level water selectivity, while minimizing solvent contamination in the aqueous phase. Improvements in solvent design for SED require a fundamental understanding of the molecular interactions governing solvation behavior. Herein, quantum chemical calculations are used to delineate different fundamental interactions in water-solvent systems involving 38 different solvents, with a particular focus on the role of the electrostatic potential characteristics. Furthermore, these interactions are rigorously analyzed using the independent gradient model based on Hirshfeld partition (IGMH), Symmetry-Adapted Perturbation Theory (SAPT), as well as the energy decomposition analysis based on force field (EDA-FF) to compare against a standard force field model. The solvation free energies and partition coefficients are also calculated to estimate water extraction performance in a bulk solvent system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigations on the hierarchy of reference frames in geodesy and geodynamics

Problems related to reference directions were investigated. Space and time variant angular parameters are illustrated in hierarchic structures or towers. Using least squares techniques, model towers of triads are presented which allow the formation of linear observation equations. Translational and rotational degrees of freedom (origin and orientation) are discussed along with and the notion of length and scale degrees of freedom. According to the notion of scale parallelism, scale factors with respect to a unit length are given. Three-dimensional geodesy was constructed from the set of three base vectors (gravity, earth-rotation and the ecliptic normal vector). Space and time variations are given with respect to a polar and singular value decomposition or in terms of changes in translation, rotation, deformation (shear, dilatation or angular and scale distortions).

Grafarend, E. W.↗

Coherent structures

In order to develop more quantitative measures of coherent structures that would have comparative value over a range of experiments, it is essential that such measures be independent of the observer. It is only through such a general framework that theories with a fundamental predictive value can be developed. The triple decomposition phi = bar-phi + phi(c) + phi(r) (where bar-phi is the mean, phi(c) is the coherent part, and phi(r) is the random part of any turbulent field phi) serves this purpose. The equations of motion for the mean and coherent flow fields, based on the triple decomposition, are presented and modeling methods for the time-averaged and phase-averaged Reynolds stress are discussed.

Speziale, Charles G.↗

Concurrent Probabilistic Simulation of High Temperature Composite Structural Response

A computational structural/material analysis and design tool which would meet industry's future demand for expedience and reduced cost is presented. This unique software 'GENOA' is dedicated to parallel and high speed analysis to perform probabilistic evaluation of high temperature composite response of aerospace systems. The development is based on detailed integration and modification of diverse fields of specialized analysis techniques and mathematical models to combine their latest innovative capabilities into a commercially viable software package. The technique is specifically designed to exploit the availability of processors to perform computationally intense probabilistic analysis assessing uncertainties in structural reliability analysis and composite micromechanics. The primary objectives which were achieved in performing the development were: (1) Utilization of the power of parallel processing and static/dynamic load balancing optimization to make the complex simulation of structure, material and processing of high temperature composite affordable; (2) Computational integration and synchronization of probabilistic mathematics, structural/material mechanics and parallel computing; (3) Implementation of an innovative multi-level domain decomposition technique to identify the inherent parallelism, and increasing convergence rates through high- and low-level processor assignment; (4) Creating the framework for Portable Paralleled architecture for the machine independent Multi Instruction Multi Data, (MIMD), Single Instruction Multi Data (SIMD), hybrid and distributed workstation type of computers; and (5) Market evaluation. The results of Phase-2 effort provides a good basis for continuation and warrants Phase-3 government, and industry partnership.

Abdi, Frank↗

Precision Calcination Mechanism of CaCO 3 to High-Porosity Nanoscale CaO CO 2 Sorbent Revealed by Direct In Situ Observations

Deploying energy storage and carbon capture at scale is hindered by the substantial endothermic penalty of decomposing CaCO 3 to CaO and CO 2 , and the rapid loss of CO 2 absorption capacity by CaO sorbent particles due to sintering at the high requisite decomposition temperatures. The decomposition reaction mechanism underlying sorbent deactivation remains unclear at the atomic level and nanoscale due to past reliance on postmortem characterization methods with insufficient spatial and temporal resolution. Thus, elucidating the important CaCO 3 decomposition reaction pathway requires direct observation by time-resolved (sub-)nanoscale methods. Here, chemical and structural dynamics during the decomposition of CaCO 3 nanoparticles to nanoporous CaO particles comprising high-surface-area CaO nanocrystallites are examined. Comparing in situ transmission electron microscopy (TEM) and synchrotron X-ray diffraction experiments gives key insights into the dynamics of nanoparticle calcination, involving anisotropic CaCO 3 thermal distortion before conversion to thermally dilated energetically stable CaO crystallites. Time-resolved TEM uncovered a novel CaO formation mechanism involving heterogeneous nucleation at extended CaCO 3 defects followed by sweeping reaction front motion across the initial CaCO 3 particle. These observations clarify longstanding, yet incomplete, reaction mechanisms and kinetic models lacking accurate information about (sub-)nanoscale dynamics, while also demonstrating calcination of CaCO 3 without sintering through rapid heating and precise temperature control.

36 MATERIALS SCIENCE↗

Mechanism and kinetic study of rare earth extraction from allanite by direct acid leaching

Direct acid leaching was conducted to recover rare earth elements (REEs) from an allanite ore assaying 9,861.4 ppm of REEs. The effects of operational variables such as acid type, sulfuric acid (H 2 SO 4 ) concentration, temperature, solid/liquid (S/L) ratio, and particle size on the total REE (TREE) recovery were systematically examined. The REE mineralogy and the involved leaching mechanisms were investigated by microscopic characterization and residue re-leaching tests. It was found that approximately 80% of REEs was extracted using 1 M H 2 SO 4 at 75 °C for 2 h. Increasing temperature significantly improved the REE extraction, while changing acid type, H 2 SO 4 concentration, and particle size did not. The TREE recovery also gradually decreased at higher S/L ratios. The low fraction of Si leached (<4.0%) suggested that the extraction did not significantly alter the crystal structure of the dominant mineral. The leaching kinetics of REEs were fast within the first 10 min, followed by a much slower stage (i.e., 10-120 min). The kinetic modeling results showed that each stage may be controlled by both chemical reaction and diffusion (i.e., mix control). The activation energies for REE leaching using 1 M H 2 SO 4 in the 0-10 min and 10-120 min were determined to be 20.3 and 10.8 kJ/mol, respectively. Finally, the residue re-leaching results suggested that the fast leaching kinetics at the beginning may be explained by the preferential decomposition of metamict allanite rather than the formation of silica gel, while the following slower leaching stage was due to the dissolution of difficult-to-leach REE-bearing minerals such as well-crystallized allanite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High Pressure X-ray Diffraction and Equation of State of Hydrazine

Synchrotron X-ray diffraction has been used to investigate the structure and equation of state (EOS) of hydrazine (N 2 H 4 ) up to 54.3 GPa at 298 K. The diffraction patterns could be fit to a monoclinic unit-cell structure and put strong constraints on previously reported phase transitions documented by vibrational spectroscopy over this pressure range. Pressure–volume ( P–V ) data were fit using a Vinet EOS, yielding parameters: V 0 = 45.2 Å 3 /molecule (fixed), K 0 = 11.8(7) GPa, and K 0 ′ = 6.5(2). Previously measured high-pressure vibrational frequency shifts were used to estimate the vibrational free energy and model P–V–T isotherms from 0 to 1200 K. The results of the P–V–T isotherms are compared to existing shock Hugoniot data on hydrazine and 298 K isotherms for assemblages of possible decomposition products. This comparison suggests dissociation at high density under shock loading. Good correspondence was found between the static lattice EOS as calculated by the model and the previously reported EOS as calculated by density functional theory. Finally, these results resolve existing uncertainties about the EOS and crystal symmetry of hydrazine at high pressure and provide valuable baseline information on this important energetic material.

diffraction↗

Aeroelastic Uncertainty Quantification Studies Using the S4T Wind Tunnel Model

This paper originates from the joint efforts of an aeroelastic study team in the Applied Vehicle Technology Panel from NATO Science and Technology Organization, with the Task Group number AVT-191, titled "Application of Sensitivity Analysis and Uncertainty Quantification to Military Vehicle Design." We present aeroelastic uncertainty quantification studies using the SemiSpan Supersonic Transport wind tunnel model at the NASA Langley Research Center. The aeroelastic study team decided treat both structural and aerodynamic input parameters as uncertain and represent them as samples drawn from statistical distributions, propagating them through aeroelastic analysis frameworks. Uncertainty quantification processes require many function evaluations to asses the impact of variations in numerous parameters on the vehicle characteristics, rapidly increasing the computational time requirement relative to that required to assess a system deterministically. The increased computational time is particularly prohibitive if high-fidelity analyses are employed. As a remedy, the Istanbul Technical University team employed an Euler solver in an aeroelastic analysis framework, and implemented reduced order modeling with Polynomial Chaos Expansion and Proper Orthogonal Decomposition to perform the uncertainty propagation. The NASA team chose to reduce the prohibitive computational time by employing linear solution processes. The NASA team also focused on determining input sample distributions.

Nikbay, Melike↗

Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds

The organic compounds found in carbonaceous chondrite meteorites provide insight into primordial solar system chemistry. Evaluating the formation and decomposition mechanisms of meteoritic amino acids may aid our understanding of the origins of life and homochirality on Earth. The amino acid glycine is widespread in meteorites and other extraterrestrial environments; other amino acids, such as isovaline, are found with enantiomeric excesses in some meteorites. The relationship between meteoritic amino acids and other compounds with similar molecular structures, such as aliphatic monoamines and monocarboxylic acids is unclear; experimental results evaluating the decomposition of amino acids have produced inconclusive results about the preferred pathways, reaction intermediates, and if the conditions applied may be compatible with those occurring inside meteoritic parent bodies. In this work, we performed extensive tandem metadynamics, umbrella sampling, and committor analysis to simulate the neutral mild hydrothermal decomposition mechanisms of glycine and isovaline and put them into context for the origins of meteoritic organic compounds. Our ab initio simulations aimed to determine free energy profiles and decomposition pathways for glycine and isovaline. We found that under our modeled conditions, methylammonium, glycolic acid, and sec-butylamine are the most likely decomposition products. These results suggest that meteoritic aliphatic monocarboxylic acids are not produced from decomposition of meteoritic amino acids. Our results also indicate that the decomposition of L-isovaline prefers an enantioselective pathway resulting in the production of (S)-sec-butylamine.

Enhanced sampling↗

Insights into Acetic Acid Binding and Ketene Formation on Anatase TiO 2 (101)

Understanding the adsorption and reactivity of carboxylic acids on oxide surfaces is of great interest in catalysis for biomass upgrading via ketonization, a carbon–carbon coupling reaction. Herein, we investigate the adsorption and reaction of acetic acid on anatase TiO 2 (101) using scanning tunneling microscopy, infrared spectroscopy, temperature programmed reaction, and density functional theory calculations. We demonstrate the adsorption of acetic acid can form two intermediates: (1) dissociated, bidentate acetate with an associated bridging hydroxyl, and (2) molecular, monodentate acetic acid. The coexistence of ordered phases with increasing monolayer (ML) saturation coverages consisting of (1) pure acetate (0.5 ML), (2) mixed acetate/acetic acid (0.67 ML), (3) mixed acetate/acetic acid (1.0 ML) and (4) pure acetic acid demonstrates similar energetics for both acetate and acetic acid species. Under ultra-high vacuum conditions, the presence of both monodentate acetic acid and bidentate acetate was observed below room temperature, while solely bidentate acetate was observed up to 575 K. The deprotonation of acetic acid produces water at 280 K, while the thermal decomposition of bidentate acetate produces ketene and acetic acid at 645 K. In conclusion, this model study provides detailed insight into the stability and reactivity of carboxylic acid surface-bound intermediates, which could participate during ketonization reactions for biomass upgrading.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hyper-Differential Sensitivity Analysis of Uncertain Parameters in PDE-Constrained Optimization

Many problems in engineering and sciences require the solution of large scale optimization constrained by partial differential equations (PDEs). Though PDE-constrained optimization is itself challenging, most applications pose additional complexity, namely, uncertain parameters in the PDEs. Uncertainty quantification (UQ) is necessary to characterize, prioritize, and study the influence of these uncertain parameters. Sensitivity analysis, a classical tool in UQ, is frequently used to study the sensitivity of a model to uncertain parameters. In this article, we introduce "hyper-differential sensitivity analysis" which considers the sensitivity of the solution of a PDE-constrained optimization problem to uncertain parameters. Our approach is a goal-oriented analysis which may be viewed as a tool to complement other UQ methods in the service of decision making and robust design. We formally define hyper-differential sensitivity indices and highlight their relationship to the existing optimization and sensitivity analysis literatures. Assuming the presence of low rank structure in the parameter space, computational efficiency is achieved by leveraging a generalized singular value decomposition in conjunction with a randomized solver which converts the computational bottleneck of the algorithm into an embarrassingly parallel loop. Two multi-physics examples, consisting of nonlinear steady state control and transient linear inversion, demonstrate efficient identification of the uncertain parameters which have the greatest influence on the optimal solution.

97 MATHEMATICS AND COMPUTING↗