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 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A Local Macroscopic Conservative (LoMaC) Low Rank Tensor Method for the Vlasov Dynamics

Abstract In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum and energy at the discrete level. This is a follow-up work of our previous development of a conservative low rank tensor approach for Vlasov dynamics ( arXiv:2201.10397 ). In that work, we applied a low rank tensor method with a conservative singular value decomposition to the high dimensional VP system to mitigate the curse of dimensionality, while maintaining the local conservation of mass and momentum. However, energy conservation is not guaranteed, which is a critical property to avoid unphysical plasma self-heating or cooling. The new ingredient in the LoMaC low rank tensor algorithm is that we simultaneously evolve the macroscopic conservation laws of mass, momentum and energy using a flux-difference form with kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables by a conservative orthogonal projection. The algorithm is extended to the high dimensional problems by hierarchical Tuck decomposition of solution tensors and a corresponding conservative projection algorithm. Extensive numerical tests on the VP system are showcased for the algorithm’s efficacy.

Guo, Wei↗

DOE OSTI2022

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A Novel Partitioned Approach for Reduced Order Model—Finite Element Model (ROM-FEM) and ROM-ROM Coupling

Partitioned methods allow one to build a simulation capability for coupled problems by reusing existing single-component codes. In so doing, partitioned methods can shorten code development and validation times for multiphysics and multiscale applications. In this work, we consider a scenario in which one or more of the “codes” being coupled are projection-based reduced order models (ROMs), introduced to lower the computational cost associated with a particular component. We simulate this scenario by considering a model interface problem that is discretized independently on two non-overlapping subdomains. Here we then formulate a partitioned scheme for this problem that allows the coupling between a ROM “code” for one of the subdomains with a finite element model (FEM) or ROM “code” for the other subdomain. The ROM “codes” are constructed by performing proper orthogonal decomposition (POD) on a snapshot ensemble to obtain a low-dimensional reduced order basis, followed by a Galerkin projection onto this basis. The ROM and/or FEM “codes” on each subdomain are then coupled using a Lagrange multiplier representing the interface flux. To partition the resulting monolithic problem, we first eliminate the flux through a dual Schur complement. Application of an explicit time integration scheme to the transformed monolithic problem decouples the subdomain equations, allowing their independent solution for the next time step. We show numerical results that demonstrate the proposed method’s efficacy in achieving both ROM-FEM and ROM-ROM coupling.

de Castro, Amy↗

DOE OSTI2022

The equation of state and shock-driven decomposition of polymethylmethacrylate (PMMA)

We present new equations of state (EOS) for polymethylmethacrylate and its shock-driven decomposition products, generated in both tabular form and as linear U S – u p fits to all previously available shock data. Different approaches are taken for materials shocked above and below the apparent threshold for chemical reactions at pressures P ≈ 30 GPa. Both EOSs are calibrated to and compared with a broad array of experimental results, as well as new ab initio molecular dynamics simulations presented here. Comparisons are good in most cases, and we focus our particular attention on those observables that test the distinction between reactants and products.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Rigorous justification for the space–split sensitivity algorithm to compute linear response in Anosov systems

Ruelle gave a formula for linear response of transitive Anosov diffeomorphisms. Recently, practically computable realizations of Ruelle's formula have emerged that potentially enable sensitivity analysis of certain high-dimensional chaotic numerical simulations encountered in the applied sciences. In this paper, we provide full mathematical justification for the convergence of one such efficient computation, the space–split sensitivity, or S3, algorithm. In S3, Ruelle's formula is computed as a sum of two terms obtained by decomposing the perturbation vector field into a coboundary and a remainder that is parallel to the unstable direction. Such a decomposition results in a splitting of Ruelle's formula that is amenable to efficient computation. Here, we prove the existence of the S3 decomposition and the convergence of the computations of both resulting components of Ruelle's formula.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Fungal community structure and function shifts with atmospheric nitrogen deposition

Fungal decomposition of soil organic matter depends on soil nitrogen (N) availability. This ecosystem process is being jeopardized by changes in N inputs that have resulted from a tripling of atmospheric N deposition in the last century. Soil fungi are impacted by atmospheric N deposition due to higher N availability, as soils are acidified, or as micronutrients become increasingly limiting. Fungal communities that persist with chronic N deposition may be enriched with traits that enable them to tolerate environmental stress, which may trade-off with traits enabling organic matter decomposition. We hypothesized that fungal communities would respond to N deposition by shifting community composition and functional gene abundances toward those that tolerate stress but are weak decomposers. We sampled soils at seven eastern US hardwood forests where ambient N deposition varied from 3.2 to 12.6 kg N ha -1 year -1 , five of which also have experimental plots where atmospheric N deposition was simulated through fertilizer application treatments (25–50 kg N ha -1 year -1 ). Fungal community and functional responses to fertilizer varied across the ambient N deposition gradient. Fungal biomass and richness increased with simulated N deposition at sites with low ambient deposition and decreased at sites with high ambient deposition. Fungal functional genes involved in hydrolysis of organic matter increased with ambient N deposition while genes involved in oxidation of organic matter decreased. One of four genes involved in generalized abiotic stress tolerance increased with ambient N deposition. In summary, we found that the divergent response to simulated N deposition depended on ambient N deposition levels. Fungal biomass, richness, and oxidative enzyme potential were reduced by N deposition where ambient N deposition was high suggesting fungal communities were pushed beyond an environmental stress threshold. Fungal community structure and function responses to N enrichment depended on ambient N deposition at a regional scale.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2012

