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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2022

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Maximum Likelihood Spectrum Decomposition for Isotope Identification and Quantification

A spectral decomposition method has been implemented to identify and quantify isotopic source terms in high-resolution gamma-ray spectroscopy in static geometry and shielding scenarios. Monte Carlo simulations were used to build the response matrix of a shielded high-purity germanium detector monitoring an effluent stream with a Marinelli configuration. The decomposition technique was applied to a series of calibration spectra taken with the detector using a multi-nuclide standard. These results are compared with decay-corrected values from the calibration certificate. For most nuclei in the standard ( 241 Am, 109 Cd, 137 Cs, and 60 Co), the deviations from the certificate values were generally no more than 6% with a few outliers as high as 10%. Furthermore, for 57 Co, the radionuclide with the lowest activity, the deviations from the standard reached as high as 25%, driven by the meager statistics in the calibration spectra. In addition, a complete treatment of error propagation for the technique is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Across‐model spread and shrinking in predicting peatland carbon dynamics under global change

Abstract Large across‐model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across‐model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO 2 levels and five temperature levels. We developed a matrix‐based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter‐model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter‐model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO 2 (eCO 2 ) concentrations in the measurements, simulated NEP responded positively to eCO 2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across‐model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.

Hou, Enqing↗

DOE OSTI2025

Water‐mediated Decomposition Pathways of ε‐Hexanitrohexaazaisowurtzitane

The dominant initial decomposition mechanism of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (HNIW) is understood from the literature to be the N─NO 2 bond scission unimolecular mechanism. Here we have found a water-mediated HONO release initial decomposition mechanism energetically competitive with the N─NO 2 bond scission mechanism modeled using ab initio nudged elastic band simulations. The activation energy of the water-mediated HONO mechanism was calculated to be 138.9 kJ/mol with the Perdew, Burke, and Ernzerhof (PBE) functional (177.8 with PBE0), while the double NO 2 release mechanism is 171.9 kJ/mol (203.3 with PBE0). Branching secondary, tertiary, and quaternary decomposition steps were also discovered including oxidation of HONO or H 2 O causing ring-opening that leads to C─N bond breaking and N 2 O or NO release. The reaction rate of HONO oxidation is much faster than HONO release, making the HONO concentration low. The Helmholtz free energy barrier of the double NO 2 release mechanism is lower than the water-mediated HONO release above 438 K (165°C) due to the vibrational contribution to the free energy. This helps explain why the literature reports more NO 2 release at higher temperatures. The present study reveals a physical mechanism for how water can catalyze the decomposition of HNIW at low enough temperatures, and reveals secondary, tertiary, and quaternary reaction mechanisms leading to NO x release. In conclusion, as water is realistically always present, the reported barriers are important for building kinetic models needed to ensure the safety and reliability of HNIW.

