Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Materials 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,711 records · Page 95

A discrete-continuum coupled finite element modelling approach for fibre reinforced concrete

Highlights: • An easy-to-implement but effective meso-scale modelling approach is developed for FRC. • It is capable of simulating and interpreting all main failure mechanisms in FRC. • It is critically validated by experiments with single or multiple random fibres under various loading conditions. • It can be used to optimise key design factors of FRC at meso/material scale by parametric studies. - Abstract: Fibre reinforced concrete (FRC) exhibits complicated failure modes such as fibre breakage, mortar cracking, crushing and spalling, fibre-mortar interfacial debonding, depending on many material properties, geometric dimensions, boundary and loading conditions. Most existing numerical models are unable to reproduce these failure modes that may occur simultaneously or sequentially in a specimen, mainly due to difficulties in generating finite element meshes with a large number of randomly-oriented fibres. Herein we develop a discrete-continuum coupled finite element modelling approach for FRC materials capable of effectively simulating all the major failure modes. The continuum damaged plasticity model is used to simulate damage and fracture behaviour of the mortar, while debonding of fibre-mortar interfaces is modelled by nonlinear cohesive interfacial elements. Unique techniques are devised to generate conforming meshes between fibres and the surrounding mortar so that the randomly-oriented fibres are easily modelled. The modelling approach is validated by simulating single fibre pullout tests with different inclination angles, notched and non-notched direct tensile tests and three-point bending beam tests with randomly-distributed multiple fibres.

36 MATERIALS SCIENCE↗

Tunable Nanoscale Evolution and Topological Phase Transitions of a Polar Vortex Supercrystal

Here, understanding the phase transitions and domain evolutions of mesoscale topological structures in ferroic materials is critical to realizing their potential applications in next-generation high-performance storage devices. Here, the behaviors of a mesoscale supercrystal are studied with 3D nanoscale periodicity and rotational topology phases in a PbTiO 3 /SrTiO 3 (PTO/STO) superlattice under thermal and electrical stimuli using a combination of phase-field simulations and X-ray diffraction experiments. A phase diagram of temperature versus polar state is constructed, showing the formation of the supercrystal from a mixed vortex and α-twin state and a temperature-dependent erasing process of a supercrystal returning to a classical α-twin structure. Under an in-plane electric field bias at room temperature, the vortex topology of the supercrystal irreversibly transforms to a new type of stripe-like supercrystal. Under an out-of-plane electric field, the vortices inside the supercrystal undergo a topological phase transition to polar skyrmions. These results demonstrate the potential for the on-demand manipulation of polar topology and transformations in supercrystals using electric fields. The findings provide a theoretical understanding that may be utilized to guide the design and control of mesoscale polar structures and to explore novel polar structures in other systems and their topological nature.

P=phase-field simulations↗

Ingrained: An Automated Framework for Fusing Atomic-Scale Image Simulations into Experiments

To fully leverage the power of image simulation to corroborate and explain patterns and structures in atomic resolution microscopy, an initial correspondence between the simulation and experimental image must be established at the outset of further high accuracy simulations or calculations. Furthermore, if simulation is to be used in context of highly automated processes or high-throughput optimization, the process of finding this correspondence itself must be automated. In this work, "ingrained," an open-source automation framework which solves for this correspondence and fuses atomic resolution image simulations into the experimental images to which they correspond, is introduced. Here, the overall "ingrained" workflow, focusing on its application to interface structure approximations, and the development of an experimentally rationalized forward model for scanning tunneling microscopy simulation are described.

36 MATERIALS SCIENCE↗

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗

Dual energy CT: How well can pseudo-monochromatic imaging reduce metal artifacts?

