Predicting Precipitation Kinetics During the Annealing of Additive Manufactured Inconel 625 Components

The prediction of solidification microstructures associated with additive manufacture of metallic components is fundamental in the identification scanning strategies, process parameters and subsequent heat treatments for optimised component properties. Interactions between the powder particles and the laser heat source result in complex thermal fields in and around the metal melt pool, which will influence the spatial distribution of chemical species as well as solid-state precipitation reactions. This paper demonstrates that a multi-component, multi-phase precipitation model can successfully predict the observed precipitation kinetics in Inconel 625, capturing the anomalous precipitation behaviour exhibited in additively manufactured components. A computer coupling of phase diagrams and thermochemistry (CALPHAD)-based approach captures the impact of dendritic segregation of alloying elements upon precipitation behaviour. The model was successful in capturing the precipitation kinetics during annealing considering the Nb-rich and Nb-depleted regions that are formed during additive manufacturing.


Introduction
The manufacture of nickel-based superalloy components through powder-bed selective laser melting (SLM) is challenging. The interaction between the laser heat source and the randomly dispersed powder particles results in complex thermal histories characterising the process. Associated with the rapid solid-liquid-vapour transitions is the development of mechanical fields connected to the development of residual stresses. Mass transport of chemical species in the melt pool will be influenced by the solidification front during dendritic growth and promote chemical segregation [1]. Variations in chemical compositions between the interdendritic zones and the truck of dendrites will lead to heterogeneous distributions of precipitates after the application of subsequent heat treatments. This has proven to be a major issue in the additive manufacture (AM) of γ strengthened superalloys where uncontrolled precipitation of the intermetallic phase will lead to reduced ductility. Precipitation during AM will influence the mechanical response and consequently the development of component distortions and residual stresses. The ability to predict location specific microstructure variations is important in both understanding and predicting component properties following post-build heat treatments and in-service response. In addition, precipitation models can also be used to guide and optimise heat treatments and predict the evolution of Zener pinning precipitates during grain growth.
The chemical composition of the powder Inconel 625 of interest is given in Table 1 [2], showing a high amount of Niobium and low concentrations of γ forming elements. The alloy is typically thought as being solid solution strengthened; however, γ and δ precipitates form during long exposures at high temperature [5]. Both the γ and δ phases are intermetallic and are formed by the ordering of niobium to create either a body centred tetragonal D0 22 structure or a orthorhombic D0a crystal structure, respectively. The γ and δ phases compete for niobium, with δ being the more thermodynamically stable phase. δ is observed to precipitate upon grain boundaries whilst γ precipitates heterogeneously within grains, with preference for nucleating upon dislocations. In conventionally processed Inconel 625, Sauve et al. [3] observed γ to form with an initial spherical morphology and grow into lens shaped precipitates during heat treatment. The δ precipitates were observed to have plate-like morphology. Although Inconel 625 is known to be a weldable alloy, the AM of the alloy poses a problem. The thermal loading and cooling rates for precipitation varies considerable considering different AM methods, such as electron beam melting (EBM), selective laser melting (SLM), or laserdirect metal deposition (L-DMD) [2]. For example, Lass et al. [2] did not observe any intermetallic precipitates in SLM builds, whilst Amato et al. [4] observed fine nanoscale γ precipitates. This suggests that the precipitation kinetics is sensitive to the process parameters for SLM powder bed processing of Inconel 625. Amato et al. [4] observed γ particles in both SLM and EBM builds, which formed columnar arrays of precipitates on low angle grain boundaries. The spacing of the columns, the precipitate size, and precipitate concentration varied for the two builds with larger plate-like precipitates in the EBM specimen. It is known that significant dendritic segregation of niobium forms during SLM, and is associated with anomalous precipitation kinetics during post-heat treatments, accelerating the precipitation of δ so that it is fully precipitated after 1 h at 870 • C [2]. For comparison, in the as-rolled condition, Sauve et al. [3] observed the onset of δ formation after 3-4 h at this temperature.
This paper presents the mathematical model and its implementation that can capture this behaviour. This involves using multiple techniques to calculate the thermo-mechanical loading to allow for the prediction of precipitation during SLM. The mean-field model is presented in [6], and describes the application of mean-field modelling of the intermetallic precipitates in Inconel 718, which is similar to Inconel 625 with the addition of γ precipitates which form due to the larger Al and Ti content in the alloy. The model presented in this work includes refinements to the description of heterogeneous nucleation and modifications to capture the sluggish nucleation kinetics of γ and δ in Inconel 625.
The paper is structured as follows. Section "Model Formulation" presents the modelling framework, outlining the mean-field model. The next section outlines the model implementation, detailing the numerical methods and input parameters. The results are then presented showing an example of the finite element analysis (FEA) predictions, and the predicted evolution of precipitates during SLM and annealing. This is followed by a discussion section, identifying strengths and weaknesses of the current model formulation. A conclusion then summarises the findings of the paper.

