Analysis of the Thermal Transfers in a VASM Crucible: Electron Beam Melting Experiment and Numerical Simulation

: A description of the Vacuum Arc Skull Melting (VASM) process is presented showing its particularly complex features because of the mixing of porous raw materials with the dense remelted metal as well as the very high temperature and the highly transient nature of the process. This paper presents a 3D transient mathematical modelling of the heat transport with the aim of bringing a better understanding of the thermal behavior of the material into the crucible during a melting cycle. The model takes into account the heat input provided by the incoming metal thanks to an adaptive meshing, as well as the latent heat of solidiﬁcation and the radiative heat transfers. An experimental validation of the model is presented where an electron beam heating source mimics the heat e ﬀ ect of the arc thanks to an excellent guidance of the beam over the melt surface. A comparison between the measured and calculated temperatures of a steel load is reported and reveals a satisfactory agreement. With very few adjustments, concerning mainly heat radiation at the top surface of metal into the crucible, the numerical model appears to be an e ﬃ cient numerical tool to simulate the VASM process at the industrial scale.


Introduction
One of the main objectives of the remelting processes of titanium alloys is to eliminate the two main types of exogenous inclusions: the low density inclusions (LDI), which include the hard-alpha (titanium enriched with interstitial elements such as nitrogen or/and oxygen), and the high density inclusions (HDI). These defects can cause premature fatigue failure of any mechanical part of plane, like critical rotating parts of aircraft engines [1,2], compromising its safety and structural integrity. Basically, melting technologies used in the titanium industry are vacuum arc remelting (VAR), plasma arc cold hearth melting (PACHM) [3], and electron beam cold hearth melting (EBCHM) [4,5]. Three successive VAR melts are requested to remove LDI and produce certified premium grade titanium alloys when only one is necessary after a prior cold hearth melting operation. In the latter case, LDI are removed more easily because of the long metal residence time within the horizontal refining hearth, enhancing the dissolution of solid inclusions. Furthermore, HDI are also trapped due to the density separation by settling in the mushy zone. Those technologies have been approved by customers for very demanding applications. However, the feedstock, which is composed of scrap and sponge, must be carefully sorted and cleaned up. In order to allow a large variety of scraps and reduce the cost, the vacuum arc skull melting (VASM) process was developed. It makes possible a massive production of titanium Metals 2020, 10 alloys based on feedstock products with an outstanding flexibility regarding shapes and dimensions of metal pieces to be recycled.

Description of the VASM Process
The furnace of the VASM process is decomposed in two parts: (i) a crucible where the melting of a consumable electrode occurs and (ii) an ingot mold in which the liquid metal is poured. The first step has the particularity to create the skull which will be used as the electrode of the next melting cycle. This is accomplished by means of a parallelepiped crucible made from water-cooled copper panels bolted together [6]. At the beginning of each cycle, scrap of any type and size (called the initial load on Figure 1) is deposited in this crucible. Above the crucible, a consumable electrode is vertically hung from a ram which can be moved up and down. As stated above, this electrode is the skull from the previous cycle. The melting sequence begins when an electric arc is ignited between the lower tip of the electrode and the initial load in the crucible, causing melting of both. Liquid metal (originating mainly from the electrode) accumulates inside the crucible where it fills in a pool at the top before partly solidifying upon contact with the water-cooled panels. The metal within the crucible, constituted of the initial load, the solidified metal and the liquid pool is called the charge in this paper. When the mass of the liquid pool reaches an expected value, typically equal to the mass of the initial load, the electric arc is switched off. The crucible is then tilted and the liquid pool is poured into the ingot mold. After cooling under vacuum or an inert gas atmosphere, the furnace is unloaded: a new ingot is stripped from the mold, and a thick crust of solid metal, hereinafter called the skull, is stripped from the crucible.
Metals 2020, 10, x FOR PEER REVIEW 2 of 21 massive production of titanium alloys based on feedstock products with an outstanding flexibility regarding shapes and dimensions of metal pieces to be recycled.

