Process-Structure-Properties-Performance Modeling for Selective Laser Melting

: Selective laser melting (SLM) is a promising manufacturing technique where the part design, from performance and properties process control and alloying, can be accelerated with integrated computational materials engineering (ICME). This paper demonstrates a process-structure-properties-performance modeling framework for SLM. For powder-bed scale melt pool modeling, we present a diffuse-interface multiphase computational ﬂuid dynamics model which couples Navier–Stokes, Cahn–Hilliard, and heat-transfer equations. A computationally efﬁcient large-scale heat-transfer model is used to describe the temperature evolution in larger volumes. Phase ﬁeld modeling is used to demonstrate how epitaxial growth of Ti-6-4 can be interrupted with inoculants to obtain an equiaxed polycrystalline structure. These structures are enriched with a synthetic lath martensite substructure, and their micromechanical response are investigated with a crystal plasticity model. The fatigue performance of these structures are analyzed, with spherical porelike defects and high-aspect-ratio cracklike defects incorporated, and a cycle-amplitude fatigue graph is produced to quantify the fatigue behavior of the structures. The simulated fatigue life presents trends consistent with the literature in terms of high cycle and low cycle fatigue, and the role of defects in dominating the respective performance of the produced SLM structures. The proposed ICME workﬂow emphasizes the possibilities arising from the vast design space exploitable with respect to manufacturing systems, powders, respective alloy chemistries, and microstructures. By digitalizing the whole workﬂow and enabling a thorough and detailed virtual evaluation of the causal relationships, the promise of product-targeted materials and solutions for metal additive manufacturing becomes closer to practical engineering application.


Introduction
Selective laser melting (SLM) is a novel technique to produce metallic parts with complex and easily customizable geometries with limited subtractive steps such as machining, typically for higher-end applications in aerospace, automotive, and biomedical industries [1]. In SLM, a metal powder layer is spread on a substrate, after which a laser melts the layer at preselected regions. New layers are repeatedly added and selectively melted, until the final part has formed. SLM is a relatively new manufacturing technology, and as a result, the current alloy compositions might need

Diffuse-Interface Melt Pool Modeling
The melt pool model can be used to accurately predict the melt pool geometry, gas pore trapping, calibrate the part scale heat conduction model, and calculate thermal gradients and cooling rates for the rapid solidification phase field model.
The model includes the following phenomena: heat-transfer, melt and gas flow, buoyancy, melt-gas interface dynamics including Marangoni effect, and solidification. This melt pool dynamics is described with a fully coupled system of equations: energy equation for temperature evolution, Navier-Stokes equations for fluid flow, Cahn-Hilliard equation for metal-gas phase tracking.
The two-phase flow with different densities and viscosity is modeled by the diffuse-interface model proposed in [16]: ∂ ∂t (ρv) + ((ρv + J) · ∇) v − ∇ · (2ηDv) + ∇p − ρg − f = 0, (1) ∇ · (v) = 0, δF with suitable initial and boundary conditions. Equations (1) and (2) are the Navier-Stokes equations, and Equations (3) and (4) are the Cahn-Hilliard equation split into two equations, where J = − ∂ρ ∂φ mg m∇µ is the correction flux that is zero in case of equal phase densities; φ mg denotes the metal-gas phase field, µ is the chemical potential; v is the volume averaged velocity; p is the pressure; ρ = ρ(φ mg ) = 1 2 ((ρ m − ρ g )φ + (ρ g + ρ m )) is the mean density, where ρ g is the gas or vapor density; and ρ m is the density of the condensed solid or liquid metal. η is the mean viscosity defined similarly to the mean density. The mobility is denoted by m = m(φ mg ). The gravitational force is denoted by g, and f gathers other possible external forces, like stress due to surface tension. Dv = 1 2 (∇v + (∇v) T ) is the symmetrized gradient (the rate of deformation tensor).
Let the scaled surface tension be denoted byσ and the interfacial width be proportional to . F := F (φ mg , ∇φ mg , φ sl , ∇φ sl , T) represents the global free energy functional, where subscripts mg and sl represent the variables associated with metal-gas and solid-liquid phase functions, respectively. Then, a suitable form of F makes it possible to model systems with multiple phases as F =σ Ω mg 2 |∇φ mg | 2 + −1 mg w(φ mg ) dΩ + Ω sl 2 |∇φ sl | 2 + (1 −φ mg ) f metal +φ mg f gas dΩ, (5) where the free energy density of the metal phase is and w(·) and g(·) are suitable double-well functions, p(·) is an interpolation function. Evolution of metal-gas phase function φ mg is governed by the Cahn-Hilliard model in Equations (3) and (4). φ sl is a phase function, which distinguishes solid and fluid domains. In general, one would need to solve the Allen-Cahn equation to simulate the dynamics of the melt interface. However, we assume the melt interface position to be at steady-state at any instant in time and use an analytical expression for the phase function φ sl as a function of temperature T: where T m is the melting temperature and ∆T is the temperature range for diffuse-interface. Then, material parameters may depend on φ sl , as well as φ mg and T.
Modeling solid-liquid phase transition and Marangoni effect requires the temperature field. Therefore, the Navier-Stokes/Cahn-Hilliard system is coupled to a variant of the convective heat-transfer equation, described in Section 2.2.
The Marangoni effect is the convection due to a gradient of the surface tension. The variations in the surface tension are caused by temperature dependence. In order to model the thermal induced Marangoni effect, the temperature-dependent surface tension coefficient σ is approximated by the Eötvös rule and it takes the form of a linear function on the temperature such that where the coefficients σ 0 > 0, a > 0, b = 0 are assumed to be constants. The constants a and b are related to the capillary and the Marangoni numbers [17].

