Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multi-material simulations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗
DOE OSTI2023

Coupling hydrodynamics to a rigid-body motion solver for fluid-structure interaction [Slides]

Quinoa is a massively parallel computational fluid dynamics with multi-material and programmed burn capabilities developed from Programmatic and LDRD funds. Overset is a well-established method of using an “inset” mesh that communicates with a background mesh. Solutions transfer freely from one to another and operate as boundary conditions on the opposite mesh. Either mesh can move at specified velocity! This project couples the existing Mesh-to-mesh transfer with Quinoa and demonstrates use of resulting Overset solver in computational simulation of blast effects on a re-entry body for survivability assessments.

42 ENGINEERING↗
DOE OSTI2026

Performance assessment of 3D printed multi-material energy absorber for automotive bumper: pedestrian lower extremity protection

Designing an energy absorber for automotive bumpers involves balancing low-speed and high-speed impacts to ensure safety, reduce repair costs, and meet regulatory standards. Here, this study explores a novel design using multi-material 3D printing and structural optimization to fabricate a lightweight and cost-efficient energy absorber. The design effectively dissipates energy in low-speed collisions and minimizes force transmission in high-speed pedestrain impacts, helping to meet both safety and performance requirements. The energy absorber design combines 20% carbon fiber-reinforced acrylonitrile butadiene styrene (CF-ABS) and thermoplastic polyurethane (TPU) for optimal stiffness and flexibility. It uses 3D-printed lattice structures optimized through finite element simulations to help meet both low-speed and high-speed impact requirements. Full-scale energy absorbers were 3D-printed using optimized CF-ABS/TPU blends and tested under high-speed impact using the Flexible Pedestrian Legform Impactor (Flex-PLI). For fair comparison, a baseline bumper with a traditional triangular lattice structure, also 3D-printed from the same CF-ABS/TPU materials, was similarly tested. Interestingly, both the optimized and baseline 3D-printed energy absorbers showed nearly identical performance, successfully meeting injury limits. Their performances were also benchmarked against an injection-molded energy absorber. While both 3D-printed and injection-molded designs met injury limits, the 3D-printed absorber exhibited a higher tibia bending moment, indicating an opportunity for further optimization. A Techno-Economic Analysis compared the costs of producing energy absorbers using traditional manufacturing and 3D printing. The analysis highlighted that 3D printing offers cost benefits for low to medium production volumes, with the total cost per energy absorber at ∼ $\$$74, compared to traditional methods that become economical beyond 2000 units.

Additive manufacturing↗
DOE OSTI2023

Reducing warpage in a hybrid large-scale additive manufacturing and compression molding process

In recent years, a hybrid manufacturing process, developed by combining extrusion-based large-scale additive manufacturing (AM) and compression molding (CM) techniques, has shown promising outcomes for producing structurally functional parts. The process can be used with both short fiber-reinforced composites and neat polymers and hence, even multi-material parts can be manufactured easily. This process offers the advantages of structural enhancement by having a desired fiber orientation using a large-scale AM process, as well as rapid manufacturing capability using a CM process. In the large-scale AM process, the alignment of fibers in the deposition direction enables significant improvement in the mechanical properties of the manufactured parts. However, the anisotropy resulting from the directional arrangement of fibers also introduces challenges related to warpage in the produced parts. This study aims to identify the causes of warpage and propose strategies to mitigate it. The research involves the use of preforms manufactured through the large-scale AM, which are then combined with the neat resin for CM manufacturing. A finite element-based numerical simulation model is developed, employing a sequentially coupled thermomechanical approach. Through a parametric study using the simulation models, optimization of printing direction and preform geometry is performed to minimize warpage. This contributes to the advancement and wider adoption of AM/CM hybrid manufacturing to produce structurally functional parts.

Jo, Eonyeon↗
DOE OSTI2024

A machine-learning-aided data recovery approach for predicting multi-material thermal behaviors in advanced test reactor capsules

