Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Instabilities and turbulence in low-$β$ guide field reconnection exhausts with kinetic Riemann simulations

The role of turbulence in low-$β$, guide-field reconnection exhausts is explored in 2D reconnection and 2D and 3D Riemann simulations. The structure of the exhaust and associated turbulence is controlled by a pair of rotational discontinuities (RDs) at the exhaust boundary and a pair of slow shocks (SSs) that are generated by counterstreaming ions beams. In 2D, the exhaust develops large-amplitude striations at the ion Larmor radius scale that are produced by electron-beam-driven ion cyclotron waves. The electron beams driving the instability are injected into the exhaust from one of the RDs. However, in 3D Riemann simulations, the additional dimension (in the out-of-plane direction) results in strong Buneman and electron-electron streaming instabilities at the RD which suppress the electron beam formation and therefore the striations in the exhaust. The strength of the streaming instabilities at the RD is controlled by the ratio of the electron thermal speed to Alfvén speed, with the lower thermal speed being more unstable. In the 3D simulations, an ion-ion streaming instability acts to partially thermalize the counterstreaming ion beams at the SSs. This instability is controlled by the ratio of the sound speed to Alfvén speed and is expected to be stable in the low $β$ solar corona. The results suggest that in a guide field reconnection exhaust with 1 $\gg$ $β >$ $m_e/m_i$, the kinetic-scale turbulence that develops will be too weak to play a significant role in energy conversion and particle acceleration. Therefore, the energy conversion will be mostly controlled by laminar physics or multi-x-line reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

An Uncertainty Structure Matrix for Models and Simulations

Software that is used for aerospace flight control and to display information to pilots and crew is expected to be correct and credible at all times. This type of software is typically developed under strict management processes, which are intended to reduce defects in the software product. However, modeling and simulation (M&S) software may exhibit varying degrees of correctness and credibility, depending on a large and complex set of factors. These factors include its intended use, the known physics and numerical approximations within the M&S, and the referent data set against which the M&S correctness is compared. The correctness and credibility of an M&S effort is closely correlated to the uncertainty management (UM) practices that are applied to the M&S effort. This paper describes an uncertainty structure matrix for M&S, which provides a set of objective descriptions for the possible states of UM practices within a given M&S effort. The columns in the uncertainty structure matrix contain UM elements or practices that are common across most M&S efforts, and the rows describe the potential levels of achievement in each of the elements. A practitioner can quickly look at the matrix to determine where an M&S effort falls based on a common set of UM practices that are described in absolute terms that can be applied to virtually any M&S effort. The matrix can also be used to plan those steps and resources that would be needed to improve the UM practices for a given M&S effort.

Green, Lawrence L.↗
DOE OSTI2025

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗
NASA NTRS2006

Porous Media Approach for Modeling Closed Cell Foam

In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is stretched to simulate the structural response of the tank during operation. The thermal expansion mismatch between the foam and the metal substrate and the thermal gradient in the foam layer causes high tensile stresses near the metal/foam interface that can lead to delamination.

Ghosn, Louis J.↗
DOE OSTI2024

Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential

Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Using Monte-Carlo Simulations to Study the Disk Structure in Cygnus X-1

As the first dynamically determined black hole X-ray binary system, Cygnus X-1 has been studied extensively. However, its broad-band spectra in hard state with BeppoSAX is still not well understood. Besides the soft excess described by the multi-color disk model (MCD), the power- law component and a broad excess feature above 10 keV (disk reflection component), there is also an additional soft component around 1 keV, whose origin is not known currently.We propose that the additional soft component is due to the thermal Comptonization process between the s oft disk photon and the warm plasma cloud just above the disk.i.e., a warm layer. We use Monte-Carlo technique t o simulate this Compton scattering process and build several table models based on our simulation results.

Yao, Y.↗
NASA NTRSDate not supplied

Experimental Observations of Mechanisms of Pattern Formation During Freeze Casting

Empirically, many different materials compositions have been freeze cast over the years. Few have been analyzed to better understand the mechanisms by which the ice crystal growth drives structure formation, particularly the substructures on the cell walls. Collating substructural features, for example, ridges, bridges, and “jellyfish-like” substructures, found on the cell walls of freeze-cast ceramics (platelets versus particles), fibrillated biopolymers (collagen, nanocellulose), and "smooth" biopolymers (chitosan, trehalose), we capture the ice crystal growth process and analyze it both in situ in 2D and "post-mortem" in 3D to gain insights into anisotropic ice crystal growth and other self-assembly mechanisms. Combining these experimental observations with 3D Phase-field simulations of structure formation in ice-templated materials, detailed in a companion presentation, it becomes possible to gain fundamental science insights and to obtain the desired structure-property-processing correlations not only at the more traditional macroscopic, but also at the meso- and microscopic scales.