Large-Scale Heat-Transfer Modeling
A computationally efficient heat conduction model is simulate temperature distribution in large volume. This model is used to roughly estimate the melt pool geometry and the associated thermal gradients and cooling rates in order to drive alloy phase field simulations, and secondary heat cycles to obtain better solid-state behavior such as back-diffusion, precipitation, and allotropic transformations.
The model interpolates the thermal conductivity and heat capacity between powder and bulk solid phases. Latent heat term is modeled as a Gaussian term. The evolution of temperature T is modeled with the following heat conduction equation: where C p (T) is the temperature-dependent heat capacity, k(T) is the temperature-dependent heat conductivity, L is the latent heat, ∆T = T l − T s is the difference between liquidus and solidus temperatures, and T mid = (T l + T s )/2 is their average value.
The effect of a powder layer is taken into account simply as the decrease in heat capacity and thermal conductivity. We consider a 20 µm powder layer, which is assumed to decrease the heat capacity by an amount corresponding to the packing fraction, assumed to be 55%. Thermal conductivity is assumed to decrease to 20% of the bulk value.
The top boundary has a fixed flux boundary condition, which is a sum of heat fluxes from the laser heat source I laser , convective cooling I convection from the inert process gas, and radiative cooling I radiation . The laser is assumed to have a Gaussian intensity profile, whose intensity at top boundary ∂Ω top is implemented as a boundary condition: η is the gray-body absorptivity, R L is the Gaussian beam width, f is the power intensity factor that describes how focused the intensity is around the beam center [1], and x center (t) is the time-dependent center of the laser beam. Convection flux is where h c is the effective convective heat-transfer coefficient and T ambient is the ambient temperature. The radiative cooling is assumed to follow a simple gray-body radiation where σ SB is the Stefan-Boltzmann constant. Note that, in Equations (9) and (11), we assume that the metal surface (of powder, liquid, and solid) is a gray body in radiative exchange equilibrium, and thereby that emissivity is equal to absorptivity. Boundaries other than the top surface are assumed to have zero heat flux boundaries. The relevant heat-transfer model parameters are shown in Table 1.