Purpose: Dual Energy CT (DECT) provides so-called monoenergetic images based on a linear combination of the original polychromatic images. At certain patient-specific energy levels, corresponding to certain patient- and slice-dependent linear combination weights, e.g., E = 160 keV corresponds to α = 1.57, a significant reduction of metal artifacts may be observed. The authors aimed at analyzing the method for its artifact reduction capabilities to identify its limitations. The results are compared with raw data-based processing. Methods: Clinical DECT uses a simplified version of monochromatic imaging by linearly combining the low and the high kV images and by assigning an energy to that linear combination. Those pseudo-monochromatic images can be used by radiologists to obtain images with reduced metal artifacts. The authors analyzed the underlying physics and carried out a series expansion of the polychromatic attenuation equations. The resulting nonlinear terms are responsible for the artifacts, but they are not linearly related between the low and the high kV scan: A linear combination of both images cannot eliminate the nonlinearities, it can only reduce their impact. Scattered radiation yields additional noncanceling nonlinearities. This method is compared to raw data-based artifact correction methods. To quantify the artifact reduction potential of pseudo-monochromatic images, they simulated the FORBILD abdomen phantom with metal implants, and they assessed patient data sets of a clinical dual source CT system (100, 140 kV Sn) containing artifacts induced by a highly concentrated contrast agent bolus and by metal. In each case, they manually selected an optimal α and compared it to a raw data-based material decomposition in case of simulation, to raw data-based material decomposition of inconsistent rays in case of the patient data set containing contrast agent, and to the frequency split normalized metal artifact reduction in case of the metal implant. For each case, the contrast-to-noise ratio (CNR) was assessed. Results: In the simulation, the pseudo-monochromatic images yielded acceptable artifact reduction results. However, the CNR in the artifact-reduced images was more than 60% lower than in the original polychromatic images. In contrast, the raw data-based material decomposition did not significantly reduce the CNR in the virtual monochromatic images. Regarding the patient data with beam hardening artifacts and with metal artifacts from small implants the pseudo-monochromatic method was able to reduce the artifacts, again with the downside of a significant CNR reduction. More intense metal artifacts, e.g., as those caused by an artificial hip joint, could not be suppressed. Conclusions: Pseudo-monochromatic imaging is able to reduce beam hardening, scatter, and metal artifacts in some cases but it cannot remove them. In all cases, the CNR is significantly reduced, thereby rendering the method questionable, unless special post processing algorithms are implemented to restore the high CNR from the original images (e.g., by using a frequency split technique). Raw data-based dual energy decomposition methods should be preferred, in particular, because the CNR penalty is almost negligible.

60 APPLIED LIFE SCIENCES↗

Quantitative measurements of dislocations in metals for advancing predictive simulations

LLNL applications require scientists to predict how materials evolve under various thermomechanical conditions. While this is achieved through physics-based simulations, uncertainty in the predictions of mechanical properties remains a serious challenge that limits the predictive capabilities of models because we lack methods to compare predictions of atomic-scale defects (dislocations) with experimental measurements. High energy X-ray diffraction (HEXRD) is the most relevant technique that can provide the necessary statistical information on dislocations. However, this technique is not yet quantitative because we lack a precise understanding of the relationship between X-ray diffraction patterns and the underlying material dislocation content and arrangements. To address this need, we used our novel computational X-ray diffraction method to simulate the effect of dislocations on the diffraction patterns. We compared virtual and experimental diffraction patterns. Results allowed us to clearly establish the relationship between X-ray diffraction patterns and the underlying dislocation structures, proving that it is feasible to quantitatively measure dislocation statistics with HEXRD. This project delivered a method that can provide the missing piece to LLNL’s mechanical property simulations in advanced metals by obtaining experimentally long-needed quantitative dislocation data, which could fully enable predictive capabilities.

36 MATERIALS SCIENCE↗

2026 Annual Molten Salt Reactor Campaign Review

This report documents the 2026 Annual Molten Salt Reactor (MSR) Campaign Review held in Albuquerque, New Mexico, from April 21–24, 2026. The review was organized by Dr. Patricia Paviet, National Technical Director of the Advanced Reactor Technology (ART) MSR program for the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), and included participation from principal investigators, federal managers, developers, regulators, and members of the broader MSR community. The 2026 review expanded beyond the ART-MSR campaign to include related DOE-NE programs supporting molten salt reactor advancement, including the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the Advanced Materials and Manufacturing Technologies (AMMT) program, Advanced Reactor Safeguards and Security (ARSS), Material Protection, Accounting and Control Technologies (MPACT), and relevant advanced fuels activities. Three days were dedicated to technical presentations and panel discussions covering thermal properties, off-gas management, modeling and simulation, safety and licensing, safeguards and security, materials and corrosion, and irradiation activities. A fourth day was dedicated to technical tours of Sandia National Laboratories and Kairos Power facilities. Attendance was strong and comparable to the prior annual review, with approximately 75 in-person participants per day, 75–80 virtual attendees per day, and approximately 30 participants in the tour day. The review fostered significant technical exchange across national laboratories, universities, industry, regulators, and international participants. This report summarizes the review structure, technical themes, participation, tours, feedback, and conclusions relevant to future planning for the MSR campaign.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Determining reference standard strength for neutron-irradiated reduced activation ferritic/martensitic steel F82H by Bayesian method

