Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Tool”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Randomized algorithms for generalized singular value decomposition with application to sensitivity analysis

The generalized singular value decomposition (GSVD) is a valuable tool that has many applications in computational science. However, computing the GSVD for large-scale problems is challenging. Motivated by applications in hyper-differential sensitivity analysis (HDSA), in this work we propose new randomized algorithms for computing the GSVD which use randomized subspace iteration and weighted QR factorization. Detailed error analysis is given which provides insight into the accuracy of the algorithms and the choice of the algorithmic parameters. We demonstrate the performance of our algorithms on test matrices and a large-scale model problem where HDSA is used to study subsurface flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Anticipating gelation and vitrification with medium amplitude parallel superposition (MAPS) rheology and artificial neural networks

Abstract Anticipating qualitative changes in the rheological response of complex fluids (e.g., a gelation or vitrification transition) is an important capability for processing operations that utilize such materials in real-world environments. One class of complex fluids that exhibits distinct rheological states are soft glassy materials such as colloidal gels and clay dispersions, which can be well characterized by the soft glassy rheology (SGR) model. We first solve the model equations for the time-dependent, weakly nonlinear response of the SGR model. With this analytical solution, we show that the weak nonlinearities measured via medium amplitude parallel superposition (MAPS) rheology can be used to anticipate the rheological aging transitions in the linear response of soft glassy materials. This is a rheological version of a technique called structural health monitoring used widely in civil and aerospace engineering. We design and train artificial neural networks (ANNs) that are capable of quickly inferring the parameters of the SGR model from the results of sequential MAPS experiments. The combination of these data-rich experiments and machine learning tools to provide a surrogate for computationally expensive viscoelastic constitutive equations allows for rapid experimental characterization of the rheological state of soft glassy materials. We apply this technique to an aging dispersion of Laponite ® clay particles approaching the gel point and demonstrate that a trained ANN can provide real-time detection of transitions in the nonlinear response well in advance of incipient changes in the linear viscoelastic response of the system.

Lennon, Kyle R. (ORCID:0000000212515461)↗
DOE OSTI2019

A Comparative Study of the Efficacy of Local/Global and Parametric/Nonparametric Machine Learning Methods for Establishing Structure–Property Linkages in High-Contrast 3D Elastic Composites

Reduced-order structure–property (S-P) linkages play a pivotal role in the tailored design of materials for advanced engineering components. There is a critical need to distill these from the simulation datasets aggregated using sophisticated, computationally expensive, physics-based numerical tools (e.g., finite element methods). The recent emergence of materials data science approaches has opened new avenues for addressing this challenge. In this paper, we critically compare the relative merits of the application of four distinct machine learning approaches for their efficacy in extracting microstructure-property linkages from the finite element simulation data aggregated on high-contrast elastic composites with different microstructures. The machine learning approaches selected for the study have included different combinations of local/global and parametric/nonparametric approaches. Furthermore, the nonparametric approaches selected for this study are based on Gaussian Process (GP) models that allow for a formal treatment of uncertainty quantification in the predicted values. The predictive performances of these different approaches have been compared against each other using rigorous cross-validation error metrics. Furthermore, their sensitivity to both the dataset size and dimensionality has been investigated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A simplified CFD approach for modeling mass transport in catalytic open-cell foams

A simplified macroscopic CFD approach is presented to model mass transport including chemical reactions in washcoated open-cell foams. The foam is treated as a porous medium. Species conversion during chemical reactions is modeled using appropriate source terms based on reaction rate expressions and modified to account for the mass transport resistances occurring at the fluid-washcoat interfaces and within the washcoat layers. As example, the catalytic CO oxidation over platinum is studied. The simulation results show good agreement with experimental data from literature. A parametric study on washcoat parameters, such as thickness, tortuosity, porosity, and size, is carried out. Increasing the washcoat thickness from 5 to 100 µm or decreasing the tortuosity to porosity ratio from 5 to 20 decreases the CO conversion by 10 %. The proposed model is found to be reliable and has the advantage of lower computational cost, making it a suitable tool for foam-based catalytic reactor design.

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

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2023

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗
DOE OSTI2019

Modeling framework for multiphysics-multiscale behavior of Si–C composite anode

To achieve the urgent requirement for high energy density in lithium-ion batteries, Si-C composite anode has been spotlighted as one of the most promising alternatives for next-generation batteries. Additionally, a comprehensive computational model is indispensable in the developing process of the excellent performance of anode material due to the low-realizability, inconvenience, and high-cost of experiments. Herein, a simultaneous multiphysics-multiscale model aiming at the Si-C composite anode was developed. This model was then used to study the effects of Si percentage, mechanical constraint and charging rate in terms of electrochemical and mechanical performances. A better design to achieve high capacity was proposed based on the computational results. Results provide a powerful tool in the design, development, and evaluation of high energy density lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2019

