Flow regime transition criteria for two-phase flow in a vertical annulus

In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug–bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been checked against air–water as well as boiling flow regime maps. In all the cases, the new developed model shows better predicting capabilities than the existing correlations most used in literature.     

Microgravity flow regime data and analysis

To utilize the advantageous properties of two-phase flow in microgravity applications, the knowledge base of two-phase flow phenomena must be extended to include the effects of gravity. In the experiment described, data regarding the behavior of two-phase flow in a conduit under microgravity conditions (essentially zero gravity) are explored. Of particular interest, knowledge of the void fraction of the gas and liquid in a conduit is necessary to develop models for heat and mass transfer, pressure drop, and wall shear. An experiment was conducted under reduced gravity conditions to collect data by means of a capacitance void fraction sensor and high speed visual imagery. Independent parameters were varied to map the flow regime regions. These independent parameters include gas and liquid volumetric flow rates and saturation pressures. Void fraction measurements were taken at a rate of 100 Hz with six sensors at two locations along the conduit. Further, statistical parameters were developed from the void fraction measurements. Statistical parameters such as variance, signal-to-noise ratio, half height value, and linear area difference were calculated and found to have characteristics allowing flow regime identification.     

Flow regime transition criteria for two-phase flow at reduced gravity conditions

Flow regime transition criteria are of practical importance for two-phase flow analyses at reduced gravity conditions. Here, flow regime transition criteria which take the frictional pressure loss effect into account were studied in detail. Criteria at reduced gravity conditions were developed by extending an existing model from normal gravity to reduced gravity conditions. A comparison of the newly developed flow regime transition criteria model with various experimental datasets taken at microgravity conditions showed satisfactory agreement. Sample computations of the model were performed at various gravity conditions, such as 0.196, 1.62, 3.71 and 9.81 m/s² corresponding to micro-gravity and lunar, Martian and Earth surface gravity, respectively. It was found that the effect of gravity on bubbly–slug and slug–annular (churn) transitions in a two-phase flow system was more pronounced at low liquid flow conditions, whereas the gravity effect could be ignored at high mixture volumetric flux conditions. While for the annular flow transitions due to flow reversal and onset of droplet entrainment, higher superficial gas velocity was obtained at higher gravity level.     

Two-phase flow operational maps for multi-microchannel evaporators

The current paper presents new operational maps for several different multi-microchannel evaporators, with and without any inlet restrictions (micro-orifices), for the two-phase flow of refrigerants R245fa, R236fa, and R1234ze(E). The test fluids flowed in 67 parallel channels, each having a cross-sectional area of 100 × 100 μm². In order to emulate the power dissipated by active components in a 3D CMOS CPU chip, two aluminium microheaters were sputtered onto the back-side of the test section providing a 0.5 cm² each. Without any inlet restrictions in the micro-evaporator, significant parallel channel flow instabilities, vapor back flow, and flow maldistribution led to high-amplitude and high-frequency temperature and pressure oscillations. Such undesired phenomena were then prevented by placing restrictions at the inlet of each channel. High-speed flow visualization distinguished eight different operating regimes of the two-phase flow depending on the tested operating conditions. Therefore, the preferred operating regimes can be easily traced. In particular, flashing two-phase flow without back flow appeared to be the best operating regime without any flow and temperature instabilities.     

Probabilistic mapping of adiabatic horizontal two-phase flow by capacitance signal feature clustering

In order to better quantify two-phase flow regime transitions, a sensor is developed which measures the electric capacitance of two-phase flows. A large number of experiments are done with air–water flow in a 9 mm ID horizontal tube. Based on a multivariate analysis, the most suitable sensor signal parameters are selected for building a flow regime classifier. This classifier is based on a fuzzy c-means clustering algorithm together with a regression technique. The output of the algorithm is used to create a probabilistic flow regime map. A comparison between a visual classification based on high speed camera images and the outcome of the flow regime classifier shows a remarkable agreement. The flow regime transitions are further quantified and discussed based on the probabilistic information and the sensor signal characterization. Probabilistic mapping makes it possible to combine flow regime dependent correlations in the two-phase flow models for heat transfer and pressure drop with smooth and appropriately quantified transitions from one flow regime to another.     

Two-phase flow distribution of air–water annular flow in a parallel flow heat exchanger