The deterministic approach widely adopted in the design of structural components relies on systematically defined design limits using empirically determined safety factors. However, this approach is not always appropriate because structures are subjected to a variety of loads in the practical environment, which may result in excessively conservative design limits. In recent years, a more rigorous probabilistic approach that incorporates material strength distributions has become an important solution. In the probabilistic approach, the probability density functions of material strength properties underpin the design criteria. Here, the objective of this study is to identify the density distribution functions that best describe tensile properties of irradiated F82H to define a reference strength for DEMO design. Due to the limited number of existing data, this study specifically employs a Bayesian prediction method based on Monte Carlo simulations to determine a material reference value with statistical reliability and to investigate its effectiveness. For example, the dependence of tensile properties of 300 °C irradiated materials on irradiation damage and the range predicted by 95% Bayesian estimation was evaluated. As a statistical model for the dose dependence of statistical parameters, the normal distribution exhibited a better fit for 0.2% proof strength and tensile strength, whereas the distribution of total elongation data gave comparable reference values for both the normal and Weibull distribution models. Both models gave comparable criteria for the distribution of total elongation data. The Weibull model also gave better results for uniform elongation. The function best describing the model was a logarithmic law for both 0.2% proof strength and tensile strength, while a power law for both total and uniform elongation, which allowed for more comprehensive data prediction of irradiation data with statistical accuracy for DEMO reactor design.

36 MATERIALS SCIENCE↗

Lost Circulation Materials in Geothermal Drilling: Thermal Degradation, Viscosity, Compression, and Flow Loop Tests

This dataset provides a set of experimental data on the performance of lost circulation materials (LCMs) used in geothermal drilling operations. It includes results from four distinct tests: thermal degradation, viscosity measurements, compression tests, and high-temperature flow loop experiments. The materials tested, such as cedar fiber, cotton seed hulls, magma fiber, microcellulose, and sawdust, were analyzed under a range of conditions to assess their effectiveness in mitigating fluid losses in geothermal wells. Data on thermal degradation explores mass loss, water loss, and gas release for each material when subjected to temperatures between 90C and 250C. Viscosity measurements capture the behavior of LCM mixtures under different temperatures, while compression tests investigate the structural integrity of both undegraded and thermally degraded materials. The high-temperature flow loop experiments simulate geothermal well conditions, analyzing how different materials and mixtures affect fluid flow and sealing performance in a fractured reservoir environment. In addition to the experimental data, the dataset includes photos documenting the materials, schematics of test setups, and several published manuscripts that provide detailed insights into the experimental methods and findings.

15 GEOTHERMAL ENERGY↗

Development of New Lunar Highland Regolith Simulant, NUW-LHT-5M

