CFD Simulation of Forced Recirculating Fired Heated Reboilers

: An advanced algorithm has been developed in order to analyze the performance of re-boiling process of crude oil flowing inside reboilers tubes. The proposed model is composed from Heptane fire heater and a tube array. The heat flux produced from burner is transferred to the crude oil flowing inside the tube. The computational model is composed of two phases—Simulation of fire by using Fire Dynamics Simulator software (FDS) version 5.0 and then a nucleate boiling computation of the crude oil. FDS code is formulated based on CFD (Computational Fluid Dynamics) of fire heater. The thermo-physical properties (such as: thermal conductivity, heat capacity, surface tension, viscosity) of the crude oil were estimated by using empirical correlations. The thermal heat transfer to evaporating two-phase crude oil mixture occur by bubble generation at the wall (nucleate boiling) has been calculated by using Chen correlation. It has been assumed that the overall convective heat transfer coefficient is composed from the nucleate boiling convective coefficient and the forced turbulent convective coefficient. The former is calculated by Forster Zuber empirical equation. The latter is computed from the Dittus-Boelter relationship. In order to validate the nucleate boiling heat transfer coefficient, a comparison has been performed to nucleate boiling convective coefficient obtained by Mostinski equation. The relative error between the nucleate boiling convective heat-transfer coefficients is 10.5%. The FDS numerical solution has been carried out by using Large Eddy Simulation (LES) method. This work has been further extended to include also the structural integrity aspects of the reboiler metal pipe by using COMSOL Multiphysics software. It was found out, that the calculated stress is less than the ultimate tensile strength of the AISI 310 Steel alloy.


Introduction
Several works have been written about crude oil boiling. Sciance [1] has conducted experiments on nucleate and film saturated pool boiling for several hydrocarbons, such as: methane, ethane at twelve reduced pressures (the ratio between the pressure and the critical pressure). These experiments have been performed between 0.02 (atmospheric pressure) and 0.9 (close to critical pressure). He has been obtained complete burnout flux (critical heat flux) vs. reduced pressure curves for each substance. The Critical Heat Flux (CHF) equation derived by Noyes [2] describes these experimental data with reasonable accuracy over the entire pressure range. Zhu [3] has applied Wall-Resolved Large Eddy Simulation (WRLES) in order to study non-reacting and reacting flows in both smooth and ribbed tubes, to show the impact of the ribs on turbulence and chemistry. The CFD (Computational Fluid Dynamics) simulations were performed with the AVBP code. This code is a parallel CFD code that solves the 3D compressible momentum equations for turbulent flows, with coupling of reduced chemistry scheme of ethane and butane cracking for reacting cases.
This work includes analysis of Fire heater and Crude boiling process. The computational process algorithm is composed of four components. Fire Dynamics Simulator (FDS) outputs a file which contains the Heptane flame simulation data, STNAJAN software which evaluates the thermophysical properties of the burned gaseous mixture (thermal conductivity, kinematic viscosity etc.), COMSOL Multi-physics transient axi-symmetrical model which performs Heat transfer and structural analyses, and the main program which calculates the nucleate boiling convective coefficient and the external flow convective coefficient of the gaseous mixture produced from Heptane burning. It also includes the radiation heat flux emitted by the soot particles. The COMSOL Multiphysics software have been loaded with the radiative and convective heat fluxes emitted by gaseous and soot mixture.