Description of the VASM Process
The furnace of the VASM process is decomposed in two parts: (i) a crucible where the melting of a consumable electrode occurs and (ii) an ingot mold in which the liquid metal is poured. The first step has the particularity to create the skull which will be used as the electrode of the next melting cycle. This is accomplished by means of a parallelepiped crucible made from water-cooled copper panels bolted together [6]. At the beginning of each cycle, scrap of any type and size (called the initial load on Figure 1) is deposited in this crucible. Above the crucible, a consumable electrode is vertically hung from a ram which can be moved up and down. As stated above, this electrode is the skull from the previous cycle. The melting sequence begins when an electric arc is ignited between the lower tip of the electrode and the initial load in the crucible, causing melting of both. Liquid metal (originating mainly from the electrode) accumulates inside the crucible where it fills in a pool at the top before partly solidifying upon contact with the water-cooled panels. The metal within the crucible, constituted of the initial load, the solidified metal and the liquid pool is called the charge in this paper. When the mass of the liquid pool reaches an expected value, typically equal to the mass of the initial load, the electric arc is switched off. The crucible is then tilted and the liquid pool is poured into the ingot mold. After cooling under vacuum or an inert gas atmosphere, the furnace is unloaded: a new ingot is stripped from the mold, and a thick crust of solid metal, hereinafter called the skull, is stripped from the crucible. The electric power is a main parameter of the process [7]. The arc must be powerful enough to maintain the molten metal in liquid state, counteracting against the radiative losses from the top surface as well as the losses to the cooled walls of the crucible. A carefully adapted electric power leads to a very limited over-heating of the melt. The heat balance evolves during the melting process, leading to an unsteady regime before the liquid pool volume remains constant.
Local intense heating of the crucible walls can lead to the formation of a crack, which might allow water from the cooling channel to penetrate the crucible and come into contact with the molten metal, hence an explosive reaction in the crucible. Ref [6] discussed the possible number of melting cycles before the next repair of the crucible, taking into account geometrical and structural criteria. The electric power is a main parameter of the process [7]. The arc must be powerful enough to maintain the molten metal in liquid state, counteracting against the radiative losses from the top surface as well as the losses to the cooled walls of the crucible. A carefully adapted electric power leads to a very limited over-heating of the melt. The heat balance evolves during the melting process, leading to an unsteady regime before the liquid pool volume remains constant.
Local intense heating of the crucible walls can lead to the formation of a crack, which might allow water from the cooling channel to penetrate the crucible and come into contact with the molten metal, hence an explosive reaction in the crucible. Ref [6] discussed the possible number of melting cycles before the next repair of the crucible, taking into account geometrical and structural criteria.
The literature on VASM process is particularly poor because of the relative novelty of the technology and of the confidentiality applied in a sensitive sector. Most of the literature states the main advantages and disadvantages of the VASM process among them the flexibility of the sources of raw material [8], the high degree of chemical homogeneity of the cast ingots [9], but with typical casting defects such as porous surface and shrinkage cavities. However, the main drawbacks can be overrun by remelting the produced ingot in a VAR furnace [9].

Objectives of the Present Study
In that context, it was decided to develop a transient 3D numerical model of the thermal transfers in the metallic load contained in the crucible, accounting for the continuous input of molten metal in order to bring a better understanding of the material behaviour during a melting cycle.
To compare the model results with experimental observation at a pilot scale, rather than using a VASM small-scale furnace, the use of an electron beam melting (EBM) device was preferred, as it allows a precisely controlled distribution of the energy input at the melt surface. The boundary conditions (see below) must obviously be derived to adapt to the specific apparatus used in the experiments.
It should also be noted that, in order to have a better control of the material homogeneity and make the experiments and temperature recordings easier, a packed bed of XC10 grade steel balls was used in this study in place of an actual load composed of Ti alloy scrap. Also, the feeding material is composed of the same steel. Therefore, the expected temperature of the liquid metal is of the order of 1550-1600 • C, slightly higher than the liquidus.
The outline of the present paper directly results from the above considerations. In Section 2, the experimental procedure is described in detail, with an emphasis on the insertion of a set of thermocouples within the initial load. The next section is devoted to the mathematical model itself. The boundary conditions correspond to the EBM device and will require an adaptation to VASM operation in a further step. Finally, the last section presents a comparison of the numerical results with the observations, as well as a simulation of the influence of an operating parameter, i.e., the total EB power.

Experiments
As stated in Section 1.2, the study involves achieving experiments using a suitable furnace in order to reproduce the behaviour of the molten metal during vacuum arc skull melting. A laboratory electron beam melting (EBM) furnace has been used for that purpose. After a first sub-section which details the experimental device, we will emphasize more particularly the operating conditions and the patterns which were effectively scanned by the electron beam.

