Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer simulation”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Effective Control of Computationally Simulated Wing Rock in Subsonic Flow

The unsteady compressible, full Navier-Stokes (NS) equations and the Euler equations of rigid-body dynamics are sequentially solved to simulate the delta wing rock phenomenon. The NS equations are solved time accurately, using the implicit, upwind, Roe flux-difference splitting, finite-volume scheme. The rigid-body dynamics equations are solved using a four-stage Runge-Kutta scheme. Once the wing reaches the limit-cycle response, an active control model using a mass injection system is applied from the wing surface to suppress the limit-cycle oscillation. The active control model is based on state feedback and the control law is established using pole placement techniques. The control law is based on the feedback of two states: the roll-angle and roll velocity. The primary model of the computational applications consists of a 80 deg swept, sharp edged, delta wing at 30 deg angle of attack in a freestream of Mach number 0.1 and Reynolds number of 0.4 x 10(exp 6). With a limit-cycle roll amplitude of 41.1 deg, the control model is applied, and the results show that within one and one half cycles of oscillation, the wing roll amplitude and velocity are brought to zero.

Kandil, Osama A.↗
NASA NTRS1997

Computational Simulation of Composite Behavior

Composite microfracture technique developed here is based on the concept of global strain energy release rate. Step-by-step procedures were outlined to evaluate composite microfracture. The technique can predict microfracture initiation/propagation, predict the "structural fracture toughness", and establish the hierarchy of fracture modes. Both unidirectional and cross-ply metal-matrix composite laminates were evaluated for microfracture, subjected to thermal and mechanical loads. Five NASA reports were published as shown in the bibliography and several presentations were made at the technical conferences/ symposia. Work was also performed to develop simplified micromechanics and macromechanics for ceramic matrix composites. A novel fiber sub-structuring technique was developed, which divides the fiber into several "slices" and applies the micromechanics equation at the "slice" level, which are then integrated upwards to get the ply and the laminate properties. It allows for a more accurate representation of interfacial conditions such as partial interface bond etc. It also predicts the stress-strain behavior of a composite up to failure by taking into account material properties dependence upon temperature as well as stress redistribution that occurs due to damage initiation and propagation. A computer code CEMCAN (Ceramic Matrix Composites Analyzer) which includes all these capabilities has been developed as a part of this research activity.

Keith, Theo G., Jr.↗
DOE OSTI2022

Concurrent Precipitation of Nb(C,N) and Metastable M 23 C 6 in Alloy 347H at 700°C and 750°C: Computer Simulations and Comparison to Experiment

Here, we present our results for the concurrent precipitation of metastable M 23 C 6 , Nb(C,N) secondary precipitates, and the Nb(C,N) primary crystals in 347H austenitic stainless steel. For precipitation modeling, we have accounted for the elastic contribution to interfacial energy, and for the Fe-spin-polarization for NbC/Fe and M 23 C 6 /Fe interfacial energy values: for NbC/Fe ~ 0.63 J/m 2 . For M 23 C 6 precipitates, an error function was used to describe the interfacial energy growth with particle size. In precipitation simulations, the average size of the primary Nb(C,N) particles remained ~ 1 μm at 700°C and ~ 0.3 μm at 750°C. The M 23 C 6 precipitates at 750°C dissolved after 120 h (our simulations) compared to 300 h (experiments). The Nb(C,N)/Fe interfacial energy was not affected by the nitrogen additions. With these modifications, reasonable agreement with the available experimental data was obtained, which allows using them in the development of the 2nd-phase particle-informed creep theory.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
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 OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computer Simulation of Complex Systems at the Extremes of pH (Final Report)

During the final funding period of the above-mentioned grant, a number of key milestones were achieved as they relate to the Separation Science program. Below are the descriptions of those efforts and the papers published in support of them (see also submitted publication list). There was a strong focus during the final funding period was on hydrated protons and the unique properties of acidic solutions, which as pertinent to many separations processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