Thermal Modeling of Temperature Distribution in Metal Additive Manufacturing Considering Effects of Build Layers, Latent Heat, and Temperature-Sensitivity of Material Properties

A physics-based analytical model is proposed in order to predict the temperature profile during metal additive manufacturing (AM) processes, by considering the effects of temperature history in each layer, temperature-sensitivity of material properties and latent heat. The moving heat source analysis is used in order to predict the temperature distribution inside a semi-infinite solid material. The laser thermal energy deposited into a control volume is absorbed by the material thermodynamic latent heat and conducted through the contacting solid boundaries. The analytical model takes in to account the typical multi-layer aspect of additive manufacturing processes for the first time. The modeling of the problem involving multiple layers is of great importance because the thermal interactions of successive layers affect the temperature gradients, which govern the heat transfer and thermal stress development mechanisms. The temperature profile is calculated for isotropic and homogeneous material. The proposed model can be used to predict the temperature in laser-based metal additive manufacturing configurations of either direct metal deposition or selective laser melting. A numerical analysis is also conducted to simulate the temperature profile in metal AM. These two models are compared with experimental results. The proposed model also well captured the melt pool geometry as it is compared to experimental values. In order to emphasize the importance of solving the problem considering multiple layers, the peak temperature considering the layer addition and peak temperature not considering the layer addition are compared. The results show that considering the layer addition aspect of metal additive manufacturing can help to better predict the surface temperature and melt pool geometry. An analysis is conducted to show the importance of considering the temperature sensitivity of material properties in predicting temperature. A comparison of the computational time is also provided for analytical and numerical modeling. Based on the obtained results, it appears that the proposed analytical method provides an effective and accurate method to predict the temperature in metal AM.