Experimental Set-Up
The experiments have been conducted with a 100-kW laboratory EBM furnace (model: ALD-Lab100. ALD vacuum technologies, Hanau, Germany) as shown schematically in Figure 2. It includes a vacuum chamber (with a pressure lower than 10 −3 mbar) and an EB gun with an accelerating voltage of 40 kV. A water-cooled copper crucible has been especially designed so as to mimic at the laboratory scale an industrial VASM crucible. The inside dimensions of the crucible are: X c × Y c × Z c = 125 × 220 × 150 mm 3 . Magnetic lenses enable a very accurate control of the beam motion, thereby of the heat power density according established pattern shapes and sizes. The automatic EB processes are operated by EB Scan and Control System (ESCOSYS ® -ALD vacuum technologies, Hanau, Germany). All the variables involved in EBM process, such as beam pattern, beam power distribution and scanning frequency are governed by this operating system.
A drip melting experiment was performed to simulate the VASM process. Hence, steel balls (d p = 13 mm) were initially loaded inside the crucible as a packed bed, leading to a porosity close to ε = 0.36. The height of the ball layer is h IL = 35 mm. The steel ball layer mimics the scraps and raw materials in the VASM process and will be referenced as the initial load in this paper. A steel bar (square section with l bar = 63 mm) was positioned above and perpendicular to the crucible (see Figure 3a). During the melting, the steel bar was fed into the EB furnace at a constant speed by a bar pusher in order to keep the bar tip location constant, as shown in Figure 3b. Ten type K thermocouples (tc 1 to tc 10 ) were inserted within the steel ball layer at different heights (see Figure 3b). The locations of the thermocouples are also reported in the Table 1. During the melting process, the EB continuously melts the tip of the bar and steel droplets fall into the crucible where they may either immediately solidify or fill in a pool depending of the beam power. The steel alloy used for both the balls and the bar is a stainless steel whose properties will be reported in Section 3.3.
Metals 2020, 10, x FOR PEER REVIEW 4 of 21 (square section with = 63 mm) was positioned above and perpendicular to the crucible (see Figure 3a). During the melting, the steel bar was fed into the EB furnace at a constant speed by a bar pusher in order to keep the bar tip location constant, as shown in Figure 3b. Ten type K thermocouples ( c 1 to c 10 ) were inserted within the steel ball layer at different heights (see Figure 3b). The locations of the thermocouples are also reported in the Table 1. During the melting process, the EB continuously melts the tip of the bar and steel droplets fall into the crucible where they may either immediately solidify or fill in a pool depending of the beam power. The steel alloy used for both the balls and the bar is a stainless steel whose properties will be reported in Section 3.3. (square section with = 63 mm) was positioned above and perpendicular to the crucible (see Figure 3a). During the melting, the steel bar was fed into the EB furnace at a constant speed by a bar pusher in order to keep the bar tip location constant, as shown in Figure 3b. Ten type K thermocouples ( c 1 to c 10 ) were inserted within the steel ball layer at different heights (see Figure 3b). The locations of the thermocouples are also reported in the Table 1. During the melting process, the EB continuously melts the tip of the bar and steel droplets fall into the crucible where they may either immediately solidify or fill in a pool depending of the beam power. The steel alloy used for both the balls and the bar is a stainless steel whose properties will be reported in Section 3.3.

EB Operating Conditions
The intensity of the EB was linearly increased up to 0.5 A during the first 20 min and then was kept constant until the end of the melt, which means an electrical power equal to 20 kW as the voltage is kept constant U = 40 kV. The melting rate was measured through the pushing speed of the bar feeder (13 mm/min once the steady state was reached) as V f us = 0.41 kg·min −1 . No experimental information is available concerning the moment t f us when the actual melting started, so we estimated it at half of the transitory phase (∆t init = 20 min) and assumed a linear variation of the melt rate from this moment until the steady state. All the operating conditions of the melting test are reported in Figure 4.

EB Operating Conditions
The intensity of the EB was linearly increased up to 0.5 A during the first 20 min and then was kept constant until the end of the melt, which means an electrical power equal to 20 kW as the voltage is kept constant = 40 . The melting rate was measured through the pushing speed of the bar feeder (13 mm/min once the steady state was reached) as = 0.41 kg. min −1 . No experimental information is available concerning the moment when the actual melting started, so we estimated it at half of the transitory phase (Δ = 20 min ) and assumed a linear variation of the melt rate from this moment until the steady state. All the operating conditions of the melting test are reported in Figure 4.

Electron Beam Pattern
During a full time period , the beam scanned three different patterns, as depicted in Figure  5: the "bar" pattern was applied for melting the bar tip, the "candle" pattern was devoted to ensure a total melting of the solid "candle" which may appear between the bar tip and the charge into the crucible, and the "pool" pattern was used to heat the surface of the metal in the crucible. The kinetic