Modelling Framework
The proposed modelling framework for the simulation of precipitation kinetics during the annealing of a SLM component is now presented. The approach involves simulation of the metal melt pool through a volume-offluid (VoF) formulation [7,8], from which the predicted melt pool geometry and the thermal loading are used to calibrate a finite element (FE) model of the heat source model when simulating the AM of a build [6]. The process-induced thermo-mechanical fields predicted by the FE analysis provide the driving forces required for solid-state precipitation reactions. Coupling of these fields to a CALPHAD thermodynamic approach, the influence of composition on phase transformations can be modelled. Statistics of the precipitation distributions are tracked through a composition dependent mean-field model of the precipitation size distribution, which accounts for nucleation, growth, coarsening and dissolution of precipitate phases. Figure 1 gives an overview of how these methods have been integrated.

Component-Scale Modelling
The simulation of an SLM build through a FE implementation is challenging due to the different length scales characterising the deposition process, heat source and size of the component. Observations of melt pool geometries suggest that large spatial gradients of temperature are generated within a depth of ≈ 80 μm. To ensure the accurate prediction of the formation of thermal strains incurred by such thermal loading, ideally the mesh size needs to be less than 20 μm. As a result, the computational cost of simulating geometries greater than 1 mm 3 with accuracy becomes prohibitive. To make simulations in a reasonable amount of time, methods for homogenising the deposition of material are needed. One approach is to calculate an eigenstrain that is characteristic of the specific SLM process and material, and then apply this eigenstrain to elements as they are activated, capturing the distortion of the geometry during the Fig. 1 The modelling framework applied to predict precipitation kinetics during annealing of a SLM component build. Alternatively, the thermal loading can be scaled to describe a larger deposition of material, simulating the SLM process and capturing the heat dissipate within the geometry of interest. The latter approach was adopted in this work as it allows for the calculation of location specific thermal loading in the component geometry, and is described in [6].
The commercial finite element package Abaqus was used to simulate the SLM process based on a decoupled thermal and mechanical analysis. The thermo-physical properties for Inconel 625 were obtained from Capriccioli and Frosi [9]. A Python script was written to automate the preprocessing operations. User subroutines were implemented, as described by Anderson et al. [6]. This includes the "property switching function", a description of the heat source, and a subroutine for updating boundary conditions. The thermal load in the FE model was calibrated so that the calculated peak temperatures aligned with predictions made using VoF approach [7,8].

