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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Node-Accessible Zirconium MOFs

High-stability, zirconium-based metal–organic frameworks are attractive as heterogeneous catalysts and as model supports for uniform arrays of subsequently constructed heterogeneous catalysts—for example, MOF-node-grafted metal–oxy and metal–sulfur clusters. For hexa-Zr(IV)-MOFs characterized by nodes that are less than 12-connected, sites not used for linkers are ideally occupied by reactive and displaceable OH/H 2 O pairs. Here, the desired pairs are ideal for grafting the aforementioned catalytic clusters, while aqua-ligand lability renders them effective for exposing highly Lewis-acidic Zr(IV) sites (catalytic sites) to candidate reactants. New single-crystal X-ray studies of an eight-connected Zr-MOF, NU-1000, reveal that conventional activation fully removes modulator ligands, but replaces them with three node-blocking formate ligands (from solvent decomposition) and only one OH/H 2 O pair, not four—a largely overlooked complication that now appears to be general for Zr-MOFs. Here we describe an alternative activation protocol that effectively removes modulators, avoids formate, and installs the full complement of terminal OH/H 2 O pairs. It does so via an unusual isolatable intermediate featuring eight aqua ligands and four non-ligated chlorides—again as supported by single-crystal X-ray data. We find that complete replacement of node-blocking modulators/formate with the originally envisioned OH/OH 2 pairs has striking consequences; here we touch upon just three. First, elimination of unrecognized formate renders aqua ligands much more thermally labile, enabling open Zr(IV) sites to be obtained at lower temperature. Second, in the absence of formate, which otherwise links and locks pairs of node Zr(IV) ions, reversible removal of aqua ligands engenders reversible contraction of MOF meso- and micropores, as evidenced by X-ray diffraction. Third, formate replacement with OH/OH 2 pairs renders NU-1000 ca.10× more active for catalytic hydrolytic degradation of a representative simulant of G-type chemical warfare agents.

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

Tully-Fisher Distances and Dynamical Mass Constraints for 24 Host Galaxies of Reverberation-Mapped AGNs

We present Tully–Fisher distances for 24 active galactic nucleus (AGN) host galaxies with black hole mass (M(sub BH)) measurements from reverberation mapping, as well as the first calibration of the V-band Tully–Fisher relation. Combining our measurements of H I 21 cm emission with Hubble Space Telescope and ground-based optical and near-infrared images allows multiple distance measurements for 19 galaxies and single measurements for the remaining 5. Separation of the nucleus from its host galaxy via surface brightness decomposition yields galaxy-only luminosities, thus allowing measurements of the distance moduli free of contamination from the AGNs. For 14 AGN hosts, these are the first reported distances independent of redshift, and hence independent of peculiar velocities. For the remaining galaxies, we show good agreement between our distances and those previously reported from surface brightness fluctuations and Cepheids. We also determine the total galaxy mass enclosed within the estimated H I radius, which when compared to the baryonic content allows for constraints on the dark matter masses. We find a typical mass fraction of M(sub DM)/M(sub DYN) = 62%, and find significant correlations between M(sub BH)–M(sub DYN) and M(sub BH)–M(sub DM). Finally, we scale our galaxy radii based on estimated relationships between visible and halo radii and assume a flat rotation curve out to the halo radius to approximate M(sub HALO). Over the range of M(sub BH) and M(sub HALO) in this sample, we find good agreement with observationally constrained relationships between M(sub BH) and M(sub HALO) and with hydrodynamical simulations.

AGN host galaxies↗
NASA NTRS1997

Reduced Navier Stokes Relaxation Procedures for Internal Flows

In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.

Rubin, Stanley G.↗
NASA NTRS1997

A Multi-Level Parallelization Concept for High-Fidelity Multi-Block Solvers

