Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD 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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Diurnal and Directional Responses of Chlorophyll Fluorescence and the PRI in a Cornfield

Determining the health and vigor of vegetation using high spectral resolution remote sensing is an important goal which has application to monitoring agriculture and ecosystem productivity and carbon exchange. Two spectral indices used to assess whether vegetation is performing near-optimally or exhibiting symptoms of environmental stress (e.g., drought or nutrient deficiency, non-optimal temperatures, etc.) are the Photochemical Reflectance Index (PRI) and solar-induced red and far-red Chlorophyll Fluorescence (Fs). Both the PRI and Fs capture the dynamics of photoprotection mechanisms within green foliage: the PRI is based on the association of the reflected radiation in the green spectrum with the xanthophyll cycle, whereas Fs measures the emitted radiation in the red and far-red spectrum. Fs was determined from retrievals in the atmospheric oxygen absorption features centered at 688 and 760 nm using a modified Fraunhofer Line Depth (FLD) method. We previously demonstrated diurnal and seasonal PRI differences for sunlit vs. shaded foliage in a conifer forest canopy, as expressed in the hotspot and darkspot of the Bidirectional Reflectance Function (BRF). In a USDA-ARS experimental field site located in Beltsville, MD, USA, measurements were acquired over a corn crop from a nadir view in 2008 with an ASD FieldSpec Pro (Analytical Spectral Devices, Inc., Boulder, CO, USA) to study the behavior of the PRI for sunlit and shaded foliage as captured in reflectance variations associated with the BRF, in a I m tall canopy in the vegetative growth stage. Those observations were compared to simulations obtained from two radiative transfer models. Measurements were then acquired to examine whether the PRI and Fs were influenced by view zenith and azimuth geometries at different times of day. Those measurements were made in 2010 with the Ocean Optics USB4000 Miniature Fiber Optic Spectrometer (Ocean Optics Inc., Dunedin, Florida, USA) at several times during the day on multiple days throughout the growing season. We found that the PRI consistently had higher values, indicating lower stress, in the BRF darkspot associated with shaded foliage than in the hotspot associated with sunlit foliage. We also found that Fs exhibited differences associated with sunlit and shaded canopy sectors, which were most pronounced for the red/far-red Fs ratio. Values indicated greater physiological stress in afternoons compared to mornings, and in the early senescent canopy as compared to the vegetative growth stage, BRFs for both the PRI and the red/far-red Fs ratio were bowl-shaped for the full azimuth sweep of the canopy. These two spectral indices (PRI, Fs ratio) provided complementary information on the photosynthetic function of the corn canopy.

Middleton, Elizabeth↗
DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
DOE OSTI2021

Experimental and computational study of the microporous layer and hydrophobic treatment in the gas diffusion layer of a proton exchange membrane fuel cell

Enabling fuel cell operation at high current density is critical for developing competitive alternative power system to replace internal combustion engine. However, liquid water management continues to be a challenge for high humidity or high current density operation. Water condensation in the porous media hinders efficient oxygen transport to the catalyst layer, which in turns, reduces fuel cell performance. To improve water management capability, the gas diffusion layer is often impregnated with Polytetrafluoroethylene (PTFE) and coated with a thin microporous layer, which have shown to improve fuel cell performance, especially under wet conditions. However, the fundamental mechanism that drives the performance enhancement is still not well understood. In this work, the effects of PTFE impregnation and MPL were studied using both experimental and computational techniques. Both limiting current and polarization tests under dry and wet operating condition are conducted to study the oxygen transport resistance and fuel cell performance. In addition, a 2-D, two-phase, multi-physics PEMFC model is developed to simulate performance and gain a fundamental understanding of local water saturation and oxygen concentration. The combined results show that 5 wt% PTFE impregnation with MPL significantly enhances liquid water management, which enables higher current density operation of a fuel cell.

25 ENERGY STORAGE↗
DOE OSTI2023

A Reconfigurable Planar Capacitor Sensor for Thru the Jacket Dielectric Constant Measurement of Cable Insulation

