Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials simulations”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Computational Fluid Dynamics Models of the NASA White Sands Test Facility Particle Impact Test Nozzles

The particle impact ignition mechanism is considered to be the most effective means for igniting metals in gaseous oxygen environments and is a concern for every oxygen system. For this reason, the National Aeronautics and Space Administration (NASA) commits a large amount of resources attempting to understand and defend against this phenomenon. The most practical means of accomplishing this is through extensive of testing. To this end, over the years NASA has developed and built a facility and various apparatuses to test material or component vulnerability to a particle impact ignition. In testing materials, NASA uses sub-sonic and supersonic nozzles to replicate the conditions that a material is expected to experience in service. However, challenges include obtaining an understanding of the gas dynamics and of the particle behavior through and up to the point of impact at the material under test. One means of addressing these challenges is the use of computational fluid dynamics (CFD). If a CFD model can be built and validated, the model can then be used to confidently interrogate changes to inlet conditions and particle geometries, as well as changes to nozzle geometries to better simulate a material’s in-situ conditions and greatly reduce or eliminate trial-and-error approaches in testing. This paper will show NASA’s approach to and results of using CFD to better understand their particle impact nozzles.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Stochastic Reconstruction of Thermal Protection Material Properties from Arc-Jet Experiments

Material response models are used to assess reliability using variances in the bond-line temperature predictions based on uncertainties in trajectory, aerothermal environment, and material properties. A key deficiency in the current approach is that input uncertainties are too often subjective, empirical, or ad-hoc, and are not rigorously linked to the arc-jet test data used to develop the TPS material model. While materials such as PICA are well understood, future missions may require more novel materials such as HEEET where unknown uncertainties have real consequences on the ability to assess reliability. A quantifiable estimate of reliability requires an iterative methodology where the parameters driving the variance in bond-line temperature (for example) are systematically identified. A test campaign to collect data or develop new models can then be identified to reduce those input uncertainties. A Bayesian inference loop defines these connections mathematically, i.e., prior knowledge about uncertainty is updated based on observation. While these concepts are well known (and often applied intuitively in a non-rigorous approach), only recent advances in reduced-order modelling have made them computationally viable methods for engineering. By replacing deterministic inverse methods with stochastic approaches, the hope is new materials proposed for future missions can more rapidly be developed with a greater understanding of the TPS material reliability. Two additional steps for the analysis of arc jet test data are discussed. The first is ability to construct a reduced-order model using material response simulations (Icarus/US3D) of the arc-jet test articles, and the second is the inclusion of this surrogate model in the Bayesian inversion process. Both capabilities will be demonstrated using prior PICA arc-jet test data. The quality of a surrogate model will be investigated and the variances on the calibrated material properties will be compared to our current understanding of the PICA material model.

Material response↗
DOE OSTI2026

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development & Experimental Validation of a Generalized Resistance-Capacitance Model for Numerical Simulation of Phase-Change Material Embedded Heat Exchangers

Latent heat thermal energy storage (LHTES) using phase change material (PCM) has attracted increased attention as a viable solution for overcoming the mismatch between energy supply and demand for renewable energy-based systems. PCM-embedded heat exchangers (PCM-HX) have the potential to significantly improve thermal performance due to high storage capacity and low temperature variation during the phase change process. Most models for simulating LHTES heat transfer use Computational Fluid Dynamics (CFD) simulations, which have high computational costs resulting from considering the complex and time-dependent physics relevant to PCM-HXs. In this paper, a Generalized Resistance Capacitance-based Model (GRCM) was developed to predict the thermal performance of arbitrary PCM-HXs in a computationally efficient manner without compromising modeling accuracy. The GRCM is exercised for three case studies: (i) verification for a single-slabbed finned PCM-HX, (ii) verification and validation for a copper foam/paraffin composite PCM-HX, and (iii) validation for a straight tube annular finned PCM-HX. The copper foam PCM-HX uses an electric heater at the top of HX, while the other two configurations utilize water as heat transfer fluid. For the single-slabbed finned PCM-HX melting case, the mean deviation in average PCM temperature predicted by the GRCM compared to the CFD model was between 0.56 – 0.73 K, with maximum temperature deviation of 2.68 K. For the HTF outlet temperature, the validation results showed that GRCM prediction matches very well with experimental data, with mean temperature deviation of 0.24 K during melting case, while for solidification case was 0.34 K. These results showcase the GRCM’s capability for accurately reproducing the thermal characteristics of PCM-HXs with considerably lower computational effort.