Rapid Microstructure Solidification Phase Field Modeling
We use phase field method to simulate the rapid solidification of Ti-6-4, where we neglect composition diffusion because the small melting range of Ti-6-4 is assumed to lead to negligible microsegregation and constitutional undercooling effects. Therefore, we model the solidification with the thermal phase field model of [18]. Heterogeneous nucleation from inoculants is modeled by incorporating thermal noise to the order parameter whose strength is determined from fluctuation-dissipation theorem. Inoculants are described with a wall field introduced in [19].
Assuming in total N unique crystalline orientations α, the evolution of order parameters φ α are governed by where τ is the interface attachment time scale, W α = W α (n α ) is the anisotropic interface width of orientation α and n α = ∇φ α /|∇φ α | is the interface normal, ω sets the strength of repulsive interaction between order parameters of different crystalline orientations [20], λ is the so-called coupling constant, η α is the thermal noise term discussed below, φ W is a static "wall field" introduced to induce heterogeneous nucleation according to [19], σ SL is the solid-liquid interface energy, h is a parameter used to set control the contact angle, and N is the total number of crystalline orientations considered. The noise term η α in Equation (12) is a random variable which is uncorrelated in space and time, and follows fluctuation-dissipation theorem: The evolution of dimensionless temperature U is where D is the heat diffusion coefficient, assumed to be constant. Matched interface asymptotics from [18] is used to relate physical material parameters (thermal capillary length d o and kinetic coefficient β) to phase field parameters (W, τ, and λ) as follows: where a 1 ≈ 0.8839 and a 2 ≈ 0.6867 are asymptotic constants. We assume a negligible capillary anisotropy due to rapid solidification conditions [21]. The kinetic coefficient β(n) is assumed to have weak cubic lattice anisotropy: where β 0 is the magnitude of the anisotropic kinetic coefficient β(n α ), and k is its anisotropy strength.
The relevant phase field model parameters are shown in Table 2. In this phase field model, we only consider the primary solidification (liquid-to-β), and neglect the martensitic transformations. However, the lamellar martensitic structure is synthesized into the polycrystalline structure given by the phase field model, and the micromechanical response of this martensitic structure is subsequently analyzed with a crystal plasticity model.