3D multiphysics simulation and analysis of a low temperature liquid metal magnetohydrodynamic power generator prototype

Magnetohydrodynamic (MHD) power generation is an interdisciplinary process, which was studied heavily in the late 1960’s but interest soon subsided. However, the current socioeconomic climate demands more diversity in renewable energy sources. MHD generators using liquid metal as the working fluid represent a viable option for this augmentation. Additionally, the steady increase in accessible computational capabilities provides a powerful new tool for interdisciplinary modeling and multiphysics simulations. These factors suggest that this is an opportune time to reevaluate MHD power generation. This paper provides mathematical background as well as a detailed multiphysics model of a low temperature liquid metal MHD power generator built in COMSOL. Simulation details are provided in addition to a comparison to experimental results. Finally, the prototype studied in this paper is the first proof of concept of a permanent magnet excited liquid metal MHD generator using Gallium as its working fluid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Toward data assimilation of ship-induced aerosol–cloud interactions

Satellite imagery can detect temporary cloud trails or ship tracks formed from aerosols emitted from large ships traversing our oceans, a phenomenon that global climate models cannot directly reproduce. Ship tracks are observable examples of marine cloud brightening, a potential solar climate intervention that shows promise in helping combat climate change. In this paper, we demonstrate a simulation-based approach in learning the behavior of ship tracks based upon a novel stochastic emulation mechanism. Our method uses wind fields to determine the movement of aerosol–cloud tracks and uses a stochastic partial differential equation (SPDE) to model their persistence behavior. This SPDE incorporates both a drift and diffusion term which describes the movement of aerosol particles via wind and their diffusivity through the atmosphere, respectively. We first present our proposed approach with examples using simulated wind fields and ship paths. We then successfully demonstrate our tool by applying the approximate Bayesian computation method-sequential Monte Carlo for data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Discovering a Transferable Charge Assignment Model Using Machine Learning

Partial atomic charge assignment is of immense practical value to force field parametrization, molecular docking, and cheminformatics. Machine learning has emerged as a powerful tool for modeling chemistry at unprecedented computational speeds given accurate reference data. However, certain tasks, such as charge assignment, do not have a unique solution. We use a machine learning algorithm to discover a new charge assignment model by learning to replicate molecular dipole moments across a large, diverse set of nonequilibrium conformations of molecules containing C, H, N, and O atoms. The new model, called Affordable Charge Assignment (ACA), is computationally inexpensive and predicts dipoles of out-of-sample molecules accurately. Furthermore, dipole-inferred ACA charges are transferable to dipole and even quadrupole moments of much larger molecules than those used for training. We apply ACA to dynamical trajectories of biomolecules and produce their infrared spectra. Additionally, we find that ACA assigns similar charges to Charge Model 5 but with greatly reduced computational cost.

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

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Carbazole-substituted benzobisoxazoles: near-UV fluorescent emitters and ambipolar hosts for organic light-emitting diodes

The design of organic light emitting diode (OLED) materials with the potential for exhibiting thermally-activated delayed fluorescence (TADF) is reported. Using computational methods (DFT/TD-DFT) as a guiding tool, six materials with a benzobisoxazole (BBO) core and donor–acceptor–donor architectures were designed by changing the conjugation position of carbazole-substituted phenyl substituents and the type of BBO isomer (cis- vs. trans-). Experimental steady-state and transient absorption spectroscopic techniques were utilized to probe the TADF activity of these molecules. Each material was then used in host–guest OLED devices as either near-UV dopants or host with low singlet-triplet energy differences. The electroluminescent properties show that when used as dopants these materials provide near-UV emission (CIE y < 0.06 and CIE x = 0.16), whereas when used as hosts, these materials show reduced operating voltages and increased performance efficiencies when compared to commercial materials.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Oxygen transport properties estimation by DSMC-CT simulations

Coupling DSMC simulations with classical trajectories calculations is emerging as a powerful tool to improve predictive capabilities of computational rarefied gas dynamics. The considerable increase of computational effort outlined in the early application of the method (Koura,1997) can be compensated by running simulations on massively parallel computers. In particular, GPU acceleration has been found quite effective in reducing computing time (Ferrigni,2012; Norman et al.,2013) of DSMC-CT simulations. The aim of the present work is to study rarefied Oxygen flows by modeling binary collisions through an accurate potential energy surface, obtained by molecular beams scattering (Aquilanti, et al.,1999). The accuracy of the method is assessed by calculating molecular Oxygen shear viscosity and heat conductivity following three different DSMC-CT simulation methods. In the first one, transport properties are obtained from DSMC-CT simulations of spontaneous fluctuation of an equilibrium state (Bruno et al, Phys. Fluids, 23, 093104, 2011). In the second method, the collision trajectory calculation is incorporated in a Monte Carlo integration procedure to evaluate the Taxman’s expressions for the transport properties of polyatomic gases (Taxman,1959). In the third, non-equilibrium zero and one-dimensional rarefied gas dynamic simulations are adopted and the transport properties are computed from the non-equilibrium fluxes of momentum and energy. The three methods provide close values of the transport properties, their estimated statistical error not exceeding 3%. The experimental values are slightly underestimated, the percentage deviation being, again, few percent.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Oxygen transport properties estimation by classical trajectory–direct simulation Monte Carlo