Introduction
Metal additive manufacturing (AM) is a "process of joining materials to make objects from 3D model data, usually layer upon layer, as opposed to subtractive manufacturing methodologies" [1]. Additive manufacturing (AM) processes have potential to be the pillar of the next industrial revolution. AM can be used to improve existing manufacturing processes and rapidly introduce new prototypes and products [2,3]. It also offers the potential to spin off entirely new industries and lead to new production methods [4]. AM offers design flexibility, the ability to produce complex parts, and lower cost due to the reduced requirement of materials and decreased lead time.
AM may also hold the potential for the repair and replacement of existing plant components [5,6]. Results from on-line monitoring and complimentary non-destructive evaluation (NDE) inspections can provide indications of component health and enable repair or replacement prior to a forced outage situation. As an example, imaging tools and software can be leveraged to create a digital image of the damaged component which can be used to 3-D print a new one. This can be especially advantageous if the component is no longer in production and/or would require a long-lead time to fabricate.
There are many challenges that necessitate being focused on this field in order to expedite the adoption of AM as an advanced manufacturing technology. The issues in this field can be classified in to: the distortion, fatigue, defects, and residual stress of the manufactured parts [7][8][9]. The modeling in additive manufacturing technology is a key to the advancement of the field due to obstacles in in-situ measurements of temperature, thermal stress, residual stress, and distortion. The available knowledge and technology to-date on the descriptions and predictions of the metal AM process have been fragmented, mostly driven by phenomenological or numerical observations [10][11][12][13], and primarily limited to macroscopic analysis in nature [14,15], thus restricting the full capability and potential of the AM process. Using numerical methods and experiments are not just expensive, but also time-consuming. On the other hand, the physics-based analytical modeling eliminates all the above-mentioned difficulties and can help to better understand the physical aspects of the metal additive manufacturing process.
The most important part of the metal AM process modeling and prediction is the prediction of the temperature induced by laser since the non-uniform temperature will cause the thermal stress to appear in the structure. As a result of thermal stress in the build material, the tensile residual stress on the surface accelerates the crack propagation and growth [16,17]. Several researchers worked on predicting the temperature profile during the additive manufacturing process. Fergani et al. introduced an analytical model to predict the temperature in the direct metal deposition process. They predict the temperature using a moving point heat source analysis. In this work, the effect of material temperature sensitivity is ignored [18]. C.Y. Yap et al. have proposed an analytical model to predict the energy input required to process different metallic materials for selective laser melting (SLM) process. The model holds many assumptions, such as a semi-circular cross-section for melt tracks, temperature-independent specific heat, no heat loss to the surroundings and absorptance of material to laser irradiation based on bulk material properties. The melting, solidification, and solid-state phase change is also not considered in their model. The simplified model is able to predict the required energy input within an order of magnitude and provide researchers with a useful model to estimate the optimal SLM parameters [19].
Predicting the temperature precisely in metal AM is the pillar for predicting the thermal stress, residual stress, and part distortion. The non-uniform heating during AM processes may lead to the thermal stress. The large thermal gradient and cooling rate during the metal AM processes can generate complex microstructures in the build material [20]. Kelly et al. used the temperature in the AM processes in order to predict the microstructure evolution in the build part. In their work, the melting/solidification phase change is not considered [21]. Hoadley and Rappaz introduced a 2D quasi-stationary model to predict the temperature in the laser cladding process. Their research focused on the influence of the laser speed and power on the layer thickness [22]. Toyserkani et al. developed a 3D model, their proposed model tried to solve the heat problem using a coupled multi-physics system. They have used thermal analysis in order to predict the melt pool shape [23]. Cao and Ayalew have developed a control-oriented multiple input multiple output modeling of the laser-aided powder deposition processes. The objective of their work is to control the height and the temperature of a layer. Their investigation described the essential role of temperature modeling to control the quality of the final part [24]. Hitzler et al. investigated the influence of scan strategy on material characteristics, such as strength, hardness, and young's modulus [25,26]. Rashid et al. worked on the effect of scan strategy on density and metallurgical properties of a build part during the selective laser melting (SLM) process. Their results showed that parts which are made using a single scan have higher levels of hardness than parts that are made by scanning each layer twice [27].
Due to the complexity of the additive manufacturing processes, such as direct metal deposition (DMD), and SLM, not only is it time-consuming to do the experiments in order to capture the physical aspects of the metal AM processes, but it is also expensive. In the past few decades, the numerical simulations appear to be the only effective way to achieve an understanding of metal additive manufacturing processes [28,29]. The numerical methods have low computational efficiency and cannot capture all the physical aspects of the metal AM processes. On the other hand, physics-based analytical models provide a deep understanding of the physical concepts of AM. The analytical solutions have the potential to predict the key AM attributes in ways significantly faster than finite element method (FEM) simulations, by two or more orders of magnitudes [30]. Efficient and accurate predictions are therefore enabled, and the optimization of metal additive manufacturing processes which would be too complicated to cope with by the majority of other studies, who have resorted to empirical and FEM attempts. It also reduces, if not completely eliminates, the need for a costly and lengthy trial and error developmental curve for new material and components [31]. A complete build analysis with high accuracy becomes computationally tractable using the analytical model.
The AM process is a coupling of many physical phenomena such as heat transfer, fluid dynamics, phase transformation and solid mechanics. Moreover, the transient nature of heat transfer phenomena and interaction of layers make it a complicated multi-physics problem. Many researchers tried to predict the temperature in metal additive manufacturing, but each of them has several limitations. For example, not considering the temperature dependent material properties, the melting/solidification phase change, and layering aspect of metal AM. The key advantage of the proposed model is the ability to capture the most physical phenomena in metal AM, which has mostly been ignored in previous works. In this work, all the above-mentioned limitations are considered in the analytical solution of temperature. It is assumed that the thermal properties of material are temperature dependent. The melting, solidification, and solid-state phase change is included by using the modified specific heat, which relates the specific heat and latent heat of fusion. As each layer is deposited, the temperature profile is predicted using the moving heat source analysis. The laser thermal energy deposited into a control volume is absorbed by the material thermodynamic latent heat and conducted through the contacting solid boundaries. The deposited energy on the first layer introduces a thermal profile. The thermal behavior in the second pass of the laser will not be the same as the first pass since the thermal interaction of the successive layers have an influence on heat transfer. The melt pool geometry is well captured based on the proposed model, since it considers most of the previous lacks.
The outline of the paper is as follows. Section 2 presents the mathematical and practical details of the proposed analytical and numerical models. Section 3 presents detail of the experimental work which is used for validation of the proposed model, the results and a detailed discussion about the obtained results. Last but not least, Section 4 presents the conclusion of this research.