Mean-Field Precipitation Modelling
As already mentioned a mean-field description is adopted to describe the evolution of the γ and δ precipitates, assuming weak interaction between particles. The precipitate morphology of γ and δ are simplified to be described by cylinders, and are characterised by the particle radius frequency density function, F(R, t), where R refers to the radius of a sphere of equivalent volume to the cylindrical particles of interest. Let F(R, t)dR describe the number of particles with a radius within the closed interval of R and R+dR within a representative unit volume. The mean particle radius R , the particle concentration N v and the volume fraction φ are obtained from moments of the F(R, t) function. The Dth moment of F(R, t) is given by the following: The particle concentration N v is given by M (0) , the mean particle radius R = M (1) /M (0) , and the volume fraction φ = (4π/3)M (3) . The evolution of the precipitate dispersion is governed by the following advection equation: where F + (R, t) and F − (R, t) are source and sink terms, describing mechanisms that affect the particle concentration, N v . In this model, they refer to nucleation and dissolution, respectively. However, they can also include precipitate coalescence phenomena [10]. The particle growth rate has the following generic form: where the term A(t) includes the diffusivity of the alloying elements, R c (t) is the critical particle radius and z(R, t) is the competitive growth correction factor. Marqusee and Ross [14] derived the following expression for the z(R, t) factor as follows: in the present work, it will be assumed that δ and γ precipitates do not directly compete for solute as δ nucleates at grain boundaries and γ nucleates within the grain interior, so that each phase has a separate z(R, t) factor. To describe the kinetics of γ and δ, the multi-component particle growth model proposed by Svoboda-Fischer-Fratzl-Kozeschnik (SFFK) [11] has been implemented. This framework was originally developed for spherical particles. To account for other precipitate morphologies, the shape factors derived by Kozeschnik et al. [12] and Svoboda et al. [13] have been used to describe the γ and δ phases as cylinders opposed to spheres, with an aspect ratio given by h = H/D, where H is the height of the cylinder and D is the diameter. Combining the SFFK model and the shape parameters of Kozeschnik et al.'s [12], the following expressions for the parameters for the particle growth rate given by Eq. 3 are obtained as follows: where γ is the interfacial energy, R g is the gas constant, and T is the absolute temperature. The terms c ki and c 0i refer to the chemical concentration of the ith alloying element in the particle and matrix phase, respectively. Similarly, μ ki and μ 0i refer to chemical potentials in the matrix and particle for the ith alloying element, and D 0i is the diffusivity of the ith alloying element. U is the misfit strain energy caused by differences in the lattice parameters between the precipitate and matrix phases, where the misfit strain is given by ε. E and v refer to the Young's modulus and Poisson's ratio. S k and O k are shape parameters and are functions of the precipitate aspect ratio, h. Nucleation is treated as outlined by Anderson et al. [6] as follows: where the Zeldovitch parameter Z, k b is the Boltzmann constant, and the form of the nuclei radius distribution function N c (R, t) are obtained from Jou et al. [15] as follows: R is approximated by the following: is the atomic volume. N 0 is the concentration of nuclei, and is estimated from the supersaturation of precipitate forming species as follows: where φ eq is the equilibrium volume fraction of precipitates.
For the δ and γ phases, the volume fraction φ(t) that enters Eq. 9 is the total volume fraction of the combined δ and γ precipitates. The term η refers to the fraction of active nucleation sites. For homogeneous nucleation, η is unity. To approximate η for heterogeneous nucleation, η is given by the ratio of the available nucleation sites divided by the total number of nucleation sites within the volume of interest. For heterogeneous nucleation upon dislocations within grains, η will be assumed to be given by the fraction of lattice sites occupied by dislocations. Let the dislocation density be defined as the total length of dislocations found within a volume, V so that ρ = L/V . The number of lattice sites along dislocations is given by L/b, where b is the Burgers vector. The total number of lattice sites within the volume can be determined by V b −3 . η is then given by the following: To complete the last expression, an estimate of the dislocation density is needed. It will be assumed that the thermo-mechanical fields induced by the deposition process will result in localised plastic distortions. Following Kocks et al. [16], the development of the dislocation density can be modelled through a continuity condition considering dislocation generation and annihilation mechanisms. Assuming that the generation rate scales with the dislocation density and dipole formation is the dominant annihilation mechanism, it can be shown that the dislocation evolution has the form as follows: where˙ p is the plastic rate, C ρ is a rate constant, and ρ s is the steady-state dislocation density. The steady-state dislocation density ρ s is given by the following [17]: where M is the Taylor factor, G is the shear modulus, v is Poisson's ratio, and σ is the effective stress. It is proposed that grain boundary dislocations act as nucleation sites for grain boundary precipitates. A lower bound estimate for the grain boundary dislocation density can be made by determining the dislocations needed to accommodate differences in lattice orientation between neighbouring grains. The dislocation density at the grain boundaries is likely to be higher, with additional dislocations formed during plastic deformation, and geometrically necessary dislocations that accommodate plastic strain. Carbides with incoherent lattice alignment with the neighbouring matrix would also provide additional dislocations where precipitates could nucleate. Let ρ gb refer to the mean dislocation density at the grain boundaries. The number of lattice sites that are found on grain boundary dislocations on a grain of diameter d is given by ρ gb π d 2 . The fraction of nucleation sites considering grain boundaries is then given by the following: Figure 2 illustrates the fraction of active nucleation sites considering heterogenous nucleation upon dislocations and upon grain boundaries, respectively.
The term β * is the atomic attachment rate, which uses the multi-component approximation derived by Svoboda et al. [11]: The energy barrier to nuclei formation G * shown in Eq. 6 is given by the following [6]: where ψ is the sphericity of the nuclei, and for cylindrical precipitates with an aspect ratio of h is given by the following: P inc is the nuclei incubation probability, and is given by the following: where τ is the incubation time. To determine the incubation probability during transient thermal loading, the incubation probability P inc is integrated numerically using the temporal derivative to Eq. 17 as follows [6]: (18) where t eq is the equivalent incubation time, t eq = −τ [ln(P inc )] −1 . The temperature dependent material coefficient C τ was introduced to delay nucleation kinetics in order to better capture the sluggish precipitation behaviour exhibited by Inconel 625.