Coupling direct simulation Monte Carlo (DSMC) simulations with classical trajectory calculations is a powerful tool to improve predictive capabilities of computational dilute gas dynamics. The considerable increase in computational effort outlined in early applications of the method can be compensated by running simulations on massively parallel computers. In particular, Graphics Processing Unit acceleration has been found quite effective in reducing computing time of classical trajectory (CT)-DSMC simulations. The aim of the present work is to study dilute molecular oxygen flows by modeling binary collisions, in the rigid rotor approximation, through an accurate Potential Energy Surface (PES), obtained by molecular beams scattering. The PES accuracy is assessed by calculating molecular oxygen transport properties by different equilibrium and non-equilibrium CT-DSMC based simulations that provide close values of the transport properties. Comparisons with available experimental data are presented and discussed in the temperature range 300–900 K, where vibrational degrees of freedom are expected to play a limited (but not always negligible) role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Fused x-ray and fast neutron CT reconstruction for imaging large and dense objects

Megavolt x-ray computed tomography (CT) is a powerful tool for three-dimensional characterization. However, its utility is limited for large objects composed of high-atomic number (Z) materials, where x rays fail to penetrate. Information from fast neutron CT (FNCT) can complement x-ray CT reconstructions since fast neutrons can more readily penetrate high-Z objects. In this work, we demonstrate a method for combining FNCT and x-ray CT data to create a single reconstruction, more accurate than could be achieved with either x rays or fast neutrons alone. The algorithm was tested on an exemplar comprising multiple concentric, nested cylinders of different materials. Simulated and empirical x-ray CT data were acquired for the exemplar using a 9 MV bremsstrahlung spectrum. Additional simulated and empirical FNCT data were acquired using an accelerator based fast neutron source. The FNCT data were used to synthesize x-ray CT data and augment the x-ray CT data missing due to lack of penetration. This approach mitigates artifacts that would otherwise negatively affect the accuracy and resolution of a single-modality reconstructed volume.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Computational quantum chemistry of metal–organic frameworks

Metal–organic frameworks (MOFs) have premium exceptional properties for a variety of functions, such as gas separation and storage and catalysis. The large variety of possible inorganometallic nodes and organic linkers provide an almost unlimited number of combinations for assembling MOFs, which makes the experimental characterization and examination of all potentially useful combinations practically impossible. Furthermore, experimental studies of MOFs typically fall short in uncovering crucial details regarding their mechanisms of action or the molecular details responsible for their functional properties, such as the nature of adsorbate binding or the structures of transition states. Computational modeling has, therefore, become an efficient and important tool for strategizing the functionalization of MOFs and explicating the mechanisms of their functions. Furthermore, we review the computational methodologies used for computational studies of MOFs, especially Kohn–Sham density functional theory and combined quantum mechanical and molecular mechanical methods for calculating their structural, electronic, and magnetic properties, as well as for understanding the mechanisms of MOFs' applications to magetic devices, thermal conduction, gas adsorption, separation, storage, and sensing, thermal catalysis, photocatalysis, and electrocatalysis.

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

Engineering an improved light-induced dimer (iLID) for controlling the localization and activity of signaling proteins

The discovery of light-inducible protein–protein interactions has allowed for the spatial and temporal control of a variety of biological processes. To be effective, a photodimerizer should have several characteristics: it should show a large change in binding affinity upon light stimulation, it should not cross-react with other molecules in the cell, and it should be easily used in a variety of organisms to recruit proteins of interest to each other. In this study, to create a switch that meets these criteria we have embedded the bacterial SsrA peptide in the C-terminal helix of a naturally occurring photoswitch, the light-oxygen-voltage 2 (LOV2) domain from Avena sativa. In the dark the SsrA peptide is sterically blocked from binding its natural binding partner, SspB. When activated with blue light, the C-terminal helix of the LOV2 domain undocks from the protein, allowing the SsrA peptide to bind SspB. Without optimization, the switch exhibited a twofold change in binding affinity for SspB with light stimulation. Here, we describe the use of computational protein design, phage display, and high-throughput binding assays to create an improved light inducible dimer (iLID) that changes its affinity for SspB by over 50-fold with light stimulation. A crystal structure of iLID shows a critical interaction between the surface of the LOV2 domain and a phenylalanine engineered to more tightly pin the SsrA peptide against the LOV2 domain in the dark. Finally, we demonstrate the functional utility of the switch through light-mediated subcellular localization in mammalian cell culture and reversible control of small GTPase signaling.

59 BASIC BIOLOGICAL SCIENCES↗