Disentangling AGN and Star Formation in Soft X-Rays

We have explored the interplay of star formation and active galactic nucleus (AGN) activity in soft X-rays (0.5-2 keV) in two samples of Seyfert 2 galaxies (Sy2s). Using a combination of low-resolution CCD spectra from Chandra and XMM-Newton, we modeled the soft emission of 34 Sy2s using power-law and thermal models. For the 11 sources with high signal-to-noise Chandra imaging of the diffuse host galaxy emission, we estimate the luminosity due to star formation by removing the AGN, fitting the residual emission. The AGN and star formation contributions to the soft X-ray luminosity (i.e., L(sub x,AGN) and L(sub x,SF)) for the remaining 24 Sy2s were estimated from the power-law and thermal luminosities derived from spectral fitting. These luminosities were scaled based on a template derived from XSINGS analysis of normal star-forming galaxies. To account for errors in the luminosities derived from spectral fitting and the spread in the scaling factor, we estimated L(sub x,AGN) and L(sub x,SF))from Monte Carlo simulations. These simulated luminosities agree with L(sub x,AGN) and L(sub x,SF) derived from Chandra imaging analysis within a 3sigma confidence level. Using the infrared [Ne ii]12.8 micron and [O iv]26 micron lines as a proxy of star formation and AGN activity, respectively, we independently disentangle the contributions of these two processes to the total soft X-ray emission. This decomposition generally agrees with L(sub x,SF) and L(sub x,AGN) at the 3 sigma level. In the absence of resolvable nuclear emission, our decomposition method provides a reasonable estimate of emission due to star formation in galaxies hosting type 2 AGNs.

Disentangling↗

NASA NTRSDate not supplied

Rolling Horizon with K-Position Search Method for Strategic Deconfliction of Package Delivery UAS

In this research, the strategic deconfliction of unmanned aircraft systems for an urban package delivery environment with two depots and multiple drop-off locations is studied. This research aims to formulate a mathematical model to compute both the departure sequence and scheduled time of departure for each unmanned aircraft system at a depot, considering temporal constraints at en-route crossing waypoints and depots for strategic deconfliction. However, the problem formulation results in an NP-hard mixed-integer nonlinear programming problem for the global optimal solution, so instead, a "rolling horizon with𝑘-position search"heuristic method is developed. The simulation studies show that an increase in the value of𝑘(the parameter used to determine the size of the local neighborhood) reduces the average ground delay at the cost of an increase in the computation time for a given problem size. The study also shows an order of magnitude increase in the maximum number of flights scheduled with the integration of rolling horizon (time decomposition) compared to those without the integration of rolling horizon in the heuristic algorithm for a given computation time cut off.

UTM↗

NASA NTRS2012

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗

DOE OSTI2022

Nanoscale Mapping and Defect-Assisted Manipulation of Surface Plasmon Resonances in 2D Bi 2 Te 3 /Sb 2 Te 3 In-Plane Heterostructures

Here, the Bi 2 Te 3 /Sb 2 Te 3 in-plane heterostructure is reported as a low-dimensional tunable chalcogenide well suited as plasmonic building block for the visible-UV spectral range. Electron-driven plasmon excitations of low-dimensional Bi 2 Te 3 /Sb 2 Te 3 are investigated by monochromated electron energy loss spectroscopy spectrum imaging. To resolve the nanoscale spatial distribution of various local plasmonic resonances, singular value decomposition is used to disentangle the spectral data and identify the individual spectral contributions of various corner, edge, and face modes. Furthermore, defect-plasmon interactions are investigated both for nanoscale intrinsic and thermally induced extrinsic polygonal defects (in situ sublimation). Signature of defect-induced red shift ranging from a several hundreds of millielectronvolts to a few electronvolts, broadening of various plasmon response, together with selective enhancement and significant variations in their intensity are detected. This study highlights the presence of a heterointerface and identifies defects as physical tuning pathways to modulate the plasmonic response over a broad spectral range. Finally, the experimental observations are compared qualitatively and validated with numerical simulations using the electron-driven discrete dipole approximation. Low-dimensional Bi 2 Te 3 /Sb 2 Te 3 as a less explored plasmonic system holds great promises as emerging platform for integrated plasmonics. Furthermore, introducing controlled structural defects can open the door for nanoengineering of plasmonic properties in such systems.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Tailoring Chemical Absorption-Precipitation to Lower the Regeneration Energy of a CO 2 Capture Solvent