Analytical Modeling
There are various engineering applications, such as turning, grinding, welding, and 3D printing in which the computation of the temperature field in the solid is modeled as a problem of heat conduction involving a moving heat source. The objective of this section is to present the mathematical formulation and the method of solution of heat conduction by considering the moving heat source, which indeed is the case in metal additive manufacturing [32].
In this study, the basic premise is that the powder is situated in the desirable location relative to the melt pool. In other words, there is no moment or mass transfer consideration in this work, and only the heat transfer is considered. Although the effect of time difference between the two consecutive irradiations on temperature profile is not considered in this work, it is worth noting that considering the existence of the time difference between two consecutive irradiations may cause an increase in predicted temperature during the metal AM processes. This is because the predicted temperature at time t + ∆t will be the materials-response-coupled superposition considering the temperature sensitivity of thermal properties at time t and t + ∆t.
By considering a line heat source of constant strength g c l (W/m) located at the x-axis and oriented parallel to the z-axis, the source releases its energy continuously over time as it moves with a constant velocity of v in the positive x-direction. The medium is initially at room temperature. It is assumed ( ∂T /∂z) = 0 everywhere in the medium. Hence, the differential equation of heat conduction in the x, y coordinates in now taken as where T ≡ T(x, y, t). k is thermal conductivity, and α is the thermal diffusivity. The line heat source g c l (W/m) is related to the equivalent volumetric source g(x, y, t) (W/m 3 ) by the delta function notation as In order to consider the moving heat source, it is assumed that the coordinate system transfers from the x, y fixed coordinate system to ζ, y coordinate moving with the line heat source by using the transformation Using the abovementioned transformation, the heat conduction equation for the moving coordinate system can be written as Equation (4) can be solved by the assumption of the quasi-stationary condition [32]. Using the separation of variables, the closed form solution of the temperature field can be obtained as where P is the laser power, K is thermal conductivity, v is scan speed (laser velocity). α is thermal diffusivity which can be calculated as K ρc . In which ρ is material density and c is material heat capacity. R = x 2 + y 2 is the radial distance from the heat source. T 0 is the initial temperature. The material is considered homogeneous and isotropic. As the laser moves along the surface it deposited some energy. Figure 1 depicts the heat transfer in metal AM. The heat loss from the surface by radiation and convection are not considered in this study. It is worth noting that the process parameters such as laser power, scanning speed, powder size, powder distribution, etc. have influence on material properties in the metal AM processes since it may change the predicted temperature [33]. As a result, the material properties are assumed to be temperature dependent as shown in Table 1 [34]. During the metal AM process such as SLM and DMD, the melting, solidification and solid-state phase transformation take place. This is considered using modified heat capacity.
In which ( ) is temperature dependent specific heat, is latent heat of fusion, and is liquid fraction which can be calculated from where, is solidus temperature and is liquidus temperature. It is worth noting that the process parameters such as laser power, scanning speed, powder size, powder distribution, etc. have influence on material properties in the metal AM processes since it may change the predicted temperature [33]. As a result, the material properties are assumed to be temperature dependent as shown in Table 1 [34]. Table 1. Material properties of Ti-6Al-4V.