Introduction: NASA has a need for large quantities of lunar simulants that closely match future manned lunar missions at landing sites near the Lunar South Pole. The current simulant was designed to approximate NU-LHT-2M, and -4M, except using natural minerals and a fully synthetic, non-basaltic, high calcium glass. Washington Mills in Niagara Falls, NY was contracted to produce this new simulant due to their electric arc furnace (EAF) fusion technology and capabilities for crushing and sizing ceramics. The primary objectives of this work were: (1) to create a commercially produced simulant that contains a high fidelity, high calcium (An90+), low Mg/Fe glass without basaltic constituents; (2) to create a simulant that uses existing, terrestrially mined, and readily available minerals blended with a synthetic glass, and (3) accomplish objectives (1) and (2) using conventional processing methods capable of production levels that can meet substantial current and future demands of NASA projects. Another consideration was that highland simulants closest to meeting these requirements, NU-LHT-2M and -4M types, are out of production and largely out of stock. Design: The design chosen for development and production for this work is NU-LHT-2M, de-signed by Doug Stoeser and Douglas Rickman, which uses Stillwater anorthosite and norite rocks and synthetic glass. The composition targeted the average composition, glass content and particle size distribution of Apollo 16 samples [1]. The design utilizes two Stillwater minerals, anorthosite (37.7 wt.%), norite (17.6 wt.%), a commercially sourced olivine (4.7 wt.%), and a synthetic glass (40 wt.%). Processing: Stillwater rocks were hand collected and initially crushed by the USGS. Further crushing and milling was performed at Washing-ton Mills. For the glass, oxides were blended in the ap-propriate ratios, blended in a v-blender, and placed in a graphite lined, water cooled pot. The mixture was fused in Washington Mills’ pilot scale, 500 kW EAF. The molten glass was then poured into water and quenched to inhibit crystallization. Several pour/quench cycles were required to produce the total quantity of glass needed. Anorthosite, norite and olivine were milled together to the target particle size in a ball mill. The glass was ball milled separately. The complete simulant was produced by blending the milled materials. Particle size distribution, phase and chemical analyses of the materials was were per-formed at Washington Mills. Results: The glass composition was found to match the Apollo 16 oxide content average [1] quite well, with the exception of SiO2 and Fe2O3. The deviation of SiO2 and Fe2O3 from batched content was believed to be due to limitations inherent in the use of the arc furnace technique for these types of materials. A particle size distribution of the full simulant formulation closely matching the average for the Apollo 16 samples was also achieved. Fig. 1 shows the morphology of the final simulant particles. The size distribution is compared to the target in Fig. 2. Conclusions: A new lunar highlands type simulant containing a high-quality synthetic glass, closely matching an average composition of the Apollo 16 regolith was produced. Processing utilized easily scalable mineral and glass making methods.

simulant↗

Molecular Simulation Results on Charged Carbon Nanotube Forest-Based Supercapacitors

In this paper, electrochemical double-layer capacitances of charged carbon nanotube (CNT) forests with tetraethyl ammonium tetrafluoro borate electrolyte in propylene carbonate are studied on the basis of molecular dynamics simulation. Direct molecular simulation of the filling of pore spaces of the forest is feasible even with realistic, small CNT spacings. The numerical solution of the Poisson equation based on the extracted average charge densities then yields a regular experimental dependence on the width of the pore spaces, in contrast to the anomalous pattern observed in experiments on other carbon materials and also in simulations on planar slot-like pores. The capacitances obtained have realistic magnitudes but are insensitive to electric potential differences between the electrodes in this model. This agrees with previous calculations on CNT forest supercapacitors, but not with experiments which have suggested electrochemical doping for these systems. Finally, those phenomena remain for further theory/modeling work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystallization of the P 3 Sn 4 Phase upon Cooling P 2 Sn 5 Liquid by Molecular Dynamics Simulation Using a Machine Learning Interatomic Potential

We performed molecular dynamics simulations to study the crystallization of the P3Sn4 phase from P 2 Sn 5 liquid using a machine learning (ML) interatomic potential with desirable efficiency and accuracy. Our results capture the liquid properties of P 2 Sn 5 at 1300 K, which is well above the melting temperature. The phase separation and crystallization are observed when P 2 Sn 5 liquid is cooled down below 832 and 505 K, respectively. The simulation results are in good agreement with the experimentally observed phase transformation behaviors and provide useful insights into the complex nucleation and crystallization process at the details of atomistic scale. Our work also demonstrated that ML interatomic potentials based on neural network deep learning are robust and capable of accurately describing the energetics and kinetics of complex materials through molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Processing of Lunar Soil Simulant for Space Exploration Applications

NASA's long-term vision for space exploration includes developing human habitats and conducting scientific investigations on planetary bodies, especially on Moon and Mars. To reduce the level of up-mass processing and utilization of planetary in-situ resources is recognized as an important element of this vision. Within this scope and context, we have undertaken a general effort aimed primarily at extracting and refining metals, developing glass, glass-ceramic, or traditional ceramic type materials using lunar soil simulants. In this paper we will present preliminary results on our effort on carbothermal reduction of oxides for elemental extraction and zone refining for obtaining high purity metals. In additions we will demonstrate the possibility of developing glasses from lunar soil simulant for fixing nuclear waste from potential nuclear power generators on planetary bodies. Compositional analysis, x-ray diffraction patterns and differential thermal analysis of processed samples will be presented.

Sen, Subhayu↗

Atomistic characterization of SAM coatings as gate insulators in Si-based FET devices

