Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational chemistry”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A new technique for low-to-high altitude predictions of ablative hypersonic flowfields

A parabolized Navier-Stokes (PNS) technique for predicting chemically reacting reentry flows over a wide range of altitudes extending from zero to 259 kft is presented. This coupled-chemistry PNS scheme can use either a strongly implicit solution approach or a two-step solution approach and is used to study a hypersonic flow over various configurations over a Mach-number range of 6.89 to 20. It is demonstrated that the two-step solution provides a combination of accuracy, stability, and computational efficiency. Nonablating low-altitude calculations show that the flowfield is in chemical equilibrium only in the near-wall region, while finite-rate chemistry effects are significant over a major portion of the shock layer, and the Blottner-type linearization of the production terms is unstable for low-altitude high-Reynolds-number conditions.

Bhutta, Bilal A.↗
NASA NTRS1993

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗
NASA NTRS2012

Compositional Effects on Nickel-Base Superalloy Single Crystal Microstructures

Fourteen nickel-base superalloy single crystals containing 0 to 5 wt% chromium (Cr), 0 to 11 wt% cobalt (Co), 6 to 12 wt% molybdenum (Mo), 0 to 4 wt% rhenium (Re), and fixed amounts of aluminum (Al) and tantalum (Ta) were examined to determine the effect of bulk composition on basic microstructural parameters, including gamma' solvus, gamma' volume fraction, volume fraction of topologically close-packed (TCP) phases, phase chemistries, and gamma - gamma'. lattice mismatch. Regression models were developed to describe the influence of bulk alloy composition on the microstructural parameters and were compared to predictions by a commercially available software tool that used computational thermodynamics. Co produced the largest change in gamma' solvus over the wide compositional range used in this study, and Mo produced the largest effect on the gamma lattice parameter and the gamma - gamma' lattice mismatch over its compositional range, although Re had a very potent influence on all microstructural parameters investigated. Changing the Cr, Co, Mo, and Re contents in the bulk alloy had a significant impact on their concentrations in the gamma matrix and, to a smaller extent, in the gamma' phase. The gamma phase chemistries exhibited strong temperature dependencies that were influenced by the gamma and gamma' volume fractions. A computational thermodynamic modeling tool significantly underpredicted gamma' solvus temperatures and grossly overpredicted the amount of TCP phase at 982 C. Furthermore, the predictions by the software tool for the gamma - gamma' lattice mismatch were typically of the wrong sign and magnitude, but predictions could be improved if TCP formation was suspended within the software program. However, the statistical regression models provided excellent estimations of the microstructural parameters based on bulk alloy composition, thereby demonstrating their usefulness.

MacKay, Rebecca A.↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
DOE OSTI2024

Unveiling the Redox Noninnocence of Metallocorroles: Exploring K-Edge X-ray Absorption Near-Edge Spectroscopy with a Multiconfigurational Wave Function Approach

X-ray absorption near-edge spectroscopy (XANES) is an advanced technique for probing the local electronic structure of catalysts, effectively identifying the noninnocent nature of ligands in transition-metal complexes. Metallocorroles with noninnocent corrole rings exhibit unusual electronic structures that challenge traditional density functional theory (DFT) methods, necessitating more rigorous approaches to describe electron correlation accurately. We explored K-edge XANES spectra of Fe, Mn, and Co metallocorroles using TDDFT and wave function-based methods. This is the first investigation employing multireference methods, specifically RASSCF, RASPT2, and MC-PDFT, to analyze the redox noninnocent nature of metallocorroles reflected in their XANES spectra. We quantified the noninnocent character of the corrole and the oxidation states of the metals, capturing more than singly excited excitations responsible for the pre-edge peak. Our findings demonstrate the importance of these advanced computational techniques for accurately predicting XANES spectra, providing a reliable understanding of the electronic properties of such complexes. In conclusion, this study offers a new strategy for investigating ligand redox noninnocence via integrated experimental and computational XANES.

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

Anisotropic interactions for continuum modeling of protein–membrane systems