Properties
Ti-6Al-4V Liquidus temperature (K) 1928 Solidus temperature (K) 1878 Thermal conductivity (W/m K) Density (Kg/m 3 ) 4420 During the metal AM process such as SLM and DMD, the melting, solidification and solid-state phase transformation take place. This is considered using modified heat capacity.
In which C p (T) is temperature dependent specific heat, L f is latent heat of fusion, and f is liquid fraction which can be calculated from where, T s is solidus temperature and T L is liquidus temperature.
The process parameters such as laser power and scan speed are defined to start the calculation of the vertical distribution of temperature during the laser-based metal AM. At first, it is assumed the powder is at room temperature. As the laser moves along the x-axis, it deposits the energy on the powder and causes the powder to melt, as the laser passes the affected region, the melt pool starts to solidify. As it creates the first layer, the temperature profile is calculated for that layer. Next, it starts the second layer with the dwell time of zero. It is possible that the first layer has not had enough time to cool down to the room temperature when the second layer is starting to build. As a result, it affects the heat transfer during the metal AM processes. Considering the layer addition also has a substantial influence on thermal stress and residual stress predictions. The fellow chart of considering the build layers is illustrated in Figure 2.
The process parameters such as laser power and scan speed are defined to start the calculation of the vertical distribution of temperature during the laser-based metal AM. At first, it is assumed the powder is at room temperature. As the laser moves along the x-axis, it deposits the energy on the powder and causes the powder to melt, as the laser passes the affected region, the melt pool starts to solidify. As it creates the first layer, the temperature profile is calculated for that layer. Next, it starts the second layer with the dwell time of zero. It is possible that the first layer has not had enough time to cool down to the room temperature when the second layer is starting to build. As a result, it affects the heat transfer during the metal AM processes. Considering the layer addition also has a substantial influence on thermal stress and residual stress predictions. The fellow chart of considering the build layers is illustrated in Figure 2.

Numerical Modeling
For further validation of this work, finite element analysis is used. The temperature profile is modeled using a moving heat source analysis. The user defined functions (UDF) code is written in ANSYS Fluent software using Equation (8) in order to run a FEA on a 2D geometry, as shown in Figure 3. The build part material is Ti-6Al-4V. The heat loss from the surface due to conduction and radiation is considered. The material properties are assumed to be temperature dependent as shown in Figure 4.