Steele, Brad A. [Lawrence Livermore National Labor↗

NASA NTRS1992

Viscous simulation method for unsteady flows past multicomponent configurations

The present algorithm for the numerical simulation of flows about complex configurations (whose multiple and nonsimilar components have arbitraty geometries) employs a hybridization of the domain decomposition techniques for grid generation as well as to reduce computer-memory requirements. A fully vectorized, finite-volume, upwind-biased, approximately factored multigrid method is used to solve 3D Reynolds-averaged unsteady and compressible Navier-Stokes equations simulating supersonic flows past an ogive-nose-cylinder near or within a cavity. The time-averaged surface pressures obtained compare favorably with the wind tunnel data.

Fouladi, Kamran↗

DOE OSTI2020

Atomic scale understanding of organic anion separations using ion-exchange resins

A combination of ab initio and classical molecular dynamic simulations was used to explore the adsorption/desorption and diffusion characteristics of ion-exchange resins for extraction of organic anions in electrically-driven separation processes. We considered two classes of carboxylate mixtures that are commonly encountered in bioprocessing separations: a short chain fatty acid mixture (acetate/butyrate) and an aromatic mixture (ferulate/coumarate). The suitability of several resin materials including PFC100E, IRC86, PFA444 and IRA67 was interrogated. The decomposition of the interaction energies by the symmetry-adapted perturbation theory and the classical molecular dynamic simulations of organic diffusion together reveal that the geometries of the organic anions and the functional groups of the resins, as well as their Columbic interactions, are the controlling factors in the diffusion process of the organic compounds in these resins. Additionally, classical simulations also show that modifying the functionality of resin beads, the magnitude of electric fields, and the ratio of organic mixtures may provide an effective way to control the diffusion rate and obtain selective separation of these organic mixtures. Finally, a general suggestion for favorable separation conditions is summarized.

36 MATERIALS SCIENCE↗

NASA NTRS1959

A Hydrogen Peroxide Hot-Jet Simulator for Wind-Tunnel Tests of Turbojet-Exit Models

A turbojet-engine-exhaust simulator which utilizes a hydrogen peroxide gas generator has been developed for powered-model testing in wind tunnels with air exchange. Catalytic decomposition of concentrated hydrogen peroxide provides a convenient and easily controlled method of providing a hot jet with characteristics that correspond closely to the jet of a gas turbine engine. The problems associated with simulation of jet exhausts in a transonic wind tunnel which led to the selection of a liquid monopropellant are discussed. The operation of the jet simulator consisting of a thrust balance, gas generator, exit nozzle, and auxiliary control system is described. Static-test data obtained with convergent nozzles are presented and shown to be in good agreement with ideal calculated values.

Runckel, Jack F.↗

NASA NTRS1959

A Hydrogen Peroxide Hot-Jet Simulator for Wind-Tunnel Tests of Turbojet-Exit Models

Runckel, Jack F.↗

DOE OSTI2020

Reshock analysis for PMMA driven above the threshold for chemical decomposition

Polymethyl methacrylate (PMMA) is a transparent thermoplastic often used in shock compression science. Previous proceedings (AIP Conf. Proc. 845, 131 (2006)) described plate impact experiments in which PMMA was reshocked to pressures of up to ~130 GPa, well above that at which it decomposes on its principal Hugoniot. Because some of the conclusions of the original analysis were surprising, the results were reanalyzed in a later proceeding (ibid. 1426, 771 (2012)). In addition to our own reanalysis of the data, we performed hydrodynamic simulations of the experiments based on a new equation of state for PMMA shocked above its threshold for decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Experimental and fuel-surrogates modeling study of the high-pressure pyrolysis of specialty cetane number fuels: implications for fall-off in ethylene unimolecular dissociation

Single pulse shock tube experiments were conducted at 50 atm nominal pressure and 4 ms nominal reaction time over a temperature range of 900–1800 K, to study the pyrolysis speciation of a multi-component jet fuel, F-24, and six cetane number (CN) specialty fuels - CN30, CN35, CN40, CN45, CN50, and CN55. Gas chromatography (GC) was used to qualitatively and quantitatively analyze the post shock gases. The relationship between the formation of key pyrolysis species and the chemically controlled combustion propensity as reflected by the cetane number of each fuel was examined. A surrogate-based mechanism from the CRECK Modelling Group and chemical-functional group based optimized surrogates (CFGO) were used to simulate the pyrolysis speciation results. The model was able to capture the chemistry of most species except two important pyrolysis intermediates – ethylene and acetylene. Chemical kinetic analyses were performed to identify the important reactions which affect the chemistry of these species; however, the rate parameters of critical reactions were found to be unsuitable for simulating the present high-pressure studies. Here, to address this unsuitability, a theory-based fall-off analysis for three reactions representing the decomposition of ethylene and subsequent formation of acetylene was performed, and these are included in an updated version of the CRECK mechanism. This update resolves discrepancies between the experimental results and simulations for ethylene and acetylene. Reaction flux analyses using the updated surrogate model were also performed to identify the important reaction pathways responsible for the formation of crucial species and to provide an analysis of the chemistry of complex multi-component fuel systems. The fundamental reactions responsible for driving pyrolysis chemistry were greatly influenced by the chemical functional groups present in these fuels. In addition to updating the rate parameters of specific reactions to improve modeling, this study also emphasizes the effectiveness of the fuel-surrogate approach, where surrogates representing the chemical functional group composition of the parent fuel serve as a valuable tool for predicting the combustion chemistry of novel fuels.

Chemical Kinetics↗

NASA NTRS1993

(abstract) 3D Electromagnetic Plasma Particle Simulations

A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.

electromagnetic plasma particles simulations 3D al↗

DOE OSTI2020

TuckerMPI: A Parallel C++/MPI Software Package for Large-scale Data Compression via the Tucker Tensor Decomposition

With this study, our goal is compression of massive-scale grid-structured data, such as the multi-terabyte output of a high-fidelity computational simulation. For such data sets, we have developed a new software package called TuckerMPI, a parallel C++/MPI software package for compressing distributed data. The approach is based on treating the data as a tensor, i.e., a multidimensional array, and computing its truncated Tucker decomposition, a higher-order analogue to the truncated singular value decomposition of a matrix. The result is a low-rank approximation of the original tensor-structured data. Compression efficiency is achieved by detecting latent global structure within the data, which we contrast to most compression methods that are focused on local structure. In this work, we describe TuckerMPI, our implementation of the truncated Tucker decomposition, including details of the data distribution and in-memory layouts, the parallel and serial implementations of the key kernels, and analysis of the storage, communication, and computational costs. We test the software on 4.5 and 6.7 terabyte data sets distributed across 100 s of nodes (1,000 s of MPI processes), achieving compression ratios between 100 and 200,000×, which equates to 99--99.999% compression (depending on the desired accuracy) in substantially less time than it would take to even read the same dataset from a parallel file system. Moreover, we show that our method also allows for reconstruction of partial or down-sampled data on a single node, without a parallel computer so long as the reconstructed portion is small enough to fit on a single machine, e.g., in the instance of reconstructing/visualizing a single down-sampled time step or computing summary statistics. The code is available at https://gitlab.com/tensors/TuckerMPI.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Stabilizing Nickel‐Rich Cathodes in Aqueous Process through Nanocellulose as Water Barrier

Nickel-rich LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NMC 811) cathode offers high voltage and high specific capacity, making it promising for high energy density batteries. However, large-scale manufacturing of aqueous-processed NMC 811 electrodes remains challenging due to proton exchange causing material decomposition and capacity loss. This work addresses this issue by constructing an in situ nanocellulose protective layer for NMC 811 particles via electrostatic interactions during the slurry preparation. For the first time, the interatomic spacing between inter-chains of nanocellulose is measured through wide-angle X-ray scattering and demonstrate the ability to effectively confine interlayer water using atomistic simulations. Moreover, this nanocellulose coverage simultaneously minimizes Li + surface segregation and mitigates water infiltration. Owing to less material decomposition during the aqueous processing, nanocellulose-protected NMC electrodes exhibit higher initial coulombic efficiency (83% vs 62% at 0.1C) and capacity (133 vs 59 mAh g −1 at 6C) than unprotected electrodes. Additionally, optimized aqueous-processed NMC electrodes offer comparable or even superior electrochemical properties compared to the electrodes fabricated using the conventional toxic organic solvent, N-methyl-2-pyrrolidone. Consequently, the developed approach enables affordable, sustainable aqueous processing for Nickel-rich NMC 811 cathodes with excellent electrochemical performances.