Solvent-based CO 2 capture consumes significant amounts of energy for solvent regeneration. To improve energy efficiency, this study investigates CO 2 fixation in a solid form through solvation, followed by ionic self-assembly-aided precipitation. Based on the hypothesis that CO 3 2- ions may bind with monovalent metal ions, we introduced Na + into an aqueous hexane-1,6-diamine solution where CO 2 forms carbamate and bicarbonate. Then, Na + ions in the solvent act as a seed for ionic self-assembly with diamine carbamate to form an intermediate ionic complex. The recurring chemical reactions lead to the formation of an ionic solid from a mixture of organic carbamate/carbonate and inorganic sodium bicarbonate (NaHCO 3 ), which can be easily removed from the aqueous solvent through sedimentation or centrifugation and heated to release the captured CO 2 . Mild-temperature heating of the solids at 80–150 °C causes decomposition of the solid CO 2 -diamine-Na molecular aggregates and discharge of CO 2 . This sorbent regeneration process requires 6.5–8.6 GJ/t CO 2 . It was also found that the organic carbamate/carbonate solid, without NaHCO 3 , contains a significant amount of CO 2 , up to 6.2 mmol CO 2 /g-sorbent, requiring as low as 2.9–5.8 GJ/t CO 2 . In conclusion, molecular dynamic simulations support the hypothesis of using Na + to form relatively less stable, yet sufficiently solid, complexes for the least energy-intensive recovery of diamine solvents compared to bivalent carbonate–forming ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effects of SO 2 poisoning and regeneration on spinel containing CH 4 oxidation catalysts

Methane oxidation under periodic conditions and the oxygen storage capacity of a bilayer Pt/Pd/Al 2 O 3 over a Mn 0.5 Fe 2.5 O 4 spinel catalyst were studied before and after SO 2 exposure, and after simulated regeneration conditions. Prior to sulfur exposure, improvement in CH 4 oxidation conversion under periodic conditions compared to steady-state conditions was observed. After sulfur exposure at 100 °C, there was a loss in CH 4 oxidation performance and a loss of oxygen storage capacity of the spinel material. Here, the extent of regeneration from sulfur poisoning depends on the ability to induce the decomposition of sulfate species, and while all regeneration methods tested in this study did improve CH 4 conversion, regeneration methods under periodic conditions induced greater sulfur species desorption from the catalyst surface leading to improved CH 4 con version. Key regeneration parameters – temperature, feed composition, modulation amplitude and frequency – were optimized to induce S species decomposition and correlated to CH 4 oxidation activity recovery.

03 NATURAL GAS↗

DOE OSTI2021

Carbon fiber synthesis from pitch: Insights from ReaxFF based molecular dynamics simulations

Despite the potential to lowering fabrication costs of pitch-based carbon fiber (CF) compared to polyacrylonitrile-based CF, large scale utilization is currently limited by the complexity in manufacturing. This is mainly due to the lack of consistency between the pitch structure/chemistry and CF properties under different processing conditions. Here, we employed full atomistic simulations based on reactive force fields to investigate the conversion of pitch precursors to CFs. Additionally, it was found that at the carbonization stage, the consumption rates of precursors, gas formations, and solid yields (carbonized CFs) are all dependent on the composition and arrangement of molecular moieties in the precursor. Specifically, precursors with sp 3 carbons in the cata-condensed core area have limited decomposition pathways (i.e., limited reaction sites), which leads to the formation of carbonized CFs with similar compositions. Furthermore, carbonized CFs possess distinct structural characteristics, and finally lead to the formation of CFs with different moduli. By analyzing CF crystallinity and elemental composition, general rationales are proposed for the selection of optimal precursors and processing conditions, and thus shed lights on parameter screening during industrial CF synthesis. For instance, inclusion of sp 3 carbons in cata- condensed rings can be preferred over the presence of sp 3 carbons on side chains.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Hydrogen Dynamics in Supercritical Water Probed by Neutron Scattering and Computer Simulations

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553–823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

NANO.PTML model for read-across prediction of nanosystems in neurosciences. computational model and experimental case of study