Electron Beam Pattern
During a full time period δt total , the beam scanned three different patterns, as depicted in Figure 5: the "bar" pattern was applied for melting the bar tip, the "candle" pattern was devoted to ensure a total melting of the solid "candle" which may appear between the bar tip and the charge into the crucible, and the "pool" pattern was used to heat the surface of the metal in the crucible. The kinetic energy of the beam is supposed to be converted to heat at the beam impingement with a yield η equal to 0.7 [10].
Metals 2020, 10, x FOR PEER REVIEW 6 of 21 energy of the beam is supposed to be converted to heat at the beam impingement with a yield equal to 0.7 [10]. For each pattern , the exposure time defines the corresponding thermal power: Table 2 shows the values of the exposure times as well as the thermal powers and heat flux densities since we have:

Mathematical Modelling
A 3D mathematical model has been developed to simulate the transient heat transfer into the metal contained in the crucible, taken particularly into account as the heat input provided by the For each pattern k, the exposure time δt k defines the corresponding thermal power: (1) Table 2 shows the values of the exposure times δt k as well as the thermal powers and heat flux densities since we have:

Mathematical Modelling
A 3D mathematical model has been developed to simulate the transient heat transfer into the metal contained in the crucible, taken particularly into account as the heat input provided by the incoming metal thanks to an adaptive meshing, as well as the latent heat of solidification. The fluid Metals 2020, 10, 1152 7 of 20 flow in the liquid pool was not explicitly calculated but the convective heat transfer was modelled through the use of an effective conductivity. This is undoubtedly the main assumption of this model.

Heat Transport Equation
The general equation of heat transport is given by: The steel alloy in the crucible experiences solidification between the liquidus and the solidus temperatures. In the mushy zone, the liquid fraction g L varies linearly from 0 to 1, which is a reasonable assumption because T solidus and T liquidus are quite close.
If < T solidus : If T solidus < T < T liquidus : If > T liquidus : The initial load is a porous media constituted of a packed bed of steel balls (d p = 13 mm). Radiative fluxes between balls are taken into account by using a radiative thermal conductivity λ ray [11].
where ε is the porosity of the packed bed, supposed to remain constant (ε = 0.36).
In order to take into account the convective and turbulent heat transfer within the liquid phase, we artificially increased the thermal conductivity in the pool by introducing a multiplying factor k λ . Such an important assumption has been used by several authors to roughly characterize the influence of turbulent convection in the liquid pool of ElectroSlag Refined (ESR) or Vacuum Arc Remelted (VAR) ingots [12][13][14]. Values for k λ typically ranged between 2 and 6. Here, a value k λ = 5 has been included in the model. With the aim of modelling all metal phases (liquid, mushy, solidified and initial load), an effective conductivity λ* is introduced and calculated as: The parameters ε, g L , and λ * take different expressions as gathered in Table 3. Table 3. Expressions of ε, g L , and λ * depending on which phase is the metal.

Metal Phase
As it will be discussed, the numerical resolution of Equation (3) is performed using a finite volume technique, hence a structured mesh composed of hexahedral cells. The melting of the steel bar leads to a rise in the height of the charge with the processing time. This growth of the metal volume is accounted for by an adaptive mesh in the Oz direction. According to the melting rate V f us , the height of the top layer of the mesh is increased at each time step and a new layer of cells is created when the top layer size becomes higher than a reference cell size. The steel droplets feeding the crucible are supposed to be at a temperature equal to: where ∆T sh is a superheat. The growth of the metal volume is taken into account by the source term S Z in the heat transfer equation as we will see in next section.

Top Surface
The boundary condition applied at the top surface of the domain is particularly complex since it is composed of three contributions as schematically illustrated Figure 6. These contributions are:

Metal Phase
As it will be discussed, the numerical resolution of Equation (3) is performed using a finite volume technique, hence a structured mesh composed of hexahedral cells. The melting of the steel bar leads to a rise in the height of the charge with the processing time. This growth of the metal volume is accounted for by an adaptive mesh in the Oz direction. According to the melting rate , the height of the top layer of the mesh is increased at each time step and a new layer of cells is created when the top layer size becomes higher than a reference cell size. The steel droplets feeding the crucible are supposed to be at a temperature equal to: where Δ ℎ is a superheat. The growth of the metal volume is taken into account by the source term in the heat transfer equation as we will see in next section.