Fired Heated Furnaces
In order to utilize the heat released by combustion, it must be set free in an enclosed space made of brick or other resistant material, known furnace. Furnaces may in general be divided into two types [4]. In the first type, the combustion is transferred by direct contact to the hot gaseous mixture produced by the combustion. In the second type, the heat transmitted to the crude oil by utilizing a heat exchanger. The combustion exhausts gases products, are completely separated from the crude oil tubes. The former type of furnace has the great advantage of cheaper construction and lower cost of maintenance and operation. It is adapted to maintain a higher temperature. The disadvantages are that the tubes may be contaminated with soot particles, smoke, or other products. The steel may also be exposed to oxidizing species such as steam, carbon monoxide, or carbon dioxide. If fuel is contaminated with Sulphur, these tubes may be exposed to Sulphur acid. On the other hand, the second type-heat exchanger has the obvious advantage of filtration of soot or smoke, more uniform control of the temperature and the chemical conditions of the existing combustion products. The disadvantages embedded in second type are heat losses and pressure losses inside the tube pipe. The maintenance is more expensive. The steel tube may be damaged by each of following mechanisms [5]: Cavitation damage, corrosion fatigue, high temperature corrosion, pitting corrosion, stress corrosion cracking, thermal creep, and fouling (water boiling surfaces usually acquire deposits, often referred to as "crud").
In the framework of this study, an analysis has been performed of the first type device. It should be noted that the algorithm described here is robust and can be adapted to simulate the second type furnace (The convective heat transfer coefficient obtained in Section 2.6 can be adapted to handle with steam or other high temperature heat carriers such as DOWTHERM A [6-9]). It is assumed that the thermal energy released by the combustion of fuel is transmitted to crude oil flowing inside the tubes. The furnace shown in Figure 1 is composed of three major components: the radiant section, the convection section, and the stack [10]. It is fueled by Heptane liquid fuel. The crude oil liquid, flowing inside the tubes, absorbs the heat mostly by radiative and convective heat transfer from the flue gases [11]. The height of the furnace is 8 m above the ground [11]. Initially, the flame velocity is normally high enough to produce large turbulent eddy motion. Friction occurs between the boundary of the flame and its surroundings, causing the surrounding air to be locally accelerated to the flame velocity. The accelerated air is then mixed with the flame and expanding it. This process is called entrainment [11]. The mixing is accomplished by turbulent eddy motion [12]. At the end of this process, the velocity of the flame and surroundings become homogenous. The analysis of such flames is based on mixing process only. In industrial furnaces, rapid fuel combustion occurs. The combustion of flame fuel in industrial furnaces as well as engines of various types is rapid process [12]. The heat released from fuel burning and carried by exhaust gaseous mixture and soot particles is transferred to reboiler tubes walls in by three main modes: forced and natural convection, thermal radiation, and conduction [13,14].

Reboiler
Reboilers are applied to provide the vapor flux to feed the bottom tray of the distillation column. The liquid exiting from the bottom of the column is partially vaporized in heated fired furnace. It is assumed that the heating medium is flue gaseous mixture produced by Heptane burning. There are three principal types of reboilers utilized in the distillation columns [15]: (1) Forced Recirculation Reboiler (2) Thermosiphon Natural Circulation Reboiler (3) Kettle type (4) Forced Recirculating Fired Heated Reboiler (see Figure 2) Distillation bottom crude oil liquids are often mixtures having substantial different boiling ranges. Thermo-physical properties of the liquid and vapor exhibit large variations throughout the reboiler. Thermodynamic calculations are required to determine the phase compositions and other properties within the reboiler. The next sections give theoretical fundamentals of the vapor-liquid equilibrium, the boiling regimes, mechanism of the nucleate boiling, and the thermo-physical properties evaluation of the crude oil.

Two-Phase Flow Heat Transfer Inside the Reboiler The Definition of Flow Regimes
The main characterizing feature of two-phase flows is the fact that an interface exists between the two phases and, in gas and liquid flows, this interface takes a wide variety of forms. In general the surface tension effects tend to create curved interfaces leading to spherical shapes (e.g., droplets or bubbles). As the discontinuities in the continuous phase are increased, the departure from a spherical shape is also increased. Thus, small droplets tend to be spherical whereas bigger ones are often deformed in the gas flow, and likewise the bubble. For vertical flows, the regimes to be defined here are illustrated in Figure 3. The various regimes are defined as follows [16]: (1) Bubble flow-the bubbles dispersed in a continuum of liquid (2) Slug or plug flow-as the concentration of bubbles increases, the bubble coalescence occurs and, progressively, the bubble diameter approaches that of the tube. Once this occurs, the bubble is bullet-shaped. (3) Churn Flow-as the gas flow is increased and the velocity of the slug-flow bubbles increased, the break-down of these bubbles occurs, leading to unstable regime. For narrow-bore tubes the oscillations may not occur and a smoother transition between the slug and annular flow may be observed. (4) Annular flow-in this regime, the liquid flows on the wall of the tubes as a film and the gas phase flows in the center. Some of the liquid phase is entrained as small droplets in the gas core. (5) Wispy annular flow-for high mass velocity flows, the droplet concentration in the gas core of annular flow increases and, ultimately droplet coalescence occurs, leading to streaks or "wisps" of liquid occurring in the gas core.

