SEARCH · Search NASA

Results for “Interactive materials model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Theory of thermal properties of magnetic materials with unknown entropy

Theoretical approaches to study thermal properties of magnetic materials typically require accurate models of magnetic interactions in order to define the entropy. Here we introduce a complementary approach for examining thermal properties in magnetic systems where an accepted model for such interactions does not exist. In place of a specific model for magnetic interactions, the approach integrates measurements of temperature dependent magnetization of the studied material into a first principles computational scheme. The approach calculates system pressure from thermally disordered microstates that properly incorporate vibrational and spin subsystems at each temperature as well as the coupling between these subsystems. We apply the approach to calculate phonon modes and to investigate the anomalously low thermal expansion of the classical Invar alloy, Fe 0.65 Ni 0.35 . Here, the calculated phonon dispersions for Invar are in excellent agreement with measured data. The Invar thermal expansion is shown to remain small between 50 K and room temperature, consistent with the experimentally observed low thermal expansion value in this same temperature range. This anomalously small thermal expansion is directly connected to a small positive contribution from lattice thermal disorder that is nearly canceled by a large negative magnetic disorder contribution. By contrast, calculations for bcc Fe show a much larger thermal expansion, consistent with experiment, which is dominated by a large contribution from lattice thermal disorder that is reduced only slightly by a small negative contribution from that of magnetism. These findings give insights into the unusual nature of magnetism and spin-lattice coupling in Invar and Fe. In addition, they give promising preliminary support to the presented new methodology as a complementary way to investigate thermal properties of magnetic materials. The success achieved on Invar and Fe motivates future testing of the approach on other magnetic materials.

36 MATERIALS SCIENCE↗

NASA NTRS2001

Nonlinear Ultrasonic Characterization of Fatigue Microstructures

Dislocation dipole substructures formed during metal fatigue are shown to produce a substantial distortion of ultrasonic waves propagating through the fatigued material. A model of ultrasonic wave-dislocation dipole interactions is developed that quantifies the wave distortion by means of a material nonlinearity parameter (beta). Application of the model to AA2024-T4 predicts a value of p approximately 300% larger in material cyclically loaded for 100 kcycles in stress-control at 276 MPa and R=0 than that measured for virgin material. Experimental measurements show a monotonic increase in p as a function of the number of fatigue cycles that closely approaches the predicted increase. The experiments also suggest that the relevant dislocation substructures are localized in the material.

Cantrell, John H.↗

NASA NTRS2011

Modeling Electrostatic Fields Generated by Internal Charging of Materials in Space Radiation Environments

Internal charging is a risk to spacecraft in energetic electron environments. DICTAT, NU MIT computational codes are the most widely used engineering tools for evaluating internal charging of insulator materials exposed to these environments. Engineering tools are designed for rapid evaluation of ESD threats, but there is a need for more physics based models for investigating the science of materials interactions with energetic electron environments. Current tools are limited by the physics included in the models and ease of user implementation .... additional development work is needed to improve models.

Minow, Joseph I.↗

DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Status of Mercury and Imp: Two Monte Carlo Transport Codes Developed Using Shared Infrastructure at Lawrence Livermore National Laboratory

The Monte Carlo Transport Project at Lawrence Livermore National Laboratory develops two Monte Carlo transport codes used in production by a sizable internal user community. Mercury is a Monte Carlo particle transport code used to model the interaction of neutrons, gammas, and light ions with a material. Imp is an implicit Monte Carlo thermal x-ray photon transport code used to model the interaction of x-ray photons with a material. This paper describes the two codes and highlights recent developments.

Pozulp, Michael↗

DOE OSTI2020

An interface-aware sub-scale dynamics multi-material cell model for solids with void closure and opening at all speeds

Here, we present a multi-material cell model (closure model) for demanding arbitrary Lagrangian-Eulerian (ALE) simulations of fluids and solids. It is based on the interface-aware sub-scale dynamics (IASSD) approach which utilizes the exact material interface geometry within the computational cell to calculate internal material interactions. Our formulation of the closure model also aims to improve the accuracy in low-speed impact events. Voids are used to represent ambient vacuum and internal free boundaries of the distinct materials. Void regions can close and open at contact surfaces, allowing a transition from contact physics to free motion in vacuum. The coupling of void closure and opening with a new formulation of the IASSD model for solids is tested on several one- and two-dimensional numerical examples, ranging from gas expansion in vacuum to planar and round object impacts at various speeds.

42 ENGINEERING↗

DOE OSTI2022

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗

NASA NTRS1990

Scattering properties of woven fibrous insulations - Effects of interference for normal incidence