The integration of high-fidelity Computational Fluid Dynamics (CFD) analysis tools with the industrial design process benefits greatly from the robust implementations that are transportable across a wide range of computer architectures. In the present work, a hybrid domain-decomposition and parallelization concept was developed and implemented into the widely-used NASA multi-block Computational Fluid Dynamics (CFD) packages implemented in ENSAERO and OVERFLOW. The new parallel solver concept, PENS (Parallel Euler Navier-Stokes Solver), employs both fine and coarse granularity in data partitioning as well as data coalescing to obtain the desired load-balance characteristics on the available computer platforms. This multi-level parallelism implementation itself introduces no changes to the numerical results, hence the original fidelity of the packages are identically preserved. The present implementation uses the Message Passing Interface (MPI) library for interprocessor message passing and memory accessing. By choosing an appropriate combination of the available partitioning and coalescing capabilities only during the execution stage, the PENS solver becomes adaptable to different computer architectures from shared-memory to distributed-memory platforms with varying degrees of parallelism. The PENS implementation on the IBM SP2 distributed memory environment at the NASA Ames Research Center obtains 85 percent scalable parallel performance using fine-grain partitioning of single-block CFD domains using up to 128 wide computational nodes. Multi-block CFD simulations of complete aircraft simulations achieve 75 percent perfect load-balanced executions using data coalescing and the two levels of parallelism. SGI PowerChallenge, SGI Origin 2000, and a cluster of workstations are the other platforms where the robustness of the implementation is tested. The performance behavior on the other computer platforms with a variety of realistic problems will be included as this on-going study progresses.

Hatay, Ferhat F.↗
DOE OSTI2024

Model associated with: "Thermodynamic control on the decomposition of organic matter across different electron acceptors"