Calculation of Nucleate Boiling Convective Coefficient
When heat is gradually transferred to the crude oil, several processes are taking place. As the temperature of the crude oil exceeds the boiling point, it begins to evaporate.
At first, individual bubbles form at the heating surface (saturation state). When they reach a particular diameter, they are breaking away and as a result, there is intensive agitation of the boundary layer and the whole mass of boiling liquid. As the heat flux increases the number of crude oil bubbles formed at the heating surface increases. This leads to a further increase in turbulence, and a rapid rise in the heat transfer coefficient. The boiling regime occurring during bubble generation in the vicinity of the tube wall surface is called "nucleate boiling".
In this boiling regime, high values of convective heat transfer are obtained with relatively small temperature gradient. Nucleate boiling is intensified as the wall temperature increases but is suppressed as the fluid velocity increases. As the forced convection component increases, the nucleate boiling component decreases (by suppression). The most widely used empirical equation which combines these two processes has been derived by Chen [17]. It takes the form: Here B h denotes the two-phase convective boiling heat transfer coefficient (W/m 2 K), FC h denotes the forced convective boiling heat transfer coefficient (W/m 2 K), and NB h is the nucleate boing convective. FZ h denotes the nucleate boiling convective heat transfer coefficient (W/m 2 K). It is calculated by using the Forster-Zuber empirical equation (see Equation (7)). This correlation gives good results when used in Equation (1). NB S is the suppression factor which account for the reduction in the coefficient due to the suppression of nucleation sites. It is related to the total Reynolds number as shown in Figure 4 [18]. The suppression factor SNB is a function of a two-phase Reynolds number, ReTP: The two-phase Reynolds number is calculated according to the following equation: 1.25 Re Re Here Re L denotes the liquid phase Reynolds number. It is defined by: c F denotes the Reynolds number factor. This factor accounts for the enhanced flow and turbulence due to the presence of vapor. This term can be approximated by [18]: c F is graphically presented in Figure 5 [18]. tt X denotes the Martinelli parameter. This term is calculated by the following equation [18]: The correlation of Forster and Zuber (1955) where sat T ∆ is the wall superheat (K). It is the wall temperature minus saturation temperature. sat p ∆ is the difference in saturation pressure (N/m 2 ) corresponding to the temperature difference sat T ∆ . L k is the thermal conductivity of the liquid (W/m K). pL c is the heat capacity of the liquid phase (J/kg K). L ρ is the density of the liquid phase (kg/m 3 ). G ρ is the density of the vapor phase (kg/m 3 ). L µ is the viscosity of the liquid phase (kg/m sec). σ is the surface tension (N/m). fg h ∆ is the enthalpy of vaporization (J/kg). The forced convective coefficient is derived from the Dittus-Boelter relationship [18]: The nucleate boiling convective coefficient hNB is related to the convective coefficient calculated for the mixture properties (hNBM) as follows [15]: y denotes the vapor mole fraction of the lighter component equilibrium with the liquid mole fraction of the component x (see Figure 6). The mole fractions x and y can be obtained by Raoult's rule. This rule states that the partial pressure of a component over a solution is the product of the vapor pressure of that component and the mole fraction of this component [20]. The parameter As is calculated by the following equation: Here p denotes the system pressure in bars. Ao denotes a constant which depends on the binary system. Stephan and Körner give a table for some combinations. They recommend a value of this term should be 1.5 if no other information is available. This method for binary mixtures has been extended to multicomponent mixtures by Stephan and Preusser [21]. The crude oil refinery is composed of several unit operations components, such as: reboilers, heat exchanger, and pumps. Analysis and design of the performance of these components require knowledge of the thermodynamic and physical properties of the crude oil fluids [22]. As described in the previous section designing a reboiler requires knowledge of thermodynamic and thermophysical properties crude oil. It will be shown later that these are: boiling temperature, pressure, heat capacity, thermal conductivity, diffusivity, and density. Heat exchanger design depends on enthalpies, thermal conductivities, and viscosity of the flowing streams (the crude oil and the steam).

