Fast three-dimensional rules-based simulation of thermal-sprayed microstructures
Not Available
SEARCH · Engineering Papers
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.
Not Available
We present metamaterial-like aerogels that provide broadband low-frequency vibration mitigation through their self-assembled locally resonant microstructures. Simulations confirm the presence of a negative dynamic stiffness for the aerogels.
FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.
FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.
Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.
The concrete biological shield of light water reactors is exposed to high neutron and gamma irradiation doses in the long term. Irradiation deteriorates the physical and mechanical properties of concrete. Such effects need to be investigated to predict the concrete’s performance in the event of a lifetime extension of a nuclear power plant. This work combines high-resolution characterization techniques with fast-Fourier transform (FFT)-based 2-D simulations to evaluate the radiation-induced volumetric expansion (RIVE) and damage in concrete microstructures under neutron irradiation. Two concrete microstructures from samples provided by the Japan Concrete Aging Management Program (JCAMP) were characterized using micro x-ray fluorescence (mXRF) to obtain elemental intensity maps, and energy-dispersive x-ray spectroscopy to complement mXRF with local elemental information for Na. Minerals and cement paste are then identified based on the elemental composition to produce high-resolution phase maps, resulting in a more accurate representation of the microstructures compared to previous work. Simulations of radiation-induced volumetric expansion (RIVE), creep, and damage in JCAMP concrete use the fast Fourier transform (FFT)-based code Microstructure Oriented Scientific Analysis of Irradiated Concrete (MOSAIC) combined with the irradiated minerals, aggregates, and concrete (IMAC) database, which contains mineral-specific RIVE models. Overall, the simulation results are in fair agreement with experimental data.
Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.
Not Available
During a loss-of-coolant accident (LOCA) in a nuclear power plant, a rapid increase in temperature is experienced in the fuel. In commercial light water reactors (LWRs), the high burnup structure (HBS) forms in regions of the UO 2 fuel where local burnup is high, characterized by smaller grain sizes and large, overpressurized bubbles. During a LOCA transient, fuel in the HBS region is susceptible to fine fragmentation, where the fuel breaks up into micron-size fragments (sometimes referred to as pulverization). To better understand the mechanisms behind this phenomenon, mesoscale computational modeling using the Idaho National Laboratory code Marmot has been employed to simulate the process of HBS formation and response to a LOCA transient. The process of complete and partial HBS formation was demonstrated in prototypical LWR conditions. The simulated HBS microstructures from these simulations were passed to a phase-field model of fracture. To provide a more direct comparison with experiments, a set of experiments conducted by Studsvik on fine fragmentation/pulverization during LOCA conditions was simulated using the BISON fuel performance code. The temperature profiles from the BISON simulations were passed to a phase-field model to determine the gas bubble pressure as a function of time. These pressure histories were used in a phase-field model of fracture to determine under what conditions and in which types of bubbles fine fragmentation was likely to occur.
This is my paper for the 2025 AIAA SciTech conference. The submission is due 12/09.
This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.
When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure–property map of the stochastic volume elements making up the statistical description of the part. In this paradigm the autoencoder can be used to directly encode the microstructure or, alternatively, its latent space can be sampled to provide likely realizations. Furthermore, we demonstrate the method on three examples using the common additively manufactured material AlSi10Mg in: (a) a comparison with direct numerical simulation of the part microstructure, (b) a push forward of microstructural uncertainty to performance quantities of interest, and (c) a simulation of functional gradation of a part with stochastic microstructure.
Concrete is a complex, heterogeneous material with a microstructure composed of aggregates, cement paste, and pores spanning multiple length scales. Understanding this microstructure is critical for advancing the performance, durability, and modeling of concrete-based systems. While experimental imaging such as X-ray computed tomography (XCT) provides valuable insights, generating large datasets with detailed ground truth annotations is both costly and labor-intensive due to challenges in segmenting similar phases, such as aggregates and cement paste, that often share similar attenuation properties. To address this, we developed a pipeline to simulate realistic 3D concrete microstructures using the open-source Python package PyCMG. This simulation effort focuses on generating high-fidelity, annotated microstructures that can serve as training or benchmarking datasets for image analysis, segmentation algorithms, and machine learning models, particularly in scenarios where experimental data is scarce.
Highlights: • The GeoMicro3D model was innovatively developed for simulating the reaction process and microstructure formation of AAS. • In GeoMicro3D, the dissolution of slag and reactions of aqueous ions are fundamentally described in the local lattice cell. • A novel approach is proposed to speed up the nucleation simulations. • GeoMicro3D was implemented and verified with the relevant experimental data and thermodynamic calculation results. For the first time, this study developed a novel model, named GeoMicro3D, to simulate the reaction process and microstructure formation of alkali-activated slag. The GeoMicro3D model consists of four modules that are designed to simulate, respectively: (i) the initial spatial distribution of real-shape slag particles in alkaline activator, (ii) the dissolution of slag and diffusion of ions via the transition state theory and lattice Boltzmann method, respectively, (iii) the spatial distribution of reaction products using a nucleation probability theory, and (iv) the chemical reactions with thermodynamic modelling. Afterwards the GeoMicro3D model was implemented and verified. The simulation results were discussed and compared with the relevant experimental data and thermodynamic calculation results using GEMS. A good agreement was found in the comparisons, showing the strong simulation capability of GeoMicro3D.
The diffraction pattern of a material contains information not only on the crystal structures of its constituting phases, but also on its mesoscale spatial distributions of phases, grains, and ferroelastic, ferroelectric, and ferromagnetic domains. While diffraction patterns from experiments such as X-ray diffraction are presented in the reciprocal or Fourier space, mesoscale microstructure models such as the phase-field method naturally produce real-space images of spatial distribution of chemical composition, structural, and ferroic domains. Although one could rather readily compute the Fourier amplitudes of chemical and structural domain distributions generated by mesoscale simulations, they only contain information about the length scale and alignment of the real-space chemical and structure domains. Therefore, a direct comparison between diffraction experiments and mesoscale microstructure simulations is not possible. Here, we develop a theoretical approach to directly compute the crystal diffraction patterns of microstructures predicted by phase-field simulations. In particular, we consider five representative examples of microstructure patterns involving purely compositional domains, a single pair of tetragonal twin structures, multiple twin variants in a hexagonal system, ferroelectric polar vortices, and polycrystalline grains. The results are compared with previous experimental observations as well as X-ray diffraction experiments performed in the present study. Furthermore, the theoretical framework allows one to directly connect material microstructures and diffraction patterns predicted from phase-field simulations and the corresponding diffraction patterns from experiments, and thus providing guidance to experimental diffraction characterization and interpretation of microstructures.
We are leveraging modern algorithms and computational science to provide a route to predictive simulation of microstructure evolution on emerging exascale architectures. We are utilizing the fastest supercomputers in the world for modeling and simulation of microstructure evolution for generation of data under AM conditions. Solidification conditions in AM can be tailored for the reliable design of materials to specific performance requirements. Developing computational tools to further characterize alloys and correlate the processing-structure-properties-performance (PSPP) relationship.
Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.
Microstructure largely dictates macroscopic material properties and is strongly affected by processing. Therefore, the simulation of microstructure evolution in response to thermal fields during processing is of significant interest within the computational materials science community. Additive manufacturing (AM) has emerged as a technique for producing complex geometries and unique microstructures. Yet, complex and rapid thermal cycles in AM pose computational challenges for existing microstructure models. This work proposes a discrete event inspired cellular automata (CA) approach, titled DECA, to accelerate simulation of grain structure evolution in AM. In contrast to conventional time-stepped CA models, this model directly solves the times capture events would take place allowing for stepping in events rather than time (a technique also found in the field of discrete-event simulation). In comparison to purely serial discrete-event models, DECA allows for temporary violation of the causality constraint, but detects and corrects these violations, leading to an emergent phenomenon dubbed causality rippling, in which previously calculated capture events are overwritten. The amount of repeated calculations, defined by the capture ratio, is taken as a measure of computational inefficiency, and the model parameters that affect this ratio are evaluated. The new DECA approach was found to be more computationally efficient than conventional time-stepped CA models while guaranteeing an accurate solution, which can only be achieved in the conventional models for vanishingly small time steps. Finally, opportunities for parallelization and scaling of the new approach are discussed.