Normalisation
The nucleation, growth, coarsening, and dissolution kinetics are calculated by solving the advection equation shown in Eq. 2 for each precipitate phase. To mitigate loss of significance errors and improve numerical stability, the problem is normalised in both space and time. Let R k and   (3). R k is either the mean precipitate radius or set to a minimum value for conditions where there are no precipitates. The spatial normalisation is updated regularly, to follow the dispersion. Let , where refers to the normalised particle radius distribution function. The normalised time and precipitate radius is given by t ≡ t/ and r ≡ R/R k , respectively. The normalised continuity equation is given by the following: (19) The normalised growth rate,ṙ(r, t ) is then given by the following: (20) where r c (t ) ≡ R c /R k and refers to the first moment of the function. The normalised nucleation rate is given by the following: (21) with the following normalisation parameters as follow: Equation 19 is solved using an explicit upwind finite difference scheme. The spatial discretisation used to generate the particle radius distribution functions are remeshed to follow the dispersions during the simulation. The values of F(R, t) which are above a certain limit are interpolated onto the new mesh using a cubic spline function, with linear interpolation applied to the remaining F(R, t) function.

CALPHAD
The multi-component, multi-phase mean-field model requires chemical potentials, phase chemistries, and diffusivities. The CALPHAD solver ThermoCalc [18] has been used to calculate these variables, using the thermodynamic database TTNi8, combined with the mobility database MOBNi1. A FORTRAN programme has been written which is coupled with ThermoCalc's [18] TQ FORTRAN interface. The phase diagram, phase chemistries, and equilibrium chemical potentials are calculated prior to the precipitation calculations. Look-up tables store this information, and are interpolated when calculating the evolution of the precipitate dispersions. During the precipitation calculation, the matrix composition is determined using mass balance, and  the matrix chemical potential and diffusivities are then calculated, allowing for the determination of the precipitate growth and nucleation rates for each phase. The simplified compositions used to model Inconel 625 in the conditions of interest are presented in Table 2. The extent of segregation in the AM condition was approximated from the energy dispersive X-ray spectrometry measurements presented by Zhang et al. [19].

Calibration
The precipitation model was calibrated using available data that describes the material in conventionally manufactured conditions. To capture the time-temperature-transition (TTT) data of Suave et al. [3], energy contributions of 500 and − 500 J/mol were added to the γ and δ phases, respectively. Figure 3 presents the calculated property diagram for the γ and δ phases for the chemistries shown in Table 2. Figure 3a and b describes the calculated phase fractions of γ and δ, respectively as a function of temperature. Three phases are considered in the thermodynamic calculations; the γ , γ , and δ phases. To evaluate the metastable γ phase, it is necessary to suspend the δ phase from the calculation. The chemistries examined show the predicted phase fraction for the bulk composition, and how this differs when considering the Nb-rich and Nb-depleted regions within the SLM component of interest.
The molar volume, misfit strains, and nucleation site fractions are given in Table 3. The coefficients for a polynomial temperature dependency for the Young's modulus and Poisson's ratio are given in Table 4, where the polynomial coefficients are defined by the following: where x is the variable of interest, the temperature, T , in Eq. 23 is given in degrees Celsius and n refers to the number of polynomial coefficients, p k . The temperature dependency of the interfacial energy and incubation probability coefficient for the γ and δ phases are presented in Tables 5 and 6, respectively. The approximation for the aspect ratio of γ " precipitates as a function of the equivalent precipitate radius was obtained from Moore et al. [20], and is shown in Fig. 4a. The aspect ratio of δ precipitates was approximated as 0.2. Figure 4b compares the predicted coarsening kinetics of γ at 650 • C with data from Moore et al. [20] and Suave et al. [3]. Figure 5 compares the calibrated TTT diagram with the measurements of Suave et al. [3]. The data refers to Inconel 625 in the as-rolled (AR) and shear-span (SS) conditions, which exhibit different precipitation kinetics. These differences may be caused by differing grain size, dislocation density, segregation of alloying elements, and residual stresses considering the two material conditions. The simulations assume that the material has been quenched from 1200 • C to room temperature at − 5 • C/s, and heated to temperature at 5 • C/s.