Instrumented experiments conducted at test reactors are essential to the deployment of new advanced reactor systems. Designing new experiments and generating data on specific reactor conditions require significant investments in terms of both time and cost. Finite element analysis software can be used to create high-fidelity models of experiment environments in order to support the actual experiments, but computation time remains a concern in terms of applying outcomes to real-time usage of data (e.g., a digital twin [DT]). Here, the present research proposes a machine-learning (ML) aided approach to making temperature and displacement predictions based on the thickness of the outer gas gap on the experimental capsule used for in-pile demonstration of a novel new thermal conductivity probe in the Advanced Test Reactor (ATR). This capsule consisted of U10Zr fuel, a rodlet, sodium, and inner and outer capsules. Gas gaps existed between the fuel and the rodlet, and between the inner and the outer capsule. The learning data pertained to an experimental capsule's radial distributions of temperature and displacement, as obtained based on Abaqus and the physical features. For the first step of ML sequence, the temperature was predicted using three positional parameters. Next, the displacement was predicted using seven additional parameters. Each physical feature was normalized in order to be both nondimensional and standardized. The temperature and displacement predictions showed good agreement with the simulation results in all cases involving interpolation and extrapolation. Furthermore, data similarity enhancement increased the similarity between the training and the target data, thereby increasing the predictive accuracy of the ML models. In certain extrapolation cases involving limited original ML model accuracy, data similarity enhancement and data recovery was able to somewhat improve this accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Multi-Objective design of interlocking metasurfaces using conditional diffusion models

Unit cell design remains a major challenge for interlocking metasurfaces, a promising joining technology for dissimilar materials, due to the complex, competing, multivariate design space and the need for rapid adaptation to varying performance requirements. This study explores Conditional Diffusion Models as a design optimization tool for interlocking metasurfaces. Given the complex, competing, multivariate design space for interlocking metasurfaces, unit cell design remains a major challenge for this joining technology. We trained a conditional diffusion model on 25,000 finite element analysis-simulated interlocking metasurface unit cells to generate designs with tailored thermo-mechanical properties (tensile strength, shear strength, and thermal conductivity) based on specified performance criteria. The model demonstrated a success rate of approximately 72 % in producing designs that met specified property bounds. The conditional diffusion model generated both thermally resistive and conductive designs, revealing clear trends in design characteristics: taller, dendritic structures were advantageous for tensile loads, while shorter, robust designs excelled in shear applications. Our findings indicate that the model's performance is more influenced by the breadth of the design space than by the quantity of training data, highlighting the importance of expansive design domains for generating innovative solutions. This work establishes conditional diffusion models as a highly efficient and adaptable tool for rapid interlocking metasurface unit cell design, paving the way for advancements in multi-material joining technologies, as well as highlighting the justification to leverage conditional diffusion models as design tools across complex design domains.

Conditional diffusion models↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2024

A Shock Stabilization of the HLLC Riemann Solver for the Carbuncle Instability

The HLLC approximate Riemann solver improves upon the HLL Riemann solver by resolving contact discontinuities. This is a particularly desirable property for multi-material codes in which problems usually contain material interfaces. However, the HLLC solver is known to suffer from the carbuncle phenomenon, a numerical instability most apparent at grid-aligned shocks in multi-dimensional simulations. Many problems of interest, including high energy-density physics applications, require the accurate resolution of both material interfaces and hydrodynamic shocks. A variety of methods have been developed to cure this instability, with varying degrees of complexity. The objective of this work is to describe a simple approach to modify the HLLC Riemann solver and prevent the carbuncle instability. Here, the method is then demonstrated for assorted two-dimensional test problems known to exhibit the shock instability. The performance of the new solver is compared with that of the standard HLL and HLLC Riemann solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A System-Level Model for Estimating Residual Strain and Life of Nuclear Reactor Coolant System Components Under Connected-System-Thermal–Mechanical Boundary Conditions