Mechanical Response with a Micromechanical Crystal Plasticity Model
The micromechanical approach relies on crystal plasticity modeling to establish the stress state as well as evolution of slip-system-specific plastic slip during cyclic loading. A multiplicative decomposition of elastic and plastic displacement gradients is used for large strains. The resolved shear stress is utilized as the criterion for plastic flow. Classical crystal plasticity framework is utilized in current work, i.e., no higher order terms inducing scale dependencies are included.
The microstructure is modeled as a BCC β-phase or martensitic HCP α -phase. The methodology enables the modeling of α + β-containing alloys or α and retained β-phase ones, but in current work, the focus is on fully α consisting microstructures, especially the martensitic α . Thus, for example, slip transmission or more complex deformation mechanisms and their interactions arising at phase interfaces are not investigated nor included. Few cases are evaluated with primary solidified β grain dominated microstructures, for reference, the α -phase, especially with respect to its acicular morphology, is of interest particularly since it is a probable outcome considering the rapid solidification conditions persistent in metal AM. With respect to the α -phase, the HCP is modeled considering the prismatic ({1010}, 1120 systems, three slip systems), basal ({0001}, 1120 systems, three slip systems), and pyramidal ({1011}, 1123}, 12 slip systems) slip systems. For the β-phase, 12 slip systems of type ({110}, 111 ) were considered, the rate of plastic slip and evolution equations for the hardening terms are described utilizing similar expressions, as presented in more detail below.
In order to capture mixed-hardening response suited for assessment of cyclic loading, both isotropic and kinematic contributions are incorporated, the model is adapted from work performed in [23]. Resolved shear stress for slip system s is obtained classically from where σ is the Cauchy true tensor, m is the vector parallel to the slip direction, and ni s the vector normal to the slip plane. The slip rate for slip system s is evaluated according to a phenomenological relation: where K and n are material parameters of the constitutive formulation associated with viscous response, γ cum is the cumulative plastic slip, and the isotropic and kinematic hardening contributions are provided by where R 0 , Q, b, and h rs are material dependent coefficients and where c and d are material related coefficients. The calibration of the various coefficients is based on numerical homogenization, following the practices presented in [13,15]. For the braces, a convention is adopted where x = max(x, 0).
Since the focus of current work is in establishing crystal plasticity-based micromechanical worklowsfor Ti-6-4 rather than quantitatively assessing specific use cases, literature data is utilized in establishing model parameterization. Different slip family critical resolved shear stresses are obtained from [24], and the viscoplastic parameterization from [23], which is adapted to match additively manufactured Ti-6-4 uniaxial tensile loading properties, as presented in [25] (by adjusting the hardening parameters to yield a comparable tensile response). Texture data is adapted from PF simulations and, for the interaction coefficients, only diagonal initial hardening modulus is used, neglecting latent hardening. Microstructural model generation, operations on the primary solidified β structure, and meshing of the resulting geometry are performed by utilizing the toolset presented in more detail in [11,26].
The micromechanical assessment of fatigue is performed utilizing the so called Fatigue Indicator Parameters (FIPs). To that effect, we utilize the Fatemi-Socie (FS) criterion, as presented in [27], and applied it to the high-cycle fatigue (HCF) regime, for example, in [28]. The FS criterion one of the most utilized FIP parameters, commonly defined as: where ∆γ p max is the change in maximum cyclic plastic strain over a finite volume (in current work interpreted as a local maximum greater than adjacent field mean value) and σ max n is the stress normal to affiliated plane. The remainder of the parameters follow a Coffin-Manson type of a representation for strain life, σ y is the cyclic yield strength and K is a parameter weighing the coupling between normal stress and cyclic plastic slip (K in current is as is common considered to have a value of unity). γ f and c are the Coffin-Manson parameters, N c are the respective cycles to fatigue failure. The parameter has been applied in different contexts commonly over cumulative plastic slip, but also over specific slip systems respective behavior. In current work the parameter is evaluated with respect to cumulative plastic slip. The fatigue damage evolution law parameters are obtained from [29], where cyclic properties and fatigue response of additively manufactured Ti-6-4 are investigated experimentally.
With respect to the fatigue analysis microstructures, the different equiaxed and columnar α and β-phase-dominated microstructures are loaded cyclically under R = 0 loading ratio in strain amplitude control to investigate the low-cycle fatigue (LCF)-and HCF-response of the respective microstructures. Additionally, since in the context of AM defects and their influence to fatigue performance are of particular interest, postulated microstructure scale defects are introduced explicitly to the microstructures and the respective effect to fatigue life prediction assessed. This consists of porelike and cracklike defects with sizes comparable to the microstructural morphological features with respect to the acicular α laths.

Diffuse-Interface Melt Pool Dynamics
The melt pool dynamics model is implemented using FEniCS Project library [30]. The system of equations is discretized in space using the mixed finite element method, and in time using an implicit method. Stable time discretizations are discussed, for example, in [31,32]. The whole domain size was 3 × 1 mm 2 , and we used triangular elements refined around the powder bed region, where the smallest element size was 6.25 µmx. The mesh consisted of 147,120 triangles in total. Quadratic elements were used for velocity vector and linear elements for pressure, temperature, phase field, and chemical potential; this results in 883,910 degrees of freedom (DOFs). The resulting nonlinear system is solved using PETSc's SNES solver [33,34] and Newton's method with line search. Additionally, we rely on FEniCS for differentiating the nonlinear residual functional to obtain the Jacobian. At each Newton iteration, the linear system is solved using MUMPS [35].
The powder bed was generated by first generating a fixed number of powder particles following a log-normal size distribution, after which the particles were iteratively packed to fill the assigned powder layer thickness.
An example of 2D melt pool dynamics is shown in Figure 2, where the voids between the powder particles coalesce into bubbles due to surface tension. These bubbles can then freeze into the structure, depending on how rapidly the solidification front advances. In the future, we will carefully consider associated model parameters, such as interface penalty term and the double well-thickness, and extend the diffuse-interface model into three dimensions.