42 ENGINEERING↗
DOE OSTI2022

Extending and Smoothing Two-dimension Equation of State Simulation Data

The Multiphase Equation of State (MEOS) project works to produce high quality equation of state tables which are used in computationally intensive simulations of materials in different conditions. An equation of state table describes the state of matter under certain physical conditions such as pressure, volume, temperature, or internal energy. MEOS uses many different models and combinations thereof to produce accurate tables consisting of continuous smooth data, derivatives, and higher-order derivatives. Accurate and smooth data are important factors in producing precise simulations. This report focuses on improving the electron tables produced from Purgatorio and Thomas-Fermi data. Purgatorio data is extremely accurate but often jagged and discontinuous in the lower temperature and density region. Thomas-Fermi is a model that provides smooth data throughout. Together, this produces an EOS data table that is both accurate and smooth. This report describes the implementation of a new feature that allows users to define more detailed regions in the Purgatorio table to be replaced with Thomas-Fermi. As a result, MEOS can generate a table that is both more smooth and still accurate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Downfolding Complex Materials Problems Onto Model Hamiltonians for Quantum Computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios Markos Alvertis↗
NASA NTRSDate not supplied

Downfolding complex materials problems onto model Hamiltonians for quantum computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios M. Alvertis↗
NASA NTRS1983

Multiaxial Cyclic Thermoplasticity Analysis with Besseling's Subvolume Method

A modification was formulated to Besseling's Subvolume Method to allow it to use multilinear stress-strain curves which are temperature dependent to perform cyclic thermoplasticity analyses. This method automotically reproduces certain aspects of real material behavior important in the analysis of Aircraft Gas Turbine Engine (AGTE) components. These include the Bauschinger effect, cross-hardening, and memory. This constitutive equation was implemented in a finite element computer program called CYANIDE. Subsequently, classical time dependent plasticity (creep) was added to the program. Since its inception, this program was assessed against laboratory and component testing and engine experience. The ability of this program to simulate AGTE material response characteristics was verified by this experience and its utility in providing data for life analyses was demonstrated. In this area of life analysis, the multiaxial thermoplasticity capabilities of the method have proved a match for the actual AGTE life experience.

Mcknight, R. L.↗
NASA NTRS1998

Damage Progression in Bolted Composites

Structural durability,damage tolerance,and progressive fracture characteristics of bolted graphite/epoxy composite laminates are evaluated via computational simulation. Constituent material properties and stress and strain limits are scaled up to the structure level to evaluate the overall damage and fracture propagation for bolted composites. Single and double bolted composite specimens with various widths and bolt spacings are evaluated. The effect of bolt spacing is investigated with regard to the structural durability of a bolted joint. Damage initiation, growth, accumulation, and propagation to fracture are included in the simulations. Results show the damage progression sequence and structural fracture resistance during different degradation stages. A procedure is outlined for the use of computational simulation data in the assessment of damage tolerance, determination of sensitive parameters affecting fracture, and interpretation of experimental results with insight for design decisions.

Minnetyan, Levon↗
NASA NTRS1998

Damage Progression in Bolted Composites

