## Abstract

Powder-bed fusion (PBF) process is a subdivision of additive manufacturing (AM) technology where a heat source at a controlled speed selectively fuses regions of a powder-bed material to form three-dimensional (3D) parts in a layer-by-layer fashion. Two of the most commercialized and powerful PBF methods for fabricating full-density metallic parts are the laser PBF (L-PBF) and electron beam PBF (E-PBF) processes. In this study, a multiphysics-based 3D numerical model is developed to compare the thermo-fluid properties of Ti-6Al-4V melt pools formed by the L-PBF and E-PBF processes. The temperature-dependent properties of Ti-6Al-4V alloy and the parameters for the laser and electron beams are incorporated in the model as the user-defined functions (UDFs). The melt-pool geometry and its thermo-fluid behavior are investigated using the finite volume (FV) method, and results for the variations of temperature, thermo-physical properties, velocity, geometry of the melt pool, and cooling rate in the two processes are compared under similar irradiation conditions. For an irradiance level of 26 J/mm3 and a beam interaction time of 1.212 ms, simulation results show that the L-PBF process gives a faster cooling rate (1. 5 K/μs) than that in the E-PBF process (0.74 K/μs). The magnitude of liquid velocity in the melt pool is also higher in L-PBF than that in E-PBF. The numerical model is validated by comparing the simulation results for the melt-pool geometry with the PBF experimental results and comparing the numerical melt-front position with the analytical solution for the classical Stephan problem of melting of a phase-change material (PCM).

## 1 Introduction

Powder-bed fusion (PBF) process is a relatively new additive manufacturing (AM) technology where the thermal energy of a computer-controlled heat source is used for selective melting and sintering of regions of a powder-bed [1]. Once a layer of an object is completed, the building platform is lowered, and more powder is spread over (usually, rolled on) the build area for a new scan. The process ends with a postprocessing step of removing all the unbound powder [2] from the fabricated object. The conventional manufacturing technologies (e.g., casting and forging) used for fabricating medical implants and components for automotive, aerospace, and space applications constrain the customization of complex geometries and consume a significant amount of material and time. The PBF process overcomes these limitations by providing the advantage of cost-effective customization with reduced assembly [2] and thus, becomes a superior AM technology in the present era. Two of the most common types of PBF processes are the laser PBF (L-PBF) and electron beam PBF (E-PBF) processes, which have brought about a revolution in the field of metal AM technology. The L-PBF process uses finely focused monochromatic coherent photons, i.e., laser, while the E-PBF process uses a beam of electrons as the heat source for melting the powder-bed. Consequently, when the laser or the electron beam scans the top surface, it melts a selective region of the powder-bed to form a liquid volume, known as the “melt pool,” which is rapidly cooled and solidified either in a vacuum (for E-PBF) or an inert gas (for L-PBF) environment [1,2]. A wide variety of materials, including the alloys of titanium, copper, nickel, aluminum, and chromium, are being used in PBF technology where the material properties significantly influence the process behavior and quality of the processed part. Titanium alloy, especially, Ti-6Al-4V, is one of the prime choices while selecting PBF raw material because of its unique properties and extensive usability in aerospace, automotive, microelectronics, and biomedical applications.

In recent years, significant growth of interest in AM technologies, especially, in the PBF processes has been marked by numerous studies on the L-PBF and E-PBF methods [39]. In most cases, the focus has been on the application-based experimental studies such as build-part microstructure, powder metallurgy, morphological features, and mechanical properties of the material [411]. Besides these experimental analyses, numerical modeling is also conducted by many to determine the thermo-physical properties and melt-pool dynamics [12,13] for investigating the process envelope, part quality, reliability, and energy efficiency. However, studies on the comparison between the two processes under equivalent process parameters are very limited [13]. Instead of studying the L-PBF and E-PBF methods separately, a comparative study of the process parameters, melt-pool geometry, and part microstructures offers more valuable information to characterize suitable applications of the processes [13]. Numerical modeling, especially, thermal modeling is, by far, one of the most convenient and cost-effective methods to conduct a robust comparative analysis between the L-PBF and E-PBF processes [13,14]. The comparison can provide an important guide to select the appropriate technology to be used in the AM industry.

Development of the thermal models for L-PBF and E-PBF requires the understanding of complex heat transport, material phase change, and intricate relations among the thermal, mechanical, and metallurgical phenomena [13,14] which make it extremely challenging to implement. The E-PBF process requires preheating of the material at a high vacuum that needs to be characterized accurately in the simulation. On the other hand, modeling of L-PBF requires the incorporation of convection with the inert gas environment at room temperature. Consequently, understanding the correlation between the process parameters and the process outcomes without costly experimentation requires comprehensive numerical modeling. While developing a robust thermal model, it is important to find a convenient numerical scheme that can accurately estimate the melt-pool geometry and determine the temperature distribution in the processed part by taking into consideration the heat source parameters and material properties. Many researchers have developed thermal models using finite difference (FD) and finite element (FE) methods at various length and time scales [1528]. As a common practice, the heat source in a PBF process is modeled as a conical volumetric heat flux under the surface of the powder-bed due to the resultant keyhole formed by an incident laser or electron beam. Qi et al. [15] used a controlled-volume FD method to develop a self-consistent model for studying the heat transfer, phase change, and fluid flow within the melt pool in the laser-based PBF process. Moraitis and Labeas [16] developed a thermo-mechanical FE model to investigate the residual stresses and distortions for aluminum lap joints in the laser beam welding process. Wang et al. [17], Roberts et al. [18], Liu et al. [19], Yang et al. [20], Lacki and Adamus [21], Shen and Chou [22], Cheng et al. [23], Chen et al. [24], and Zäh and Lutzmann [25] developed FE models incorporating the Gaussian heat source, powder porosity, and thermal properties to simulate the transient heat transfer in the PBF processes. Recent progress on FE analyses to investigate the effect of process parameters in PBF can be attributed to the studies of Andreotta et al. [26], Sadowski et al. [27], and Ladani et al. [28]. All these studies included numerical and/or experimental analyses of either the L-PBF or E-PBF process, but the comparison between the two processes, which could facilitate the selection of the suitable one in the industry, was not outlined.

Thermo-fluid models based on the computational fluid dynamics (CFD) and finite volume (FV) method become more effective than the FE thermal models when fluid flow and heat convection in the melt pool are dominant factors in the process outcomes [29]. Studies show that thermo-fluid models using CFD can effectively provide quantitative information about the part geometry, thermal cycle, cooling rate, and solidification process with the same accuracy as the FE models [29,30]. Wang et al. [30] developed a three-dimensional (3D) volume-of-fluid method to measure the real-time melt-pool shape and obtained the distribution of temperature in laser keyhole welding. Cho et al. [31] studied the effect of fluid motion in the melt pool using thermo-fluid simulation. Rai et al. [32,33] and Li et al. [34] showed that fluid convection inside the melt pool resulted in an increase in heat transfer and gave a better correlation between numerical and experimental results of the melt-pool geometry. Chahine [35] studied the effects of the current and exposure time of an electron beam on the temperature distribution and fluid flow of a melt spot using the CFD technique. Yuan and Gu [36] used FV simulation and laser experiments to investigate the melt-pool evolution and thermal behavior of TiC/AlSi10Mg powder-bed in the L-PBF process. Rahman et al. [3742] conducted CFD-based thermo-fluid modeling of the Ti-6Al-4V [3741] and Cu-Cr-Zr [42] melt pools to study the thermal features and melt-pool dynamics in the PBF processes. Jamshidinia et al. [43,44] developed 3D thermal and fluid flow models of E-PBF, where the influence of fluid convection on the melt-pool geometry was investigated and the effects of changing process parameters were studied numerically and experimentally. However, a comparison between the L-PBF and E-PBF processes based on thermo-fluid modeling is yet to be studied rigorously.

The study aims to compare the thermo-fluid properties of Ti-6Al-4V alloy while undergoing the L-PBF and E-PBF processes by developing a 3D multiphysics CFD model for each method. The effects of the laser and electron beam parameters on the temperature distribution, melt-pool geometry, melting, and solidification criteria are studied and compared under similar heat source specifications. The laser and the electron beam are considered as the Gaussian heat sources to perform simulations for the L-PBF and E-PBF models separately. The CFD simulations are conducted in ansys fluent 18.2 by setting appropriate solver specifications and utilizing the user-defined functions (UDFs) for the heat source and the material properties of Ti-6Al-4V. Results obtained from the CFD simulations are validated by comparing with analytical and experimental results to corroborate the comparative study of the L-PBF and E-PBF processes.