In this work, a model for anisotropic interactions between proteins and cellular membranes is proposed for large-scale continuum simulations. The framework of the model is based on dynamic density functional theory, which provides a formalism to describe the lipid densities within the membrane as continuum fields while still maintaining the fidelity of the underlying molecular interactions. Within this framework, we extend recent results to include the anisotropic effects of protein–lipid interactions. As applications, we consider two membrane proteins of biological interest: a RAS–RAF complex tethered to the membrane and a membrane embedded G protein-coupled receptor. A strong qualitative and quantitative agreement is found between the numerical results and the corresponding molecular dynamics simulations. Combining the scope of continuum level simulations with the details from molecular level particle simulations enables research into protein–membrane behaviors at a more biologically relevant scale, which crucially can also be accessed via experiment.

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

Integration of LIBS with Machine Learning for Real-Time Monitoring of Feedstock in H 2 Gasification Applications

This project, funded by the U.S. Department of Energy (DOE) – Office of Fossil Energy under Award Number DE-FE0032177, aimed to assess the feasibility of an integrated Laser-Induced Breakdown Spectroscopy (LIBS) system with advanced machine learning (ML) models for real-time characterization and potential control of hydrogen gasifiers running on waste materials as feedstocks. This was a multidisciplinary effort that encompassed the acquisition and standardized analysis of individual and blended feedstocks—comprising biomass, coal waste, and plastic waste, followed by the development of a dynamic LIBS bench system for material sample analysis and development of predictive ML models. Comprehensive laboratory testing enabled the creation of a robust elemental dataset that served as the foundation for ML model training. Techniques such as Random Forest, Gradient Boosting, Support Vector Regression, and Neural Networks were employed to predict key feedstock properties, including higher heating value (HHV), moisture content, thermal conductivity, and ash composition with high accuracy. The results were validated against experimental data and demonstrated strong potential for real-time application in gasifier control systems. The project concluded with a study on the integration of the LIBS+ML approach for gasifier control and a techno-economic analysis of the implementation of the approach into hydrogen (H 2 ) gasification systems. Dissemination of results was carried out at a DOE meeting. This work establishes a scalable framework for automated, in-line feedstock quality assessment, offering significant implications for process optimization and emissions reduction in hydrogen production.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Subsurface mafic and ultramafic rock mapping and analysis for carbon mineralization in the United States (SubMAP-CO 2 )

Achieving durable, gigaton-scale carbon dioxide (CO₂) removal will require the development of storage strategies that are both permanent and scalable. Among the available options, carbon mineralization in mafic and ultramafic rocks offers a uniquely robust pathway for long-term CO₂ sequestration, converting CO₂ into thermodynamically stable carbonate minerals. While surface exposures of these lithologies have been extensively studied, particularly in ophiolitic and volcanic settings, the potential for carbon mineralization in the subsurface remains comparatively underexplored at national scale. This gap limits the ability to assess storage capacity, accessibility, and deployment pathways for mineralization-based carbon management in the United States.

08 HYDROGEN↗
DOE OSTI2026

Electronic and Structural Properties of the Polymer-Electrolyte Interphase in Electrochemically Doped Polymers

The advance of soft, polymer-based (photo)electrochemical energy transformation and storage applications requires a framework for polymer and electrolyte design that includes a deep understanding of the polymer–electrolyte interphase. Here, we report on the investigations of two napthalenediimide (NDI)–bithiophene (T2)-based semiconducting copolymers using computational modeling and in situ, ex situ, and operando techniques to reveal how changes in electrolyte and polymer chemistry modulate the electronic and structural properties of polymer electrodes during electrochemical (de)doping. These systems are shown to host an ensemble of polarons, in contrast with the single polaron-like character often reported, whose properties vary with the nature of the local environment. Importantly, these polarons serve as reporters of the nanoscale environments in which they reside. We demonstrate that controlling the polymer and electrolyte chemistry regulates the nature of the charge carriers generated upon electrochemical doping and/or exciton dissociation in a photoelectrochemical solar cell: For instance, divalent counterions enable polaron and bipolaron formation at lower reducing potentials, while supporting more bipolaron formation than monovalent counterions. A novel application of NEXAFS reveals insights into charge (de)localization, providing a pathway for future investigation of electron transport mechanisms. Finally, simulations of polymer swelling of an amorphous interphase show that charge formation has a large impact on polymer swelling and ion penetration. These studies deliver insights to enable the control of charge-carrier and ion transport, the rates of electron transfer and catalytic efficiency, device stability, and overall device performance.

14 SOLAR ENERGY↗
NASA NTRS2002

An Improved Retrieval of Tropospheric Nitrogen Dioxide from GOME