Numerical Modeling
For further validation of this work, finite element analysis is used. The temperature profile is modeled using a moving heat source analysis. The user defined functions (UDF) code is written in ANSYS Fluent software using Equation (8) in order to run a FEA on a 2D geometry, as shown in Figure 3. The build part material is Ti-6Al-4V. The heat loss from the surface due to conduction and radiation is considered. The material properties are assumed to be temperature dependent as shown in  The laser power distribution on the laser beam focus plane is described by the Gaussian equation [32], as where, is the total laser power input, is laser spot radius, is scanning speed is gaussian shape factor, and is a numerical parameter used to fit the experimental data. It accounts for the absorptivity of the material, the heat lost to the metal powder before it falls into the melt pool and the angle of the surface with the laser beam. The values of the process parameters are listed in Table 2. The melting temperature of Ti-6Al-4V is in the range of 1538-1649 °C. In this study, 1620 is selected as the melting temperature of Ti-6Al-4V. The two-dimensional heat transfer in a rectangular surface could be described by where is material density, is specific heat, is thermal conductivity, and is the heat sink. The boundary condition on the laser heating surface is defined as where ( , ) is laser power input, ℎ is the heat transfer coefficient, is the thermal radiation coefficient, is the material emissivity, is the ambient temperature. The initial condition could be as  The laser power distribution on the laser beam focus plane is described by the Gaussian equation [32], as where, is the total laser power input, is laser spot radius, is scanning speed is gaussian shape factor, and is a numerical parameter used to fit the experimental data. It accounts for the absorptivity of the material, the heat lost to the metal powder before it falls into the melt pool and the angle of the surface with the laser beam. The values of the process parameters are listed in Table 2. The melting temperature of Ti-6Al-4V is in the range of 1538-1649 °C. In this study, 1620 is selected as the melting temperature of Ti-6Al-4V. The two-dimensional heat transfer in a rectangular surface could be described by where is material density, is specific heat, is thermal conductivity, and is the heat sink. The boundary condition on the laser heating surface is defined as where ( , ) is laser power input, ℎ is the heat transfer coefficient, is the thermal radiation coefficient, is the material emissivity, is the ambient temperature. The initial condition could be as The laser power distribution on the laser beam focus plane is described by the Gaussian equation [32], as where, P is the total laser power input, r is laser spot radius, v is scanning speed B is gaussian shape factor, and D is a numerical parameter used to fit the experimental data. It accounts for the absorptivity of the material, the heat lost to the metal powder before it falls into the melt pool and the angle of the surface with the laser beam. The values of the process parameters are listed in Table 2.
The melting temperature of Ti-6Al-4V is in the range of 1538-1649 • C. In this study, 1620 is selected as the melting temperature of Ti-6Al-4V. The two-dimensional heat transfer in a rectangular surface could be described by where ρ is material density, C is specific heat, k is thermal conductivity, and S is the heat sink. The boundary condition on the laser heating surface is defined as where q(x,y) is laser power input, h is the heat transfer coefficient, σ is the thermal radiation coefficient, ε is the material emissivity, T 0 is the ambient temperature. The initial condition could be as In this section, the temperature profile, maximum temperature and surface temperature are predicted and compared to the experimental results. A moving heat source analysis is used in order to predict the temperature distribution associated with the dynamic heat deposition. The explicit and closed-form temperature solutions are calculated in Section 2.1. The general differential equation of heat conduction in the 2D plane is used. In order to consider the moving heat source, it is assumed that the coordinate system moves with the heat source by using a transformation as shown in Equation (3). Finally, using the separation of variables, the closed-form solution of temperature is obtained in Equation (5). The material properties are assumed to be temperature dependent. The melting/solidification phase change is also considered. The analytical and numerical analysis are conducted in this work.
In order to validate the proposed model, the experimental temperature data are used from the work of Pauzet [35]. The Ti-6Al-4V samples are manufactured using the DMD machine. The dimensions of the samples are 2 mm in width, 70 mm in depth and 80 mm in length. The temperature on the build part surface is measured using the thermocouple of type K. In order to control the experimental setup, the authors used a thermal-camera and a high-speed camera to provide comparison bases for the temperature and the melt-pool size. The DMD machine has used the laser with the wavelength of 1030 nm. The scanning speed of 0.2 m/min and 0.4 m/min and the laser power of 400 W and 600 W are studied. The initial temperature of each layer depends on the final temperature of the previous layer, as the process is multi-layered. Figure 5 shows the temperature profile of the build part. The temperature is predicted using both the analytical model and the numerical model. The laser moves along the x-axis from left to right. The small red spot on top shows the laser location. The layer thickness is chosen to be 80 µm. The distance of the laser from the powder is 0.4 mm. For the same power, as the velocity is increased the maximum temperature is decreased since the powder has less time to absorb the energy. Different combinations of the process parameters are presented in Figure 5, specifically scanning speed and laser power. The evolution of the surface temperature is plotted as a function of time for each case as shown in Figure 6. A study point will be chosen from the 2D geometry. When the laser is far away from the study point, the powder is at room temperature. As the laser approaches the study point, the temperature increases continuously. The maximum temperature on the curve corresponds to the moment that the laser is above the study point. After the laser passes the point, the temperature is decreased which shows that the material is cooling down. As shown in these plots, the cooling rate in the AM process is substantially high. The evolution of the surface temperature is plotted as a function of time for each case as shown in Figure 6. A study point will be chosen from the 2D geometry. When the laser is far away from the study point, the powder is at room temperature. As the laser approaches the study point, the temperature increases continuously. The maximum temperature on the curve corresponds to the moment that the laser is above the study point. After the laser passes the point, the temperature is decreased which shows that the material is cooling down. As shown in these plots, the cooling rate in the AM process is substantially high. In order to understand the influence of the process parameters on the maximum temperature, and surface temperature, a sensitivity study is designed to investigate both the scan speed and laser power. The short computational time associated with the analytical modeling approach allows for a better understanding of the influence of the process parameters as discussed previously. Figure 7 depicts the influence of the scan speed and laser power on temperature, as predicted by the analytical model and compared to the experimental results.
The results of the simulations from the analytical model illustrates that the maximum temperature decreases linearly as the scan speed increases since the material has less time to absorb the energy. On the other hand, for the fixed scanning speed, as the power increases the maximum temperature increases. The four experimental data are also pointed in Figure 7. The predicted temperature from the analytical model is slightly higher than the experimental values. This error is mainly because the temperature is measured using thermocouples which are a little below the surface. In order to understand the influence of the process parameters on the maximum temperature, and surface temperature, a sensitivity study is designed to investigate both the scan speed and laser power. The short computational time associated with the analytical modeling approach allows for a better understanding of the influence of the process parameters as discussed previously. Figure 7 depicts the influence of the scan speed and laser power on temperature, as predicted by the analytical model and compared to the experimental results.
The results of the simulations from the analytical model illustrates that the maximum temperature decreases linearly as the scan speed increases since the material has less time to absorb the energy. On the other hand, for the fixed scanning speed, as the power increases the maximum temperature increases. The four experimental data are also pointed in Figure 7. The predicted temperature from the analytical model is slightly higher than the experimental values. This error is mainly because the temperature is measured using thermocouples which are a little below the surface.  Figure 8 represents the influence of the laser on the surface temperature. As the power increases from 200 W to 600 W, the surface temperature increases for a fix scanning speed. On the other hand, the surface temperature will decrease as the scanning velocity increases from 0.1 m/min to 0.6 m/min for a fix laser power as shown in Figure 9.   Figure 8 represents the influence of the laser on the surface temperature. As the power increases from 200 W to 600 W, the surface temperature increases for a fix scanning speed. On the other hand, the surface temperature will decrease as the scanning velocity increases from 0.1 m/min to 0.6 m/min for a fix laser power as shown in Figure 9.  Figure 8 represents the influence of the laser on the surface temperature. As the power increases from 200 W to 600 W, the surface temperature increases for a fix scanning speed. On the other hand, the surface temperature will decrease as the scanning velocity increases from 0.1 m/min to 0.6 m/min for a fix laser power as shown in Figure 9.   As explained, the proposed model considers the multi-layer aspects of metal additive manufacturing. The effect of considering the layer addition on peak temperature is compared to the obtained peak temperature without considering the layer addition, and also compared to the experimental results.
To further validate the proposed model, the peak temperature is plotted as a function of scanning speed for different laser powers. Two different values of laser power (400 W and 600 W) and scanning speed (0.2 m/min and 0.4 m/min) are chosen. The temperature considering the layer addition, the temperature not considering the layer addition, and also experimental values are compared. The values are listed in Table 3. The observations show that considering layer addition improves the prediction of temperature, as shown in Figure 10. For example, the predicted temperature for scanning velocity of 0.2 m/min and laser power of 400 W without considering the layer addition is 2042 °C, but when considering the layering aspect of AM, the predicted temperature reduces to 1802.8 °C which shows that it affects the heat transfer mechanisms.  As explained, the proposed model considers the multi-layer aspects of metal additive manufacturing. The effect of considering the layer addition on peak temperature is compared to the obtained peak temperature without considering the layer addition, and also compared to the experimental results.
To further validate the proposed model, the peak temperature is plotted as a function of scanning speed for different laser powers. Two different values of laser power (400 W and 600 W) and scanning speed (0.2 m/min and 0.4 m/min) are chosen. The temperature considering the layer addition, the temperature not considering the layer addition, and also experimental values are compared. The values are listed in Table 3. The observations show that considering layer addition improves the prediction of temperature, as shown in Figure 10. For example, the predicted temperature for scanning velocity of 0.2 m/min and laser power of 400 W without considering the layer addition is 2042 • C, but when considering the layering aspect of AM, the predicted temperature reduces to 1802.8 • C which shows that it affects the heat transfer mechanisms.  A comparison is also conducted among the analytical model considering the layer addition and dwell time, numerical model and experimental values as shown in Figure 11. Overall, the temperature on the surface in terms of magnitude is well captured by both analytical and numerical approaches. The analytical model better approached the experimental measurements. This comparison shows the capability to accurately predict the temperature profile on the surface using the analytical modeling. The analytical approach also provides the power of a short computational time.
In order to illustrate the importance of considering the temperature dependent material properties, a sensitivity analysis is conducted to compare the predicted surface temperature with and A comparison is also conducted among the analytical model considering the layer addition and dwell time, numerical model and experimental values as shown in Figure 11. A comparison is also conducted among the analytical model considering the layer addition and dwell time, numerical model and experimental values as shown in Figure 11. Overall, the temperature on the surface in terms of magnitude is well captured by both analytical and numerical approaches. The analytical model better approached the experimental measurements. This comparison shows the capability to accurately predict the temperature profile on the surface using the analytical modeling. The analytical approach also provides the power of a short computational time.
In order to illustrate the importance of considering the temperature dependent material properties, a sensitivity analysis is conducted to compare the predicted surface temperature with and Overall, the temperature on the surface in terms of magnitude is well captured by both analytical and numerical approaches. The analytical model better approached the experimental measurements. This comparison shows the capability to accurately predict the temperature profile on the surface using the analytical modeling. The analytical approach also provides the power of a short computational time.
In order to illustrate the importance of considering the temperature dependent material properties, a sensitivity analysis is conducted to compare the predicted surface temperature with and without considering the property's temperature-sensitivity. The obtained results demonstrate a significant difference between them as shown in Figure 12. The thermal conductivity of the Ti-6Al-4V is 6.7 W/m· • C which results in a low rate of heat transfer in the build part. However, the thermal conductivity of Ti-6Al-4V varies from 6 to 35 W/m· • C with respect to temperature. The increase in heat transfer rate induced by the increase in thermal conductivity, causes the predicted surface temperature decrease. In the cases that the temperature sensitivity of the material properties is considered, as the velocity increases from 0.2 m/min to 0.4 m/min, the variation of predicted surface temperature is less than 100 • C. However, when the temperature sensitivity of material properties is not considered, the variation of temperature is more than 1000 • C. As it is shown in Figure 12 the predicted temperature can be quite unrealistic without considering the material properties sensitivity to temperature. without considering the property's temperature-sensitivity. The obtained results demonstrate a significant difference between them as shown in Figure 12. The thermal conductivity of the Ti-6Al-4V is 6.7 W/m·°C which results in a low rate of heat transfer in the build part. However, the thermal conductivity of Ti-6Al-4V varies from 6 to 35 W/m·°C with respect to temperature. The increase in heat transfer rate induced by the increase in thermal conductivity, causes the predicted surface temperature decrease. In the cases that the temperature sensitivity of the material properties is considered, as the velocity increases from 0.2 m/min to 0.4 m/min, the variation of predicted surface temperature is less than 100 °C. However, when the temperature sensitivity of material properties is not considered, the variation of temperature is more than 1000 °C. As it is shown in Figure 12 the predicted temperature can be quite unrealistic without considering the material properties sensitivity to temperature.