This model data package is associated with the publication “Thermodynamic control on the decomposition of organic matter across different electron acceptors” submitted to Soil Biology and Biochemistry (Zheng et al., 2023; https://doi.org/10.1016/j.soilbio.2024.109364).In this research, a thermodynamic modeling framework is built to flexibly incorporate both organic matter (OM) molecules and electron acceptors for estimating potential free energy release from various redox reactions and to further predict reaction rates based on Microbial Transition State Theory. The model package includes scripts for thermodynamic modeling and postprocessing. Input Fourier-transform ion cyclotron resonance (FTICR) data are from a previous experimental study (Boye et al., 2018), and model outputs are free energy predictions and stoichiometric coefficients associated with all possible redox reactions.This data package is associated with the project GitHub repository found at MM_bioenergetic_modeling.This data package contains four folders (Input_FTICR, Model, Output, and Output_processing), a file-level metadata (FLMD) csv, and a data dictionary (dd) csv. Please see Zheng_bioenergetic_modeling_flmd.csv for a list of all files contained in this data package and descriptions for each. The Zheng_bioenergetic_modeling_dd.csv file describes the csv column headers. The “Model” folder contains scripts to run energy balance calculations for each electron acceptor. The “Output” folder contains csv files with stoichiometric information from model simulations. And the "Output_processing" folder contains scripts for reaction rate calculations and to generate plots.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multi-modelling predictions show high uncertainty of required carbon input changes to reach a 4‰ target

Soils store vast amounts of carbon (C) on land, and increasing soil organic carbon (SOC) stocks in already managed soils such as croplands may be one way to remove C from the atmosphere, thereby limiting subsequent warming. The main objective of this study was to estimate the amount of additional C input needed to annually increase SOC stocks by 4‰ at 16 long-term agricultural experiments in Europe, including exogenous organic matter (EOM) additions. We used an ensemble of six SOC models and ran them under two configurations: (1) with default parametrization and (2) with parameters calibrated site-by-site to fit the evolution of SOC stocks in the control treatments (without EOM). We compared model simulations and analysed the factors generating variability across models. The calibrated ensemble was able to reproduce the SOC stock evolution in the unfertilised control treatments. We found that, on average, the experimental sites needed an additional 1.5 ± 1.2 Mg C ha -1 year -1 to increase SOC stocks by 4‰ per year over 30 years, compared to the C input in the control treatments (multi-model median ± median standard deviation across sites). That is, a 119% increase compared to the control. While mean annual temperature, initial SOC stocks and initial C input had a significant effect on the variability of the predicted C input in the default configuration (i.e., the relative standard deviation of the predicted C input from the mean), only water-related variables (i.e., mean annual precipitation and potential evapotranspiration) explained the divergence between models when calibrated. Here, our work highlights the challenge of increasing SOC stocks in agriculture and accentuates the need to increasingly lean on multi-model ensembles when predicting SOC stock trends and related processes. To increase the reliability of SOC models under future climate change, we suggest model developers to better constrain the effect of water-related variables on SOC decomposition.

4 per 1000 initiative↗
DOE OSTI2025

Beyond real: alternative unitary cluster Jastrow models for molecular electronic structure calculations on near-term quantum computers

Near-term quantum devices require wavefunction ansätze that are expressive while also of shallow circuit depth in order to both accurately and efficiently simulate molecular electronic structure. While the unitary coupled cluster ansatz (e.g., UCCSD) has become a standard, the high gate count associated with the implementation of this limits its feasibility on noisy intermediate-scale quantum (NISQ) hardware. k -Fold unitary cluster Jastrow (uCJ) ansätze mitigate this challenge by providing O( kN 2 ) circuit scaling and favorable linear depth circuit implementation. Previous work has focused on the real orbitalrotation (Re-uCJ) variant of uCJ, which allows an exact (Trotter-free) implementation. Here we extend and generalize the k -fold uCJ framework by introducing two new variants, Im-uCJ and g-uCJ, which incorporate imaginary and fully complex orbital rotation operators, respectively. Similar to Re-uCJ, both of the new variants achieve quadratic gate-count scaling. Our results focus on the simplest k = 1 model, and show that the uCJ models frequently maintain energy errors within chemical accuracy (∼1 kcal mol −1 ). Both g-uCJ and Im-uCJ are more expressive in terms of capturing electron correlation and are also more accurate than the earlier Re-uCJ ansatz. We further show that Im-uCJ and g-uCJ circuits can also be implemented exactly, without any Trotter decomposition. Numerical tests using k = 1 on H 2 , H 3 + , Be 2 , C 2 H 4 , C 2 H 6 and C 6 H 6 in various basis sets confirm the practical feasibility of these shallow Jastrow-based ansätze for applications on near-term quantum hardware.

Tkachenko, Nikolay V. [University of California, B↗
DOE OSTI2025

Risk Assessment in a Chemical Laboratory Following an Explosive Incident Involving a Novel Diazonium Compound: Retrospective Analysis and Lessons Learned

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. Furthermore, the steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Computer simulations↗
DOE OSTI2020

Attack-resilient algorithms and testbed federation for wide-area protection and control in smart grid

Today’s energy infrastructure is undergoing a massive transformation across all generation, transmission, and distribution systems to provide reliability, efficiency, and sustainability to the power system network. During the past 10 years, there has been a remarkable enhancement in developing close-loop wide-area protection and control (WAPAC) applications, such as automatic generation control (AGC), wide-area protection scheme (WAPS), synchrophasor-based wide-area voltage control system (WAVCS), etc., in the energy management system (EMS) to provide an adequate situational awareness to control center operators during physical disturbances and natural events. The rapid digitalization, unencrypted communication, and numerous data-sharing devices have increased digital access points that have exposed the grid network to countless cyberattacks pertinent to the WAPAC cybersecurity. Given the amount of complex grid network coupled with legacy infrastructures, it is untrustworthy to completely rely on information technology (IT)-based conventional security measures to defend against evolving cyberattacks. Hence, there is a compelling need to go beyond the traditional paradigm of ‘security by obscurity’ and ‘bolt-on’ security measures, and develop a suite of innovative security solutions at multiple layers to transform the current “fault-resilient” grid into “fault and attack-resilient” grid of the future. This dissertation specifically focuses on developing application-specific attack-resilient solutions for critical WAPAC applications in the smart grid. The first part of the study, Attack-Resilient Wide-Area Protection Scheme, presents a cyber-physical attack-resilient system (CPARS) for WAPS, where two anomaly detection methods are proposed to detect multi-level cyberattacks followed by the rules-based intrusion mitigation system (RIMS) to initiate effective mitigation strategies, tailored to these attacks, to restore the normal grid operation after disturbances. The first anomaly detection component, Machine learning-based Anomaly Detection for Centralized WAPS, discusses a novel architecture and methodology for developing a machine-learning-based anomaly detection system (ADS) for the centralized wide-area protection system (CWAPS) by applying the variational mode decomposition (VMD) technique and decision tree (DT) algorithms using cyber logs and synchrophasor measurements. The second component, Multi-Agent-based Anomaly Detection for Decentralized WAPS, presents a two-level hierarchical multi-agent-based decentralized WAPS against the system-aware stealthy cyber-attacks. Finally, for the RIMS, the state transition diagram-based mitigation measures are defined based on the grid operation states and experimental testing and evaluation are presented in a cyber-physical testbed environment. The second part of the study, Attack-Resilient Wide-Area Voltage Control System, addresses the existing vulnerabilities in the wide-area voltage control system (WAVCS), presents the methodology for performing an impact analysis during data integrity attacks, and proposes a data-driven attackresilient system (DARS) by incorporating machine learning-based ADS and rules-based anomaly mitigation system (RAMS) to detect these attacks and provide necessary mitigation actions to maintain the transient voltage stability after disturbances. Finally, the third part, Networked Federation Testbed (NEFTSec) for Smart Grid Cybersecurity, focuses on developing an interconnected federated testbed through a common network like the internet, describes its conceptual architectures, and also presents its design components for experimental testing, validation, and evaluation. Further, it presents the co-simulation interface algorithm (CIA) to facilitate the geographically-dispersed real-time distributed simulation (GD-RTDS) during the cyber-physical federation. In general, this dissertation presents innovative attack-resilient solutions for WAPAC applications that are experimentally tested and evaluated using the PowerCyber testbed and NEFTSec platform to accelerate the transition of state-of-the-art research works to the real-field deployments.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Elliptically-Contoured Tensor-variate Distributions with Application to Image Learning

Statistical analysis of tensor-valued data has largely used the tensor-variate normal (TVN) distribution that may be inadequate for data arising from distributions with heavier or lighter tails. We study a general family of elliptically contoured (EC) TV distributions and derive its characterizations, moments, marginal, and conditional distributions. We describe procedures for maximum likelihood estimation from data that are (1) uncorrelated draws from an EC distribution, (2) from a scale mixture of the TVN distribution, and (3) from an underlying but unknown EC distribution, for which we extend Tyler’s robust estimator. A detailed simulation study highlights the benefits of choosing an EC distribution over the TVN for heavier-tailed data. We develop TV classification rules using discriminant analysis and EC errors and show that they better predict cats and dogs from images in the Animal Faces-HQ dataset than the TVN-based rules. A novel tensor-on-tensor regression and TV analysis of variance (TANOVA) framework under EC errors is also demonstrated to better characterize gender, age, and ethnic origin than the usual TVN-based TANOVA in the celebrated labeled faces of the wild dataset.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials (Second Annual Progress Report)

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of nuclear reactor parts. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which can introduce internal material flaws, such as pores and anisotropy. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to exposure to high temperature, radiation and corrosive environment in a nuclear reactor. Thermal tomography (TT) provides a capability for non-destructive evaluation of sub-surface defects in arbitrary size structures. Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer, and through reconstruction of data obtained from experimental measurements. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and Inconel 718 plate produced with laser powder bed fusion (LPBF) method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBH’s have the same sizes as the imprinted defects in the LPBF specimens, but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects, and allow for estimation of their sizes and locations. Detection sensitivity of TT is limited by noises. We investigate separation of signal from noise in thermography images using several machine learning (ML) methods, including new spatio-temporal blind source separation (STBSS) and spatio-temporal sparse dictionary learning (STSDL) methods. Performance of the ML methods is benchmarked using thermography data obtained from imaging stainless steel 316L and Inconel 718 specimens produced LPBF method with imprinted calibrated porosity defects. The ML methods are ranked by F-score and execution runtime. Finally, we investigate TT of AM stainless steel 316L specimen with imprinted internal porosity defects using relatively low-cost, small form factor infrared (IR) camera based on uncooled micro bolometer detector. Sparse coding related K-means singular value decomposition (SVD) machine learning, image processing algorithms are developed to improve quality of TT images through removal of Additive white Gaussian noise without blurring the images. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact TT can also be used for in-service nondestructive evaluation (NDE). With capability to perform in-service NDE of the AM component lifecycle, TT data can be used for development of a component digital twin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2024

Sound propagation in realistic interactive 3D scenes with parameterized sources using deep neural operators

We address the challenge of acoustic simulations in three-dimensional (3D) virtual rooms with parametric source positions, which have applications in virtual/augmented reality, game audio, and spatial computing. The wave equation can fully describe wave phenomena such as diffraction and interference. However, conventional numerical discretization methods are computationally expensive when simulating hundreds of source and receiver positions, making simulations with parametric source positions impractical. To overcome this limitation, we propose using deep operator networks to approximate linear wave-equation operators. This enables the rapid prediction of sound propagation in realistic 3D acoustic scenes with parametric source positions, achieving millisecond-scale computations. By learning a compact surrogate model, we avoid the offline calculation and storage of impulse responses for all relevant source/listener pairs. Our experiments, including various complex scene geometries, show good agreement with reference solutions, with root mean squared errors ranging from 0.02 to 0.10 Pa. Notably, our method signifies a paradigm shift as—to our knowledge—no prior machine learning approach has achieved precise predictions of complete wave fields within realistic domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Lepidocrocite to Maghemite to Hematite: A way to have Magnetic and Hematitic Martian Soil

We examined decomposition products of lepidocrocite, which were produced by heating the phase in air at temperatures up to 525 C for 3 and 300 hr, by XRD, TEM, magnetic methods, and reflectance spectroscopy (visible and near-IR). Single-crystal lepidocrocite particles dehydroxilated to polycrystalline particles of disordered maghemite which subsequently transformed to polycrystalline particles of hematite. Essentially pure maghemite was obtained at 265 and 223 C for the 3 and 300 hr heating experiments, respectively. Its saturation magnetization (J(sub s)) and mass specific susceptibility are approximately 50 A(sq m)/kg and approximately 40 cubic micrometers/kg, respectively. Because hematite is spectrally dominant, spectrally-hematitic samples (i.e., characterized bv a minimum near 860 nm and a maximum near 750 nm) could also be strongly magnetic (J(sub s) up to approximately 30 A(sq m)/kg) from the masked maghemite component. TEM analyses showed that individual particles are polycrystalline with respect to both maghemite and hematite. The spectrally-hematitic and magnetic Mh+Hm particles can satisfy the spectral and magnetic constraints for Martian surface materials over a wide range of values of Mh/(Mh+Hm) and as either pure oxide powders or (within limits) as components of multiphase particles. These experiments are consistent with lepidocrocite as the precursor of Mh+Hm assemblages on Mars, but other phases (e.g., magnetite) that decompose to Mh and Hm are also possible precursors. Simulations done with a copy of the Mars Pathfinder Magnet Array showed that spectrally hematitic Mh+Hm powders having J(sub s) equal to 20.6 A(sq m)/kg adhered to all five magnets.

Morris, Richard V.↗
NASA NTRS1998

Parallel 3D Multi-Stage Simulation of a Turbofan Engine

A 3D multistage simulation of each component of a modern GE Turbofan engine has been made. An axisymmetric view of this engine is presented in the document. This includes a fan, booster rig, high pressure compressor rig, high pressure turbine rig and a low pressure turbine rig. In the near future, all components will be run in a single calculation for a solution of 49 blade rows. The simulation exploits the use of parallel computations by using two levels of parallelism. Each blade row is run in parallel and each blade row grid is decomposed into several domains and run in parallel. 20 processors are used for the 4 blade row analysis. The average passage approach developed by John Adamczyk at NASA Lewis Research Center has been further developed and parallelized. This is APNASA Version A. It is a Navier-Stokes solver using a 4-stage explicit Runge-Kutta time marching scheme with variable time steps and residual smoothing for convergence acceleration. It has an implicit K-E turbulence model which uses an ADI solver to factor the matrix. Between 50 and 100 explicit time steps are solved before a blade row body force is calculated and exchanged with the other blade rows. This outer iteration has been coined a "flip." Efforts have been made to make the solver linearly scaleable with the number of blade rows. Enough flips are run (between 50 and 200) so the solution in the entire machine is not changing. The K-E equations are generally solved every other explicit time step. One of the key requirements in the development of the parallel code was to make the parallel solution exactly (bit for bit) match the serial solution. This has helped isolate many small parallel bugs and guarantee the parallelization was done correctly. The domain decomposition is done only in the axial direction since the number of points axially is much larger than the other two directions. This code uses MPI for message passing. The parallel speed up of the solver portion (no 1/0 or body force calculation) for a grid which has 227 points axially.

Turner, Mark G.↗
DOE OSTI2022

Decomposition of high-rank factorized unitary coupled-cluster operators using ancilla and multiqubit controlled low-rank counterparts

The factorized form of the unitary coupled-cluster (UCC) approximation is one of the most promising methodologies to prepare trial states for strongly correlated systems within the variational quantum eigensolver framework. The factorized form of the UCC Ansatz can be systematically applied to a reference state to generate the desired entanglement. The difficulty associated with such an approach is the requirement of simultaneously entangling a growing number of qubits, which quickly exceeds the hardware limitations of today's quantum machines. In particular, while circuits for singles and double excitations can be performed on current hardware, higher-rank excitations require too many gate operations. In this work, we propose a set of schemes that trade off using extra qubits for a reduced gate depth to decompose high-rank UCC excitation operators into significantly lower depth circuits. These results will remain useful even when fault-tolerant machines are available to reduce the overall state-preparation circuit depth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Mechanical Testing of PMCs under Simulated Rapid Heat-Up Propulsion Environments: In-Plane Compressive Behavior - II

Carbon fiber thermoset polymer matrix composites (PMC) with high temperature polyimide based in-situ polymerized monomer reactant (PMR) resin has been used for some time in applications which can see temperatures up to 550 F. Currently, graphite fiber PMR based composites are used in several aircraft engine components including the outer bypass duct for the GE F-404, exit flaps for the P&W F-100-229, and the core cowl for the GE/Snecma CF6-80A3. Newer formulations, including PMR-II-50 are being investigated as potential weight reduction replacements of various metallic components in next generation high performance propulsion rocket engines that can see temperatures which exceed 550 F. Extensive FEM thermal modeling indicates that these components are exposed to rapid heat-up rates (up to -200 F/sec) and to a maximum temperature of around 600 F. Even though the predicted maximum part temperatures were within the capability of PW-II-50, the rapid heat-up causes significant through-thickness thermal gradients in the composite part and even more unstable states when combined with moisture. Designing composite parts for such extreme service environments will require accurate measurement of intrinsic and transient mechanical properties and the hygrothermal performance of these materials under more realistic use conditions. The mechanical properties of polymers degrade when exposed to elevated temperatures even in the absence of gaseous oxygen. Accurate mechanical characterization of the material is necessary in order to reduce system weight while providing sufficient factors of safety. Historically, the testing of PMCs at elevated temperatures has been plagued by the antagonism between two factors. First, moisture has been shown to profoundly affect the mechanical response of these materials at temperatures above their glass transition temperature while concurrently lowering the material's Tg. Moisture phenomena is due to one or a combination of three effects, i.e., plastization of polymeric material by water, the internal pressure generated by the volatilization of water at elevated temperatures, and hydrolytic chemical decomposition. However, moisture is lost from the material at increasing rates as temperature increases. Second, because PMCs are good thermal insulators, when they are externally heated at even mild rates large thermal gradients can develop within the material. At temperatures where a material property changes rapidly with temperature the presence of a large thermal gradient is unacceptable for intrinsic property characterization purposes. Therefore, long hold times are required to establish isothermal conditions. However, in the service environments high-heating-rates, high temperatures, high-loading rates are simultaneous present along with residual moisture. In order to capture the effects of moisture on the material, holding at- temperature until isothermal conditions are reached is unacceptable particularly in materials with small physical dimensions. Thus, the effects due to moisture on the composite's mechanical characteristics, ie., their so-called analog response, may be instructive. One approach employed in this program was rapid heat-up (approx. 200 F/sec.) and loading of both dry and wet in-plane compressive specimens to examine the effects of moisture on this resin dominated mechanical property of the material.

Stokes, Eric H.↗
DOE OSTI2020

Learning in Modal Space: Solving Time-Dependent Stochastic PDEs Using Physics-Informed Neural Networks

One of the open problems in scientific computing is the long-time integration of nonlinear stochastic partial differential equations (SPDEs), especially with arbitrary initial data. We address this problem by taking advantage of recent advances in scientific machine learning and the spectral dynamically orthogonal (DO) and borthogonal (BO) methods for representing stochastic processes. The recently introduced DO/BO methods reduce the SPDE to solving a system of deterministic PDEs and a system of stochastic ordinary differential equations. Specifically, we propose two new physics-informed neural networks (PINNs) for solving time-dependent SPDEs, namely the neural network (NN)-DO/BO methods. The proposed methods incorporate the DO/BO constraints into the loss function (along with the modal decomposition of the SPDE) with an implicit form instead of generating explicit expressions for the temporal derivatives of the DO/BO modes. Hence, the NN-DO/BO methods can overcome some of the drawbacks of the original DO/BO methods. For example, we do not need the assumption that the covariance matrix of the random coefficients is invertible as in the original DO method, and we can remove the assumption of no eigenvalue crossing as in the original BO method. Moreover, the NN-DO/BO methods can be used to solve time-dependent stochastic inverse problems with the same formulation and same computational complexity as for forward problems. Furthermore, we demonstrate the capability of the proposed methods via several numerical examples, namely: (1) A linear stochastic advection equation with deterministic initial condition: we obtain good results with the proposed methods, while the original DO/BO methods cannot be applied directly in this case. (2) Long-time integration of the stochastic Burgers' equation: we show the good performance of NN-DO/BO methods, especially the effectiveness of the NN-BO approach for such problems with many eigenvalue crossings during the whole time evolution, while the original BO method fails. (3) Nonlinear reaction diffusion equation: we consider both the forward problem and the inverse problems, including very noisy initial point values, to investigate the flexibility of the NN-DO/BO methods in handling inverse and mixed type problems. Taken together, these simulation results demonstrate that the NN-DO/BO methods can be employed to effectively quantify uncertainty propagation in a wide range of physical problems, but future work should address the efficiency issue of PINNs for forward problems.

97 MATHEMATICS AND COMPUTING↗