Background: Environmental-assisted fatigue (EAF) is a major issue for the long-term survival of nuclear power plant fleets in the U.S. and worldwide. Multi-material welded regions (e.g., nozzles) and other high-stress regions of reactor coolant system (RCS) components are prone to EAF-related damage. Objective: The discussed work describes a system-level finite element (FE) model of RCS components of a pressurized water reactor (PWR). This is with the goal of predicting the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Methods: Herein, the FE model was developed considering system-level loading conditions (under connected system thermal-mechanical boundary conditions). Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the FE model results, the strain residuals, strain amplitudes and resulting fatigue lives of RCS components were predicted. Results: The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. Higher residual strain can lead to accelerated cyclic hardening of material and the associated effect of EAF. The simulated component-level strain profile (under realistic multi-axial-multi-physics loading cycle) can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests, which is a focus of future works. Conclusions: Despite similar geometry and material the RCS component can have significantly different strain profiles and resulting fatigue lives. The discussed approach can help to identify and prioritize the RCS components for conducting expensive nondestructive evaluation (NDE) inspections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

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

Conservative remapping of material-dependent fields between possibly misaligned material regions

In this work, we propose an interpolation or remapping algorithm of material-dependent fields on polyhedral meshes where any source or target cell contains only one material. It is conservative and it preserves sharp material boundaries on the target mesh, even if the source and target regions delineating the same material are slightly misaligned. If those material regions are aligned, then the algorithm is also linearity-preserving and bounds-preserving. For a given material, it consists of a conservative field reconstruction on a target mesh part from a source mesh part associated with that material, followed by a repair step in case of misaligned boundaries. No assumption is made regarding the topology of the input meshes

36 MATERIALS SCIENCE↗
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 OSTI2022

Multi-material swept face remapping on polyhedral meshes

Remapping is a conservative interpolation of a discretized intensive quantity between two meshes. In this article, we propose a novel multi-material flux remapping method that avoids the geometric computation of mesh-mesh intersections needed for an accurate intersection based remap. The flux remap is applicable to scalar quantities such as material density describing the multi-material flow between meshes with the same connectivity but small mesh displacements. Herein, the method is described for two- and three-dimensional polygonal/polyhedral meshes as it is implemented in Portage. Another open source library, Tangram, is used to calculate material interfaces in cells containing more than one material. Performance and accuracy of the flux remap are discussed with respect to Arbitrary Lagrangian-Eulerian simulations and compared to an accurate intersection based remap. In particular, cyclic remapping shows that the accuracy of the flux remap is limited to first order on material boundaries while maintaining second order accuracy in pure material regions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Property optimized energy absorber for automotive bumpers utilizing multi-material and structural design strategies

This study proposes a novel design for automotive bumper using optimized lattice structures and multi-materials to balance low-speed collision and high-speed pedestrian impact performance. Different blends of 20 % carbon fiber-reinforced acrylonitrile butadiene styrene with thermoplastic polyurethane were used to tailor material properties. The energy absorber features lattice structures with customized mechanical responses, created by varying the incline angle θ from 0 to 180°. We conducted 576 finite element simulations on a half-scale model to optimize energy absorption and stiffness, leading to 66 optimized designs that met both low-speed and high-speed impact criteria. Two sub-scale optimized energy absorbers with different peak forces—both meeting low-speed impact requirements—were 3D printed and validated through drop-weight testing. The one with lower peak stress demonstrated a more compliant response, exhibiting approximately 90 % lower initial peak force and an increase in energy absorption of around 33 % (from 24 J to 32 J). Compared to the baseline triangular lattice, the optimized absorber increased energy absorption by 68 % from (19 J to 32 J) and reduced peak stress by 70 %. It also showed near-complete recovery with minimal fractures, making it suitable for repeated use. This design improves safety while offering a lightweight, durable, and cost-effective bumper system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ultra-light antennas via charge programmed deposition additive manufacturing