We present a retrieval of tropospheric nitrogen dioxide (NO2) columns from the Global Ozone Monitoring Experiment (GOME) satellite instrument that improves in several ways over previous retrievals, especially in the accounting of Rayleigh and cloud scattering. Slant columns, which are directly fitted without low-pass filtering or spectral smoothing, are corrected for an artificial offset likely induced by spectral structure on the diffuser plate of the GOME instrument. The stratospheric column is determined from NO2 columns over the remote Pacific Ocean to minimize contamination from tropospheric NO2. The air mass factor (AMF) used to convert slant columns to vertical columns is calculated from the integral of the relative vertical NO2 distribution from a global 3-D model of tropospheric chemistry driven by assimilated meteorological data (Global Earth Observing System (GEOS)-CHEM), weighted by altitude dependent scattering weights computed with a radiative transfer model (Linearized Discrete Ordinate Radiative Transfer), using local surface albedos determined from GOME observations at NO2 wavelengths. The AMF calculation accounts for cloud scattering using cloud fraction, cloud top pressure, and cloud optical thickness from a cloud retrieval algorithm (GOME Cloud Retrieval Algorithm). Over continental regions with high surface emissions, clouds decrease the AMT by 20- 30% relative to clear sky. GOME is almost twice as sensitive to tropospheric NO2 columns over ocean than over land. Comparison of the retrieved tropospheric NO2 columns for July 1996 with GEOS-CHEM values tests both the retrieval and the nitrogen oxide radical (NOx) emissions inventories used in GEOS-CHEM. Retrieved tropospheric NO2 columns over the United States, where NOx emissions are particularly well known, are within 18% of GEOS-CHEM columns and are strongly spatially correlated (r = 0.78, n = 288, p less than 0.005). Retrieved columns show more NO2 than GEOS-CHEM columns over the Transvaal region of South Africa and industrial regions of the northeast United States and Europe. They are lower over Houston, India, eastern Asia, and the biomass burning region of central Africa, possibly because of biases from absorbing aerosols.

Martin, Randall V.↗
NASA NTRS2005

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing During Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on simulations using the NASA Global Modeling and Assimilation Office (GMAO) finite-volume general circulation model (fvGCM) with Microphyics of clouds in Relaxed Arakawa Schubert Scheme (McRAS). The aerosol loading are prescribed from three-dimensional monthly distribution of tropospheric aerosols viz., sulfate, black carbon, organic carbon, soil dust, and sea salt from output of the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The aerosol extinction coefficient, single scattering albedo, and asymmetric factor are computed as wavelength-dependent radiative forcing in the radiative transfer scheme of the fvGCM, and as a function of the aerosol loading and ambient relative humidity. We find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excites a planetary scale teleconnection pattern in sea level pressure, temperature and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, i.e., South Asia, East Asia, and northern and western Africa. Additionally, atmospheric heating is found in regions with large loading of dust (over Northern Africa, and Middle East), and black carbon (over South-East Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east-west dipole anomaly with strong cooling over the Caspian Sea, and warming over central and northeastern Asia, where aerosol concentration are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations.

Kim, Maeng-Ki↗
NASA NTRS2005

Rare Gases Inserted into Biological Building Blocks: A Theoretical Study of Glycine - Rg Compounds (Rg-Xe, Kr, Ar)

Compounds formed by insertion of rare-gas atoms (Xe, Kr, and Ar) into glycine molecule are investigated using accurate ab initio computational methods. Identification of such insertion compounds may open new frontiers in the field of rare-gas chemistry, such as possible existence of biological molecules that include chemically bound rare gas atoms. The most stable glycine-Rg configuration is found to correspond to insertion of Rg atoms into the 0-H bond of glycine. These NH2CH2COORgH compounds are metastable , but separated by sizable potential barriers from the Rg + glycine dissociation products. Preliminary calculations show that NH2CH2COOXeH compound is energetically stable with respect to another (3-body) dissociation channel (NH2CH2COO + Rg + H), while the corresponding Ar species is not stable in this respect. The compound with the inserted Kr is a borderline case, with the 3-body dissociation products being close in energy to the NH2CH2COOKrH minimum.

Chaban, Galina M.↗
NASA NTRSDate not supplied

Simulation of the Aerosol Size Distribution Using a Neural Network Surrogate for the Modal Aerosol Module (MAM7)

One objective of atmospheric simulations is to quantify the distribution of aerosols and their properties. Accurate parameterizations of the processes governing aerosol mass, particle number, and particle size distribution are important for predicting the Earth’s net radiative balance and aerosol-cloud interactions. The Modal Aerosol Module (MAM7) is a two-moment aerosol model that simulates mass, number, and size distribution of seven modes comprised of internally mixed aerosol species. The two-moment scheme adds significant computational expense but allows for the prediction of varying particle size distribution relative to the bulk method which predicts only total mass. In this work, we developed a neural network surrogate model for MAM7 (MAMnet) to predict the aerosol number concentration in NASA’s Global Earth Observing System (GEOS) without adding prohibitive computational expense. MAMnet, can be driven by output from a single moment, mass-based, aerosol scheme (Goddard Chemistry Aerosol and Radiation model (GOCART)) or from reanalysis products (Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2)). MAMnet was trained using number concentrations from a 5-year GEOS/MAM7 simulation at 1-degree horizontal resolution and using the total mass calculated across modes as inputs, as well as temperature and air density. The model architecture for MAMnet was based on AlexNet, the 2012 winner of the ImageNet Large Scale Visual Recognition Challenge. While some modifications were necessary to accommodate our problem, important aspects of the network were preserved. MAMnet was able to reproduce zonal dynamics and spatial distributions of the aerosol number concentration however predictability in the upper troposphere was poor.