Planar capacitor sensors have been proposed to measure the aging-related dielectric constant variation of the jacket or insulation materials of cables. Although capacitor sensors have been used to monitor the aging of cable insulation by removing the jacket, through the jacket measurement of cable insulation dielectric constant variation is difficult because once the sensor is placed on the jacket it would either measure the aging of the jacket or that of the jacket and insulation together. A reconfigurable sensor along with innovative cable-specific models could be used to determine the jacket and insulation aging using the same sensor. For the proposed novel reconfigurable sensor, the electric field penetration depth is changed by activating and deactivating PIN diode switches. In one instance, the sensor measures the dielectric constant variation of the jacket, while in the next, it measures the dielectric constant variation of both the jacket and the insulation. Here, by leveraging previously developed capacitance estimation models from large-scale finite element simulations for a specific cable, these two sets of measurement data are used to evaluate the aging-related dielectric constant variation of the insulation. Experimental data from aged cable specimens are presented to demonstrate the operation of the sensor.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Machine learning based prediction of airflow maldistribution in air-to-refrigerant heat exchangers

Flow maldistribution is a common challenge in heat exchanger (HX) design and particularly important for air-to-refrigerant geometries where capacity losses can approach 65%. This has a major impact on central air conditioning systems, as compact duct design motivates the use of A-type HXs which are known to be affected by airflow maldistribution. Because velocity profiles are difficult to predict, components are often oversized leading to increased material cost, system footprint, and refrigerant charge. Several studies detail airflow maldistribution for individual HXs and packages, but findings cannot always be extrapolated to new designs. In this work, a machine learning (ML) based flow profile prediction framework is developed and applied to two common package configurations: (i) A-type and (ii) U-type HXs, across a broad range of HX geometries and flow rates. Porous media CFD simulations are validated against independent data for both package types as well as comprehensive in house measurements for a finless geometry with shape optimized non-round tubes, which validates the framework for new heat transfer surfaces. The ML models are trained on the porous media CFD simulations, predicting volumetric flow rate (VFR) within 1.1% and 1.9% with maximum relative L 2 norm errors of 0.48 and 0.65, respectively, while also delivering 10 5 speed up factor compared to full porous media CFD. HX level simulations show an up to 9% reduction in heat transfer from flow maldistribution, with greater losses occurring at smaller half apex angles. This framework enables rapid and highly accurate prediction of airflow maldistribution induced capacity degradation.

42 ENGINEERING↗
DOE OSTI2025

Investigation of residual stress distribution in wire-arc directed energy deposited refractory molybdenum alloy utilizing numerical thermo-mechanical analysis and neutron diffraction method

Directed energy deposition (DED), a metal additive manufacturing (AM) technique, offers higher deposition rates and energy efficiency, making it suitable for fabricating components from refractory molybdenum alloys, such as molybdenum-titanium-zirconium (TZM). However, large thermal gradients and non-equilibrium thermal cycles in DED could generate high residual stress in the component, potentially deteriorating quality and performance. Thus, this study aims to investigate residual stress generation and its distribution in wire-arc DED of TZM thin-wall, utilizing thermo-mechanical analysis and high-fidelity neutron diffraction (ND) method. Two interpass temperatures (50°C and 200°C) have been considered to investigate their impact on residual stress formation. During experiments, in-situ thermal data has been recorded using thermocouples, which have been utilized for calibrating the thermal model. Thermocouple data shows a good agreement with the simulation results, having a difference of less than 10 %. Post-deposition part deformation has been observed, which is measured using a coordinate measuring machine, showing maximum values of 0.93 mm and 0.78 mm for interpass temperatures of 50°C and 200°C, respectively. Numerical predictions of distortion deviated by less than 15% from the experimental results. ND measurement and simulation results indicate that residual stress magnitude and evolution vary across the TZM deposits, revealing microstructural anisotropy in both conditions. Notably, lower interpass temperatures resulted in higher residual stresses, confirmed by experimental and simulation data. Further, this study demonstrated that an integrated experimental and thermo-mechanical analysis can potentially reveal the temperature history, part deformation, and residual stress formation in wire-arc DED TZM alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Powder metallurgy hot isostatic pressing - workshop report

The findings of the powder metallurgy (PM) hot isostatic pressing (HIP) workshop confirmed that the accelerated adoption of this technology requires the following action items: 1. By 2035, the HIP infrastructure must be scaled to sizes larger than 4 .0 meters for manufacturing components relevant to energy and defense industries at large. 2. The community must develop cost-effective and agile standards, qualifications, and certification pathways for these HIP parts. 3. At the same time, domestic powder production should be expanded to a wide range of metallic alloys, including Cr-Mo steels, stainless steels, and titanium- and nickel-based superalloys, as well as to emerging alloy classes in the refractory space. 4. The above initiatives can be enhanced by implementing digital modeling and controls, including AI-driven simulations, in situ monitoring, and digital twins.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

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

All-Atom Molecular Dynamics Simulations of the Temperature Response of Densely Grafted Polyelectrolyte Brushes