Experimental Validation Based on Melt Pool Geometry Measurement
In order to predict the morphology of the manufactured part, most of the researchers have used FEA or empirical models [36,37]. The proposed analytical model is used to predict the melt pool size. A comparison between the model and experimental results are conducted. The different process parameters such as laser power and scanning speed are used to predict the melt pool geometry. Figure 13 shows the experimental measurement of melt pool size from Peyre [38]. In this experiment, a high-speed C-Mos camera (Fastcam Photron) is used to measure the melt pool size which is generated by the DMD process. Figure 14 demonstrates the predicted melt pool size and geometry for different process parameters in metal AM. The laser distance from powder is 1 μm.

Experimental Validation Based on Melt Pool Geometry Measurement
In order to predict the morphology of the manufactured part, most of the researchers have used FEA or empirical models [36,37]. The proposed analytical model is used to predict the melt pool size. A comparison between the model and experimental results are conducted. The different process parameters such as laser power and scanning speed are used to predict the melt pool geometry. Figure 13 shows the experimental measurement of melt pool size from Peyre [38]. In this experiment, a high-speed C-Mos camera (Fastcam Photron) is used to measure the melt pool size which is generated by the DMD process. Figure 14 demonstrates the predicted melt pool size and geometry for different process parameters in metal AM. The laser distance from powder is 1 µm. As shown in Figure 14, the melt pool depth and length are obtained using the analytical solution of temperature that is given in Section 2.1. The maximum error in length and depth is 7.6% and 3.7%, respectively. Table 4 listed the process parameters, predicted melt pool size, the experimental values, and also the corresponding error. Based on the calculated error, it is shown that the proposed 2D model can accurately capture the melt pool size. As a result, it eliminates the needs for doing costly experiments and also time-consuming FEM.   As shown in Figure 14, the melt pool depth and length are obtained using the analytical solution of temperature that is given in Section 2.1. The maximum error in length and depth is 7.6% and 3.7%, respectively. Table 4 listed the process parameters, predicted melt pool size, the experimental values, and also the corresponding error. Based on the calculated error, it is shown that the proposed 2D model can accurately capture the melt pool size. As a result, it eliminates the needs for doing costly experiments and also time-consuming FEM.   Figure 14, the melt pool depth and length are obtained using the analytical solution of temperature that is given in Section 2.1. The maximum error in length and depth is 7.6% and 3.7%, respectively. Table 4 listed the process parameters, predicted melt pool size, the experimental values, and also the corresponding error. Based on the calculated error, it is shown that the proposed 2D model can accurately capture the melt pool size. As a result, it eliminates the needs for doing costly experiments and also time-consuming FEM.

