Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials modelling”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The role of alkali metal exchange in zeolite-templated carbon synthesis

Despite great efforts to achieve ideal atomistic packing of carbon in the pore networks of even the largest pore zeolites, templating inaccuracies are ubiquitous in all presently synthesized zeolite-templated carbon (ZTC) materials. Idealized models suggest that the long-sought schwarzite-like structures could in principle be obtained by zeolite templating if the appropriate zeolite template were chosen (e.g., faujasite for D surface schwarzites) and if perfect template fidelity (insertion of a pristine layer of pure carbon directly on the surface of the zeolite) could be achieved. A requirement to achieve such structures is increased carbon density within the zeolite. We report the investigation of a series of alkali metal cation-exchanged zeolites to determine how the periodic trends in the group 1 elements influence zeolite templating, with a specific focus on the metric of structural packing density (SPD) as resolved by ex situ thermogravimetry. In a survey based on controlled synthesis temperature, time, and flow conditions, an increasing SPD was observed with decreasing cation size, an effect that is consistent with the increasing strength of cation-π interactions. Furthermore, this effect could be promising for future work to increase the SPD of ZTCs for the synthesis of closed-tube, schwarzite-like carbonaceous solids.

03 NATURAL GAS↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
NASA NTRS2012

Lattice Thermal Conductivity from Atomistic Simulations: ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
NASA NTRS2012

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
DOE OSTI2025

Possible Spin-Triplet Excitonic Insulator in the Ultraquantum Limit of HfTe 5

More than 50 years ago, excitonic insulators formed by the pairing of electrons and holes due to Coulomb interactions were first predicted [A. N. Kozlov and L. A. Maksimov, Sov. J. Exp. Theor. Phys. 21, 790 (1965); L. V. Keldysh and Y. V. Kopaev, Sov. Phys. Solid State 6, 2219 (1965); D. Jérome, T. M. Rice, and W. Kohn, Phys. Rev. 158, 462 (1967)]. Since then, excitonic insulators have been observed in various classes of materials, including quantum Hall bilayers, graphite, transition metal chalcogenides, and more recently in moiré superlattices. In these excitonic insulators, an electron and a hole with the same spin bind together, and the resulting exciton is a spin singlet. Here, we report the experimental observation of a spin-triplet excitonic insulator in the ultra-quantum limit of a three-dimensional topological material HfTe 5 . We observe that the spin-polarized zeroth Landau bands dispersing along the field direction cross each other beyond a characteristic magnetic field in HfTe 5 , forming the one-dimensional Weyl mode. Transport measurements reveal the emergence of a gap of about 2⁢5⁢0 μ⁢eV when the field surpasses a critical threshold. By performing the material-specific modeling, we identify this gap as a consequence of a spin-triplet exciton formation, where electrons and holes with opposite spin form bound states, and the translational symmetry is preserved. The system reaches charge neutrality following the gap opening, as evidenced by the zero Hall conductivity over a wide magnetic field range (10–72 T). In conclusion, our finding of the spin-triplet excitonic insulator paves the way for studying novel spin transport including spin superfluidity, spin Josephson currents, and Coulomb drag of spin currents in analogy to the transport properties associated with the layer pseudospin in quantum Hall bilayers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling Wavelength Dependent Mid-infrared (5.5–25 μm) Optical Constants of Silicate Glasses: A Genetic Algorithm Approach

Wavelength-dependent mid-infrared (400–1,800 cm -1 ; 5.56–25 μm) optical constants (real and imaginary indices of refraction; n and k) are determined using reflectance spectra at a spectral sampling of 4 cm -1 for several silicate glasses of varying SiO 2 wt% which include (a) basaltic volcanic glass from Kilauea, Hawaii, (b) synthetic andesite, (c) two synthetic dacites, (d) obsidian volcanic glass from Mount Lassen, California, (e) synthetic rhyodacite, and (f) rhyolitic volcanic glass from Mexico. Because glasses are optically isotropic, no specific orientation was required for spectral measurements, and polished glass samples were measured at random orientations using micro-FTIR spectrometer. Lorentz-Lorenz dispersion theory and Fresnel reflectance model for high symmetry materials were used to model the optical constants by optimizing oscillator parameters in modeled spectra to match laboratory spectra. A genetic algorithm (GA) approach automatically finds the natural oscillators and their parameters for each spectrum, and then uses these parameters in a non-linear least squares optimization routine. The study compared spectral parameters such as Christiansen feature, reststrahlen bands (RBs), and the peaks centered around 860–1,100 (peak 1) and 400–480 cm -1 (peak 2) of n and k to their respective SiO 2 wt% of the glasses. CF, RBs, peak 1 of n and k, and peak 2 of n shift to higher wavenumbers with increased SiO 2 wt%, whereas peak 2 of k shifts to lower wavenumbers with increased SiO 2 wt%. Derived optical constants of these glasses will improve quantitative abundance mapping of volcanic materials on the surfaces of silicate targets in the Solar System.

58 GEOSCIENCES↗
DOE OSTI2023

Hydrogen Storage System Modeling: Public Access, Maintenance, and Enhancements