Structural durability, damage tolerance, and progressive fracture characteristics of bolted graphite/epoxy composite laminates are evaluated via computational simulation. Constituent material properties and stress and strain limits are scaled up to the structure level to evaluate the overall damage and fracture propagation for bolted composites. Single and double bolted composite specimens with various widths and bolt spacings are evaluated. The effect of bolt spacing is investigated with regard to the structural durability of a bolted joint. Damage initiation, growth, accumulation, and propagation to fracture are included in the simulations. Results show the damage progression sequence and structural fracture resistance during different degradation stages. A procedure is outlined for the use of computational simulation data in the assessment of damage tolerance, determination of sensitive parameters affecting fracture, and interpretation of experimental results with insight for design decisions.

Minnetyan, Levon↗
DOE OSTI2021

A fourth-order Lagrangian discontinuous Galerkin method using a hierarchical orthogonal basis on curvilinear grids

The existing high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic methods are restricted to using quadratic meshes with quadratic polynomials (P2), which in turn, yield up to third-order accuracy. Additionally, these existing DG hydrodynamic schemes, when extended to work with cubic meshes and cubic polynomials (P3), can be unstable on strong-shock problems. Therefore, this paper presents a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate compressible material dynamics ( e.g., gasses, fluids, and solids) with strong-shocks using cubic meshes and cubic polynomials, and delivers up to fourth-order accuracy on smooth flows. The stability on shock problems is achieved using new hierarchical orthogonal basis functions and a new subcell mesh stabilization (SMS) scheme for cubic meshes. The accuracy and robustness of the new high-order accurate Lagrangian DG hydrodynamic method is demonstrated by simulating a diverse suite of challenging test problems covering gas and solid dynamic problems on curvilinear meshes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Liquid to crystal Si growth simulation using machine learning force field

Machine learning force field (ML-FF) has emerged as a potential promising approach to simulate various material phenomena for large systems with ab initio accuracy. However, most ML-FFs have been used to study the phenomena relatively close to the equilibrium ground states. In this work, we have studied a far from equilibrium system of liquid to crystal Si growth using ML-FF. Here, we found that our ML-FF based on ab initio decomposed atomic energy can reproduce all the aspects of ab initio simulated growth, from local energy fluctuations to transition temperatures, to diffusion constant, and growth rates. We have also compared the growth simulation with the Stillinger-Weber classical force field and found significant differences. A procedure is also provided to correct a systematic fitting bias in the ML-FF training process, which exists in all training models, otherwise critical results like transition temperature will be wrong.

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

On the measurement of shape: With applications to lunar regolith

With the renewed commitment from NASA and other commercial entities for a presence on the Moon, the importance of understanding the characteristics of lunar regolith and how to utilize it have become the target of increasing scrutiny. Much of what is known about lunar regolith was collected during and immediately after the Apollo program, however, analytical techniques and instrumentation have advanced in leaps and bounds in the subsequent decades. Specifically, dynamic image analysis systems have advanced to the point that millions of particles can have morphological characteristics automatically determined in relatively short time frames. Particle morphological data was collected on several lunar samples and a pair of widely used regolith simulants to ascertain the accuracy of these simulants and to explore statistical analysis methods of these large datasets. It is found that these morphology datasets can vary widely depending on the particle size of the particles, and simple averaging of the data skews the results heavily towards the numerically abundant size ranges, the fines. Different reporting methods are suggested to ameliorate these problems. When applied to the lunar regolith, the particles are noted to be less morphologically complex than initially suspected. Compared to the lunar material, the simulants are found to contain some more morphological variability and angular grains. Such difference is likely due to the wildly different comminution processes that these different powder systems are subjected to.

2D shape↗
NASA NTRS2005

JSC-1: Lunar Simulant of Choice for Geotechnical Applications and Oxygen Production