Many nano-material systems are currently under consideration as possible candidates for gate dielectric insulators in both metal-oxide-semiconductor (MOSFET) and organic (OFET) field-effect transistors. In this contribution, the possibility of employing self-assembled monolayers (SAMs) of hydroxylated octadecyltrichlorosilane (OTS) chains on a (111) Si substrate as gate dielectrics is discussed; in particular ab initio theoretical simulations have been employed to study the structural properties, work function modifications, and the insulating properties of OTS thin film coatings on Si substrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sequential lateral solidification of silicon thin films on low-k dielectrics for low temperature integration

We present the excimer laser crystallization of amorphous silicon on a low dielectric constant (low-k) insulator for very large scale integration monolithic 3D integration and demonstrate that low dielectric constant materials are suitable substrates for 3D integration through laser crystallization of silicon thin films. We crystallized 100 nm amorphous silicon on top of SiO{sub 2} and SiCOH (low-k) dielectrics, at different material thicknesses (1 μm, 0.75 μm, and 0.5 μm). The amorphous silicon crystallization on low-k dielectric requires 35% less laser energy than on an SiO{sub 2} dielectric. This difference is related to the thermal conductivity of the two materials, in agreement with one dimensional simulations of the crystallization process. We analyzed the morphology of the material through defect-enhanced microscopy, Raman spectroscopy, and X-ray diffraction analysis. SEM micrographs show that polycrystalline silicon is characterized by micron-long grains with an average width of 543 nm for the SiO{sub 2} sample and 570 nm for the low-k samples. Comparison of the Raman spectra does not show any major difference in film quality for the two different dielectrics, and polycrystalline silicon peaks are closely placed around 517 cm{sup −1}. From X-ray diffraction analysis, the material crystallized on SiO{sub 2} shows a preferential (111) crystal orientation. In the SiCOH case, the 111 peak strength decreases dramatically and samples do not show preferential crystal orientation. A 1D finite element method simulation of the crystallization process on a back end of line structure shows that copper (Cu) damascene interconnects reach a temperature of 70 °C or lower with a 0.5 μm dielectric layer between the Cu and the molten Si layer, a favorable condition for monolithic 3D integration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling of Dual Gate Material Hetero-dielectric Strained PNPN TFET for Improved ON Current

The tunnel field effect transistor (TFET) is considered to be a promising alternative device for future low-power VLSI circuits due to its steep subthreshold slope, low leakage current and its efficient performance at low supply voltage. However, the main challenging issue associated with realizing TFET for wide scale applications is its low ON current. To overcome this, a dual gate material with the concept of dielectric engineering has been incorporated into conventional TFET structure to tune the tunneling width at source–channel interface allowing significant flow of carriers. In addition to this, N+ pocket is implanted at source–channel junction of the proposed structure and the effect of strain is added for exploring the performance of the model in nanoscale regime. All these added features upgrade the device characteristics leading to higher ON current, low leakage and low threshold voltage. The present work derives the surface potential, electric field expression and drain current by solving 2D Poisson’s equation at different boundary conditions. A comparative analysis of proposed model with conventional TFET has been done to establish the superiority of the proposed structure. All analytical results have been compared with the results obtained in SILVACO ATLAS device simulator to establish the accuracy of the derived analytical model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simulation of Creep Deformation and Failure in Graded AM Microstructures

This report describes modeling tools and techniques developed to simulate the long-term material performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF). A physics-based Crystal Plasticity Finite Element model is used to simulate creep in microstructures and to study the roles of grain morphology, porosity, and texture. We describe our modeling methodology, including an orientation-mapping technique to capture the spatially varying crystallographic orientation that results from the build conditions. Our study of microstructural features shows that AM microstructures produced by LPBF tend to creep faster in the build direction, while texture and grain boundaries strengthen the transverse directions. However, when grain-boundary porosity and the consequent cavity growth are included in the model, the transverse directions begin to creep faster. In examining texture, the results indicate that spatially varying orientation arising from the build conditions increases anisotropy in the material, making it critical to account for orientation gradients in the material to accurately model its mechanical behavior. We also describe a material-model calibration campaign in which we calibrated the constitutive model specifically for LPBF 316H stainless steel at 725℃ for both solution-annealed and as-built conditions. Finally, these tools and techniques are used to model creep in microstructures representing different regions of an LPBF material with graded microstructure, owing to intentional variation in processing conditions. The creep simulation results show good agreement with experimental data across all three microstructures, with future work planned to study rupture in the material.

36 MATERIALS SCIENCE↗