Results
The finite element model described in [6] was applied to predict the thermo-mechanical loading during AM of the geometry illustrated in Fig. 6a. The component geometry contains several arches. The supporting material on either side of the arches are referred to as legs. Precipitation simulations have been performed for the first two legs, showing similar behaviour. Figure 6b provides an example of the thermal loading predicted within the centre of the vertical cross section taken through the first leg.
The precipitation kinetics of γ and δ was calculated accounting for the segregation of alloying elements, shown in Table 2. The dislocation density in the AM condition was approximated by the steady-state dislocation density  1.00E-04 1.00E-05 1.00E-05 described in Eq. 12. The predicted von Mises stress and dislocation site fraction in leg 1 is shown in Fig. 7. The δ nucleation site fraction and the incubation coefficient of δ were both assumed to be 10× larger than in the conventional processed condition due to changes in grain boundary dislocation content and grain size. No significant precipitation was predicted in the AM condition in either the Nb-rich or Nb-depleted regions. Neither γ or δ was predicted to form during annealing within the Nb-depleted regions. The predicted size and volume fraction of the cylindrical γ and δ precipitates in the Nb-rich annealed condition are shown in Fig. 8. The contour map of γ follows the same shape as the contour of the dislocation nucleation site fraction η shown in Fig. 7; however, the differences in statistics across the geometry is minimal. The temporal evolution of the size and volume fraction of γ and δ at the centre of the vertical cross section of leg1 is given in Fig. 9. The predicted TTT kinetics for the AM condition is presented in Figs. 10 and 11, for the Nb-rich and Nb-depleted regions respectively. The model predictions are in reasonable agreement with Stoudt et al. [21], however over-predicts kinetics at high temperatures.