Large-Scale Heat Transfer
The large-scale heat-transfer model implemented in MOOSE framework [36]. The mesh consists of first-order hexahedral lagrangian elements with adaptive mesh refinement, which are set to resolve solid-liquid boundary region. On highest refinement level, the element size is 9 µm. We use a two-step backward differentation formula (BDF2) for time integration, with adaptive time stepping targeted to maintain six nonlinear iterations per time step. The system of equations uses Single Matrix Preconditioning (SMP), which is solved with Newton's method.
The resulting temperature distribution in a single line scan is shown in Figure 3 (top left). The model gives an estimate of the overall melt pool geometry (top right). The model can also be used to estimate the local thermal conditions that drive the solidification-the thermal gradients are estimated along certain directions around the trailing edge of the melt pool (bottom-left). In addition, thermal cycles can be estimated in different scanning patterns: Figure 3 (bottom-right) shows the temperature distribution in a zig-zag scan pattern. Regarding the large-scale heat-transfer model, we believe that, specifically, the following model parameters are both hard to estimate from literature or theory and affect the heat transfer significantly: gray-body absorptivity η and powder conductivity fraction f powder cond . This model assumes that convective heat transfer in the melt pool is negligible. In the future, we plan to coarse-grain convection effects from the melt pool level model of Section 3.1 into this simple and efficient heat conduction model for more realistic large-scale heat-transfer modeling.

Rapid Solidification Microstructures
The phase field simulations use the finite difference method with explicit Euler forward time stepping, with the time step size set to 0.8 of the linear stability limit for the heat diffusion equation. The mesh was adaptively refined to capture gradients of order parameter and dimensionless temperature with the software platform introduced in [37]. The used phase field model parameters are shown in Table 2.
We perform polycrystalline Ti-6-4 solidification simulations under a constant initial dimensionless undercooling corresponding to ∆ = 0.5. Figure 4 (left column) shows the competitive growth of nuclei initiated at the bottom of the undercooled melt pool; dark red corresponds to the grain orientation with fastest growth direction closest to the vertical axis, and thereby, it eventually dominates the solidification, leading to an epitaxial polycrystalline structure. The right column of Figure 4 shows how the epitaxial growth is interrupted with randomly placed inoculants in the melt pool; the bottom of the melt pool is initialized with a number of epitaxial fast-growing grains. The inoculants are large enough to nucleate multiple grains of different orientations simultaneously, leading to several high-angle triple junctions. The heterogeneously nucleated grains grow rapidly and radially, giving rise to straight grain boundaries. In contrast, the inoculant-free system (Figure 4, left) has smoother grain boundaries due to the gradual spreading of the fast-growing (dark red) orientation.
Note that, in this simulation, the constant initial undercooling leads to rapid nucleation from all inoculants added to the melt pool. For more quantitative results, we are extending the analysis to a more realistic temperature distribution, coupled to the melt pool model of Section 3.1 and the large-scale heat-transfer model of Section 3.2. The repulsive grain-grain term in Equation (12), whose strength is determined by interaction parameter Ω, gives qualitatively the right interaction between grains, but the solid-solid interface energy cannot be directly controlled with this approach. We are extending this approach in the future to a more quantitative grain-grain repulsion.