Lunar simulant JSC-1 was produced as the result of a workshop held in 1991 to evaluate the status of simulated lunar material and to make recommendations on future requirements and production of such material (McKay et al., 1991). JSC-1 was prepared from a welded tuff that was mined, crushed, and sized from the Pleistocene San Francisco volcanic field, northern Arizona. As the initial production of approxiamtely 12,300kgs is nearly depleted, new production has commenced. The mineralogy and chemical properties of JSC-1 are described in McKay et al. (1994) and Hill et al. (this volume); description of its geotechnical properties appears in Klosky et al. (1996). Although other lunar-soil simulants have been produced (e.g., MLS-1: Weiblen et al., 1990; Desai et al., 1992; Chua et al., 1994), they have not been as well standardized as JSC-I; this makes it difficult to standardize results from tests performed on these simulants. Here, we provide an overview of the composition, mineralogy, strength and deformation properties, and potential uses of JSC-1 and outline why it is presently the 'lunar simulant of choice' for geotechnical applications and as a proxy for lunar-oxygen production.

Taylor, Lawrence A.↗
NASA NTRS1987

Simulation studies for surfaces and materials strength

A realistic potential energy function comprising angle dependent terms was employed to describe the potential surface of the N+O2 system. The potential energy parameters were obtained from high level ab-initio results using a nonlinear fitting procedure. It was shown that the potential function is able to reproduce a large number of points on the potential surface with a small rms deviation. A literature survey was conducted to analyze exclusively the status of current small cluster research. This survey turned out to be quite useful in understanding and finding out the existing relationship between theoretical as well as experimental investigative techniques employed by different researchers. Additionally, the importance of the role played by computer simulation in small cluster research, was documented.

Halicioglu, Timur↗
NASA NTRS2002

Mixed-Mode Decohesion Finite Elements for the Simulation of Delamination in Composite Materials

A new decohesion element with mixed-mode capability is proposed and demonstrated. The element is used at the interface between solid finite elements to model the initiation and non-self-similar growth of delaminations. A single relative displacement-based damage parameter is applied in a softening law to track the damage state of the interface and to prevent the restoration of the cohesive state during unloading. The softening law for mixed-mode delamination propagation can be applied to any mode interaction criterion such as the two-parameter power law or the three-parameter Benzeggagh-Kenane criterion. To demonstrate the accuracy of the predictions and the irreversibility capability of the constitutive law, steady-state delamination growth is simulated for quasistatic loading-unloading cycles of various single mode and mixed-mode delamination test specimens.

Camanho, Pedro P.↗
DOE OSTI2025

Effect of pore fluid chemistry on the mechanical behavior of a divalent compacted bentonite, an experimental and constitutive study

Ongoing research in isolating high-level nuclear waste and spent fuel has highlighted compacted bentonite as a suitable material for engineered barrier systems in deep geological repositories due to its extraordinary swelling and retention properties. This research focuses on the chemo-mechanical behavior of compacted bentonite exposed to different pore fluids with different concentrations and loading conditions. The study involves swelling pressure and compressibility experiments along with mineralogy analysis employing X-ray diffraction (XRD) and Cation exchange. The tests were conducted on BCV (a Mg/Ca- bentonite) compacted at a dry density of 1.48 ± .02 Mg/m 3 . An advanced chemical-mechanical constitutive model for unsaturated highly expansive clays was adopted to simulate the material response and better understand its behavior. The model is able to account for the main phenomena at both macro and microstructural levels and the interactions between them. The model successfully replicated experimental observations. The XRD analyses support the macroscopic observation, indicating that salinity impacts crystalline swelling as demonstrated by the reduction of basal spacing from 19.27 Å to 15.68 Å when the osmotic suction increases from 0 MPa to 33 MPa. The results suggested that the osmotic pressure generated by the concentration in the pore fluids promotes a reduction in swelling pressures, swelling strains, and crystalline swelling of clay minerals. Also, it affects the pre-consolidation stress and the compressibility of the compacted samples. In conclusion, it was also observed that both solution type and solution concentration impact the clay swelling pressure.

Chemo-mechanical constitutive model↗