Estimation of the Thermo-Physical Properties of the Crude Oil
Due to the complexity of the composition crude oils, it is not possible to measure or calculate accurately all of these properties [22]. Furthermore, calculation methods developed for pure hydrocarbons are not always applicable. Over the years, the petroleum refining engineers have developed empirical correlations in order to estimate the properties of crude oil fraction by measuring properties such as specific gravity (SG) and normal boiling point (NBP). Another important term connected with SG is called API gravity. This term is a measure of the relative density of a crude oil liquid and the density of water (i.e., how heavy or light a crude oil liquid is compared to water) [23]. These correlations require minimum input data, and are usually presented in graphical form. These quantities are described in the following section.
In order to calculate the nucleate boiling heat flux, it is essential to estimate the thermo-physical properties of the Crude oil. These properties are described in this section. The thermal conductivity of the crude oil as a function of the temperature is shown in Figure 7   The density of liquid crude oil is shown in Figure 8. The viscosity of liquid crude oil is shown in Figure 9. The latent heat of the crude oil is shown In Figure 10   °C as a function of API (see Figure 11) [26]. It can be seen from this Figure that the surface tension decreases with the API and the temperature. The viscosity of vapor crude oil is shown in Figure 12. Riazi and Al-Sahhaf [27] evaluated different crude oil properties such as the boiling point, density, refractive index, critical temperature, pressure, and density as a function of the molecular weight by using the following general equation: Here θ can be any one of the properties mentioned above. M is the molecular weight and θ ∞ is the limiting value for any property as M → ∞ . This empirical equation can be applied in order to calculate the following properties. b T is the mean average boiling point in K, SG is the specific gravity and c p is the critical pressure in bar. The constants a, b, and c for each property are given in Table 1 [22,27].

Fire Dynamic Simulation (FDS) Modeling of the Fired Heater
The FDS has been developed at the National Institutes of Standards and Technology (NIST) [28,29]. It solves simultaneously the classical conservation equations, which includes the momentum, energy, mass species equations and thermodynamics-based state equation of a perfect gas. The chemical combustion reaction data are taken from a library of different fuel sources. By solving these equations, it calculates the velocity, pressure by solving the momentum equation), temperature (by solving the Energy equation), density (by solving the Equation of state), and chemical composition (by solving the Diffusion equation) within each numerical grid cell as a function of time. It is capable to calculate the heat flux. FDS is composed of three major components: Hydrodynamic Model-it is formulated based on of fire-driven fluid flow (i.e., natural convection). The numerical solution is carried out using Large Eddy Simulation (LES). This technique is applied for relatively low speeds (or low Reynolds number).
Combustion Model-two types of combustion models are employed in FDS software. The first (regarded as the default) employs the mixture fraction, a quantity representing the fuel and the products of combustion. The second model makes use of the specified Arrhenius reaction parameters for individual gas species. The mass fractions of all of the major reactants and products can be calculated from the mixture fraction by using means of "state relations", empirical correlations obtained by a combination of simplified analysis and measurements [30].
Radiation Transport-radiative heat transfer is obtained by solving of the radiation transport equation (RTE) for a non-scattering gray gas. In a several number of cases, a wide band model can be applied in place of the gray gas model. A Finite Volume Method (FVM) technique is applied in FDS software in order to solve this equation [29].
The postprocessor of this software is named "smoke-view". As we shall see in Section 3, it is capable to perform post-processing images of the temperature, velocity, and soot fields. This code has been employed in this algorithm in order to simulate the Heptane combustion performance.

Governing Transport Equations of FDS Software
The following subsections describe the conservation momentum, mass (continuity), diffusion, and energy transport equations for the Newtonian fluid.

Mass and Species Transport Equations
Mass conservation is expressed as a function of the density, ρ [28]: where D α represents the diffusion coefficient of α component of the mixture [m 2 /s].

Momentum Transport Equation
The momentum equation is described in the following equation [28]: ( ) The term Sij represents the symmetric rate of strain tensor. μ denotes the dynamic viscosity of the fluid. In this work the Heptane fire heater simulation has been carried out by applying LES method, in which the large-scale eddies are calculated directly, and the subgrid-scale dissipative processes are modeled. FDS uses LES in most fire simulation applications. The following subsection contains a description of how the transport dissipative terms are modeled in FDS software.