Top Surface
The boundary condition applied at the top surface of the domain is particularly complex since it is composed of three contributions as schematically illustrated  The beam scans only the central part of the crucible S pool in order to avoid any beam-to-wall impingement. This pool surface is itself divided into four subzones as shown in Figure 7. Because of the shadowing influence of the steel bar, the subzone 3 is not heated by the beam. For the three remaining subzones, the heat density supplied by the EB is calculated as: where δt pool,k is the exposure time of the EB on S pool k . All the operating data are reported in the Table 4.
impingement. This pool surface is itself divided into four subzones as shown in Figure 7. Because of the shadowing influence of the steel bar, the subzone 3 is not heated by the beam. For the three remaining subzones, the heat density supplied by the EB is calculated as: where , is the exposure time of the EB on . All the operating data are reported in the Table 4.  7. Subzones of the "pool" surface.

Enthalpy Flux Density Provided by the Droplets
In this section, we will make use of the following notations: Figure 7. Subzones of the "pool" surface. In this section, we will make use of the following notations: • h 0 is the specific enthalpy of the top layer at any time t 0 : • h 1 is the specific enthalpy of the top layer at any time t 1 = t 0 + ∆t: • h drop is the specific enthalpy of the liquid drops falling during ∆t: In the preceding, h re f is the reference enthalpy at a standard state. The enthalpy flux density provided by the droplets is given by: As stated previously, this flux is incorporated in the model through an alike source term for all the cells in the top layer. At any time t 0 , the volume of the top layer is equal to V 0 . At t 1 = t 0 + ∆t, this volume becomes: Between t 0 and t 1 , if heat transfer within the crucible is neglected, the conservation of energy implies: By combining Equations (15) and (16), we get: ρ(h 1 − h 0 )V 1 represents the rate of heat accumulation in the volume V 1 during ∆t. The source term S Z in Equation (3) must therefore be: As the melting rate is V f us = ρ.V drop ∆t , we finally get: Heat Loss by Radiation The heat transfer by radiation between the steel bar, the top surface of the charge and the inner crucible wall plays an important role on the global heat balance. This is the reason why a precise calculation based on the determination of surface radiosities has been set up to get an accurate prediction of the net flux.
The space between the top of the charge and the tip of the steel bar is divided into 12 volumes having a rectangular parallelepiped shape (in purple in Figure 8), leading to a system of 52 surfaces. All solid surfaces are supposed to behave as grey bodies (i.e., their emissivity is independent on the temperature). All the shape factors between surfaces j and k can be derived from geometrical considerations. The radiosity emitted by each surface k is readily obtained by the inversion of the radiosity matrix [15]: where is the emissivity of surface k. The net radiative flux density lost by the surface k is finally: where σ is the Stefan-Boltzmann constant.

Crucibles Walls
The walls of the hearth are water-cooled copper panels. If a determination of the heat transfer coefficient for a casting mold is readily available [16,17], there is no accurate analysis in the literature dedicated to water-cooled refining hearths. A Fourier condition is applied to model the heat resistance between the solid steel, copper and the cooling water as: Figure 8. Meshing of the domain above the crucible top including the bar tip for radiative calculation purpose, the meshing of the charge is addressed in Section 3.3.
All the shape factors F jk between surfaces j and k can be derived from geometrical considerations. The radiosity J k emitted by each surface k is readily obtained by the inversion of the radiosity matrix [15]: where ε k is the emissivity of surface k. The net radiative flux density lost by the surface k is finally: where σ is the Stefan-Boltzmann constant.

Crucibles Walls
The walls of the hearth are water-cooled copper panels. If a determination of the heat transfer coefficient for a casting mold is readily available [16,17], there is no accurate analysis in the literature dedicated to water-cooled refining hearths. A Fourier condition is applied to model the heat resistance between the solid steel, copper and the cooling water as: where • h eq : equivalent heat transfer coefficient W·m 2 ·K −1 • T water : temperature of the cooling water (K) The water temperature is kept constant (T water = 20 • C). The heat transfer coefficient should account for the solidification shrinkage and gap formation at the metal/crucible interface. In such conditions, the maximum value of h eq can be estimated at 200 W·m 2 ·K −1 [18]. Furthermore, since the initial load is composed of steel balls, the contact between it and the crucible is very bad, leading to a weak heat transfer. Therefore, we have chosen an average value of the heat transfer coefficient equal to 50 W·m 2 ·K −1 at all the crucible walls.