25 ENERGY STORAGE↗

DOE OSTI2024

Initiation and Carbene Induced Radical Chain Reactions in CH 2 F 2 Pyrolysis

High temperature dissociations of organic molecules typically involve a competition between radical and molecular processes. In this work, we use a modeling, experiment, theory (MET) framework to characterize the high temperature thermal dissociation of CH 2 F 2 , a flammable hydrofluorocarbon (HFC) that finds widespread use as a refrigerant. Initiation in CH 2 F 2 proceeds via a molecular elimination channel; CH 2 F 2 →CHF+HF. Here we show that the subsequent self-reactions of the singlet carbene, CHF, are fast multichannel processes and a facile source of radicals that initiate rapid chain propagation reactions. These have a marked influence on the decomposition kinetics of CH 2 F 2 . The inclusion of these reactions brings the simulations into better agreement with the present and literature experiments. Additionally, flame simulations indicate that inclusion of the CHF+CHF multichannel reaction leads to a noticeable enhancement in predictions of laminar flame speeds, a key parameter that is used to determine the flammability of a refrigerant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Structure–Property Effects in the Generation of Transient Aqueous Benzoic Acid Anhydrides by Carbodiimide Fuels

The design of dissipative systems, which operate out-of-equilibrium by consuming chemical fuels, is difficult. As yet, there are a few examples of privileged fuel chemistries that can be broadly applied in abiotic systems in the same way that ATP hydrolysis is exploited throughout biochemistry. The key issue is that designing nonequilibrium systems is inherently about balancing the relative rates of coupled reactions. The use of carbodiimides as fuels to generate transient aqueous carboxylic anhydrides has recently been used in examples of new nonequilibrium materials and supramolecular assemblies. Here, we explore the kinetics of formation and decomposition of a series of benzoic anhydrides generated from the corresponding acids and EDC under typical conditions (EDC = N-(3-dimethylaminopropyl)-N'-ethylcarbodiimide hydrochloride). The reactions can be described by a simple mechanism that merges known behavior for the two processes independently. Structure–property effects in these systems are dominated by differences in the anhydride decomposition rate. The kinetic parameters allow trends in concentration-dependent properties to be simulated, such as reaction lifetimes, peak anhydride concentrations, and yields. For key properties, there are diminishing returns with the addition of increasing amounts of fuel. These results should provide useful guidelines for the design of functional systems making use of this chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Spatiotemporal Energy Cascade in Three-Dimensional Magnetohydrodynamic Turbulence

We present a new scale decomposition method to investigate turbulence in wavenumber-frequency space. Using 3D magnetohydrodynamic turbulence simulations, we show that magnetic fluctuations with time scales longer than the nonlinear time exhibit an inverse cascade toward even smaller frequencies. Low frequency magnetic fluctuations support turbulence, acting as an energy reservoir that is converted into plasma kinetic energy, the latter cascading toward large wavenumbers and frequencies, where it is dissipated. Our results shed new light on the spatiotemporal properties of turbulence, potentially explaining the origin and role of low frequency turbulent fluctuations in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