Water acts a “good” solvent for most polyelectrolyte (PE) systems due to its high dielectric strength. Theoretical and experimental studies have highlighted some remarkably unconventional properties of strongly water-swollen PE brushes such as a reduction in brush height with increase in temperature. However, an understanding of the atomistic-scale response of the grafted PE chains as well as the local arrangement and organization of the brush-supported counterions and water molecules to changes in temperature is missing. Here, we conduct an all-atom molecular dynamics (MD) study to probe the influence of temperature on the “microstructure” of densely grafted, fully ionized polyacrylic acid (PAA) brushes. The atomistic insights obtained from our study shed light on prior experimental observations and elucidate the key role played by changes in hydrogen bonding network in dictating the hydrophilicity of the brushes. In addition to the temperature-mediated reduction in height of the PE brush layer, we study the influence of temperature on several properties of the brush-supported water molecules such as the number distribution and kinetics of water-water and water-PE hydrogen bonds (HBs), mass density, orientational tetrahedral order parameter, and the activation energy associated with the self-diffusion of water. The effect of temperature on the properties of the brush-trapped counterions is also quantified via changes in their solvation structure, counterion-water radial distribution function (RDF), and the translational mobility. Our findings unravel several hitherto unknown phenomena such as the enhanced propensity for attaining PE brush-induced “water-in-salt”-like scenarios at elevated temperatures, and an increase in the activation energy for the self-diffusion of water with increase in the degree of PE brush-induced nanoconfinement (from ‘bulk water’ to ‘sparsely grafted PE brushes’ to ‘densely grafted PE brushes’). Further, we anticipate that this work will be helpful in improving our predictive capabilities of the temperature response of water-swollen PE brushes, which would be instrumental for better design of several PE brush systems used in a multitude of applications.

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

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA’s diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Diamond under extremes

Diamond is, by virtue of the covalent bonding between atoms and the very strong carbon to carbon bonds, the hardest natural material. It has been a fascinating material since its discovery, first as a decorative gem and more recently, for its numerous industrial uses because of its extreme hardness, elastic modulus, and optical transparency. In recent years, it has become a preferred ablator for laser shock experiments, and this has led to its choice as the capsule material for fusion experiments at the National Ignition Facility. Further, this review covers both experimental and computational (including machine learning) advancements in research on diamond subjected extreme conditions of temperature and pressure. The synergy between shock and ramp loading experiments and atomic level simulations is proving to be powerful in advancing our understanding of diamond under extremes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Janssen

Differentiable and Modular Optical Microscopy Simulations and Inversions

Mukherjee, Debangshu [Oak Ridge National Laborator↗
DOE OSTI2022

Building energy efficiency and load flexibility optimization using phase change materials under futuristic grid scenario

This work investigates the energy savings and load flexibility capacity of phase change material (PCM) integrated into building envelope under a future energy generation scenario, where 80% of the energy load comes from renewable sources. Using community-scale modeling, we first determine the net thermal storage required to fully manage the demand variability at a utility grid, and then using whole-building simulations we optimize operating parameters - such as PCM thickness, latent heat, interior setpoint profile, PCM distribution, and heat transfer coefficient - to determine the optimal conditions required for maximum energy reduction and load flexibility without compromising occupants' thermal comfort. The optimal PCM-integrated envelope proposed in this study can provide annual load flexibility up to 33.6% and annual energy savings up to 10.8% in a lightweight residential building located in Baltimore, MD. This study is relevant given the increasing contributions of renewable energy in the total energy generation mix, leading to a significant time-imbalance between peak energy demand and peak energy production. While thermal energy storage using PCM is a recognized technique, there is no known prior study dedicated to examining the thermal performance of PCM-integrated envelope under future energy generation scenarios. This study bridges the research gap by investigating a practical way to implement PCMs in the buildings and maximize the energy efficiency as well as load flexibility related benefits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Computational modeling of graphite degradation in molten salt reactors: Role of infiltration

Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Outstanding Questions and Future Research on Magnetic Reconnection

This short article highlights unsolved problems of magnetic reconnection in collisionless plasma. Advanced in-situ plasma measurements and simulations have enabled scientists to gain a novel understanding of magnetic reconnection. Nevertheless, outstanding questions remain concerning the complex dynamics and structures in the diffusion region, cross-scale and regional couplings, the onset of magnetic reconnection, and the details of particle energization. We discuss future directions for magnetic reconnection research, including new observations, new simulations, and interdisciplinary approaches.

79 ASTRONOMY AND ASTROPHYSICS↗