Kaiyang Yin↗
NASA NTRSDate not supplied

Experimental Observations of Mechanisms of Pattern Formation During Freeze Casting

Empirically, many different materials compositions have been freeze cast over the years. Few have been analyzed to better understand the mechanisms by which the ice crystal growth drives structure formation, particularly the substructures on the cell walls. Collating substructural features, for example, ridges, bridges, and “jellyfish-like” substructures, found on the cell walls of freeze-cast ceramics (platelets versus particles), fibrillated biopolymers (collagen, nanocellulose), and "smooth" biopolymers (chitosan, trehalose), we capture the ice crystal growth process and analyze it both in situ in 2D and "post-mortem" in 3D to gain insights into anisotropic ice crystal growth and other self-assembly mechanisms. Combining these experimental observations with 3D Phase-field simulations of structure formation in ice-templated materials, detailed in a companion presentation, it becomes possible to gain fundamental science insights and to obtain the desired structure-property-processing correlations not only at the more traditional macroscopic, but also at the meso- and microscopic scales.

Kaiyan Yin↗
DOE OSTI2020

Predictive multiscale modeling for Unidirectional Carbon Fiber Reinforced Polymers

This paper presents a predictive multiscale modeling scheme for Unidirectional (UD) Carbon Fiber Reinforced Polymers (CFRP). A bottom-up modeling procedure is discussed for predicting the performance of UD structures. UD material responses are computed from high fidelity Representative Volume Elements (RVEs). A data-driven Reduced Order Modeling approach compresses RVEs into Reduced Order Models so material responses can be computed in a concurrent fashion along with the structural level simulation. The approach presented in this paper is validated against experimental data and is believed to provide design guidance for future fiber reinforced polymers development.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Contribution of Subtropical Moisture Within an Atmospheric River on Moisture Flux, Cloud Structure, and Precipitation Over the Salmon River Mountains of Idaho Using Moisture Tracers

The impact of an atmospheric river (AR) on the flux of subtropical moisture across Idaho's Salmon River Mountains and precipitation over the mountains is evaluated using the Weather, Research, and Forecasting model with water vapor tracers (WRF-WVT). The AR impacted Idaho between 17 and 19 January 2017 during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE) campaign. WRF-WVT is configured to isolate the subtropical moisture contribution to the AR, the moisture flux, and precipitation. Subtropical water vapor advected by the AR into Idaho is tagged and tracked in three-dimensional space throughout the run. This allows the contribution of the subtropical moisture to the vertical distribution of water vapor and the precipitation to be directly calculated. The simulated cloud structure is compared with airborne radar data collected during two SNOWIE intensive operation periods. This study found that more than 70% of the moisture flux and more than 80% of the precipitation across the Idaho Mountains during SNOWIE IOP 4 could be attributed to subtropical moisture within the AR. Nearly all of the moisture flux in the upper cloud and 50% of the moisture in the lower cloud was attributable to the subtropical moisture. The subtropical moisture contribution within the AR to precipitation ranged from 35% in northern Idaho to more than 90% in southern Idaho. Across the entire period of impact of the AR, more than 60% of precipitation in Idaho was attributable to the subtropical moisture within the AR, with this percentage increasing toward the south across the state.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Molecular determinants of pH sensing in the proton-activated chloride channel

In response to acidic pH, the widely expressed proton-activated chloride (PAC) channel opens and conducts anions across cellular membranes. By doing so, PAC plays an important role in both cellular physiology (endosome acidification) and diseases associated with tissue acidosis (acid-induced cell death). Despite the available structural information, how proton binding in the extracellular domain (ECD) leads to PAC channel opening remains largely unknown. Here, through comprehensive mutagenesis and electrophysiological studies, we identified several critical titratable residues, including two histidine residues (H130 and H131) and an aspartic acid residue (D269) at the distal end of the ECD, together with the previously characterized H98 at the transmembrane domain–ECD interface, as potential pH sensors for human PAC. Mutations of these residues resulted in significant changes in pH sensitivity. Some combined mutants also exhibited large basal PAC channel activities at neutral pH. By combining molecular dynamics simulations with structural and functional analysis, we further found that the β12 strand at the intersubunit interface and the associated “joint region” connecting the upper and lower ECDs allosterically regulate the proton-dependent PAC activation. Our studies suggest a distinct pH-sensing and gating mechanism of this new family of ion channels sensitive to acidic environment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗
NASA NTRS2016

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗
NASA NTRS2006