Micromechanical Properties and Performance Modeling
The phase-field-simulated epitaxial and equiaxed Ti-6-4 polycrystalline structures are used as direct input microstructures to the crystal plasticity modeling of monotonic stress-strain response as well as micromechanical fatigue performance modeling. The polycrystalline structure is assumed to be completely transformed to a simplified description of lathlikemartensite, as shown in Figure 5. First, a simple tensile loading was performed for these structures, as shown in Figure 6. The equiaxed structure results in a smoother distribution of stresses and strains, whereas the epitaxial structure shows anisotropic stress and strain peaks. In addition, the equiaxed α -dominated structures in general consist of a larger number of misoriented interfaces, which, in the crystal plasticity modeling, manifests itself as increased hardening and greater strength of the resulting microstructure. Results for the cyclic fatigue loading of these structures are shown in Figure 7. The upper row shows the cumulative plastic slip at different strain amplitudes, for a constant number of loading cycles computed such that a stabilized condition has been obtained. The bottom row introduces two different postulated defect types at a high-angle grain boundary (HAGB): a spherical porelike defect due to excessive laser key-holing, and a high-aspect-ratio cracklike defect due to solidification cracking or incomplete powder melting. The fatigue failure of these structures was estimated by utilizing the FS FIP to evaluate their significance with respect to cycles to fatigue failure. By considering different amplitudes with different types of defects and the micromechanical measures regarding when they are expected to lead to the nucleation of a macroscale crack identifiable as structural failure (or short fatigue crack nucleation which can already be considered as descriptive), these cyclic loading simulations can be used to generate cycles to failure vs. amplitude fatigue graphs, as shown in Figure 8. For high-strain amplitudes, the nominally-similar different defect-free microstructures show a very similar behavior and break down at similar cycle counts, although the differences are easily several tens of percents. As the amplitude is decreased, there is more scatter in the number of cycles to fatigue failure for the considered microstructures. The introduction of microstructural defects influences the predicted outcome significantly. The simulations show that, compared to spherical pores, high-aspect ratio cracklike defects have a poorer fatigue performance, which is consistent with the literature [1]. Overall, the detrimental effect to fatigue performance of metal AM microstructures is reproduced, and the defects that have been introduced could act as nuclei for, for example, HCF failure. This adds to the confidence that the proposed analysis approaches are beneficial in assessing the micromechanical performance and further developing fatigue-resistant metal AM solutions.