Thermophysical Properties and Numerical Solution
The density and specific heat capacity of the stainless steel are kept constant with the temperature. According to the thermal properties reported in [19], the thermal conductivity linearly decreases from room temperature up to 800 • C. Hence, the conductivity halves over the temperature range as shown in Table 5, which gathers the set of data used in the simulation. Table 5. Thermodynamic values of steel alloy, data from [19].
The heat transfer Equation (3) is solved using a home-made finite volume code written in Fortran 90 using an implicit scheme with a time step equal to 1 s. The structured mesh is composed of N x × N y × N z (t) rectangular parallelepipedic cells. The total number of cells is equal to 256,000 at the start of the melting operation and reaches 768,000 cells at the end.
A global thermal balance was carried out to control the accuracy of the numerical procedure. All the components of this thermal balance are therefore calculated at each time step and the deviation of the balance is expressed as: error (%) = P electron + P drop + P rad + P cooling − ∆H acc ∆t P electron + P drop + P rad + P cooling × 100 (23) In the simulation reported below, this error remains lower than 2% at every moment.

Results
This chapter is divided into three sub-sections. The first examines the experimental results and the comparison with the numerical results is presented. Because the numerical simulation allows a deeper thermal analysis, a detailed investigation is reported in the second sub-section with a particular attention to the process heat balance. In order to guarantee the formation of a liquid pool at the top of the charge, a simulation with a higher EB power is eventually presented and analyzed.

Experimental Results and Comparison with the Numerical Simulation
The skull was easily removed from the copper crucible after cooling. Figure 9 reveals that the steel ball layer was not remelted and that the infiltration of liquid metal through the porous bed was hardly noticeable. Moreover, the top surface of the skull was far from flat and horizontal, and the presence of a candle of solid metal reveals that the liquid steel heated at the tip of the bar was partly solidified before being in contact with the top of the charge.
presence of a candle of solid metal reveals that the liquid steel heated at the tip of the bar was partly solidified before being in contact with the top of the charge.
Based on that observation, we can conclude that the remelted steel at the tip of the bar spread on the surface of the charge but also quickly solidified. According the EB thermal power condition prevailing in this experiment, the fill level was not uniform (as it is clearly observed in Figure 9b) which is not well modelled through the adaptive mesh. The calculated height of the charge at the end of the remelting test (z = 8.25 cm) is in good agreement with the average value measured (see Figure 9a). This clearly indicates that, as assumed in the simulation, the liquid steel did not infiltrate the raw charge. Figure 10 confirms that the maximum temperature calculated by the numerical model remains lower than the liquidus temperature, i.e., the EB power applied at the top of the charge did not allow for the formation of an actual liquid pool. Based on that observation, we can conclude that the remelted steel at the tip of the bar spread on the surface of the charge but also quickly solidified. According the EB thermal power condition prevailing in this experiment, the fill level was not uniform (as it is clearly observed in Figure 9b) which is not well modelled through the adaptive mesh.
The calculated height of the charge at the end of the remelting test (z = 8.25 cm) is in good agreement with the average value measured (see Figure 9a). This clearly indicates that, as assumed in the simulation, the liquid steel did not infiltrate the raw charge. Figure 10 confirms that the maximum temperature calculated by the numerical model remains lower than the liquidus temperature, i.e., the EB power applied at the top of the charge did not allow for the formation of an actual liquid pool. As stated previously, ten thermocouples have been installed within the steel ball layer (their location is shown in Figure 3). The actual and simulated thermograms are reported in Figures 11 and  12 respectively, for the 40 min of the melting trial. Again, all temperatures are much lower than the steel liquidus temperature, which implies that no liquid pool was formed during this melting operation. At the beginning of the experiment, the thermocouples c 6 and c 8 reached higher temperatures than the other ones. It could be explained by some rare infiltrations of metal within the initial bed of steel balls. Obviously because the infiltrations are not modelled, their thermal effect cannot be obviously observed in the numerical results. The other thermocouples display a similar dynamic behavior. We can observe a significant time-lag (around 10 min) required for the heat to diffuse from the top surface to the location of the thermocouples, followed by a temperature increase at a rate lying between 35 and 50 °C/min. Finally, the maximum temperature (800-1050 °C) was reached at t = 25 min. The calculated curves (Figure 12) match relatively well with these dynamics, showing a characteristic S-shape for a diffusive transfer. This is particularly true for the group of thermocouples (tc3, tc7, tc9 and tc10) where the dynamic behaviour (heating rate around 50 °C/min) and the maximum temperature (around 1000 °C) are well simulated. The agreement appears however As stated previously, ten thermocouples have been installed within the steel ball layer (their location is shown in Figure 3). The actual and simulated thermograms are reported in Figures 11 and 12 respectively, for the 40 min of the melting trial. Again, all temperatures are much lower than the steel liquidus temperature, which implies that no liquid pool was formed during this melting operation. At the beginning of the experiment, the thermocouples tc 6 and tc 8 reached higher temperatures than the other ones. It could be explained by some rare infiltrations of metal within the initial bed of steel balls. Obviously because the infiltrations are not modelled, their thermal effect cannot be obviously observed in the numerical results. The other thermocouples display a similar dynamic behavior. We can observe a significant time-lag (around 10 min) required for the heat to diffuse from the top surface to the location of the thermocouples, followed by a temperature increase at a rate lying between 35 and 50 • C/min. Finally, the maximum temperature (800-1050 • C) was reached at t = 25 min. The calculated curves ( Figure 12) match relatively well with these dynamics, showing a characteristic S-shape for a diffusive transfer. This is particularly true for the group of thermocouples (tc 3 , tc 7 , tc 9 and tc 10 ) where the dynamic behaviour (heating rate around 50 • C/min) and the maximum temperature (around 1000 • C) are well simulated. The agreement appears however weaker for another group of thermocouples (tc 1 , tc 2 , tc 4 and tc 5 ) since the time-lag and maximum temperatures do not match properly. This could be explained by the mediocre filling of the crucible with the incoming partly liquid steel.