Katherine H Breen↗
DOE OSTI2024

Effect of Si Content in Electrode and SiO2 Additions to the Slag during Electroslag Remelting

Evolution of Si concentration in 316 stainless steel electrodes was observed during a melting campaign consisting of recycling electroslag remelted (ESR) ingots to make new electrodes using vacuum induction melting (VIM). This campaign consisted of iterations of VIM + ESR operations to optimize melting parameters. The effect of Si content on the melt parameters and ingot quality was further evaluated and additions of SiO2 to the slag chemistry were studied using research-scale experiments, x-ray diffraction (XRF), combustion analysis, visual inspections, and computational tools. The Si concentration was found to decrease by approximately 600 ppm following ESR of 150 lb. research-scale electrodes. Eventually, this led to failure of the slag skin and direct ingot/crucible contact. Additions of SiO2 to the slag at levels matching the original calculated Si concentration in the electrode did not eliminate the slag skin failure and the current during steady state increased to maintain a constant melt rate. In this investigation, we propose a mechanism of slag skin failure consisting of local concentration of current density due to absence, or breakage, of the SiO2 layer around the molten metal drop during ESR. This theory was reinforced by additional experiments in which Nb was added to the electrode to change the structure of the oxide layer around the drops.

Jablonski, Paul↗
NASA NTRS1987

Insights into the general circulation of the lower stratosphere from TOMS

Total ozone is controlled by dynamical advection as well as chemistry. At least for day-to-day variations, dynamical processes appear to be in control of total ozone. There also appears to be good evidence that seasonal and secular changes in total ozone are also dynamically controlled. For example, the zonal mean changes in total ozone in the two hemispheres in spring appear to be quite different. The TOMS total ozone data suggest a south polar spring upwelling while the Northern Hemisphere shows a clear downwelling in the same period. Radiative transfer computations support this conclusion. The secular changes in total ozone over the South Pole in spring indicate a change in dynamics rather than chemistry.

Schoeberl, Mark R.↗
NASA NTRS1979

A sensitivity analysis of volcanic aerosol dispersion in the stratosphere

A computer sensitivity analysis was performed to determine the uncertainties involved in the calculation of volcanic aerosol dispersion in the stratosphere using a 2 dimensional model. The Fuego volcanic event of 1974 was used. Aerosol dispersion processes that were included are: transport, sedimentation, gas phase sulfur chemistry, and aerosol growth. Calculated uncertainties are established from variations in the stratospheric aerosol layer decay times at 37 latitude for each dispersion process. Model profiles are also compared with lidar measurements. Results of the computer study are quite sensitive (factor of 2) to the assumed volcanic aerosol source function and the large variations in the parameterized transport between 15 and 20 km at subtropical latitudes. Sedimentation effects are uncertain by up to a factor of 1.5 because of the lack of aerosol size distribution data. The aerosol chemistry and growth, assuming that the stated mechanisms are correct, are essentially complete in several months after the eruption and cannot explain the differences between measured and modeled results.

Butler, C. F.↗
NASA NTRS1990

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.

Yungster, Shaye↗