Development of Design Analysis Methods for C/SiC Composite Structures

The stress-strain behavior at room temperature and at 1100 C (2000 F) was measured for two carbon-fiber-reinforced silicon carbide (C/SiC) composite materials: a two-dimensional plain-weave quasi-isotropic laminate and a three-dimensional angle-interlock woven composite. Micromechanics-based material models were developed for predicting the response properties of these two materials. The micromechanics based material models were calibrated by correlating the predicted material property values with the measured values. Four-point beam bending sub-element specimens were fabricated with these two fiber architectures and four-point bending tests were performed at room temperature and at 1100 C. Displacements and strains were measured at various locations along the beam and recorded as a function of load magnitude. The calibrated material models were used in concert with a nonlinear finite element solution to simulate the structural response of these two materials in the four-point beam bending tests. The structural response predicted by the nonlinear analysis method compares favorably with the measured response for both materials and for both test temperatures. Results show that the material models scale up fairly well from coupon to subcomponent level.

Sullivan, Roy M.↗
NASA NTRS1996

Structure and dynamics of the coronal magnetic field

The last few years have seen a marked increase in the sophistication of models of the solar corona. This has been brought about by a confluence of three key elements. First, the collection of high-resolution observations of the Sun, both in space and time, has grown tremendously. The SOHO (Solar Heliospheric Observatory) mission is providing additional correlated high-resolution magnetic, white-light and spectroscopic observations. Second, the power and availability of supercomputers has made two- and three-dimensional modeling routine. Third, the sophistication of the models themselves, both in their geometrical realism and in the detailed physics that has been included, has improved significantly. The support from our current Space Physics Theory grant has allowed us to exploit this confluence of capabilities. We have carried out direct comparisons between observations and models of the solar corona. The agreement between simulated coronal structure and observations has verified that the models are mature enough for detailed analysis, as we will describe. The development of this capability is especially timely, since observations obtained from three space missions that are underway (Ulysses, WIND and SOHO) offer an opportunity for significant advances in our understanding of the corona and heliosphere. Through this interplay of observations and theory we can improve our understanding of the Sun. Our achievements thus far include progress modeling the large-scale structure of the solar corona, three-dimensional models of active region fields, development of emerging flux and current, formation and evolution of coronal loops, and coronal heating by current filaments.

VanHoven, Gerard↗
DOE OSTI2023

Constructing the Mechanism of Dinoflagellate Luciferin Bioluminescence Using Computation

Dinoflagellate luciferin bioluminescence is unique since it does not rely on decarboxylation but is poorly understood compared to that of firefly, bacteria, and coelenterata luciferins. Here, in this work, we computationally investigate possible protonation states, stereoisomers, a chemical mechanism, and the dynamics of the bioluminescence intermediate that is responsible for chemiexcitation. Using semiempirical dynamics, time-dependent density functional theory static calculations, and a correlation diagram, we find that the intermediate’s functional group that is likely responsible for chemiexcitation is a 4-member ring, a dioxetanol, that undergoes [2π + 2π] cycloreversion and the biolumiphore is the cleaved structure. The simulated emission spectra and luciferase-dependent absorbance spectra agree with the experimental data, giving support to our proposed mechanism and biolumiphore. We also compute circular dichroism spectra of the intermediate’s four stereoisomers to guide future experiments in differentiating them.

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

Human protein synthesis requires aminoacyl-tRNA pivoting during proofreading

Abstract Rigorous studies have characterized the aa-tRNA selection mechanism in bacteria, which is essential for maintaining translational fidelity. Recent investigations have identified critical distinctions in humans, such as the requirement of subunit rolling and a tenfold slower proofreading step. Although these studies captured key intermediates involved in tRNA selection, they did not elucidate the transitions of aa-tRNA between intermediates. Through diverse structure-based simulations, we simulated 1856 aa-tRNA accommodation events into the human ribosomal A site. Here we show the requirement for a distinct ~30° pivoting of aa-tRNA about the anticodon stem within the accommodation corridor. This pivoting is crucial for navigating the crowded accommodation corridor, which becomes more constrained due to subunit rolling. Subunit rolling-dependent crowding increases the steric contributions of the accommodation corridor during aa-tRNA accommodation, consistent with the tenfold reduction in the rate of proofreading. Furthermore, we show that eEF1A interacts with the accommodating aa-tRNA through conserved basic residues, limiting premature aa-tRNA dissociation from the A site. These findings provide a structural description of the human aa-tRNA selection process and demonstrate that the aa-tRNA alignment relative to the ribosomal catalytic sites is a critical determinant of translational fidelity.

Science & Technology - Other Topics↗