Analysis of the Numerical Results
The numerical simulation calculates the 3D temperature field at any time of the melting process as shown in Figure 13. An increase in the surface temperature with time is noticed, and the maximum surface temperature is located just below the edge of the steel bar (i.e., approximately above the thermocouple tc 8 ).

Analysis of the Numerical Results
The numerical simulation calculates the 3D temperature field at any time of the melting process as shown in Figure 13. An increase in the surface temperature with time is noticed, and the maximum surface temperature is located just below the edge of the steel bar (i.e., approximately above the thermocouple tc8). Associated with this thermal field, Figure 14 reports the global heat balance of the steel charge into the crucible, as assessed by the numerical simulation. The dotted green line and dotted purple line draw the thermal heat provided by the electron beam on the metal surface when the stationary melting rate is attained (4.01 kW) and on the bar (9.99 kW) respectively. Obviously, their sum is equal to the total beam thermal power: Associated with this thermal field, Figure 14 reports the global heat balance of the steel charge into the crucible, as assessed by the numerical simulation. The dotted green line and dotted purple line draw the thermal heat provided by the electron beam on the metal surface when the stationary melting rate is attained (4.01 kW) and on the bar (9.99 kW) respectively. Obviously, their sum is equal to the total beam thermal power: P max thermal = ηP elec = 0.7 × 20 = 14 kW (24) Thermal power of the EB (kW) 14.0 100% Thermal power lost by the steel bar (radiation) (kW) 0.9 6.4% Thermal power lost by the charge (water cooling and radiation) (kW) 8.0 57.1% Heating power of the charge (kW) 5.1 36.5%

Simulation with a Higher EB Power
Since the obvious purpose of a process such as VASM is to actually remelt the raw material, a numerical simulation with a higher value for the EB power has been investigated in order to get an actual liquid pool. The steel drops cast into the crucible brought an enthalpy flux (continuous red line) according to the equation: When the stationary melting rate is attained this enthalpy flux reaches P drop = 9.06 kW, so that the heat lost by radiation from the steel bar to the EB furnace corresponds to the difference between P max thermal,bar and P max drop , i.e., ∆P loss 0.9 kW. The heat lost by radiation at the surface of the charge (blue curve) and to the cooling system (turquoise curve) evolves following the same trend and are approximately of the same order of magnitude, with a maximum value of 4 kW each. The thermal accumulation, which is not reported in the graph, compensates this heat balance, since the thermal capacity of the charge continuously rises with the processing time (growth of the mass of the steel block and increase in its temperature). Table 6 gathers the items of the heat balance near the end of the melting trial, and shows that more than 50% of the thermal EB power is lost through the cooling circuit of the crucible and by radiation at the top surface, while less than 40% remains in the growing charge. Finally, this lab-scale melting test consumed 0.81 kWh per kg of remelted steel (P elec /V f us ) whereas the typical yield is around 1 kWh per kg for a single VAR [10].