Abstract Neurodegenerative diseases involve progressive neuronal death. Traditional treatments often struggle due to solubility, bioavailability, and crossing the Blood-Brain Barrier (BBB). Nanoparticles (NPs) in biomedical field are garnering growing attention as neurodegenerative disease drugs (NDDs) carrier to the central nervous system. Here, we introduced computational and experimental analysis. In the computational study, a specific IFPTML technique was used, which combined Information Fusion (IF) + Perturbation Theory (PT) + Machine Learning (ML) to select the most promising Nanoparticle Neuronal Disease Drug Delivery (N2D3) systems. For the application of IFPTML model in the nanoscience, NANO.PTML is used. IF-process was carried out between 4403 NDDs assays and 260 cytotoxicity NP assays conducting a dataset of 500,000 cases. The optimal IFPTML was the Decision Tree (DT) algorithm which shown satisfactory performance with specificity values of 96.4% and 96.2%, and sensitivity values of 79.3% and 75.7% in the training (375k/75%) and validation (125k/25%) set. Moreover, the DT model obtained Area Under Receiver Operating Characteristic (AUROC) scores of 0.97 and 0.96 in the training and validation series, highlighting its effectiveness in classification tasks. In the experimental part, two samples of NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) were synthesized by thermal decomposition of an iron(III) oleate (FeOl) precursor and structurally characterized by different methods. Additionally, in order to make the as-synthesized hydrophobic NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) soluble in water the amphiphilic CTAB (Cetyl Trimethyl Ammonium Bromide) molecule was employed. Therefore, to conduct a study with a wider range of NP system variants, an experimental illustrative simulation experiment was performed using the IFPTML-DT model. For this, a set of 500,000 prediction dataset was created. The outcome of this experiment highlighted certain NANO.PTML systems as promising candidates for further investigation. The NANO.PTML approach holds potential to accelerate experimental investigations and offer initial insights into various NP and NDDs compounds, serving as an efficient alternative to time-consuming trial-and-error procedures.

60 APPLIED LIFE SCIENCES↗

NASA NTRS1977

Characterization of secondary ignition sources in unattended compartments and full-scale baseline test

The characteristics of five fuel loads burned within a metal lavatory were identified. In 15 of the tests the lavatory door remained closed for the 30-minute test period while in 15 additional tests the door was opened after the fire had developed. Upon completion of these tests the most severe source was selected for use in the baseline test. In the baseline test, the lavatory and adjacent panels, all of which were constructed of contemporary materials, were tested for a period of 1 hour. Thermal, environmental, and biological data were obtained for all fuel loads, door conditions, and the baseline test. All tests were conducted in a cabin fire simulator with separate ventilation of the cabin and lavatory representative of an inflight condition. The baseline test established that by using the most severe fuel source: (1) the exposed animal subject survived without complications; (2) no toxic levels of gas within the cabin were detected; (3) a propagating fire did not develop in adjacent structures; (4) the lavatory containing the fire remained structurally intact; (5) decomposition of portions of the lavatory did occur; and (6) cabin visibility would have presented a problem after 5 minutes.

Klink, D. M.↗

NASA NTRS1992

On testing models for the pressure-strain correlation of turbulence using direct simulations

Direct simulations of homogeneous turbulence have, in recent years, come into widespread use for the evaluation of models for the pressure-strain correlation of turbulence. While work in this area has been beneficial, the increasingly common practice of testing the slow and rapid parts of these models separately in uniformly strained turbulent flows is shown in this paper to be unsound. For such flows, the decomposition of models for the total pressure-strain correlation into slow and rapid parts is ambiguous. Consequently, when tested in this manner, misleading conclusions can be drawn about the performance of pressure-strain models. This point is amplified by illustrative calculations of homogeneous shear flow where other pitfalls in the evaluation of models are also uncovered. More meaningful measures for testing the performance of pressure-strain models in uniformly strained turbulent flows are proposed and the implications for turbulence modeling are discussed.

Speziale, Charles G.↗

NASA NTRS2007

POD Model Reconstruction for Gray-Box Fault Detection

Proper orthogonal decomposition (POD) is the mathematical basis of a method of constructing low-order mathematical models for the "gray-box" fault-detection algorithm that is a component of a diagnostic system known as beacon-based exception analysis for multi-missions (BEAM). POD has been successfully applied in reducing computational complexity by generating simple models that can be used for control and simulation for complex systems such as fluid flows. In the present application to BEAM, POD brings the same benefits to automated diagnosis. BEAM is a method of real-time or offline, automated diagnosis of a complex dynamic system.The gray-box approach makes it possible to utilize incomplete or approximate knowledge of the dynamics of the system that one seeks to diagnose. In the gray-box approach, a deterministic model of the system is used to filter a time series of system sensor data to remove the deterministic components of the time series from further examination. What is left after the filtering operation is a time series of residual quantities that represent the unknown (or at least unmodeled) aspects of the behavior of the system. Stochastic modeling techniques are then applied to the residual time series. The procedure for detecting abnormal behavior of the system then becomes one of looking for statistical differences between the residual time series and the predictions of the stochastic model.

Park, Han↗