The air and water flow distribution are experimentally studied for a round header – flat tube geometry simulating a parallel flow heat exchanger. The number of branch flat tube is 30. The effects of tube outlet direction, tube protrusion depth as well as mass flux, and quality are investigated. The flow at the header inlet is identified as annular. For the downward flow configuration, the water flow distribution is significantly affected by the tube protrusion depth. For flush-mounted configuration, most of the water flows through frontal part of the header. As the protrusion depth increases, more water is forced to the rear part of the header. The effect of mass flux or quality is qualitatively the same as that of the protrusion depth. Increase of the mass flux or quality forces the water to rear part of the header. For the upward flow configuration, however, most of the water flows through rear part of the header. The protrusion depth, mass flux, or quality does not significantly alter the flow pattern. Possible explanations are provided based on the flow visualization results. Negligible difference on the water flow distribution was observed between the parallel and the reverse flow configuration.     

Flow regime development analysis in adiabatic upward two-phase flow in a vertical annulus

In this work radial and axial flow regime development in adiabatic upward air-water two-phase flow in a vertical annulus has been investigated. Local flow regimes have been identified using conductivity probes and neural networks techniques. The inner and outer diameters of the annulus are 19.1 mm and 38.1 mm, respectively. The equivalent hydraulic diameter of the flow channel, D_H, is 19.0 mm and the total length is 4.37 m. The flow regime map includes 1080 local flow regimes identifications in 72 flow conditions within a range of 0.01 m/s < 〈j_g〉 < 30 m/s and 0.2 m/s < 〈j_f〉 < 3.5 m/s where 〈j_g〉 and 〈j_f〉 are, respectively, superficial gas and liquid velocities. The local flow regime has been classified into four categories: bubbly, cap-slug, churn-turbulent and annular flows. In order to study the radial and axial development of flow regime the measurements have been performed at five radial locations. The three axial positions correspond to z/D_H = 52, 149 and 230, where z represents the axial position. The flow regime indicator has been chosen as some statistical parameters of local bubble chord length distributions and self-organized neural networks have been used as mapping system. This information has been also used to compare the results given by the existing flow regime transition models. The local flow regime is characterized basically by the void fraction and bubble chord length. The radial development of flow regime shows partial and complete local flow regime combinations. The radial development is controlled by axial location and superficial liquid velocity. The radial flow regime transition is always initiated in the center of the flow channel and it is propagated towards the channel boundaries. The axial development of flow regime is observed in all the flow maps and it is governed by superficial liquid velocity and radial location. The prediction results of the models are compared for each flow regime transition.     

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included ground-based tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.     

Nitrogen/non-Newtonian fluid two-phase upward flow in non-circular microchannels

This paper presents experimental investigations on nitrogen/non-Newtonian fluid two-phase flow in vertical noncircular microchannels, which have square or triangular cross-section with the hydraulic diameters being D_h = 2.5, 2.886 and 0.866 mm, respectively, by visualization method. Three non-Newtonian aqueous solutions with typical rheological properties, i.e., 0.4% carboxymethyl cellulose (CMC), 0.2% polyacrylamide (PAM) and 0.2% xanthan gum (XG) are chosen as the working fluids. The common flow patterns are identified as slug flow, churn flow and annular flow. The dispersed bubble flow is only found in the case with nitrogen/CMC solution two-phase flow in the largest channel. A new flow pattern of nitrogen/PAM solution two-phase flow, named chained bubble/slug flow, is observed in all the test channels. The flow regime maps are also developed and the results show that the rheological properties of the non-Newtonian fluid have remarkable influence on the flow pattern transitions. The geometrical factors of the microchannel such as the cross-section shape and hydraulic diameter of the channel can also affect the flow regime map. Finally, the results obtained in this work are compared with the available flow pattern transitions.     

Two-phase flow patterns and size distribution of droplets in a microfluidic T-junction: Experimental observations in the squeezing regime

Generating micrometer sized droplets has been studied in a microfluidic system with T-junction geometry 250 μm in internal diameter and with PTFE capillary tubing. Several experiments were conducted by varying the flow rate of the dispersed phase from $2.78\cdot {10}^{-11}{\text{ m}}^3/\text{s}$ to $5.28\cdot {10}^{-9}{\text{ m}}^3/\text{s}$ and that of the continuous phase from $2.78\cdot {10}^{-10}{\text{ m}}^3/\text{s}$ to $1.94\cdot {10}^{-9}{\text{ m}}^3/\text{s}$. The visualization of different flow regimes (drop, plug, and annular) was carried out for three configurations (not inverted in a horizontal position, inverted in a horizontal position, and inverted in a vertical position) for low capillary numbers. The model of Gauss was also chosen for a droplet size distribution in the dispersed phase, with the flow quality x varying from 0.016 to 0.44. The evolution of the drop size distribution as a function of the flow quality in the dispersed phase shows that the variation coefficient of the droplet's diameter is inversely proportional to the flow quality.     