Discussion
The model predictions presented in this work show reasonable agreement with the experimental data, however Fig. 5 The simulated TTT diagram for the bulk composition compared with the TTT measurements from Suave et al. [3] for Inconel 625 in the as-rolled (AR) and shear-span (SS) conditions deviate at high temperature considering the TTT kinetics in the SLM condition as shown in Fig. 10. The δ precipitates are predicted to be stable at higher temperatures than Stoudt et al. [21] observed experimentally. Lindwall et al. [22] have also modelled precipitation kinetics during the heat treatment of AM Inconel 625, and have captured this behaviour. The approaches are similar in the use of FEA to predict the thermal history during AM and application of CALPHAD-based mean-field modelling of precipitation kinetics. They have demonstrated the use of the commercial software DICTRA [18] to simulate the formation of dendritic segregation of alloying elements during SLM, and have used TC-PRISMA [18] to simulate precipitation kinetics for different chemistries predicted by the DICTRA calculation.
Both the calculations presented in this paper and those by Lindwall et al.'s [22] share a similar problem in that the simulations do not account for any homogenisation Lindwall et al. [22] have captured the high temperature TTT kinetics shown in Fig. 10 through calibration of the interfacial energy and nucleation site concentration during conditions where the homogenisation of the dendritic segregation of Nb would impact precipitation behaviour. One concern is that such a model calibration may adversely impact the predictive capability when applied to other AM conditions where the degree of segregation of alloying elements may vary.
This work aimed at identifying model parameters that can capture precipitation behaviour in both conventionally processed and SLM conditions to give greater confidence when applying the model to new conditions. It was assumed that the interfacial energy did not change significantly for different processing conditions and that the differences in precipitation behaviour are caused by the extent of the segregation of alloying elements and differences in nucleation site densities. No attempts were made to calibrate the model to capture the high temperature TTT behaviour in the AM condition as the model does not yet account for the homogenisation of the dendritic segregation which would impact the property diagram and the chemical driving force for the precipitate phase transformations.
To improve upon the accuracy of the precipitation kinetics at long ageing times and during solid solution treatments, it is necessary to couple the precipitation calculation within the multi-component diffusion model simulating dendritic segregation. This would allow for the simulation of the formation of dendritic segregation during SLM followed by the homogenisation of the diffusion fields during subsequent heat treatments.
The precipitation kinetics in Inconel 625 are slow compared to Inconel 718, which contains higher amounts of Nb-and γ '-forming elements. A phenomenological parameter was introduced to delay nucleation, and better capture the precipitate behaviour. Brooks and Bridges [23] proposed that the slower kinetics of γ and δ in Inconel 718 is a result of the slower diffusivity of Nb. However, the diffusivity of Nb obtained from the MOBNi1 and MOBNi2 mobility databases are not significantly slower than γforming species. It is possible that the Nb diffusivities obtained from these databases are not accurate, however attempts to artificially reduce the diffusivity of Nb did not  have a significant impact upon the calculated incubation time or the growth rates of γ .
The particle growth rates of Svoboda et al. [13] and Moore et al. [20] go into greater detail regarding the description of the interfacial energy, precipitate misfit strain energy and the evolution of the aspect ratio. The evolution of the aspect ratio of the precipitates can be calculated as a function of the interfacial energy on the planar (γ p ) and curved surfaces (γ c ) of the cylinder and consider the misfit strain energy as a function of aspect ratio (U(h)). The particle growth rate of Moore et al. [20] may be expressed in the generic form given in Eq. 3: where ζ is a shape factor, and the terms G U , G γ , and G h describe the contributions to the critical particle radius Fig. 9 The temporal evolution of statistical information regarding the dispersions at a point in leg 1. Panels a and c present the evolution of the mean height and diameter of the cylindrical γ and δ precipitates, respectively. Panel figures b and d present the evolution of the precipitate volume fractions arising from changes to the misfit strain energy, interfacial energy, and aspect ratio, respectively. These terms are given below as follows: The impact of these factors upon precipitation kinetics can be discussed in terms of the classical regimes of precipitation: nucleation, growth, coarsening, and dissolution. For precipitates to nucleate, the critical particle radius must be smaller or equivalent to the critical nuclei radius. For a unimodal precipitate dispersion during the regime of growth, R c < R . During steady-state coarsening, R c ≈ R , and during dissolution R c > R . Once the precipitate Fig. 10 The predicted time-temperature-transition behaviour within the Nb-rich region of a AM component. The data for Inconel 625 in the AM condition is measured by Stoudt et al. [21], and refers to the time needed for the formation of ≈ 1% volume fraction of δ reaches a dynamic equilibrium shape, misfit strain energy and interfacial energy contributions balance, and dh dR → 0. The contributions to R c shown in Equation set (25) would diminish once the precipitates morphology changes sufficiently to reach a dynamic equilibrium shape, with reduced impact during coarsening kinetics compared to the regime of growth. The inclusion of these details may improve the ability to capture the rate at which the volume fraction of γ and δ are predicted to reach a dynamic equilibrium. The difficulties encountered when pursuing this approach is the determination of the planar and curved interfacial energies, and the misfit strains of the γ and δ phases. These vary as a function of composition and temperature. Atomistic calculations can provide an estimate for these parameters; however, such methods are difficult to extend to the high number of alloying elements commonly used nickel-based superalloys. Fig. 11 The predicted time-temperature-transition behaviour within the Nb-depleted region of a AM component

Conclusions
• A mean-field theoretical framework has been developed for precipitation in nickel-based superalloys during AM, which captures nucleation, growth, coarsening and dissolution. • An extension to classical nucleation theory is presented that accounts for the influence of grain size and deformation upon precipitation. • Evolution of the precipitate size distributions is handled through a multi-component, multi-phase description of particle growth rate as proposed by the SFFK model. The chemical driving force of stable and metastable phases can be calculated, capturing the competition between the γ and δ phases. The CALPHAD approach can link the impact of variations in chemistry upon precipitation behaviour. • The precipitation model suggests that the anomalous precipitation behaviour exhibited by AM Inconel 625 during annealing is largely caused by the dendritic segregation of Nb. No precipitation is predicted to form within the Nb-depleted regions. Differences in dislocation density and grain size are predicted to have a smaller impact upon precipitation behaviour in comparison to the role of dendritic segregation of Nb. • The low volume fraction of precipitates in Inconel 625 means that the mean-field assumption of weak interactions between precipitates is reasonable when modelling this alloy. • To avoid the excessive precipitation of δ, the anneal temperature needs to be high enough to homogenize Nb whilst minimising grain growth. If possible, AM process parameters and build strategies which minimise the segregation of Nb would help mitigate issues arising from the rapid precipitation kinetics of δ and γ .