Abstract The demand for lightweight antennas in 5 G/6 G communication, wearables, and aerospace applications is rapidly growing. However, standard manufacturing techniques are limited in structural complexity and easy integration of multiple material classes. Here we introduce charge programmed multi-material additive manufacturing platform, offering unparalleled flexibility in antenna design and the capability for rapid printing of intricate antenna structures that are unprecedented or necessitate a series of fabrication routes. Demonstrating its potential, we present a transmitarray antenna composed of an interconnected, multi-layered array of dielectric/conductive S-ring unit cells, reducing 94% mass of conventional antenna configurations. A fully printed circular polarized transmitarray system fed by a source and a Risley prism antenna system operating at 19 GHz both show close alignment between testing results and numerical simulations. This printing method establishes a universal platform, propelling discovery of new antenna designs and enabling data-driven design and optimizations where rapid production of antenna designs is crucial.

Science & Technology - Other Topics↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Role of mechanical stress localizations on the radiation hardness of AlGaN/GaN high electron mobility transistors

Multi-material, multi-layered systems such as AlGaN/GaN high electron mobility transistors (HEMTs) contain residual mechanical stresses that arise from sharp contrasts in device geometry and materials parameters. These stresses, which can be either tensile or compressive, are difficult to detect and eliminate because of their highly localized nature. We propose that their high-stored internal energy makes potential sites for defect nucleation sites under radiation, particularly if their locations coincide with the electrically sensitive regions of a transistor. In this study, we validate this hypothesis with molecular dynamic simulation and experiments exposing both pristine and annealed HEMTS to 2.8 MeV Au +3 irradiation. Our unique annealing process uses mechanical momentum of electrons, also known as the electron wind force (EWF) to mitigate the residual stress at room temperature. High-resolution transmission electron microscopy and cathodoluminescence spectra reveal the reduction of point defects and dislocations near the two-dimensional electron gas region of EWF-treated devices compared to pristine devices. The EWF-treated HEMTs showed relatively higher resilience with approximately 10% less degradation of drain saturation current and ON-resistance and 5% less degradation of peak transconductance. Both mobility and carrier concentration of the EWF-treated devices were less impacted compared to the pristine devices. Our results suggest that the lower density of nanoscale stress localization contributed to the improved radiation tolerance of the EWF-treated devices. Intriguingly, the EWF is found to modulate the defect distribution by moving the defects to electrically less sensitive regions in the form of dislocation networks, which act as sinks for the radiation induced defects and this assisted faster dynamic annealing.

AlGaN/GaN HEMTs↗
DOE OSTI2022

Revisiting Multi-Material Composite Structures with Homogenized Composite Properties

Composite structures inherently develop residual stresses during their curing process. Driven predominately by mismatched thermal strains between differing materials or ply orientations, but also affected by curing process phenomena like polymer shrinkage, these residual stresses can lead to failure within composite structures. There are several methods varying in complexity that can be used to model the development of residual stresses, all of which are capable of capturing sufficient detail to understand the residual stress state at the ply level. However, explicitly modeling all plies of a layup in a composite structure can be prohibitively expensive based on the number of plies, structure size, and required element size. The computational cost can be reduced through the homogenization of the composite layup without losing much fidelity of the overall response of the structure. The homogenization process reduces the many plies of a laminate to a single lamina that reduces complexity and increases the mesh size where a single element can span multiple plies. This report focuses on verification and validation efforts for a homogenization process using a suite of finite element simulations rather than an analytic solution derived from classical laminate theory. Initial verification using representative element volumes indicated there was minimal error in the homogenization process; however, this compounded to a small, but acceptable error in strip and split ring experimental composite structures. The error does under predict the residual stress state in the strip and split ring and should be accounted for when simulating composite structures with homogenized properties.

36 MATERIALS SCIENCE↗