Experimental investigation of forced flow regime transition in a dual bell nozzle by secondary fluidic injection

Experiments on an axisymmetric dual-bell nozzle were performed at EDITH nozzle test facility of CNRS in Orléans, France. The main purpose of the study was to explore the possibility of controlling the flow regime transition by a secondary fluidic injection in the dual bell nozzle. The main focus of the present paper is to investigate the impact of the secondary injection parameters on the flow regimes transition in such nozzles. Secondary injection has been found to effectively control the flow regime transition and consequently to increase the propulsive performance of the device. It has also been pointed out that even a very low injected secondary mass flow rate leads to the control of the transition and contributes to reducing the lateral loads which can exist, moreover, when transitions are operated without injection.     

Two-phase flow instabilities in a silicon microchannels heat sink

Two-phase flow instabilities are highly undesirable in microchannels-based heat sinks as they can lead to temperature oscillations with high amplitudes, premature critical heat flux and mechanical vibrations. This work is an experimental study of boiling instabilities in a microchannel silicon heat sink with 40 parallel rectangular microchannels, having a length of 15 mm and a hydraulic diameter of 194 μm. A series of experiments have been carried out to investigate pressure and temperature oscillations during the flow boiling instabilities under uniform heating, using water as a cooling liquid. Thin nickel film thermometers, integrated on the back side of a heat sink with microchannels, were used in order to obtain a better insight related to temperature fluctuations caused by two-phase flow instabilities. Flow regime maps are presented for two inlet water temperatures, showing stable and unstable flow regimes. It was observed that boiling leads to asymmetrical flow distribution within microchannels that result in high temperature non-uniformity and the simultaneously existence of different flow regimes along the transverse direction. Two types of two-phase flow instabilities with appreciable pressure and temperature fluctuations were observed, that depended on the heat to mass flux ratio and inlet water temperature. These were high amplitude/low frequency and low amplitude/high frequency instabilities. High speed camera imaging, performed simultaneously with pressure and temperature measurements, showed that inlet/outlet pressure and the temperature fluctuations existed due to alternation between liquid/two-phase/vapour flows. It was also determined that the inlet water subcooling condition affects the magnitudes of the temperature oscillations in two-phase flow instabilities and flow distribution within the microchannels.     

Two-phase damping for internal flow: Physical mechanism and effect of excitation parameters

Two-phase flow induced-vibration is a major concern for the nuclear industry. This paper provides experimental data on two-phase damping that is crucial to predict vibration effects in steam generators. An original test section consisting of a tube subjected to internal two-phase flow was built. The tube is supported by linear bearings and compression springs allowing it to slide in the direction transverse to the flow. An excitation system provides external sinusoidal force. The frequency and magnitude of the force are controlled through extension springs. Damping is extracted from the frequency response function of the system. It is found that two-phase damping depends on flow pattern and is fairly proportional to volumetric fraction for bubbly flow. Measurements are completed by the processing of high-speed videos which allow to characterize the transverse relative motion of the gas phase with respect to the tube for bubbly flow. It is shown that the bubble drag forces play a significant role in the dissipation mechanism of two-phase damping.     

An application of the wavelet analysis technique for the objective discrimination of two-phase flow patterns

This paper presents an application of the wavelet analysis technique for two-phase flow pattern identification by using the void fraction signals obtained from a multi-channel Impedance Void Meter (IVM) in a vertical-upward air–water flow. A new method for the objective discrimination of the two-phase flow pattern has been developed to provide information regarding the local energy of void fraction signals at a given scale on the joint time–frequency diagram. The void signals are processed with Continuous Wavelet Transform (CWT) to get the local wavelet energy coefficients map on the time–frequency diagram. The effective local wavelet energy and the effective scale are then calculated. Then the criteria for flow pattern identification are, finally, obtained. A series of void fraction measurements were conducted over a wide range of air–water vertical-upward flow condition to provide an extensive database to cover several types of flow patterns. The results show that the proposed method has a high precision for characterizing different flow regimes in two-phase flow, and is considerably more promising for the online recognition of two-phase flow patterns due to the short time of data processing.     