## 2 Material and Methods

### 2.1 Material Modeling.

Ti-6Al-4V provides a unique combination of physical and mechanical properties including lightweight, high strength-to-weight ratio, and excellent resistance to corrosion and fatigue. However, the thermo-physical properties of metallic powder materials are significantly different from those of the corresponding solid bulk material [45,46], especially, the thermal conductivity, melting point, specific heat capacity, and density. When the phase change occurs from solid to liquid and vice versa, properties of the liquid material also have significant effects on the process outcomes. Table 1 shows the temperature-dependent properties of Ti-6Al-4V in solid, powder, and liquid states [39,40,46] that are used in both the L-PBF and E-PBF models as the UDFs.

Table 1

Thermo-physical properties of Ti-6Al-4V alloy used as the user-defined functions [39,40,46]

PropertiesMaterial stateTemperature range (°C)UDFs
Specific heat capacity, Cp (J/kg K)Powder23 <T <1650Cp = [0.52036–(8.34 × 10−6) T (°C)–(4.46 × 10−7) T2 (°C) + (5.44 × 10−10) T3 (°C)] × 1000
Solid23 <T <1650Cp = [0.54058 + (1.02 × 10−4) T (°C) + (1.35 × 10−7) T2 (°C) –(6.50 × 10−11) T3 (°C)] × 1000
Liquid1650 <T <2700830
Thermal conductivity, k (W/m K)Powder23 <T <1650k = 0.9315 – 0.00339 T (°C) + (6.55 × 10−6) T2 (°C) – (1.41 × 10−9) T3(°C)
Solid23 <T <1650k = 6.95757 + 0.00224 T (°C) + (1.69 × 10−5) T2 (°C)—(7.58 × 10−9) T3(°C)
Liquid1650 <T <2700k = −1.6614 + 0.0183 T (°C)
Emissivity, ɛMelt-pool front23 <T <1650ɛ = 0.43356 + (2.94 × 10−4T) (°C) + (5.48 × 10−7) T2 (°C)–(5.53 × 10−10) T3 (°C)
Powder23 <T <16500.6
Melt-pool top1650 <T <27000.4
Density, ρ (kg/m3)Solid23 <T <1650ρsolid = 4420–0.154 (T–25 °C)
Powder23 <T <1650ρpowder = (1–porosity) × ρsolid
Liquid1650 <T <2700ρliq = 3920–0.68 (T–1650 °C)
Viscosity, μ (Pa · s)LiquidT >1605μ(T) = C eE/RT; with $E=0.431Tl1.348$ where C is a constant, E is the activation energy, Tl is the liquidus temperature, and R is the molar gas constant [47]
PropertiesMaterial stateTemperature range (°C)UDFs
Specific heat capacity, Cp (J/kg K)Powder23 <T <1650Cp = [0.52036–(8.34 × 10−6) T (°C)–(4.46 × 10−7) T2 (°C) + (5.44 × 10−10) T3 (°C)] × 1000
Solid23 <T <1650Cp = [0.54058 + (1.02 × 10−4) T (°C) + (1.35 × 10−7) T2 (°C) –(6.50 × 10−11) T3 (°C)] × 1000
Liquid1650 <T <2700830
Thermal conductivity, k (W/m K)Powder23 <T <1650k = 0.9315 – 0.00339 T (°C) + (6.55 × 10−6) T2 (°C) – (1.41 × 10−9) T3(°C)
Solid23 <T <1650k = 6.95757 + 0.00224 T (°C) + (1.69 × 10−5) T2 (°C)—(7.58 × 10−9) T3(°C)
Liquid1650 <T <2700k = −1.6614 + 0.0183 T (°C)
Emissivity, ɛMelt-pool front23 <T <1650ɛ = 0.43356 + (2.94 × 10−4T) (°C) + (5.48 × 10−7) T2 (°C)–(5.53 × 10−10) T3 (°C)
Powder23 <T <16500.6
Melt-pool top1650 <T <27000.4
Density, ρ (kg/m3)Solid23 <T <1650ρsolid = 4420–0.154 (T–25 °C)
Powder23 <T <1650ρpowder = (1–porosity) × ρsolid
Liquid1650 <T <2700ρliq = 3920–0.68 (T–1650 °C)
Viscosity, μ (Pa · s)LiquidT >1605μ(T) = C eE/RT; with $E=0.431Tl1.348$ where C is a constant, E is the activation energy, Tl is the liquidus temperature, and R is the molar gas constant [47]

### 2.2 Physical Model.

A finite volume based 3D model is developed in ansys fluent for the transient thermo-fluid simulation of the PBF process. The geometry of the 3D PBF model is shown in Fig. 1, where the physical domain consists of a solid rectangular block of Ti-6Al-4V, known as the “substrate” and a layer of Ti-6Al-4V powder on top of the substrate. The dimensions of the substrate and the powder layer are considered as 14 mm × 4 mm × 4 mm and 14 mm × 4 mm × 0.07 mm, respectively [13], which implies that the powder layer thickness is 0.07 mm. The laser, or the electron beam, scans the top surface of the powder-bed in the y-direction. A single-track scan is considered for comparing the L-PBF and E-PBF models. The PBF model is first used to simulate the L-PBF process by assigning the laser as the heat source and setting the boundary and initial conditions for the L-PBF process. The same PBF model is then converted to simulate the E-PBF process by assigning the electron beam as the heat source and setting the appropriate conditions for the E-PBF process. The spot size and scanning speed for both the electron beam and the laser are kept the same (0.4 mm and 330 mm/s, respectively). They scan the top surface of the domain starting from (0, 2 mm, 0) to the endpoint at (0, 12 mm, 0) as shown in Fig. 1.

Fig. 1
Fig. 1
Close modal

The maximum heat of the laser or electron beam is located at the center of a target surface, and the intensity varies radially from the center of the heat source. The powder-bed medium is homogeneous and continuous. The heat source is applied at the top of the powder-bed in the form of heat flux density obeying the Gaussian distribution. The top surface of the melt pool is assumed to be flat and all the nodes remained in both vertical and horizontal positions. During the melting and solidification processes, the thermo-physical properties including density, specific heat, thermal conductivity, and viscosity are considered as temperature-dependent parameters. Heat transfer by convection on the top surface is neglected in E-PBF due to the vacuum environment. The thermal boundary conditions applied in the simulations of L-PBF and E-PBF are similar, but unlike E-PBF, L-PBF included convection heat transfer on the top surface. The convective heat transfer coefficient between the powder-bed and surrounding gas in L-PBF is taken as a constant. The top surface is exposed to radiation heat transfer with an ambient temperature of 298 K for both the L-PBF and E-PBF models. The sidewalls and bottom of the domain in E-PBF are adiabatic at a temperature of 1003 K. As L-PBF does not require preheating, the side walls and bottom of the domain are considered adiabatic at a temperature of 298 K.

## 3 Mathematical Formulation

### 3.1 Governing Equations.