Simulation with a Higher EB Power
Since the obvious purpose of a process such as VASM is to actually remelt the raw material, a numerical simulation with a higher value for the EB power has been investigated in order to get an actual liquid pool.
The new operating conditions called thereafter as "power+" differ from the experimental operating conditions labelled as "reference" case only by the total EB power (see Table 7), all the other conditions and properties remaining the same. The numerical results of this "power+" case simulation show that the initial load composed of steel balls is still not remelted since the top surface temperature of this ball bed remains below the liquidus temperature. The first drops of liquid steel (t = 10 min) solidify when contacting the steel balls but soon after the model predicts a liquid pool to be formed, as presented in Figure 15. The mass of the liquid bath increases slowly and continuously until reaching 1.5 kg, when the EB is switched off. The new operating conditions called thereafter as "power+" differ from the experimental operating conditions labelled as "reference" case only by the total EB power (see Table 7), all the other conditions and properties remaining the same. The numerical results of this "power+" case simulation show that the initial load composed of steel balls is still not remelted since the top surface temperature of this ball bed remains below the liquidus temperature. The first drops of liquid steel (t = 10 min) solidify when contacting the steel balls but soon after the model predicts a liquid pool to be formed, as presented in Figure 15. The mass of the liquid bath increases slowly and continuously until reaching 1.5 kg, when the EB is switched off. Figure 15. Predicted mass of the initial load (in black), of the liquid pool (in red) and of the solidified metal including the mushy zone (in blue) vs time. "power+" case.
At the end of the melting operation, the charge can be divided into three zones which are observed in the computed cross-section on Figure 16: the initial load made of steel balls, the solidified metal and eventually the liquid pool. The latter only involves only 6.4% of the total mass of the charge (see Table 8). This value is still far from the VASM industrial conditions where a final mass fraction of liquid pool of a third is expected (since the new consumable electrode is the skull of the previous melting cycle [8]). At the end of the melting operation, the charge can be divided into three zones which are observed in the computed cross-section on Figure 16: the initial load made of steel balls, the solidified metal and eventually the liquid pool. The latter only involves only 6.4% of the total mass of the charge (see Table 8). This value is still far from the VASM industrial conditions where a final mass fraction of liquid pool of a third is expected (since the new consumable electrode is the skull of the previous melting cycle [8]). Metals 2020, 10, x FOR PEER REVIEW 19 of 21 Figure 16. Computed cross-section of the crucible ( = 6.25 cm) at the end of melting (t = 40 min), showing the initial load, solidified steel and liquid pool predicted by the numerical model. "power+" case. Such a difference is readily explained by the small dimensions of the lab-scale crucible, which enhance cooling at the crucible walls. Indeed, the ratio of the wall surface and the volume yields: which decreases when scale-up the crucible size. The corresponding thermal balance is reported in Table 9. Not surprisingly, comparison of Tables 6 and 9 indicates that an increase in the EB power results in an increase in the accumulation of enthalpy, correlated to the formation of the liquid pool. Table 9. Heat thermal balance at t = 40 min (power+ case).
Thermal power of the EB (kW) 25.2 100% Thermal power lost by the steel bar (radiation) (kW) 1.6 6.3% Thermal power lost by the charge (water cooling and radiation) (kW) 13.2 52.4% Heating power of the charge (kW) 10.4 41.3% The ratio 41.3%/52.4% reported in Table 9 would then grow by increasing crucible size because of the aspect ratio mentioned earlier. The numerical modelling sounds a good approach for the assessment of the heat efficiency and liquid pool volume at the industrial scale.

Conclusions
In order to represent the thermal transfers within the liquid and solid parts of a metallic piece contained in a rectangular crucible continuously fed with liquid metal, a 3D numerical model has been set-up and applied to simulate a remelting test of stainless steel using a lab-scale electron beam melting (EBM) furnace. In addition to the evolution of the 3D temperature field with time, the model  Such a difference is readily explained by the small dimensions of the lab-scale crucible, which enhance cooling at the crucible walls. Indeed, the ratio of the wall surface and the volume yields: which decreases when scale-up the crucible size. The corresponding thermal balance is reported in Table 9. Not surprisingly, comparison of Tables 6  and 9 indicates that an increase in the EB power results in an increase in the accumulation of enthalpy, correlated to the formation of the liquid pool. The ratio 41.3%/52.4% reported in Table 9 would then grow by increasing crucible size because of the aspect ratio mentioned earlier. The numerical modelling sounds a good approach for the assessment of the heat efficiency and liquid pool volume at the industrial scale.

Conclusions
In order to represent the thermal transfers within the liquid and solid parts of a metallic piece contained in a rectangular crucible continuously fed with liquid metal, a 3D numerical model has been set-up and applied to simulate a remelting test of stainless steel using a lab-scale electron beam melting (EBM) furnace. In addition to the evolution of the 3D temperature field with time, the model calculates