The hydrogen storage system modeling project aims to develop new and enhance existing material-based hydrogen storage system models and make them accessible to the research community through a public web page. Updating existing storage systems models supports material developers in evaluating the performance of their new materials in mobile and stationary applications relative to the available DOE Technical Targets. Developing new tools helps researchers evaluate the performance of hydrogen storage materials developed under HyMARC activities or other fundamental hydrogen storage materials discovery research. New models will also work for validating alternatives to material-based systems (liquefied & gaseous hydrogen), comparing various mobile and stationary use cases, and expand to include medium-, heavy-duty, and mining vehicles and stationary application(s).

AMR↗
DOE OSTI2024

Leveraging Curvature on N–Doped Carbon Materials for Hydrogen Storage

Carbon sorbent materials have shown great promise for solid-state hydrogen (H 2 ) storage. Modification of these materials with nitrogen (N) dopants has been undertaken to develop materials that can store H 2 at ambient temperatures. In this work density functional theory (DFT) calculations are used to systematically probe the influence of curvature on the stability and activity of undoped and N-doped carbon materials toward H binding. Specifically, four models of carbon materials are used: graphene, [5,5] carbon nanotube, [5,5] D 5d -C 120 , and C 60 , to extract and correlate the thermodynamic properties of active sites with varying degrees of sp 2 hybridization (curvature). From the calculations and analysis, it is found that graphitic N-doping is thermodynamically favored on more pyramidal sites with increased curvature. In contrast, it is found that the hydrogen binding energy is weakly affected by curvature and is dominated by electronic effects induced by N-doping. These findings highlight the importance of modulating the heteroatom doping configuration and the lattice topology when developing materials for H 2 storage.

08 HYDROGEN↗
DOE OSTI2022

The effect of Cr alloying on defect migration at Ni grain boundaries

Mass transport along grain boundaries in alloys depends not only on the atomic structure of the boundary, but also its chemical make-up. In this work, we use molecular dynamics to examine the effect of Cr alloying on interstitial and vacancy-mediated transport at a variety of grain boundaries in Ni. We find that, in general, Cr tends to reduce the rate of mass transport, an effect which is greatest for interstitials at pure tilt boundaries. However, there are special scenarios in which it can greatly enhance atomic mobility. Cr tends to migrate faster than Ni, though again this depends on the structure of the grain boundary. Further, grain boundary mobility, which is sometimes pronounced for pure Ni grain boundaries, is eliminated on the time scales of our simulations when Cr is present. We conclude that the enhanced transport and grain boundary mobility often seen in this system in experimental studies is the result of non-equilibrium effects and is not intrinsic to the alloyed grain boundary. These results provide new insight into the role of grain boundary alloying on transport that can help in the interpretation of experimental results and the development of predictive models of materials evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2021

High-efficiency perovskite photovoltaic modules achieved via cesium doping

Perovskite solar modules have been attracting increasing attention due to their market potential, yet publications concerned with theintrinsic scale-up potential of different perovskite compositions remain relatively scarce. On the other hand, while great success is being made towards improving the power conversion efficiency (PCE) of perovskite solar cells (PSCs) by cesium cation (Cs + ) doping of the perovskite, more attention is being paid to the perovskite phase stabilization effect of Cs + doping, and less to other properties that are critical to understand and futher improve the PSC's. In this work, moderately-Cs-doped MAPbI 3 was employed as a model perovskite material in order to exclude the phase stabilization effect. Our systematic study revealed the influence of Cs + in organic-inorganic hybrid perovskites on the crystal structure, crystallization process, trap state density, band structure and charge (i.e., ions or photo-carriers) transport. Markedly, it has been observed that Cs + doping can greatly increase the carrier diffusion length in the perovskite films, thus improving the potential to scale-up PSC's.The PCE of small area devices (0.09 cm 2 ) was increased to 21.72% from 19.73%, with decreased hysteresis behavior and increased operational stability (T85 = 1000 h) after Cs + doped, where T85 refers to the retention of 85% of the initial PCE. Moreover, a PCE of 21.08% was obtained for a Cs + -containing perovskite module with an active area > 30 cm 2 , which demonstrates a better "reproducibility" than the reference sample (MAPbI3-based perovskite modules, PCE = 18.26%).

14 SOLAR ENERGY↗
DOE OSTI2024

Delivering laser performance conditions to enable fusion ignition, and beyond at the National Ignition Facility

On December 5th, 2022, controlled fusion ignition was demonstrated for the first time at the National Ignition Facility (NIF), a major achievement in the field of Inertial Confinement Fusion (ICF) requiring a multi-decadal effort involving broad national and international collaborations. To drive the fusion ignition reaction with the compressed fuel capsule, that yielded 3.15 MJ of nuclear energy [1], the NIF laser delivered a high-precision pulse shape with 2.05 MJ of ultra-violet (UV) laser energy and a peak power of 440 TW. This laser energy was an increase of ~8 % compared to that delivered on the previous “threshold of ignition” record yield experiment (1.37 MJ of yield for 1.89 MJ of laser energy) on August 8th, 2021. Here we explain how the results of our extensive research in laser technology and UV optics damage mitigation led to major improvements in the NIF laser, enabling this energy increase along with additional accuracy, precision, and power balance enhancements. Furthermore, we will discuss on-going efforts that have enabled operations at 2.2 MJ of UV energy as well as potential new initiatives to push the laser performance –accuracy and delivered energy– to even higher levels in the future as previously demonstrated on a small subset of NIF beams.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