Discussion
We presented a modeling framework to investigate the process-structure-properties performance chain for selective laser melting (SLM), where the presented models can be readily coupled in a sequential manner.
The diffuse-interface melt pool model can be used to assess the feasibility of process parameters-primarily, scanning speed, laser power, beam intensity intensity profile, and layer thickness. It gives a high-fidelity estimate of the temperature distribution around the melt pool, which can be used to parameterize a larger-scale heat conduction model, or to drive solidification on a smaller scale via phase field method, or polycrystalline structure formation via cellular automaton or kinetic Monte Carlo type models.
In particular, the model can predict the formation of trapped gas pores, overall melt pool instability due to balling or spattering, and irregularity of the solidified-free surface. The model gives an accurate description of the heat transfer inside the melt pool, where convection, Marangoni effect, and evaporation are crucial features for accurate description of heat transfer in addition to simple conduction. A realistic ray-tracing model should be implemented, such as one used in [10]. We employed a simple gray-body assumption to describe the laser absorption. In reality, for example, the vapor phase can cause significant deviation from this gray-body behavior. Absorptivity can also have a significant temperature dependence. Using a more complicated laser absorption model would require careful experimental validation.
In the future, the melt pool model requires quantitative parameterization and validation, by using experimental data on temperature distribution around the melt pool, based on in situ temperature measurements via thermocouples or thermal cameras, and estimation of melt pool dimensions based on post mortem microstructure characterization. The model includes multiple parameters, which are difficult to estimate or find in the existing literature. Parameter sensitivity analysis of the model will help in the task of testing robustness of the simulation results for a range of different physical parameters, representing different materials and manufacturing conditions. Large computational cost of the model due to the multiphysics coupling prohibits the use of the model in the design optimization loop. One way to mitigate the cost is to construct the surrogate model. In this context, the calibrated large-scale heat-transfer model can be used for the coarse grained description of costly melt pool diffuse-interface simulation.
The large-scale heat-transfer model can be used for relatively fast estimation of overall temperature distribution and secondary heat cycles using a specific scanning strategy, and to assess the risk of forming lack-of-fusion type defects. This will be useful for assessment of scanning strategies, heat accumulation across layers, formation of residual stresses, and microstructure annealing in secondary heat cycles. Efforts to improve the computational efficiency of these models will be of primary interest, in order to provide realistic overall temperature distribution on a level of a full part. Innately, multiscale and adaptive finite element methods and on-the-fly coarse graining between solvers is expected to play an ever increasing role towards the future, in order to push forward both the increase in computational resolution and simultaneously target the reduction of the computational cost.
The presented rapid solidification phase field model can be used to investigate the formation of polycrystalline microstructures as a function of process conditions, alloying, and, e.g., inoculants. In particular, an accurate description of solid-state transformations, such as α martensite formation, will be of great interest. Quantitative usage of lower scale descriptions, such as phase field crystal type models, will be of great interest to better understand the formation of nonequilibrium structures, and, e.g., dislocation formation.
Finally, the mechanical response of these rapid solidification microstructures can be predicted with a micromechanical crystal plasticity model. In addition, these models can assess the importance of different defect types and morphologies, hard oxide inclusions, or the martensite lath spacing. Among other microstructural heterogeneities, defects become particularly important in considering material failure and structural performance, for example, in high-cycle-fatigue conditions. Close coupling between phase field and crystal plasticity models will become desirable to describe, for example, solid-state transformations or microstructural damage. This can be achieved either through a priori phase field modeling studies for crystal plasticity model parameter calibration, concurrent two-way coupled descriptions, or, as would be preferred, fully coupled models.
The presented physics-based models can be readily coupled sequentially, and they provide a basis for investigating causalities within the composition-process-structure-properties performance chain for SLM. However, via brute force computation, it is impractical to cover the relevant materials engineering design space in terms of an optimal alloy composition and optimal processing parameters, and assess whether the resulting structure produces feasible properties and performance behavior. Taking full advantage of the presented modeling framework, one needs high-throughput [38], in situ experiments, and materials informatics methods [39].
Overall, the simulations lead to vast amounts of data, from which relevant information needs to be extracted. These issues can be tackled with concepts from materials informatics [39], where the idea is to make use of present and emerging methods from data science to analyze and drive the modeling chains.
New physics can be extracted, by either choosing between competing models or by determination of hard-to-measure material properties; for example, anisotropic kinetic coefficient or solute-trapping parameters in a phase field model could be extracted from so-called dynamic transmission electron microscopy (DTEM) experiments [40][41][42][43][44]. Expensive, high-fidelity experiments (or simulations) can be better designed and targeted to give maximal useful information through so-called adaptive learning schemes. Relevant microstructural features can be identified, reduced descriptions of microstructures to generate surrogate models for process-structure or structure-property-performance mapping. Inverse design schemes can be used to turn process-structure-properties performance modeling chain around. The uncertainty and sensitivity of an optimal design can be better assessed.
The presented modeling chain will be used together with high-throughput screening of feasible alloy compositions, quick-and-dirty based on CALPHAD-thermodynamic databases, and completely new compositions can be explored with density functional theory (DFT)-based approaches. This can be complemented with high-throughput experimentation, including rapid combinatorial materials synthesis and characterization schemes [38,45].
Future efforts should include in-depth analyses of each presented model separately, with explicit comparisons to experiments. This includes validation of phase field simulations of rapid solidification via so-called Dynamic TEM (DTEM) experiments [43], and in situ SEM tensile testing and nanoindentations to calibrate the crystal plasticity model for an SLM material.

Conclusions
We presented an integrated computational materials engineering framework (ICME), demonstrating the behavior of a diffuse-interfaced melt pool model for detailed powder-bed scale process modeling, complemented by an efficient large-scale heat-transfer model. The process models can be used to estimate the risk of forming different type of defects, temperature distribution in solidification, and solid-state secondary heat cycles. Phase field modeling was used to predict the formation of rapid solidification microstructures for Ti-6-4 alloy in epitaxial and equiaxed growth regimes. These solidification microstructures, with potential defects, were analyzed with respect to their micromechanical response and fatigue performance using a crystal plasticity model. Future work will improve and expand the existing models, develop experimental calibration and validation procedures, and form closer couplings between the presented models.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:

AM
Additive Manufacturing CP Crystal Plasticity FIP Fatigue Indicator Parameters SLM Selective Laser Melting PSPP Process-structure-properties-performance ICME Integrated Computational Materials Engineering