Conclusions
Analytical models and numerical model are used to predict the temperature in laser-based metal additive manufacturing configurations of either direct metal deposition or selective laser melting.
In the past few decades, many researchers have been trying to understand the relationships between the process parameters and temperature using FEM. The numerical methods have low computational efficiency, and it cannot capture all the physical aspects of the metal AM processes. The lack of a physics-based analytical model that captures all the physical phenomena of the AM processes is sensible. The physics-based analytical modeling provides accurate results. The high computational efficiency and easy implementation are the other advantages of the analytical model for the additive manufacturing modeling.
In this work, an analytical model is proposed to predict the distribution of the temperature profile by considering the interaction of the layers during the laser metal additive manufacturing process. The material properties are assumed to be temperature dependent, and also the melting/solidification phase change is considered in this work. The temperature profile, the peak temperature, and the evolution of surface temperature are obtained from the proposed model. The analytical model of the temperature is based on the moving heat source assumption, as described in Section 2. The general differential equation of heat conduction is used to obtain the closed-form temperature solution using the separation of variables in a semi-infinite medium. The material is assumed to be homogeneous and isotropic. The predicted temperature from the analytical model are compared with the experimental values and FEM results. For further validation, a comparison of peak temperature considering the layer addition and without considering the layer addition is conducted and compared with experimental values.
The results of the temperature distribution considering the layering aspects of metal additive manufacturing showed better agreement with experimental values in comparison with the predicted temperature not considering the layer addition. The observations suggest that for a fixed laser power, the laser speed increases as the temperature decreases, since the material has less time to absorb the energy. Also, for a given scanning speed, the laser power increases as the maximum temperature increases.
A numerical model is also used to predict the temperature in the metal additive manufacturing process. The material properties are assumed to be temperature dependent. In the numerical model, the heat loss due to convection and radiation is considered. The temperature is well obtained using numerical models.
A comparison is conducted in order to capture the effect of considering the temperature sensitivity of material properties. A significant difference is observed between them. The main reason is that the thermal conductivity of Ti-6Al-4V is low, so the heat transfer rate decreases and causes the surface temperature to increase. However, by considering the temperature dependent material properties, the thermal conductivity increases by increasing the temperature. As a result, the heat transfer rate increases and causes the obtained surface temperature to be less compared to the case in which the thermal conductivity is constant.
The proposed model can also predict the melt pool size with the error margin being less than 7.6%. Hence, it eliminates the costly experiments and time-consuming FEM for predicting the melt pool size. This 2D model also shows that there is no need for doing 3D simulations in order to predict the melt pool size and geometry. The proposed analytical model shows a good agreement with the experimental values. The proposed analytical model reduces the computational time to a fraction when compared to finite element analysis. The analytical model has also eliminated the costly experiments in order to understand the physical concepts of laser metal additive manufacturing. The influence of scanning speed and laser power on the temperature profile, surface temperature, and also peak temperature are investigated and the relations between them are established.