This paper examines the scattering properties of woven fibrous materials in order to model the interactions between the scattered radiation from individual fibers. A normally incident plane electromagnetic wave is considered and different representative geometries are analyzed theoretically. Experimental results are presented for specific woven thermal insulation fabrics. Previous studies in the literature have mostly ignored the effect of constructive or destructive addition of scattered waves from individual fibers with the exception of the case of parallel fibers. Whereas the light-scattering characteristics of a single long fiber depend only on the size parameter and on the material optical constants, the additional mechanism of interference between the scattered waves from fibers of parallel or perpendicular orientation introduces an additional dependence on fiber locations. In this study, analytical models are developed for obtaining the radiative scattering characteristics of woven fibrous media taking into account the mechanism of interference. Experimental measurements of forward- and backward-scattered light from a woven silica cloth over a broad range of wavelengths in the infrared are obtained and compared with the predicted results for a system of two orthogonal families of infinite parallel fibers.

Kumar, Sunil↗

DOE OSTI2025

Higher-order factorization machine for accurate surrogate modeling in material design

Efficient and robust optimization is important in material science for identifying optimal structural parameters and enhancing material performance. Surrogate-based active learning algorithms have recently gained great attention for their ability to efficiently navigate large, high-dimensional design spaces. Among surrogate models, 2 nd -order factorization machine (FM) models are widely employed as the surrogate model in active learning algorithms due to their balance between simplicity and effectiveness. However, their quadratic nature limits their capacity to capture complex, higher-order interactions among variables, often leading to suboptimal solutions. To overcome this limitation, we propose an active learning scheme integrating a 3 rd -order FM model, capable of modeling three-variable interactions and more intricate relationships in material systems. We comprehensively evaluate the surrogate modeling performance of the 3 rd -order FM case using various objective functions. Furthermore, we examine the optimization reliability and efficiency of the 3 rd -order FM-based active learning in a real-world material design task (e.g., nanophotonic structures for transparent radiative cooling). Our study shows that the 3 rd -order FM outperforms the 2 nd -order model in both surrogate accuracy and optimization performance, highlighting higher-order models’ promises for material design and optimization problems.

Factorization machine↗

DOE OSTI2025

A GPU ‐Accelerated 3D Unstructured Mesh Based Particle Tracking Code for Multi‐Species Impurity Transport Simulation in Fusion Tokamaks

ABSTRACT This paper presents the multi‐species global impurity transport capability developed in a GPU‐accelerated fully 3D unstructured mesh‐based code, GITRm, to simultaneously track multiple impurity species and handle interactions of these impurities with mixed‐material surfaces. Different computational approaches to model particle‐surface interaction or surface response have been developed and compared. Sheath electric field is taken into account by employing a fast distance‐to‐boundary calculation, which is carried out in parallel on distributed or partitioned meshes on multiple GPUs without the need for any inter‐process communication during the simulation. Several example cases, including two for the DIII‐D tokamak, that is, one with the SAS‐V divertor and the other with the collector probes, are used to demonstrate the utility of the current multi‐species capability. For the DIII‐D probe case, the capability of GITRm to resolve the spatial distribution of particles in localized regions, such as diagnostic probes, within non‐axisymmetric tokamak geometries is demonstrated. These simulations involve up to 320 million particles and utilize up to 48 GPUs.