Large Eddy Simulation (LES) Turbulence Method
Large Eddy Simulation (LES) is applied in FDS software in order to model the dissipative terms (such as: viscosity, thermal conductivity, and material diffusivity) which occur at length scales smaller than those that are explicitly resolved on the numerical computational grid cell. They are replaced by surrogate expressions that "model" their impact on the approximate form of the governing equations (It means that these properties cannot be applied directly in most practical simulations). The dissipation rate, ε , [Pa/s] is the rate at which kinetic energy is converted to thermal energy by the viscosity of the gaseous mixture [28]: The viscosity μ is computed according to the following equation Here ∆ is a length on the order of the size of a grid cell and Cs represents the empirical constant. The bar above the quantities appearing in Equation (18) It is assumed that the two terms, turbulent Schmidt number t Sc (defined as the ratio of momentum diffusivity to mass diffusivity) and turbulent Prandtl number t Pr (defined as the ratio of momentum diffusivity to thermal diffusivity), are constant for a given scenario. The model for the viscosity, LES µ , serves two roles: first, it has the appropriate mathematical form to describe the dissipation of kinetic energy from the flow. Second, it provides a stabilizing effect in the numerical algorithm, damping out numerical instabilities as they arise in the flow field, especially where vorticity is generated.

Energy Transport Equation
The energy transport equation is described as a function of the sensible enthalpy, s h [J/kg] [28]: This term depends on the temperature (according to the ideal law assumption) [

Equation of State
The pressure is computed by applying the ideal gas equation of state [28]: where T is the temperature in [K], W is the molecular mass of the gaseous mixture in [J/mole], R is the gas constant, and p is the pressure in [Pa].

Fire Dynamics Simulation (FDS) Modeling of the Heptane Burner
The geometric model of the Heptane burner is shown in Figure 13. The height of the burner model is 2.4 m. It should be noted that most tube failures in fired heaters are caused by high tube metal temperatures (TMTs). Flame impingement on the tube external surface caused by long flames results in higher tube metal temperatures and coking that affect the performance of the heater and decrease the strength of the tubes (such as thermal creep) [31]. The length and the width of the burner are 1.2 m. The mesh size of the burner is: 27,468 cells. At the bottom of the burner, liquid Heptane is injected and ignited. The green points indicate the positions of the thermocouples and the soot particle sensors. The heat of combustion of Heptane is 44,600 (kJ/kg) [32]. The thermal conductivity of the liquid Heptane is 0.14 (W/(m·K)) [33]. Its specific heat is 2.24 (kJ/(kg·K)) [34], and its density is 679.5 (kg/m 3 ) [35]. The soot yield of Heptane is: 0.037 [36].
Initial condition-it is assumed that the initial temperature, the component concentration in the air and pressure are:

Turbulent Flow through Reboiler tube bundle
Schematics of the reboiler tube bundle system are shown at Figure 14. The hot gas mixture enters from above. The crude oil flows through the pipe and enters from the side. The tubes are assumed to be made of Stainless Steel. The arrows in Figure 14 indicate the flows directions of the crude oil and hot combustion products. The corresponding 2D domain appears in Figure 15. The arrows presented in Figure 16 indicate the flow direction the Heptane combustion products.

Convective Heat Transfer across Banks of Tubes
This section deals with heat transfer to bundle of tubes inside the reboiler. For flue gaseous mixture flowing across tube bundles composed of 10 or more rows, Grimison [37] has obtained a correlation of the form:  (27) According to Incropera et al. [38], the values of the constants C1 and m are: ,max Re D denotes the Reynolds number of the external flow. This term is defined in Equation (29):  (30): The value of the flue gas velocity, V, is obtained from FDS calculation of flame. The required transport properties of the flue gaseous mixture were calculated by using STANJAN software [39]. The required mole fraction of the products species of the oxidation reaction were calculated from the following reaction: The flame temperature of the Heptane flame was obtained by using the FDS software. The transport thermo-physical properties are listed in Table 2.

Calcaulation of the Radiative Heat Flux Emitted by the Mixture Flue Gaseous and the Soot Particles
The of radiative heat transfer within the fire heater enclosure are composed of two components: the high temperature burned gaseous mixture and the soot particles within the Heptane flame [13]. The burned gaseous mixture is composed mainly of the carbon dioxide and water vapor. They are not grey but emit and absorb radiation only in the discrete bands in the spectrum and therefore, have mean absorption coefficients which vary with gas temperature and the radiation source. The emissivity of the gaseous mixture can be obtained by adding the absorption coefficient for the gas and soot in each spectral line. The overall emissivity of the gas-soot mixture is given in the following equation [40]: The weighting factors and absorption coefficients for a three-term equation have been evaluated for the case of a mixture of the products of oil combustion. soot C is the soot concentration. It is obtained by applying FDS software (see Figure 23). The weighting factors,  Table 3. With these values, Equations (32-34) define the gas and flame emissivity in an oil fired furnace where the radiation path length is in the range of 0.2 m to 6 m, and the temperature is in the range of 1100 to 1800 K.

Thermal Structural Analysis-Reboiler Steel Tube
It is assumed that the reboiler tubes are made from Stainless Steel AISI 310. Figure 16 shows the geometry of the steel tube.The dimensions of the tube are given in Table 4. Table 4. Dimensions of the tube.

Value [m] D (outer diameter)
0.121 w (tube thickness) 0.006 AISI 310 stainless steel is a heat resistant alloy. It has moderately high temperature strength and oxidation resistance. The COMSOL software has been loaded with the thermomechanical properties of Stainless steel of AISI 3010 (thermal conductivity, heat capacity, density, Young modulus, coefficient of thermal expansion, and Poisson's ratio). The thermo-physical and thermomechanical properties of the AISI 310 stainless steel alloy are listed in Table 5.

Structure of the Compoutation Process
The computational process algorithm discussed in the previous section is summarized in Figure 16.
The algorithm is composed of four modules. FDS output file which contains the Heptane flame simulation data, Stanjan software which evaluates the thermo-physical properties of the burned mixture (thermal conductivity, kinematic viscosity etc.), COMSOL Multi-physics transient axisymmetrical model which performs Heat transfer and structural analyses [41], and the main program which calculate the nucleate boiling convective coefficient and the convective coefficient of the external flow (the gaseous mixture produced from Heptane burning). It provides the required heat loads for COMSOL Multiphysics computation. The schematics of the COMSOL computational process algorithm is shown in Figure 17.

Results.
This section is composed of three parts. The first part (Section 3.1) presents the thermal results of Fire Dynamics Simulation (FDS) software. The second part (Section 3.2) presents the calculation and validation results for the nucleate boiling convective coefficient of the crude oil. The last part (Section 3.3) composed of the COMSOL thermal structural analyses results for the reboiler tubes.

Fire Dynamics Simulator Software Results for Heptane Burner
FDS software has provided detailed information results regarding the Heptane burner performance. The results include the local and transient gas velocity, gaseous mixture temperature, soot concentrations, radiative heat flux, and convective heat flux. The temperature field at t = 80.5 s is shown in Figure 18 (It will be shown later that in this time interval, the radiative heat flux is maximal). It can be seen from this figure that the maximal temperature at time = 80.5 s reaches to 915 °C. The peak temperature results obtained by FDS software (915 °C) has been compared to the literature. It is claimed that typical post-flashover room fire will have more commonly a peak temperature of 900 to 1000 °C (1652 to 1832 °F) [31]. At the flame center the temperature is higher. The velocity field of flue gases at t = 80.5 s is presented in Figure 19 (It will be shown later that in this time interval, the radiative heat flux is almost maximal). From Figure 19 it can be seen that gaseous mixture flows upwards. The flow movement is controlled by buoyancy. The soot mass fraction field is shown in Figure 20.  Figure 20 shows that at the flame core the soot concentration is higher. The soot flows upward because of the buoyancy forces. Figure 21 shows the radiative and total heat flux emitted by the heptane burner. From Figure 21 it can be seen that that the maximal radiation heat flux approaches to 120,000 (W/m 2 ). This model which is described in this work is robust. The Heat Release Rate (HRR) can be increased by replacing the fuel considered here with higher heating value. Five thermocouples and soot concentration sensors have been installed inside the burner model. They were installed at heights: 0.5 m, 1 m, 1.5 m, 2 m, and 2.5 m above the burner surface at the burner axis. Based on the collected data gathered from these thermocouples and soot sensor readings, two 3D plots have been constructed. The three dimensional plot of the temperature field as a function of the height and the time is shown in Figure 22. The three dimensional plot of the soot concentration field as a function of the height and the time is shown in Figure 23. The soot concentration and the temperature profiles as a function of the time are similar to the radiation heat flux profile (see Figure 21). In Section 2.7, the soot concentration field is applied in order to calculate the radiative heat transfer emitted by the soot particles.

Calculation and Validation Results for the Nucleate Boiling Convective Coefficient
After substitution of the thermodynamic values described in Section 2.2 in Equations (1) and (7) The reduced pressure correlation is simple to use and gives values that are as reliable as those given by more complex equations. Here: pc denotes the liquid critical pressure in [bar]. The value of the critical pressure has been taken from the

Thermal Structural Analyses Results for the Reboiler Tube
A thermal structural analyses study has been performed in order to analyze the influence of the radiative and convective heat flux on the reboiler tube structural integrity. The calculated external convective coefficient of gaseous mixture, obtained by the algorithm described in section 2.6 is 78.3 W/(m 2 ·K). The overall emissivity obtained by Equation (32) is 0.42 for distance of 6 m. The COMSOL software has been loaded with the thermal loads (convective and radiative heat fluxes). The numerical results for the temperature Von-Mises stress are presented in this section. Figure 24 shows the 3D temperature field inside the steel tube. From Figure 24 it can be seen that temperature gradient inside the steel is small. This is because the thermal conductivity of the steel is high. The 3D Von-Mises stress distribution of reboiler tube is shown in Figure 25. According to Figure 25, the maximal stress reaches to 1.1 GPa at the tube top corner. The stress at the middle section of the tube is 391 MPa. The tensile ultimate strength of several stainless steels alloys as a function of the temperature is shown in Figure 26 [43]. Figure 26. Plot of ultimate tensile strengths (UTS) of several stainless steels compared to the stress obtained by COMSOL software (horizontal line) shown as a function of temperature (The thermomechanical data were taken from Reference [43]).
This figure shows that the UTS decrease with the temperature. As can be seen from this figure, the calculated stress obtained in this work is less than the ultimate tensile strength of the AISI 310 Steel alloy. The structural integrity of the tube is maintained.

Discussion
An advanced simulation algorithm has been developed in order to analyze the performance of re-boiling process of crude oil flowing inside reboilers tubes. The computational process algorithm is composed of four components: FDS output file which contains the Heptane flame simulation data, Stanjan software which evaluates the thermo-physical properties of the burned mixture (thermal conductivity, kinematic viscosity etc.), COMSOL Multi-physics transient axi-symmetrical model which performs Heat transfer and structural analyses, and the main program which calculates the nucleate boiling convective coefficient and the convective coefficient of the external flow (the gaseous mixture produced from Heptane burning). It provides the required heat loads COMSOL Multiphysics. It also integrates the various components of the algorithm. The proposed model is composed of Heptane fire heater and the tube array. The heat flux produced from burner is transferred to the crude oil flowing inside the pipe. The computational model is composed of two phases-Simulation of fire by using Fire Dynamics Simulator software (FDS) version 5.0 and then a nucleate boiling computation of the crude oil. FDS code is formulated based on CFD of Heptane fire. The FDS has Large Eddy Simulation (LES) method. It calculates simultaneously the temperature, density, pressure, velocity, soot concentration, chemical composition within each numerical grid cell as a function of the time. It is capable to calculate the convective and radiative heat fluxes.
The thermo-physical properties (such as: thermal conductivity, heat capacity, surface tension, viscosity) of the crude oil were estimated by using empirical correlations. The thermal heat transfer to evaporating two-phase crude oil mixture occur by bubble generation at the wall (nucleate boiling) has been calculated by using Chen empirical correlation. It has been assumed that overall convective heat transfer coefficient is composed from the nucleate boiling convective coefficient and the forced turbulent convective coefficient. The former is calculated by Forster Zuber empirical equation. The latter is computed from the Dittus-Boelter relationship. The results obtained for the nucleate boiling convective coefficient has been compared to nucleate boiling convective coefficient obtained by Mostinski equation. It has been found out that the relative error between the nucleate boiling convective heat-transfer coefficients is: 10.5%. This work has been further extended to include also the structural integrity aspects of the reboiler metal pipe by using COMSOL Multiphysics software. It was found out, that the calculated stress is less than the ultimate tensile strength of the AISI 310 Steel alloy. The simulation described in this work revealed that the Heat Release Rate produced by Heptane burner can provide the necessary radiative and convective heat fluxes required for the evaporation of the crude oil (the calculated maximal temperature of the steel reached to 523 °C. This temperature is greater than the boiling temperature of the lube oil, motor oil, and grease). It has been found that the calculated stress obtained in this work is less than the ultimate tensile strength of the AISI 310 Steel alloy. The structural integrity of the tube is maintained.

Conclusions
A new tool for simulating crude oil boiling flowing inside reboiler tube has been proposed in this paper. This heat supplied to the crude oil is produced by using other heat source heptane furnace burner.
Funding: This research hasn't received external funding. No funding sponsors had any role in the numerical analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Conflicts of Interest:
The author declares no conflict of interest.