The governing equations of fluid motion and heat transfer generally follow three basic physical conservation laws—the conservations of mass, momentum, and energy. Based on the assumptions stated in Physical Model section, the mass, momentum, and energy conservation equations for a 3D Cartesian coordinate system can be expressed by Eqs. (1) through (3). The continuity equation is given by
$∂ρ∂t+∂(ρui)∂xi=0$
(1)
The conservation of momentum equation is given by
$∂(ρuj)∂t+∂(ρuiuj)∂xi=−∂P∂xj+∂∂xi(μ∂uj∂xi)+∂∂xi(μ∂ui∂xj)−ρgzβ(T−Tref)+ρgz−CM((1−fL)2fL3+B)uj−vs∂(ρuj)∂xi$
(2)
In Eq. (2), the scanning speed vs is in the y-direction and gravitational acceleration gz is in the z-direction. Therefore, only the y-momentum equation contains the last term associated with vs which is the relative motion between the heat source and the work piece. Here, β is the coefficient of volume expansion, fL is the liquid fraction, CM is a constant that accounts for the mushy zone morphology, and B is a very small computational constant introduced to avoid division by zero [33]. The fourth, fifth, and sixth terms in the right side of Eq. (2) represent the buoyancy force, gravity, and the frictional drag in the mushy zone during the solid–liquid–solid transition, respectively. Considering the motion of the heat source vs in y-direction, the transient conservation of energy equation (i.e., the heat equation) in terms of temperature [13,22,23,25] is given by
$∂T∂t+vs∂T∂y=∇⋅(k(T)ρ(T)cp(T)∇⋅T)+Q˙(x,y,z,t)ρ(T)cp(T)+Sgzρ(T)cp(T)+μ(T)ρ(T)cp(T)ΦV$
(3)
where k is the effective thermal conductivity, $Q˙(x,y,z,t)$ is the heat source, $Sgz$ is the source term due to gravity, and ΦV is the viscous dissipation term which is defined by Rahman et al. [39] as follows:
$ΦV=2[(∂u∂x)2+(∂v∂y)2+(∂w∂z)2]+(∂u∂y+∂v∂x)2+(∂v∂z+∂w∂y)2+(∂w∂x+∂u∂z)2−23(∂u∂x+∂v∂y+∂w∂z)2$
(4)
The solid–liquid–solid phase change is encountered by an enthalpy method [39,44,48] where the temperature-dependent enthalpy H [48] is a function of the specific heat capacity cp, liquidus temperature TL, solidus temperature TS, latent heat of fusion Lf, and the liquid fraction fL, and can be defined as
$H(T)=∫0TcpdT+LffL$
(5)
where the liquid fraction fL varies from 0 to 1 according to the following relation:
$fL={0,TTL$
(6)
Equation (6) can be rearranged with the information of Eq. (5) to yield the relation between the enthalpy and temperature as given by Eq. (7) [39]
$T={Hcp,HcpTL+Lf$
(7)

### 3.2 Modeling of the Heat Source.

The laser and electron beams are modeled as conical volumetric heat sources with the Gaussian distribution as shown in Fig. 2, where the maximum power intensity is at the center and the intensity decreases with the increase of the depth and width.

Fig. 2
Fig. 2
Close modal
The 3D conical volumetric Gaussian heat source model considered for this analysis is expressed by Eq. (8) [13,37,38]
$Q˙(x,y,z)=η×Ixy×IZS$
(8)
where the Gaussian surface intensity profile Ixy and the beam penetration function, i.e., vertical intensity profile IZ are given by Eqs. (9) and (10), respectively
$Ixy=2PHπΦ2exp{−2[(x−xs)2+(y−ys)2]Φ2}$
(9)
$Iz=10.75(−2.25(zS)2+1.5(zS)+0.75)$
(10)
The power of the electron beam is given by PH = VIb, where V is the acceleration voltage and Ib is the beam current. The values of efficiency η for the laser and electron beams are usually different and depend on a number of factors including the beam control [25], focusing [49], vacuum [49] or convective environment, inclination angle [50], and the energy absorption by the target material [49,50]. The values of the penetration depth S (also known as absorption depth) for laser and electron beam are also different as their wavelengths are not the same [51]. The value of S for the electron beam can be determined by the following relation:
$S=SE=2.1×10−5V2ρ$
(11)
where SE is the penetration depth of electron beam in μm, V is the electron beam potential in V, and ρ is density of powder-bed in kg/m3 [43]. Using a voltage of 60 kV and a powder density of 2150 kg/m3, the value of SE is found to be 35.16 μm. In case of laser, the parameter S is set to be the optical penetration depth SL which is defined as the depth at which the intensity of the laser drops to 1/e of its initial value at the interface [52], and can be determined by
$S=SL=1a=12.303×Alt$
(12)
where a is the absorption coefficient [5053], lt is the powder layer thickness, and A is the optical absorbance [53] of the laser beam while penetrating the Ti-6Al-4V powder-bed. The absorbance of electron beam in Ti-6Al-4V powder-bed is higher than that of the laser beam because the photons are mostly deflected rather than absorbed into the material [51]. Taking lt = 70 μm and S = 35.16 μm in Eq. (12), the value of A for electron beam is found 0.8645. Considering a solid state yttrium-aluminum garnet doped with neodymium ions (Nd:YAG) laser with a wavelength of 1060 nm, the absorbance of laser beam in Ti-6Al-4V alloy is considered as 0.49 [54] which gives SL = 62 μm for the L-PBF simulations. Figure 3 shows the values of beam penetration function for the static laser and electron beams along the vertical coordinate of the domain where the z values are taken such that 0 ≤ zS. The higher value of S for the laser beam results in deeper distribution of its intensity as compared with the intensity of the electron beam within the specified range of vertical coordinate z.
Fig. 3
Fig. 3
Close modal

### 3.3 Initial and Boundary Conditions.

Since the E-PBF process includes preheating of the entire domain to a temperature of 1003 K, the initial conditions are at t = 0, u = v = w = 0, and T = Tpreheat = 1003 K, in the entire domain. The top surface in the E-PBF model is exposed to radiation only at 298 K while the side walls and the bottom of the substrate are considered as adiabatic surfaces. For the L-PBF model, the initial temperature is 298 K and boundary conditions are same as the E-PBF model except for the top surface, which is exposed to radiation and convection at 298 K. The temperature coefficient of surface tension (Marangoni coefficient) is set as −2.6 × 10−4 N/m K at the top surface for both cases [40,44].

## 4 Numerical Simulation

### 4.1 Simulation Procedure.

Transient thermo-fluid analyses of the 3D configuration are performed numerically using ansys fluent 18.2. The thermal properties and the specifications of the moving heat source are assigned as UDFs to simulate the transient melting and solidification for both the L-PBF and E-PBF models. The mass, momentum, and energy conservation equations are discretized and solved using the control volume method with appropriate boundary conditions. ansys Design Modeler is used to create the geometry, Mesh tool is used to generate the structured mesh, and the mathematical model is followed to define the boundary types of the 3D computational domain. Under the pressure-based solver, the “Coupled” algorithm was selected for solving the conservation equations. The Pressure Staggering Option (PRESTO) is selected for pressure discretization, while density, momentum, and energy are discretized with the second-order upwind method. The transient formulation is done by choosing the first-order upwind method. The Courant number is set equal to 1 and all the residual criteria are set to be 10−5. The time-step size is restricted to be 0.001 s or lower to satisfy the convergence criteria while the maximum iterations per time-step are taken as 50/s for running the calculation. After the solution, the postprocessing steps begin to observe and save the results.

### 4.2 Mesh Convergence Study.

The 3D computational domain considered for the analysis is discretized using a structured mesh with hexahedral cells as shown in Fig. 4. The structured mesh is formed by biasing the grid in the powder layer region and around the scanning path of the moving heat source to have a very fine mesh in the target zone.

Fig. 4
Fig. 4
Close modal

A mesh convergence (i.e., the mesh sensitivity) study is conducted for the structured mesh of the 3D domain considering the variation of melt-pool temperature with the increase of the number of nodes. The temperature at location $(x,y,z)$ = (0, 5 mm, 0.03 mm) is monitored for several different mesh densities at 0.009 s when the beam spot size is 0.4 mm, scanning speed is 330 mm/s, and the effective power is 216 W. The value of temperature inside the melt pool converges to 2571 K with the increase of the number of nodes in the domain for the E-PBF model. Figure 5 depicts the results for the mesh convergence study where the temperature at the fixed point remains unchanged after 200,889 nodes corresponding to 190,040 hexahedral cells. Results for both the L-PBF and E-PBF models are obtained for the converged mesh.

Fig. 5
Fig. 5
Close modal

## 5 Model Validation

### 5.1 Comparison of Modeling Results with Analytical Results.

The results obtained from the Fluent simulation for the melting of pure titanium (Ti) is compared with the analytical solution of the classical Stefan problem [42,55] of the melting of a phase-change material (PCM) with pure conduction. Figure 6 shows the standard geometry for the Stefan problem, where the PCM is semi-infinite and initially (at t = 0) solid at its melting temperature Tm [42,55]. The wall temperature Tw is raised to Tw ˃ Tm for melting the PCM in a linear fashion starting at x = 0. According to the Stefan condition, the solution for the transient temperature distribution in the liquid is given by the following relation [42,55]:
$Tl(x,t)−TwTm−Tw=erf(x/2αlt)erf(λ)$
(13)
Fig. 6
Fig. 6
Close modal
The melting front moves forward in the x-direction as time increases. The position of the melting front s(t), measured from x = 0, is given by
$s(t)=2λαlt$
(14)
where αl = kl/ρlcp,l is the thermal diffusivity of the liquid PCM. The parameter λ is calculated from the interfacial melt-front equation and the Stefan number (Ste) as defined below [42,55]:
$λeλ2erf(λ)=Steπ$
(15)
$Ste=cp,l(Tw−Tm)L$
(16)

The validation of the fluent result is conducted by predicting the motion of the liquid–solid interface during the melting of pure Ti. The parameters shown in Table 2 are used in the ansys fluent simulation.

Table 2

List of the simulation parameters for Ti melting

ParametersValues
Density of liquid Ti, ρl (kg/m3)4500
Specific heat capacity of liquid Ti, cp,l (J/kg K)528
Effective viscosity, µ (kg/m s)4.3 × 10−3
Thermal conductivity of liquid Ti, kl (W/m K)17
Latent heat of fusion, Lf (kJ/kg)435.4
Melting temperature, Tm (K)1923
Wall temperature, Tw (K)2073
Solidus temperature, TS (K)1923
Liquidus temperature, TL (K)1943
ParametersValues
Density of liquid Ti, ρl (kg/m3)4500
Specific heat capacity of liquid Ti, cp,l (J/kg K)528
Effective viscosity, µ (kg/m s)4.3 × 10−3
Thermal conductivity of liquid Ti, kl (W/m K)17
Latent heat of fusion, Lf (kJ/kg)435.4
Melting temperature, Tm (K)1923
Wall temperature, Tw (K)2073
Solidus temperature, TS (K)1923
Liquidus temperature, TL (K)1943

The simulation results of the change in interface position with respect to time during the melting of pure Ti show a good agreement with the analytical results. The liquid fraction contours obtained from the simulation at t = 0.37 s and t = 1.8 s are shown in Fig. 7, where the melt front moves forward with the increase of time.

Fig. 7
Fig. 7
Close modal

Figure 8 shows the comparison between the analytical and simulation results for the melt-front position with respect to time considering x = 1 mm.

Fig. 8
Fig. 8
Close modal

Results for temperature distribution also show a good match between the analytical and simulation results. At x = 1 mm, the comparison between the analytical and simulation results for centerline temperature at three different times is shown in Table 3.

Table 3

Analytical versus numerical results for temperature

Time (s)Temperature (K)
AnalyticalNumericalDeviation (%)
0.51936.971940.500.182
1.01975.701979.450.189
102041.912048.070.300
Time (s)Temperature (K)
AnalyticalNumericalDeviation (%)
0.51936.971940.500.182
1.01975.701979.450.189
102041.912048.070.300

### 5.2 Experimental Validation.

The numerical results for the melt-pool geometry are also validated by comparing with the experimental results. The experimental procedure for the E-PBF process, conducted by Jamshidinia et al. [44] with Ti-6Al-4V, is followed to validate the proposed multiphysics model. Jamshidinia et al. [44] compared the results for the variation of the average melt-pool width and depth with the change in scanning speed. Using a constant electron beam spot size of 0.4 mm, a beam current of 14 mA, and a voltage of 60 kV, they applied three levels of scanning speed, namely, 100 mm/s, 300 mm/s, and 500 mm/s to measure the average melt-pool width and depth. The differences between their modeling results and experimental results ranged from −3.5% to +3% for the melt-pool width, and from +2.1% to +3.5% for the melt-pool depth.

Following a similar approach, the simulated results for the melt-pool geometry from the proposed multiphysics CFD model are compared with the experimental results presented by Jamshidinia et al. [44] using the converged mesh of 190,040 hexahedral cells connected by 200,889 nodes. The comparison gave a good agreement between the simulation results and the experimental results of Jamshidinia et al. with a maximum deviation of 3.73%, which indicates a good accuracy of the proposed multiphysics CFD model. The comparison of results for the melt-pool width and depth is illustrated in Fig. 9.

Fig. 9
Fig. 9
Close modal

## 6 Results and Discussion

Results for the thermo-fluid properties and melt-pool geometry are generated at similar irradiation conditions to make a valid comparison. The simulations for the L-PBF and E-PBF models are conducted using the UDFs and the parameters shown in Table 4.

Table 4

List of the simulation parameters

ParametersValues
Solidus temperature, TS (K)1878
Liquidus temperature, TL (K)1938
Latent heat of fusion, Lf (kJ/kg)440
Spot size of laser or electron beam, Φ (mm)0.4
Scanning speed, vs (mm/s)330
Acceleration voltage, V (kV)60
Electron beam current, Ib (mA)4
Laser power, P (W)240
Preheat temperature in E-PBF, Tpreheat (K)1003
Initial temperature in L-PBF, TL-PBF (K)298
Electron beam efficiency, ηe0.9
Laser absorption efficiency, ηl0.865
Powder porosity (%)50
Powder layer thickness, lt (mm)0.07
Electron beam penetration depth, SE (µm)35.16
Laser beam penetration depth, SL (µm)62
Convective heat transfer coefficient, h (W/m2 K)10
Effective viscosity of liquid, µ (kg/m s)UDF
Specific heat, cp (J/kg K)UDF
Thermal conductivity, k (W/m K)UDF
Emissivity, ɛUDF
Density, ρ (kg/m3)UDF
ParametersValues
Solidus temperature, TS (K)1878
Liquidus temperature, TL (K)1938
Latent heat of fusion, Lf (kJ/kg)440
Spot size of laser or electron beam, Φ (mm)0.4
Scanning speed, vs (mm/s)330
Acceleration voltage, V (kV)60
Electron beam current, Ib (mA)4
Laser power, P (W)240
Preheat temperature in E-PBF, Tpreheat (K)1003
Initial temperature in L-PBF, TL-PBF (K)298
Electron beam efficiency, ηe0.9
Laser absorption efficiency, ηl0.865
Powder porosity (%)50
Powder layer thickness, lt (mm)0.07
Electron beam penetration depth, SE (µm)35.16
Laser beam penetration depth, SL (µm)62
Convective heat transfer coefficient, h (W/m2 K)10
Effective viscosity of liquid, µ (kg/m s)UDF
Specific heat, cp (J/kg K)UDF
Thermal conductivity, k (W/m K)UDF
Emissivity, ɛUDF
Density, ρ (kg/m3)UDF
Results for both the L-PBF and E-PBF models are obtained for the converged mesh having 190,040 hexahedral cells with 200,889 nodes and considering a single scan by the heat source at the top surface of the powder-bed. Both the laser and the electron beam scanning simulations are performed keeping the same spot size Φ = 0.4 mm and the same scanning speed vs = 330 mm/s. The same energy density of ED = 26 J/mm3 and interaction time ti = 1.212 ms are used for both the L-PBF and E-PBF models which are determined by the following relations [56,57]:
$ED=PHvs×lt×Φ$
(17)
$ti=Φvs$
(18)

### 6.1 Temperature and Thermo-Physical Properties.

The heat source scans the top surface in the y-direction in both the L-PBF and E-PBF models. The temperature contours at the top surface for L-PBF and E-PBF are shown in Fig. 10.

Fig. 10
Fig. 10
Close modal

The images for the temperature contours are captured when both the laser and electron beams are at y = 8.0 mm (t = 0.018 s). The melt region is longer in the contour for E-PBF than that in L-PBF. The position y = 8.0 mm at the top surface is chosen to show the complete tailing effect during melting and consolidation in the two processes. A cross section is considered at y = 8.0 mm (t = 0.018 s) along the xz-plane to show the comparison of the results in the form of a two-dimensional representation as shown in Figs. 11 and 12. The temperature contours at the cross section for L-PBF and E-PBF are shown in Fig. 11.

Fig. 11
Fig. 11
Close modal
Fig. 12
Fig. 12
Close modal

As the temperature in the domain increases, density and viscosity decrease following the temperature-dependent equations given in Table 1. The results for density and viscosity are not represented here but they show similar patterns as presented by Rahman et al. in Refs. [39] and [42], respectively. The density of the powder-bed is less than that of the solid substrate due to its porosity. The values of other thermal properties including thermal conductivity, specific heat capacity, and enthalpy increase with the increase of temperature in the domain. Before scanning by the heat source, the thermal conductivity of the powder material is lower than that of the solid substrate as the porosity of the powder reduces the thermal conductivity. When the heat source is applied and the powder material is melted, the thermal conductivity of the liquid melt pool becomes very high due to the high temperature. Contour plots for thermal conductivity in L-PBF and E-PBF are shown in Fig. 12. All the results are obtained at an energy density of 26 J/mm3 and an interaction time of 1.212 ms for both the laser and electron beams. Any increase in energy density due to the increase of power causes a temperature rise in the domain, and eventually affects the temperature-dependent properties significantly.

### 6.2 Melt-Pool Geometry.

The evolution of the melt pool depends on several factors including the material properties, processing parameters, energy absorption, and thermo-fluid interactions. Under the same energy density of 26 J/mm3, the results for the maximum length, width, and depth of penetration of the melt pool at y = 8.0 mm for L-PBF and E-PBF are shown in Table 5.

Table 5

Comparison of the melt-pool dimensions

ProcessED = 26 J/mm3 and ti = 1.212 ms
Length (mm)Width (mm)Depth (mm)
L-PBF1.20.60.08
E-PBF2.10.6050.12
ProcessED = 26 J/mm3 and ti = 1.212 ms
Length (mm)Width (mm)Depth (mm)
L-PBF1.20.60.08
E-PBF2.10.6050.12

The melt-pool volume is larger in E-PBF than that in L-PBF as obtained from the simulation results. The high absorption rate of electron beam, preheating condition, and lack of heat dissipation due to convection result in the formation of a larger volume of the melt pool in E-PBF. Figure 13 depicts the comparison of the length and width of the melt pool when the laser and the electron beam are at y = 8.0 mm.

Fig. 13
Fig. 13
Close modal

A parametric study on the effects of processing parameters on the evolution of the melt pool is also conducted with the numerical simulations of the L-PBF and E-PBF processes. First, the effect of increasing the power of the laser and electron beams on the depth of penetration of the melt pool is investigated while keeping the same spot size of 0.4 mm and scanning speed of 330 mm/s. As expected, the depth of the melt pool increases with the increase of beam power. The comparison of the melt-pool depth in the L-PBF and E-PBF processes are shown in Fig. 14.

Fig. 14
Fig. 14
Close modal

In contrast to the beam power, the depth of the melt pool decreases as the scanning speed increases at a given power and a spot size of the laser or electron beam. The simulation results for melt-pool depth versus scanning speed for L-PBF and E-PBF at a power of 240 W and a spot size of 0.4 mm are shown in Fig. 15. Results show that the melt-pool depth is more sensitive to the change in beam power as compared with the change in scanning speed at a given spot size. For instance, due to an increase of 127.27% of the scanning speed from 330 mm/s to 750 mm/s in the L-PBF model, the percentage of decrease in the melt-pool depth is 55% (as shown in Fig. 15). However, a power increase of 87.5% from 240 W to 450 W in the L-PBF model results in a 433.33% increase in the melt-pool depth, which is calculated from the values shown in Fig. 14. Therefore, variation in the melt-pool depth due to the change in beam power is significantly higher than the variation caused by the change in scanning speed.

Fig. 15
Fig. 15
Close modal

### 6.3 Liquid Flow Inside the Melt Pool.

The simulation results for the velocity distribution inside the melt pool obtained from the L-PBF and E-PBF models are also compared under the same energy density of 26 J/mm3 and heat source interaction time of 1.212 ms. The velocity of liquid inside the melt pool in L-PBF is higher than that in E-PBF due to the greater convection in L-PBF. Along the yz-plane corresponding to the origin (i.e., the longitudinal section), the velocity contours inside the melt pool when the laser and electron beam are at y = 8.0 mm are shown in Fig. 16. The maximum velocity for the L-PBF model is found in the middle of the melt pool, whereas the maximum melt-pool velocity in the E-PBF model is detected toward the tail end from the center.

Fig. 16
Fig. 16
Close modal

The temperature variation leads to a surface tension gradient, which causes the Marangoni flow from low surface tension area to high surface tension area of the melt pool as described by Yuan and Gu [36]. The cooler liquid near the edge of the melt pool having higher surface tension tends to pull the liquid away from the melt-pool center. However, the magnitudes of the maximum velocity in the melt pool for the L-PBF and the E-PBF models are about 18.6 mm/s and 15.4 mm/s, respectively, which confirm that the values of the Reynolds number Re(= ρuiΦ/μ) are very low and the flow is laminar in both the L-PBF and E-PBF cases.

### 6.4 Cooling Rate.

The cooling rates for the Ti-6Al-4V melt pool in L-PBF and E-PBF are found very fast due to the combined heat transfer. For the given specifications of the laser and electron beam, a point on the top surface at y = 5 mm along the scan is considered to observe the temperature variation with respect to time. Any temperature above the liquidus temperature (1938 K) indicates a complete liquid state of that point. If the temperature of that point is below the solidus temperature (1878 K), it specifies the point to be in the solid state [39,41]. The material is in the mushy zone when it remains between the liquidus and solidus temperatures. The heating and cooling of the point at $(x,y,z)$ = (0, 5 mm, 0) with respect to time for L-PBF and E-PBF are shown in Figs. 17(a) and 17(b), respectively. The variation of heating and cooling rates (in K/µs) with respect to time [22] for the L-PBF and E-PBF models are shown in Figs. 18(a) and 18(b), respectively, where the vertical axes values are multiplied by (−1). In both cases, the time count starts when the heat source strikes the point. The liquid melt pool takes 4.5 ms to cool down from the maximum temperature to the solidus temperature in L-PBF, yielding an average cooling rate of 1.5 K/µs. Under the same irradiation condition, the liquid melt pool takes about 8 ms for cooling in E-PBF, which gives an average cooling rate of 0.74 K/µs. These results are obtained for an energy density of 26 J/mm3 and the interaction time of 1.212 ms of the laser and electron beams.

Fig. 17
Fig. 17
Close modal
Fig. 18
Fig. 18
Close modal

## 7 Conclusions

Numerical modeling for the transient thermo-fluid properties of a Ti-6Al-4V powder-bed in the PBF process is conducted to determine the differences between the L-PBF and electron beam PBF (L-PBF) processes. The unique features of the two processes are characterized by implementing physical assumptions. A moving heat source with the Gaussian distribution is applied as a user-defined function and the thermo-fluid properties of the melt pool are obtained using the finite volume method. Modeling results for the temperature distributions, thermo-physical properties, melt-pool geometries, and solid–liquid–solid phase changes are compared at similar operating conditions. Based upon the comparison, the following conclusions can be drawn:

• A similar irradiance condition is considered for both the laser and the electron beam to make a valid comparison. The absorption efficiency in the Ti-6Al-4V powder-bed is determined by adjusting the powder porosity, layer thickness, optical absorbance, optical penetration depth, and efficiency of the heat source.

• For an irradiance level of 26 J/mm3 and interaction time of 1.212 ms, simulation results show that the liquid melt pool in E-PBF cools down in 8 ms whereas the melt pool in L-PBF cools down 4.5 ms. Thus, the cooling rate in L-PBF is faster than that in E-PBF. The lower volume of melt pool, absence of preheating, and convection heat transfer at the top surface make the cooling process faster in L-PBF.

• At the same scanning speed and beam spot size, the depth of the melt pool in E-PBF is higher than that in L-PBF. The lack of penetration of laser causes shallow melt-pool depth in L-PBF. Although the length and depth are significantly different, the width of the melt pool is quite similar in both L-PBF and E-PBF under similar operating conditions.

• The melt-pool volume is larger in E-PBF than that in L-PBF under the same energy density of 26 J/mm3 and interaction time of 1.212 ms. The melt-pool volume increases with the increase of beam power but decreases with the increase of scanning speed as evinced by the simulation results and experimental validation.

• The porosity of Ti-6Al-4V powder is taken as 50% for both L-PBF and E-PBF to keep the consistency of comparison. However, if the powder porosity is increased, the maximum temperature in the melt pool becomes higher due to the lower density of powder. This eventually makes the cooling rate slower than the current cases.

• Due to greater convection, the melt-pool velocity in the L-PBF process is higher than that in the E-PBF process. The fluid flow is laminar in both cases.

• A comparative study, differentiating the effects of the laser and electron beams under similar irradiation conditions, provides a thorough understanding of the physics involved in the two processes. The study facilitates the design for correct experiments prior to the actual production by giving room to optimize the process parameters and control the energy transfers in the L-PBF and E-PBF processes.

Numerical simulation of thermal behavior and the melt-pool dynamics as a result of the interaction between the moving heat source and powder zone for a single scan is the foundation for obtaining feedback of laser or electron beam processing parameters in the PBF process. The residual stress analysis of the processed part and shape optimization of the melt zone also depend on the thermal history and melt-pool evolution during the process. Therefore, the thermo-fluid models presented in this study are focused on characterizing the thermal history and melt-pool dynamics of the two significant PBF processes and offer a wide range of results with good accuracy. The model can incorporate various materials and operating conditions for further analyses of the material and process behavior in the powder-bed fusion additive manufacturing process.

## Funding Data

• This research is funded by the National Science Foundation (award number OIA-1541079) and the Louisiana Board of Regents.

## Conflict of Interest

There are no conflicts of interest.

## Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

## Nomenclature

• a =

absorption coefficient of laser (1/μm)

•
• h =

convective heat transfer coefficient (W/m2 K)

•
• k =

effective thermal conductivity (W/m K)

•
• s =

melt-front position with respect to time

•
• t =

time (s)

•
• z =

distance in the direction of penetration (mm)

•
• A =

optical absorbance

•
• B =

a computational constant to avoid division by zero

•
• C =

a constant for viscosity calculation

•
• H =

total enthalpy (J/kg)

•
• P =

pressure (Pa)

•
• R =

molar gas constant (kcal/mol K)

•
• S =

penetration depth (µm)

•
• T =

temperature (K)

•
• V =

acceleration voltage of the electron beam (kV)

•
• cp =

specific heat capacity (J/kg K)

•
• cp,l =

specific heat capacity of liquid material (J/kg K)

•
• fL =

liquid fraction

•
• gz =

gravitational acceleration (m/s2)

•
• kl =

thermal conductivity of liquid Ti (W/m K)

•
• lt =

powder layer thickness (mm)

•
• ti =

heat source interaction time (s) or (ms)

•
• vs =

beam scanning speed (mm/s)

•
• xi =

distance along the Cartesian coordinates (mm)

•
• xs =

instantaneous position of heat source in the x-direction (mm)

•
• ys =

instantaneous position of heat source in the y-direction (mm)

•
• Aw =

atomic weight (g/mol)

•
• CM =

a constant regarding mushy zone morphology

•
• ED =

energy density (J/mm3)

•
• Ib =

electron beam current (mA)

•
• Ixy =

the Gaussian surface intensity profile

•
• IZ =

beam penetration function

•
• Lf =

latent heat of fusion (kJ/K)

•
• PH =

power of the heat source (W)

•
• $Q˙(x,y,z)$ =

absorbed heat flux (W/m2)

•
• $Sgz$ =

source term due to gravity

•
• TL =

liquidus temperature (K)

•
• Tm =

melting temperature (K)

•
• Tpreheat =

preheat temperature (K)

•
• Tref =

reference temperature (K)

•
• TS =

solidus temperature (K)

•
• Tw =

wall temperature (K)

•
• u, v, w =

velocity components in the x-, y-, and z-directions, respectively (m/s)

•
• ui, uj =

velocity along the Cartesian coordinates (m/s)

•
• Ste =

Stefan number

•
• αl =

thermal diffusivity of liquid (m2/s)

•
• β =

coefficient of volume expansion

•
• ɛ =

emissivity

•
• η =

efficiency of laser or electron beam

•
• λ =

parameter in interfacial melt-front equation

•
• μ =

absolute viscosity (N·s/m2) or (Pa·s)

•
• μm =

viscosity of liquid metal (N·s/m2) or (Pa·s)

•
• ρ =

density of the material (kg/m3)

•
• ρl =

density of liquid Ti (kg/m3)

•
• σ =

Stefan–Boltzmann constant (W/m2 K4)

•
• Φ =

laser or electron beam spot size (mm)

•
• ΦV =

viscous dissipation term

## References

1.
Sames
,
W. J.
,
List
,
F. A.
,
Pannala
,
S.
,
Dehoff
,
R. R.
, and
Babu
,
S. S.
,
2016
, “
The Metallurgy and Processing Science of Metal Additive Manufacturing
,”
Int. Mater. Rev.
,
61
(
5
), pp.
315
360
. 10.1080/09506608.2015.1116649
2.
Goodridge
,
R.
, and
Ziegelmeier
,
S.
,
2017
, “7—Powder Bed Fusion of Polymers,”
,
,
Duxford, UK
, pp.
181
204
. 10.1016/C2014-0-03891-9
3.
Sing
,
S. L.
,
An
,
J.
,
Yeong
,
W. Y.
, and
Wiria
,
F. E.
,
2015
, “
Laser and Electron-Beam Powder-Bed Additive Manufacturing of Metallic Implants: A Review on Processes, Materials and Designs
,”
J. Orthop. Res.
,
34
(
3
), pp.
369
385
. 10.1002/jor.23075
4.
Zhao
,
X.
,
Li
,
S.
,
Zhang
,
M.
,
Liu
,
Y.
,
Sercombe
,
T. B.
, and
Wang
,
S.
,
2016
, “
Comparison of the Microstructures and Mechanical Properties of Ti–6Al–4 V Fabricated by Selective Laser Melting and Electron Beam Melting
,”
Mater. Des.
,
95
, pp.
21
31
. 10.1016/j.matdes.2015.12.135
5.
Gong
,
H.
,
Rafi
,
K.
,
Gu
,
H.
,
Janaki Ram
,
G. D.
,
Starr
,
T.
, and
Stucker
,
B.
,
2015
, “
Influence of Defects on Mechanical Properties of Ti–6Al–4 V Components Produced by Selective Laser Melting and Electron Beam Melting
,”
Mater. Des.
,
86
, pp.
545
554
. 10.1016/j.matdes.2015.07.147
6.
Chastand
,
V.
,
Quaegebeur
,
P.
,
Maia
,
W.
, and
Charkaluk
,
E.
,
2018
, “
Comparative Study of Fatigue Properties of Ti-6Al-4V Specimens Built by Electron Beam Melting (EBM) and Selective Laser Melting (SLM)
,”
Mater. Charact.
,
143
, pp.
76
81
. 10.1016/j.matchar.2018.03.028
7.
Wysocki
,
B.
,
Maj
,
P.
,
Sitek
,
R.
,
Buhagiar
,
J.
,
Kurzydłowski
,
K. J.
, and
Święszkowski
,
W.
,
2017
, “
Laser and Electron Beam Additive Manufacturing Methods of Fabricating Titanium Bone Implants
,”
Appl. Sci.
,
7
(
7
), pp.
1
20
. 10.3390/app7070657
8.
Gokuldoss
,
P. K.
,
Kolla
,
S.
, and
Eckert
,
J.
,
2017
, “
Additive Manufacturing Processes: Selective Laser Melting, Electron Beam Melting and Binder Jetting-Selection Guidelines
,”
Materials (Basel, Switzerland)
,
10
(
6
), pp.
1
12
. 10.3390/ma10060672
9.
Raplee
,
J.
,
Plotkowski
,
A.
,
Kirka
,
M. M.
,
Dinwiddie
,
R.
,
Okello
,
A.
,
Dehoff
,
R. R.
, and
Babu
,
S. S.
,
2017
, “
Thermographic Microstructure Monitoring in Electron Beam Additive Manufacturing
,”
Sci. Rep.
,
7
(
43554
), pp.
1
16
. 10.1038/srep43554
10.
Siddiqui
,
S. F.
,
Fasoro
,
A. A.
,
Cole
,
C.
, and
Gordon
,
A. P.
,
2019
, “
Mechanical Characterization and Modeling of Direct Metal Laser Sintered Stainless Steel GP1
,”
ASME J. Eng. Mater. Technol.
,
141
(
3
), p.
031009
. 10.1115/1.4042867
11.
Riedlbauer
,
D.
,
Scharowsky
,
T.
,
Singer
,
R. F.
,
Steinmann
,
P.
,
Körner
,
C.
, and
Mergheim
,
J.
,
2017
, “
Macroscopic Simulation and Experimental Measurement of Melt Pool Characteristics in Selective Electron Beam Melting of Ti-6Al-4V
,”
,
88
(
5
), pp.
1309
1317
. 10.1007/s00170-016-8819-6
12.
Cook
,
P. S.
, and
Murphy
,
A. B.
,
2020
, “
Simulation of Melt Pool Behaviour During Additive Manufacturing: Underlying Physics and Progress
,”
,
31
(
100909
), pp.
1
23
13.
Rahman
,
M. S.
,
Schilling
,
P. J.
,
Herrington
,
P. D.
, and
Chakravarty
,
U. K.
,
2018
, “
A Comparative Study Between Selective Laser Melting and Electron Beam Additive Manufacturing Based on Thermal Modeling
,”
Proceedings of the ASME International Mechanical Engineering Congress and Exposition
,
Pittsburgh, PA
,
Nov. 9−15, 2018
, IMECE2018-86428, Vol.
1
,
, pp.
1
13
.
14.
Bikas
,
H.
,
Stavropoulos
,
P.
, and
Chryssolouris
,
G.
,
2015
, “
Additive Manufacturing Methods and Modelling Approaches: A Critical Review
,”
,
83
(
1
), pp.
389
405
. 10.1007/s00170-015-7576-2
15.
Qi
,
H.
,
Mazumder
,
J.
, and
Ki
,
H.
,
2006
, “
Numerical Simulation of Heat Transfer and Fluid Flow in Coaxial Laser Cladding Process for Direct Metal Deposition
,”
J. Appl. Phys.
,
100
(
2
), p.
024903
. 10.1063/1.2209807
16.
Moraitis
,
G. A.
, and
Labeas
,
G. N.
,
2008
, “
Residual Stress and Distortion Calculation of Laser Beam Welding for Aluminum Lap Joints
,”
J. Mater. Process. Technol.
,
198
(
1–3
), pp.
260
269
. 10.1016/j.jmatprotec.2007.07.013
17.
Wang
,
L.
,
Felicelli
,
S.
,
Gooroochurn
,
Y.
,
Wang
,
P. T.
, and
Horstemeyer
,
M. F.
,
2008
, “
Optimization of the LENS® Process for Steady Molten Pool Size
,”
Mater. Sci. Eng., A
,
474
(
1–2
), pp.
148
156
. 10.1016/j.msea.2007.04.119
18.
Roberts
,
I. A.
,
Wang
,
C. J.
,
Esterlein
,
R.
,
Stanford
,
M.
, and
Mynors
,
D. J.
,
2009
, “
A Three-Dimensional Finite Element Analysis of the Temperature Field During Laser Melting of Metal Powders in Additive Layer Manufacturing
,”
Int. J. Mach. Tools Manuf.
,
49
(
12–13
), pp.
916
923
. 10.1016/j.ijmachtools.2009.07.004
19.
Liu
,
C.
,
Wu
,
B.
, and
Zhang
,
J.
,
2010
, “
Numerical Investigation of Residual Stress in Thick Titanium Alloy Plate Joined With Electron Beam Welding
,”
Metall. Mater. Trans. B
,
41
(
5
), pp.
1129
1138
. 10.1007/s11663-010-9408-y
20.
Yang
,
J.
,
Sun
,
S.
,
Brandt
,
M.
, and
Yan
,
W.
,
2010
, “
Experimental Investigation and 3D Finite Element Prediction of the Heat Affected Zone During Laser Assisted Machining of Ti–6Al–4 V Alloy
,”
J. Mater. Process. Technol.
,
210
(
15
), pp.
2215
2222
. 10.1016/j.jmatprotec.2010.08.007
21.
Lacki
,
P.
, and
,
K.
,
2011
, “
Numerical Simulation of the Electron Beam Welding Process
,”
Comput. Struct.
,
89
(
11–12
), pp.
977
985
. 10.1016/j.compstruc.2011.01.016
22.
Shen
,
N.
, and
Chou
,
K.
,
2012
, “
Thermal Modeling of Electron Beam Additive Manufacturing Process: Powder Sintering Effects
,”
Proceedings of the ASME Manufacturing Science and Engineering Conference
,
Notre Dame, IN
,
June 4–8, 2012
, MSEC2012-7253, pp.
287
295
.
23.
Cheng
,
B.
,
Price
,
S.
,
Lydon
,
J.
,
Cooper
,
K.
, and
Chou
,
K.
,
2014
, “
On Process Temperature in Powder-Bed Electron Beam Additive Manufacturing: Model Development and Validation
,”
ASME J. Manuf. Sci. Eng.
,
136
(
6
), p.
061018
. 10.1115/1.4028484
24.
Chen
,
Y. X.
,
Wang
,
X. J.
, and
Chen
,
S. B.
,
2014
, “
The Effect of Electron Beam Energy Density on Temperature Field for Electron Beam Melting
,”
,
900
, pp.
631
638
. 10.4028/www.scientific.net/AMR.900.631
25.
Zäh
,
M. F.
, and
Lutzmann
,
S.
,
2010
, “
Modelling and Simulation of Electron Beam Melting
,”
Product. Eng. Res. Develop.
,
4
(
1
), pp.
15
23
. 10.1007/s11740-009-0197-6
26.
Andreotta
,
R.
,
,
L.
, and
Brindley
,
W.
,
2017
, “
Finite Element Simulation of Laser Additive Melting and Solidification of Inconel 718 With Experimentally Tested Thermal Properties
,”
Finite Elem. Anal. Des.
,
135
, pp.
36
43
. 10.1016/j.finel.2017.07.002
27.
,
M.
,
,
L.
,
Brindley
,
W.
, and
Romano
,
J.
,
2017
, “
Optimizing Quality of Additively Manufactured Inconel 718 Using Powder Bed Laser Melting Process
,”
,
11
, pp.
60
70
28.
,
L.
,
Romano
,
J.
,
Brindley
,
W.
, and
Burlatsky
,
S.
,
2017
, “
Effective Liquid Conductivity for Improved Simulation of Thermal Transport in Laser Beam Melting Powder Bed Technology
,”
,
14
, pp.
13
23
29.
Wen
,
S.
, and
Shin
,
Y. C.
,
2010
, “
Modeling of Transport Phenomena During the Coaxial Laser Direct Deposition Process
,”
J. Appl. Phys.
,
108
(
4
), p.
044908
. 10.1063/1.3474655
30.
Wang
,
R.
,
Lei
,
Y.
, and
Shi
,
Y.
,
2011
, “
Numerical Simulation of Transient Temperature Field During Laser Keyhole Welding of 304 Stainless Steel Sheet
,”
Opt. Laser Technol.
,
43
(
4
), pp.
870
873
. 10.1016/j.optlastec.2010.10.007
31.
Cho
,
W.
,
Na
,
S.
,
Thomy
,
C.
, and
Vollertsen
,
F.
,
2012
, “
Numerical Simulation of Molten Pool Dynamics in High Power Disk Laser Welding
,”
J. Mater. Process. Technol.
,
212
(
1
), pp.
262
275
. 10.1016/j.jmatprotec.2011.09.011
32.
Rai
,
R.
,
Burgardt
,
P.
,
Milewski
,
J.
,
Lienert
,
T.
, and
DebRoy
,
T.
,
2009
, “
Heat Transfer and Fluid Flow During Electron Beam Welding of 21Cr–6Ni–9Mn Steel and Ti–6Al–4 V Alloy
,”
J. Phys. D: Appl. Phys.
,
42
(
2
), p.
02550
. 10.1088/0022-3727/42/2/025503
33.
Rai
,
R.
,
Palmer
,
T. A.
,
Elmer
,
J. W.
, and
Debroy
,
T.
,
2009
, “
Heat Transfer and Fluid Flow During Electron Beam Welding of 304L Stainless Steel Alloy
,”
Weld. J.
,
88
(
3
), pp.
54
61
. 10.1088/0022-3727/42/2/025503
34.
Li
,
J. F.
,
Li
,
L.
, and
Stott
,
F. H.
,
2004
, “
A Three-Dimensional Numerical Model for a Convection-Diffusion Phase Change Process During Laser Melting of Ceramic Materials
,”
Int. J. Heat Mass Trans.
,
47
(
25
), pp.
5523
5539
. 10.1016/j.ijheatmasstransfer.2004.07.024
35.
Chahine
,
G.
,
2011
, “
Application of Digital Engineering in the Development of a Bio-Adaptable Dental Implant
,”
Ph.D. thesis
,
Southern Methodist University
,
Dallas, TX
.
36.
Yuan
,
P.
, and
Gu
,
D.
,
2015
, “
Molten Pool Behaviour and Its Physical Mechanism During Laser PBF of TiC/AlSi10Mg Nanocomposites: Simulation and Experiments
,”
J. Phys. D: Appl. Phys.
,
48
(
3
), p.
035303
. 10.1088/0022-3727/48/3/035303
37.
Rahman
,
M. S.
,
Schilling
,
P. J.
,
Herrington
,
P. D.
, and
Chakravarty
,
U. K.
,
2016
, “
Thermo-Fluid Characterizations of Ti-6Al-4V Melt Pool in Powder-Bed Electron Beam Additive Manufacturing
,”
Proceedings of the ASME 2016 International Mechanical Engineering Congress and Exposition
,
Phoenix, AZ
,
Nov. 11–17, 2016
, IMECE2016-65854, Vol.
1
,
, pp.
1
9
.
38.
Rahman
,
M. S.
,
Schilling
,
P. J.
,
Herrington
,
P. D.
, and
Chakravarty
,
U. K.
,
2017
, “
Thermal Analysis of Electron Beam PBF Using Ti-6Al-4V Powder-Bed
,”
Proceedings of the ASME 2017 International Mechanical Engineering Congress and Exposition
,
Tampa, FL
,
Nov. 3–9, 2017
, IMECE2017-71663, Vol.
1
,
, pp.
1
13
.
39.
Rahman
,
M. S.
,
Schilling
,
P. J.
,
Herrington
,
P. D.
, and
Chakravarty
,
U. K.
,
2019
, “
Thermofluid Properties of Ti-6Al-4V Melt Pool in Powder-Bed Electron Beam Additive Manufacturing
,”
ASME J. Eng. Mater. Technol.
,
141
(
4
), p.
041006
. 10.1115/1.4043342
40.
Rahman
,
M. S.
,
Schilling
,
P. J.
,
Herrington
,
P. D.
, and
Chakravarty
,
U. K.
,
2020
, “
Heat Transfer and Melt-Pool Evolution During Powder-Bed Fusion of Ti-6Al-4V Parts Under Various Laser Irradiation Conditions
,”
Proceedings of the ASME International Mechanical Engineering Congress and Exposition
,
Portland, OR
,
Nov. 16–19, 2020
, IMECE2020-23838, Vol.
1
,
, pp.
1
10
.
41.
Rahman
,
M. S.
,
Zeng
,
C.
,
Wen
,
H.
,
Guo
,
S.
,
Schilling
,
P. J.
,
Herrington
,
P. D.
, and
Chakravarty
,
U. K.
,
2020
, “
An Investigation on the Melting and Solidification Behavior of Ti64 Powder in the Laser Powder-Bed Fusion Process
,”
Proceedings of the Louisiana EPSCoR RII CIMM 2020 Symposium
,
Baton Rouge, LA
,
July 20, 2020
, pp.
9
12
.
42.
Rahman
,
M. S.
,
Schilling
,
P. J.
,
Herrington
,
P. D.
, and
Chakravarty
,
U. K.
,
2019
, “
Thermal Behavior and Melt-Pool Dynamics of Cu-Cr-Zr Alloy in Powder-Bed Selective Laser Melting Process
,”
Proceedings of the ASME International Mechanical Engineering Congress and Exposition
, Vol.
1
,
,
Salt Lake City, UT
,
Nov. 11–14, 2019
, IMECE2019-11581, pp.
1
9
.
43.
Jamshidinia
,
M.
,
Kong
,
F.
, and
Kovacevic
,
R.
,
2012
, “
Temperature Distribution and Fluid Flow Modeling of Electron Beam Melting® (EBM)
,”
Proceedings of the ASME 2012 International Mechanical Engineering Congress and Exposition
,
Houston, TX
,
Nov. 9–15, 2012
, IMECE2012-89440, Vol.
7
,
Part D, Fluid and Heat Transfer
, pp.
3089
3102
.
44.
Jamshidinia
,
M.
,
Kong
,
F.
, and
Kovacevic
,
R.
,
2013
, “
Numerical Modeling of Heat Distribution in the Electron Beam Melting of Ti–6Al–4 V
,”
ASME J. Manuf. Sci. Eng.
,
135
(
6
), p.
061010
. 10.1115/1.4025746
45.
Sih
,
S. S.
, and
Barlow
,
J. W.
,
2004
, “
The Prediction of the Emissivity and Thermal Conductivity of Powder Beds
,”
Part. Sci. Technol.
,
22
(
3
), pp.
291
304
. 10.1080/02726350490501682a
46.
Arce
,
A. N.
,
2012
, “
Thermal Modeling and Simulation of Electron Beam Melting for Rapid Prototyping on Ti6Al4 V Alloys
,”
Ph.D. thesis
,
North Carolina State University
,
Raleigh, NC
.
47.
Pei
,
Y. T.
,
Ocelik
,
V.
, and
Hosson
,
J. T. M. D.
,
2002
, “
SiCp/Ti6Al4V Functionally Graded Materials Produced by Laser Melt Injection
,”
Acta Mater.
,
50
(
8
), pp.
2035
51
. 10.1016/S1359-6454(02)00049-6
48.
Esen
,
A.
, and
Kutluay
,
S.
,
2002
, “
A Numerical Solution of the Stefan Problem With a Neumann-Type Boundary Condition by Enthalpy Method
,”
Appl. Math. Comput.
,
148
(
2004
), pp.
321
329
. 10.1016/S0096-3003(02)00846-9
49.
Yan
,
W.
,
Ge
,
W.
,
Smith
,
J.
,
Lin
,
S.
,
Kafka
,
O. L.
,
Lin
,
F.
, and
Liu
,
W. K.
,
2016
, “
Multi-scale Modeling of Electron Beam Melting of Functionally Graded Materials
,”
Acta Mater.
,
115
, pp.
403
412
. 10.1016/j.actamat.2016.06.022
50.
Carriere
,
P. R.
, and
Yue
,
S.
,
2017
, “
Energy Absorption During Pulsed Electron Beam Spot Melting of 304 Stainless Steel: Monte-Carlo Simulations and In-Situ Temperature Measurements
,”
Vacuum
,
142
, pp.
114
122
. 10.1016/j.vacuum.2017.04.039
51.
Yan
,
W.
,
Smith
,
J.
,
Ge
,
W.
,
Lin
,
F.
, and
Liu
,
W. K.
,
2015
, “
Multiscale Modeling of Electron Beam and Substrate Interaction: A New Heat Source Model
,”
Comput. Mech.
,
56
(
2
), pp.
265
276
. 10.1007/s00466-015-1170-1
52.
Brown
,
M. S.
, and
Arnold
,
C. B.
,
2010
,
Fundamentals of Laser-Material Interaction and Application to Multiscale Surface Modification, Laser Precision Microfabrication
, Vol.
135
,
Springer Series in Materials Science Springer
,
Berlin, Germany
, pp.
91
120
.
53.
Ilican
,
S.
,
Caglar
,
M.
, and
Caglar
,
Y.
,
2007
, “
Determination of the Thickness and Optical Constants of Transparent Indium-Doped ZnO Thin Films by the Envelope Method
,”
Materials Science-Poland
,
25
(
3
), pp.
709
718
.
54.
Lisiecki
,
A.
,
2019
, “
Study of Optical Properties of Surface Layers Produced by Laser Surface Melting and Laser Surface Nitriding of Titanium Alloy
,”
Materials
,
12
(
3112
), pp.
1
14
. 10.3390/ma12193112
55.
Ogoh
,
W.
, and
Groulx
,
D.
,
2010
, “
Stefan’s Problem: Validation of a One-Dimensional Solid–Liquid Phase Change Heat Transfer Process
,”
Proceedings of the COMSOL Conference
,
Boston, MA
,
Oct. 7–9, 2010
, pp.
1
6
.
56.
Dilip
,
J. J. S.
,
Zhang
,
S.
,
Teng
,
C.
,
Zeng
,
K.
,
Robinson
,
C.
,
Pal
,
D.
, and
Stucker
,
B.
,
2017
, “
Influence of Processing Parameters on the Evolution of Melt Pool, Porosity, and Microstructures in Ti-6Al-4V Alloy Parts Fabricated by Selective Laser Melting
,”
,
2
(
3
), pp.
157
167
. 10.1007/s40964-017-0030-2
57.
Matilainen
,
V.
,
Piili
,
H.
,
Salminen
,
A.
, and
Nyrhilä
,
O.
,
2015
, “
Preliminary Investigation of Keyhole Phenomena During Single Layer Fabrication in Laser Additive Manufacturing of Stainless Steel
,”
Phys. Procedia
,
78
, pp.
377
387
. 10.1016/j.phpro.2015.11.052