Two-phase gas-liquid flow regimes for smooth and ribbed rectangular ducts

In this research, adiabatic two-phase air-water flow was investigated, and results for smooth and ribbed rectangular ducts are presented here. The test fluids were air and water at approximately atmospheric conditions. Three ribs of different heights were used; the rib width and pitch were held constant. The ribs were positioned in the duct at three different locations to establish three different conditions: on the bottom wall (water side), on the top wall (air side) and on both the top and bottom walls. The flow regimes in the smooth and ribbed ducts, which were recorded with a video camera, were classified as plug, stratified, slug and wavy flow. The location of the ribs in the duct did not alter the shape of the flow regimes, but the regime boundaries were considerably changed (repositioned). The effects of using ribs of different heights on regime boundaries are represented with flow map diagrams and discussed in detail. Compared to the smooth duct, the ribbed duct had different regime boundary positions. Increasing the rib height initiated hydrodynamical instability at lower fluid velocities. These findings are relevant for the operation and design of pipes, boilers and heat exchangers.     

Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

Heat transfer measurements and correlations for air–water flow of different flow patterns in a horizontal pipe

Heat transfer coefficients were measured and new correlations were developed for two-phase, two-component (air and water) heat transfer in a horizontal pipe for different flow patterns. Flow patterns were observed in a transparent circular pipe using an air–water mixture. Visual identification of the flow patterns was supplemented with photographic data, and the results were plotted on the flow regime map proposed by Taitel and Dukler and agreed quite well with each other. A two-phase heat transfer experimental setup was built for this study and a total of 150 two-phase heat transfer data with different flow patterns were obtained under a uniform wall heat flux boundary condition. For these data, the superficial Reynolds number ranged from 640 to 35,500 for the liquid and from 540 to 21,200 for the gas. Our previously developed robust two-phase heat transfer correlation for a vertical pipe with modified constants predicted the horizontal pipe air–water heat transfer experimental data with very good accuracy. Overall the proposed correlations predicted the data with a mean deviation of 1.0% and an rms deviation of 12%.     

On two-phase flow patterns and transition criteria in aqueous methanol and CO2 mixtures in adiabatic,rectangular microchannels

This paper presents flow map investigations of adiabatic two-phase flow in square cross-sectioned, 200 μm deep microchannels fabricated in silicon, employing laser induced fluorescence microscopy. The influence of surface tension and nozzle geometry on the flow pattern transition was investigated using two nozzle widths (orifices of 30 μm and 50 μm, respectively) and methanol–water solutions with CO₂ as the gas phase. It was found and quantified that smaller nozzle geometries and smaller liquid surface tension promote the propagation of capillary gas bubbles at lower superficial gas and liquid velocities. Within the measurement domain of superficial gas (0.01–0.625 m/s) and liquid (0.0005–0.5000 m/s) velocities, we observed dispersed bubbly, regularly ordered bubbly, wedging, slug and annular flows, thus extending the experimental knowledge base to smaller superficial liquid velocities by almost two orders of magnitude. With the help of the flow maps presented herein, we were able to characterize the observed regularly ordered bubbly flow as the transition regime between dispersed bubbly and wedging flow. The results of the present investigation are of direct relevance to the operation of small-scale direct methanol fuel cells.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

Numerical study of flow and heat transfer during a high-speed micro-drop impact on thin liquid films

Numerical simulation of high-speed micro-droplet impingement on thin liquid film covering a heated solid surface has been carried out. Effect of droplet Weber number and liquid film thickness on the characteristics of flow and heat transfer has been investigated using the coupled level set and volume of fluid method. The code is validated against both the experimental and numerical results from the literature. Results show that the crown dynamics is mostly affected by variations in the initial film thickness but is weakly influenced by changes in the Weber number. The liquid within the film can be categorized as three regions based on the heat transfer distribution: the static film region, the transition region, and the impact region. The transient local wall temperature shows three stages: first stage when the temperature decreases rapidly, followed by a second stage in which the temperature starts to rise and then becomes almost constant in the third stage. After drop impact, the local Nusselt number continuously increases until reaching a maximum value, and then decreases approaching the initial impact stage. Our analysis of the change in Weber number shows that larger Weber number contributes to intense temperature variation at the crater core relative to other radial locations. Lastly, the results reveal that the thinner liquid film leads to lower wall temperature and hence, higher average Nusselt number.