Nath, Dhyanjyoti D. [Scientific Computation Resear↗

DOE OSTI2021

Tuning magnetic antiskyrmion stability in tetragonal inverse Heusler alloys

The identification of materials supporting complex, tunable magnetic order at ambient temperatures is foundational to the development of new magnetic device architectures. In this work, we report the design of Mn 2 XY tetragonal inverse Heusler alloys that are capable of hosting magnetic antiskyrmions whose stability is sensitive to elastic strain. Here We first construct a universal magnetic Hamiltonian capturing the short- and long-range magnetic order which can be expected in these materials. This model reveals critical combinations of magnetic interactions that are necessary to approach a magnetic phase boundary, where the magnetic structure is highly susceptible to small perturbations such as elastic strain. We then computationally search for quaternary Mn 2 (X 1 , X 2 ) Y alloys where these critical interactions may be realized and which are likely to be synthesizable in the inverse Heusler structure. We identify the Mn 2 Pt 1–z X z Ga family of materials with X = Au, Ir, Ni as an ideal system for accessing all possible magnetic phases, with several critical compositions where magnetic phase transitions may be actuated mechanically.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Gutzwiller hybrid quantum-classical computing approach for correlated materials

Rapid progress in noisy intermediate-scale quantum (NISQ) computing technology has led to the development of novel resource-efficient hybrid quantum-classical algorithms, such as the variational quantum eigensolver (VQE), that can address open challenges in quantum chemistry, physics, and material science. Proof-of-principle quantum chemistry simulations for small molecules have been demonstrated on NISQ devices. While several approaches have been theoretically proposed for correlated materials, NISQ simulations of interacting periodic models on current quantum devices have not yet been demonstrated. Here, we develop a hybrid quantum-classical simulation framework for correlated electron systems based on the Gutzwiller variational embedding approach. We implement this framework on Rigetti quantum processing units (QPUs) and apply it to the periodic Anderson model, which describes a correlated heavy electron band hybridizing with noninteracting conduction electrons. Our simulation results quantitatively reproduce the known ground state quantum phase diagram including metallic, Kondo and Mott insulating phases. This is the first fully self-consistent hybrid quantum-classical simulation of an infinite correlated lattice model executed on QPUs, demonstrating that the Gutzwiller hybrid quantum-classical embedding framework is a powerful approach to simulate correlated materials on NISQ hardware. This benchmark study also puts forth a concrete pathway towards practical quantum advantage on NISQ devices.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Cascade of phase transitions in a planar Dirac material

Here we investigate a model of interacting Dirac fermions in 2 + 1 dimensions with M flavors and N colors having the U(M)×SU(N ) symmetry. In the large-N limit, we find that the U(M) symmetry is spontaneously broken in a variety of ways. In the vacuum, when the parity-breaking flavor-singlet mass is varied, the ground state undergoes a sequence of M first-order phase transitions, experiencing M + 1 phases characterized by symmetry breaking U(M)→U(M – k)×U(k) with k ϵ {0, 1, 2, · · · , M}, bearing a close resemblance to the vacuum structure of three-dimensional QCD. At finite temperature and chemical potential, a rich phase diagram with first and second-order phase transitions and tricritical points is observed. Also exotic phases with spontaneous symmetry breaking of the form as U(3)→U(1) 3 , U(4)→U(2)×U(1) 2 , and U(5)→U(2)2×U(1) exist. For a large flavor-singlet mass, the increase of the chemical potential μ brings about M consecutive first-order transitions that separate the low-μ phase diagram with vanishing fermion density from the high-μ region with a high fermion density.

1/N expansion↗

DOE OSTI2025

Characterization of ablator dynamics initiated by picket-pulse conditions

Direct-drive laser-fusion targets may include low-cost plastic ablator material. Here, this work describes a tabletop experimental system designed to enable the investigation of the interaction of such ablator material as a function of the laser wavelength and the associated excitation mechanism. This experimental system offers time-resolved imaging of the material response with subpicosecond temporal resolution and ≈1-μm spatial resolution. The system involves synchronized femtosecond and picosecond lasers to initiate and characterize the dynamics of polystyrene planar targets irradiated by conditions similar to a “picket” prepulse from a direct-drive pulse shape. Ultrafast imaging of the initial excitation process followed by the resulting critical plasma formation, shock-wave propagation, and blowoff plasma expansion is used to compare the impact of excitation wavelength, 355 nm or 266 nm, on the ensuing dynamics. This work aims to investigate fundamental laser material interactions to help advance physics models that can inform the development of laser driver and target designs.

42 ENGINEERING↗

NASA NTRS2009

Simulation of RCC Crack Growth Due to Carbon Oxidation in High-Temperature Gas Environments

The carbon wall oxidation technique coupled with a CFD technique was employed to study the flow in the expanding crack channel caused by the oxidation of the channel carbon walls. The recessing 3D surface morphing procedure was developed and tested in comparison with the arcjet experimental results. The multi-block structured adaptive meshing was used to model the computational domain changes due to the wall recession. Wall regression rates for a reinforced carbon-carbon (RCC) samples, that were tested in a high enthalpy arcjet environment, were computationally obtained and used to assess the channel expansion. The test geometry and flow conditions render the flow regime as the transitional to continuum, therefore Navier-Stokes gas dynamic approach with the temperature jump and velocity slip correction to the boundary conditions was used. The modeled mechanism for wall material loss was atomic oxygen reaction with bare carbon. The predicted channel growth was found to agree with arcjet observations. Local gas flow field results were found to affect the oxidation rate in a manner that cannot be predicted by previous mass loss correlations. The method holds promise for future modeling of materials gas-dynamic interactions for hypersonic flight.

Titov, E. V.↗

DOE OSTI2021

Progressive growth of the solid–electrolyte interphase towards the Si anode interior causes capacity fading

The solid–electrolyte interphase (SEI), a layer formed on the electrode surface, is essential for electrochemical reactions in batteries and critically governs the battery stability. Active materials, especially those with extremely high energy density, such as silicon (Si), often inevitably undergo a large volume swing upon ion insertion and extraction, raising a critical question as to how the SEI interactively responds to and evolves with the material and consequently controls the cycling stability of the battery. Here, by integrating sensitive elemental tomography, an advanced algorithm and cryogenic scanning transmission electron microscopy, we unveil, in three dimensions, a correlated structural and chemical evolution of Si and SEI. Furthermore, corroborated with a chemomechanical model, we demonstrate progressive electrolyte permeation and SEI growth along the percolation channel of the nanovoids due to vacancy injection and condensation during the delithiation process. Consequently, the Si–SEI spatial configuration evolves from the classic ‘core–shell’ structure in the first few cycles to a ‘plum-pudding’ structure following extended cycling, featuring the engulfing of Si domains by the SEI, which leads to the disruption of electron conduction pathways and formation of dead Si, contributing to capacity loss. The spatially coupled interactive evolution model of SEI and active materials, in principle, applies to a broad class of high-capacity electrode materials, leading to a critical insight for remedying the fading of high-capacity electrodes.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Energy-dependent implementation of secondary electron emission models in continuum kinetic sheath simulations

The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