Next Article in Journal
Nonlinear Methods Most Applied to Heart-Rate Time Series: A Review
Next Article in Special Issue
Lognormality in Turbulence Energy Spectra
Previous Article in Journal
Using the Maximal Entropy Modeling Approach to Analyze the Evolution of Sedentary Agricultural Societies in Northeast China
Previous Article in Special Issue
Effect of Incidence Angle on Entropy Generation in the Boundary Layers on the Blade Suction Surface in a Compressor Cascade
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation on the Thermodynamic Characteristics of a Liquid Film upon Spray Cooling Using an Air-Blast Atomization Nozzle

1
School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, China
2
Institute of Engineering Thermophysics, North China University of Water Resources and Electric Power, Zhengzhou 450045, Henan, China
3
Advanced Research Center of Thermal and New Energy Technologies, Xingtai Polytechnic College, Xingtai 054035, Hebei, China
4
Chengyi Academy of PKUHS, Peking Univeristy, Beijing 100080, China
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(3), 308; https://doi.org/10.3390/e22030308
Submission received: 3 February 2020 / Revised: 26 February 2020 / Accepted: 8 March 2020 / Published: 9 March 2020
(This article belongs to the Special Issue Entropy Production in Turbulent Flow II)

Abstract

:
This paper developed a three-dimensional model to simulate the process of atomization and liquid film formation during the air-blast spray cooling technological process. The model was solved using the discrete phase model method. Several factors including the thermodynamic characteristics of the liquid film as well as the spray quality with different spray mass flow rates under different spray heights were numerically investigated and discussed. The results show that the varied spray height has little effect on the Sauter Mean Diameter (d32) of the spray droplet, while the thermodynamic characteristics of liquid film including the liquid film height, the liquid film velocity, and the liquid film generation rate are sensitive to the change of the spray height. With the growth of spray mass flow rates, d32, the liquid film generation rate and liquid film height become larger, while the liquid film velocity with different spray mass flow rates has a similar velocity distribution, indicating that the spray mass flow rate has little effect on the liquid film velocity. The average d32 of droplet size shows a sharp drop when sprayed from the nozzle in a short period of time (<1.5 ms), then approaching smoothness, below a value of 40 μ m , the spray status tends to be stable.

1. Introduction

During the last several decades, the miniaturization, as well as the highly integrated, high- density device packaging, and inefficiency of electronics have evoked great interest in modern industry and a wide variety of aids from industrial and environmental fields have been applied [1,2], for instance, ground-based devices like the lateral diffused metal oxide semiconductor field effect transistor (LD-MOSFET) [3] and space devices like the NASA Goddard mars orbiter laser altimeter [4]. However, large amounts of waste heat generated at high heating power, cause a super high temperature so as to decrease the service lifetime and efficiency of the electronic devices. Therefore, a high-efficiency thermal management approach, whose heat dissipation requirement would be over an order of 102 W/cm2 [5], is increasingly needed to remove the waste heat due to the dramatic growth in both electrical power and local heat flux of the electronics devices.
Tremendous studies [5,6,7] regarding the heat transfer attainability and discovering that spray cooling has long enjoyed a reputation as a cooling method being one of the most effective methods for removing high heat-waste flux. It has been proved that its utility in dissipating high heat fluxes could be up to 1200 W/cm2 using water as the working fluid [7]. Besides, spray cooling has several other remarkable advantages including its latent heat of vaporization which absorbs a great amount of heat with a minimal temperature variation and thermal uniformity with low coolant mass flux in a given condition. The spray cooling phenomenon can be described as follows: working fluid is forced through the atomization nozzle where the liquid is broken into numerous tiny droplets due to the pressure. Subsequently, droplets impinge onto a heated surface, forming a thin liquid film upon the surface and taking away heat. The waste heat can be dissipated by the heat transfer mechanism of forced convection, film evaporation, nucleate boiling, and second nucleation [8,9].
For enhancing the cooling capacity of spray cooling, numerous researches regarding the complicated heat transfer mechanisms have been carried out through both numerical simulation investigations and experimental studies under different conditions. In the existing literature, it has been found that droplet-impacting and the thin liquid film are the critical factors affecting the performance of the spray cooling system when using the same coolant. In general, a better spray cooling performance can be acquired by decreasing the Sauter Mean Diameter (SMD) of the spray droplet size [10] and increasing the velocity of droplet-impacting [11]. It has been found that the twin- fluid atomization method can obtain a higher quality of spray, and the air-blast nozzle is one representative example [12,13,14]. On the other side, the cooling ability of the spray shows that the heat transfer characteristics are influenced significantly by low gravity, which may result in a lower cooling effectiveness than higher gravity and normal gravity due to the increasing thickness of the average liquid film [15].
Experimental research is the major method to investigate the heat transfer performance of spray cooling technology. However, the conclusions achieved by experimental conditions have been usually case- specific and could not obtain a mechanical arrangement. Spray cooling involves many complicated processes, such as liquid atomization, droplet interaction, and liquid film formation, thus it is very difficult to complete research of the entire mechanism and model the full process. Therefore, many researchers have attempted to build a model to investigate the mechanism of the spray cooling technology for achieving an ideal heat transfer capacity. Some researchers set up a three-dimensional model to investigate the heat transfer mechanism of spray cooling technology. In these studies, a multiphase flow model was used to simulate the surface temperature, the liquid film thickness, and the film velocity distribution [16,17,18,19,20]. The two-dimensional model also has been established for modelling a single droplet impacting on a flat surface with a constant surface temperature [2,21,22,23]. The volume of fluid method was used in these researches. The droplet diameter was 50–150 μ m and the droplet velocity was 1.0–10.0 m/s. The simulation results indicated that the decrease of the droplet diameter and the increase of the velocity can promote the performance of spray cooling technology. Some researchers studied the cooling performance of spray cooling in harmful environments [24,25], such as the over-load environment.
Looking through the previous researches, most of the investigations have focused on the average liquid film thickness, not on the thermodynamic characteristics of the liquid film. In the present study, the Commercial CFD software ANSYS (15.0) (ICEM CFD and FLUENT) was used to carry out the numerical investigation with the discrete phase model (DPM) method, simulating the process of atomization and liquid film formation and summarizing the SMD of the droplet size and the thermodynamic characteristics. This included the distribution of liquid film height (LFH), liquid film velocity (LFV), and liquid film generation rate (LFGR) when using the air-blast nozzle under different spray height (SH) and varied spray mass flow rates (SMFR).

2. Computational Model and Numerical Analysis

2.1. Physical Model and Its Finite Element Mesh Generation

A 3-D CFD model of a cylindrical region equivalent to the working spray chamber with a diameter of D (100 mm) and a height of H (50 mm) had been established by the commercial 3D modeling software SOLIDWORKS (2017), which schematic view is shown in Figure 1 and the detailed structure of the air-blast nozzle is also displayed. Since the FLUENT’s air-blast atomizer model does not include the internal gas flows, the inlets of air streams were created on the top the model including the mass-flow-inlet of central-flow, velocity-inlet of swirling-flow and pressure- inlet. The bottom is for the outlet of pressure and the other surface of the cylindrical region is for the numerical escape walls with no slip. The target surface on the top of the internal cylinder, has height and area of 8 mm and 1.131 cm2 (diameter, d = 1.2 cm) respectively.
Figure 1b demonstrates the mesh generation of the computational model. The structured hexahedron grid combined with the unstructured tetrahedron grid were applied to generate the finite element meshes of the model. Additionally, the advanced O-type grid was used to promote the quality of mesh. In order to improve the accuracy of calculation and simulation results, mesh densification was carried out for the jet center area and the target surface region. The solution domain was divided into 2.5 × 10 6 meshes and this grid refinement could be shown to lead to grid independence solutions.

2.2. Governing Equations and Boundary Conditions

In the present study, the Commercial CFD software ANSYS (15.0) (ICEM CFD and FLUENT) was used to carry out the finite element numerical investigation. A Discrete-Phase particle trajectory model based on the Euler–Lagrangian method was adopted, in which the gas-phase is regarded as a continuous phase, and the Euler method was used to solve the Navier–Stokes conservation equations for the continuous phase; the particle phase is treated as the discrete phase, and the particle trajectory is solved by the Lagrangian method. Except for the heat and mass transfer caused by particles, it is assumed that the movements of droplet particles along their respective orbits do not interfere with each other, and there is no turbulent diffusion of droplet particles.

2.2.1. Continuous Phase Model

The continuous phase is the air, and the mass equation, momentum equation, and energy equation can be described as Equations (1) to (3) respectively [2]:
ρ t + x ( ρ u x ) + y ( ρ u y ) + z ( ρ u z ) = 0
[ ( ρ u x ) t + ( ρ u x U ) = ( μ u x ) ρ x + F x ( ρ u y ) t + ( ρ u y U ) = ( μ u y ) ρ y + F y ( ρ u z ) t + ( ρ u z U ) = ( μ u z ) ρ z + F z ]
( ρ E ) t + [ U ( ρ E + p ) ] = [ λ e f f T j J j T r e f T c p , j d T + ( τ e f f U ) ] + S h
where E is the specific energy including kinetic energy, intrinsic energy, and potential energy; λ e f f is the effective thermal conductivity; J j is the component diffusion flux, and S h is the volume heat source due to the heat transfer between droplet and the continuous phase.
The standard k ε model is used for the turbulence model based on the turbulent flow energy and the diffusivity. The turbulent kinetic energy equation is an exact equation while the diffusivity equation is derived from empirical formulas. The turbulent kinetic energy k-equation and the diffusivity ε -equation are described as Equations (4) and (5) respectively [2]:
ρ D k D t = x i [ ( μ l + μ t σ k ) k x i ] + G k + G b ρ ε
ρ D ε D t = x i [ ( μ l + μ t σ ε ) ε x i ] + C 1 ε ε k ( G k + G 3 ε G b ) C 2 ε ρ ε 2 k
where μ l and μ t are the viscous coefficients of laminar and turbulence respectively; G k and G b are the turbulence energies due to the laminar velocity gradient and buoyant force respectively; C μ is the turbulent constant; C 1 ε , C 2 ε , C 3 ε , σ k and σ ε are all empirical constants whose values are 1.44, 1.92, 0.09, 1.0, and 1.3 respectively.

2.2.2. Discrete Phase Model

The air-blast nozzle, which existed in the model library of the FLUENT, was used to investigate the spray characteristic. The air-blast atomization process is to accelerate the breakup of liquid sheets from the atomizer by means of an additional air stream passing through the atomizer. ANSYS FLUENT’s air-blast atomization model assumes that the sheet breakup is always due to short waves. This assumption is a consequence of the greater sheet thickness commonly found in air-blast atomizers. Hence the ligament diameter is assumed to be linearly proportional to the wavelength of the fastest-growing wave on the sheet, and which can be calculated by Equation (6).
d L = 2 π C L K s
where C L is the ligament constant which is equal to 0.5 by default; where K s is the wave number corresponding to the maximum growth rate.
In the atomization field, KHRT (Kelvin–Helmholt–Rayleigh–Taylor) is used to solve the two- phase flow. The breakup process is the result of the interaction of the KH model and the RT model, and the formation liquid radius and breakup time can be calculated by Equations (7) and (8) respectively:
r c = π C R T K R T
τ RT = C τ Ω R T
Ω RT = 2 ( g t ( ρ p ρ g ) ) 3 2 / 3 3 σ ( ρ p + ρ g )
K R T = g t ( ρ p ρ g ) / 3 σ
where Ω RT is the fastest growing unstable wave frequency of the jet; K R T is the wave number; ρ p and ρ g are density of liquid droplet and gas-phase respectively; σ is the surface tension.
The particle motion equation in Lagrange coordinates can be obtained by Newton’s second law:
[ d u p , x d t = F D ( u x u p , x ) + g x ( ρ p ρ v ) ρ p + F x d u p , y d t = F D ( u y u p , y ) + g y ( ρ p ρ v ) ρ p + F y d u p , z d t = F D ( u z u p , z ) + g z ( ρ p ρ v ) ρ p + F z ]
where u is the velocity of the liquid-phase, m/s; u p is the velocity of the liquid droplet particle, m/s; F D is the drag force on the droplet of unit mass, m/s2.

2.2.3. Coupling Calculation of Continuous Phase and Discrete Phase

The numerical investigation of spray in ANSYS FLUENT belongs to a two-phase simulation of air and liquid coupling calculation. As a continuous phase, air has an effect on the orbital distribution as well as the heat and mass transfer of droplets in the discrete term; in turn, the discrete phase will also affect the continuous phase. The change values of momentum, mass and heat of droplets, which are achieved by the source term, can be calculated by Equations (12) to (14):
[ F x = ( 3 μ C D Re 4 ρ d d p 4 ( u v , x u p , x ) + g x ( ρ p ρ v ) ρ p D T , P 1 m p T T x ) m p Δ t F y = ( 3 μ C D Re 4 ρ d d p 4 ( u v , y u p , y ) + g y ( ρ p ρ v ) ρ p D T , P 1 m p T T y ) m p Δ t F z = ( 3 μ C D Re 4 ρ d d p 4 ( u v , z u p , z ) + g z ( ρ p ρ v ) ρ p D T , P 1 m p T T z ) m p Δ t ]
M = Δ m p m p , 0 m p , 0
Q = ( m p ¯ m p , 0 c p Δ T p + Δ m p m p , 0 ( h f g + T r e f T p c p , h 2 o d T ) ) m p , 0

2.3. Boundary Conditions and Coupling Calculation Procedure

Reasonable setting of the boundary conditions of the computational model according to the physical environment of the experiment can not only accelerate the convergence speed of the simulation process, but also improve the accuracy of the calculation. The present study set the corresponding initial state and boundary conditions according to the actual working conditions of the spray experiment, specifically as follows:
(1)
The pressure, temperature, and gravity acceleration of the initial states were 1 atm, 293 K, and 9.8 m/s2 respectively. The injection direction was Y-axis;
(2)
The velocity of mass-flow-inlet of central-flow, velocity-inlet of swirling-flow and the pressure-inlet were 9.167 × 10−5 kg/s, 20 m/s, and 1m/s respectively, whose values were all obtained from the experimental data. The rotating tangential vector of the swirling-flow was 0.7071 and the gas-phase was air;
(3)
The vertical wall surface was defined as an escape and no-slip boundary condition and the bottom side was the pressure outlet;
(4)
The target surface was no-slip condition and the formation of liquid film due to atomized droplets impinging on the surface was simulated by the FLUENT’s Eulerian-wall film model.
The present study assumed that the thickness of liquid film was less than 500 μ m and the simulation process was transient with an unsteady solution method. The coupled calculation of continuous phase and discrete phase could be carried out in the following step, which is shown in Figure 2 and described as follows: (1) Without adding the discrete phase, calculate the convergent or partially convergent continuous phase field; (2) Create a nozzle model for adding the discrete phase to simulate the orbital equation of the droplet particles and obtain the trajectory, velocity, mass and temperature of the particle; (3) Exchange the heat, momentum, and mass terms in each fluid unit and apply to the continuous phase calculation to resolve the continuous phase flow field; (4) Repeat the above steps and iterate repeatedly until the continuous phase flow field does not change with the calculation process.

2.4. Simulating Condition Arrangement

The target of the current investigation is to study the effects of SH and SMFR on the thermodynamic characteristics of the liquid film as well as the average droplets SMD, which all influence the heat transfer performance of the spray cooling technology using the air-blast nozzle. The related simulating operation conditions are listed in Table 1. Transparently, cases 1 to 4 were arranged to study the influence of various SHs, while the cases 2, 5–7 were arranged to research the effect of SMFRs on the thermodynamic characteristics of the liquid film and spray droplet size. It should be noted that the heat flux was 100 W/cm2 in the present study.

3. Results and Discussion

The thermodynamic characteristics of the air-blast atomization spray cooling technology were simulated with an explored 3-D model using the commercial software of ANSYS 15.0 (ICEM CFD and FLUENT), exploring the effects of SH (15 mm, 20 mm, 25 mm, and 30 mm) and SMFR (0.020 kg/s, 0.025 kg/s,0.030 kg/s, and 0.035 kg/s) on SMD of the droplet and the thermodynamic characteristics of the liquid film (LFH, LFV, and LFGR). Figure 3 shows the cloud picture of the spray process under an SH of 30 mm. It took 0.6 ms for the liquid droplets from the nozzle to reach the surface and the spray completely covered the surface in 0.9 ms.

3.1. Effects of the SH

The distance between the nozzle outlet and the target surface has a great influence on the thermodynamic property of the generated liquid film upon the target surface. The aim in the current section is to explore the effects of SH on the air-blast atomization spray cooling process, seeking efficient thermodynamic characteristics of the obtained droplets and thin liquid film.

3.1.1. SMD of the Spray Droplet Size

The SMD is a critical factor to measure the atomization quality and affect the performance of spray cooling. In general, a better spray cooling performance can be acquired by decreasing the SMD of the spray droplet size [10]. Figure 4 shows the variation of the SMD with time under different SHs of 15 mm, 20 mm, 25 mm, and 30 mm.
In Figure 4, it is obvious that the average SMD of the droplet size has a sharp drop being sprayed from the nozzle in a short period of time (<1.5 ms), then approaching smoothness where it is below a value of 40 μ m where the spray status tends to be stable. The three curves show the same trend, indicating that the SH has little effect on the SMD because of the same initial conditions with a constant spray flow rate of 0.025 kg/s and three invariable inlet conditions as previously mentioned.

3.1.2. The Hydraulic Characteristics of the Liquid Film

The main hydraulic characteristics of the liquid film are all discussed in detail as follows, including the liquid film height (LFH), the liquid film velocity (LFV), and the liquid film generation rate (LFGR) under different SH. Especially, the cloud pictures and the corresponding curves were all obtained under the condition in which the liquid film completely covered the target surface with a time of 3.5 ms.
1. The LFH and the LFV
The cloud pictures including the LFH and LFV are displayed in Figure 5a,b respectively. It is obvious that the variation of the LFH and the LFV present a circular distribution with the center below the nozzle. Near the central area, the atomization droplets directly impinged on the surface and were just below the jet center. With the rebound and splash of a great amount of droplets, the LFH in this area is low and the LFV has a large speed fluctuation. After the liquid film spreads rapidly along the radius, the LFH presents a trend of first increasing and then decreasing, while the LFV shows the opposite trend. This is due to the extrusion and flow of a large number of droplets around the central area, resulting in the growth of the LFH and the LFV. On continuing outward, the liquid film obtains an adequate diffusion space so that the LFV continues to increase, while the LFH gradually decreases due to the momentum conservation law. In addition, comparing with the LFH under the four SHs (15 mm, 20 mm, 25 mm, and 30 mm), it is found that the average LFH at the nozzle height of 15 mm is higher than that of the other three heights as well as the LFV being the largest. However, the distribution of liquid film under the highest SH condition of 30 mm is the most uniform on the whole target surface.
For the purpose of better studying the flow characteristics of the liquid film, several points along the radius of the X-axis on the circular surface were taken to quantitatively analyze the LFH and the LFV under different SH conditions. The specific hydraulic characteristic curves distribution is shown in Figure 6. The maximum LFHs of 28.1 μ m , 26.8 μ m , 23.8 μ m , and 22.4 μ m appear at X = 2.1 mm, X = 3 mm, X = 3.3 mm, and X = 3.9 mm respectively. The lower SH has a higher maximum LFH due to the smaller projection circle of the spray on the surface, resulting in a reduction in the radius of the central area directly affected by the droplets of the atomized main flow region. It can be seen from Figure 6. Figure 6b that at the area near the center, the LFV fluctuates greatly under different SH conditions especially for a nozzle height of 30 mm. After leaving the central area, the LFV shows an increasing trend, and the lower the SH, the larger the LFV. The reason being the LFV is affected by the droplet impact, splash, and collision in the central area. With the growing amount of spray, the liquid film flows rapidly around and the LFV is gradually increased, reaching a maximum value at the edge of the target surface.
2. The LFGR
When the atomization droplets impinge on the target surface, a liquid film is formed, whose quality is also a critical factor influencing the ability of the spray cooling technology. The film generation rate ( ) is used for comparing the ratio of the atomization droplets impinging and forming liquid film on the surface ( m l f ) to total amount of the working fluid ( m ) under different conditions. The calculation can be expressed by Equation (15) [18].
= m l f m
Figure 7 shows the variation curves of the LFGR with time under different SHs. It is obvious that all four curves show an upward trend first then gradually flatten, indicating that the formation of liquid film gradually completes to an invariable shape meaning that the LFH and the LFV are fixed in a range. The LFGR under an SH of 15 mm is higher than that of the other three SHs, because the decrease of SH narrows down the atomization range and less droplets collide and splash during the movement, resulting in the supplementation of the quality of liquid film formation as well as an increase of the LFGR.

3.2. Effects of the SMFR

The SMFR of the input coolant that are charged into the air-blast nozzle will have a major effect on the property of the thermodynamic characteristics. In this section, four different SMFRs of 0.020 kg/s, 0.025 kg/s, 0.030 kg/s, and 0.035 kg/s are changed to investigate the influence on the thermodynamic properties of the generated droplets and the characteristics of the liquid film on the target surface.

3.2.1. SMD of the Spray Droplet Size

Figure 8 shows the variation of the average droplets SMD with time under different SMFRs. It can be observed obviously that all four curves show similar trends, the average SMD first has a sharp drop with a short period time of 1.5 ms, then approaching smoothness below the value of 40 μ m . With the growth of SMFRs from 0.020 kg/s to 0.035 kg/s, the relatively positions of the four curves increase. The green curve representing the maximum SMFR of 0.035 kg/s stays at the highest position, indicating that the average droplets SMD generated by the maximum SMFR have the largest value.
Generally speaking, liquid coolant comes into the nozzle, where it is mixed with high-pressure air then atomized into thin liquid droplets. The SMD of the generated droplets depends on the state of the nozzle inlet including SMFR and the high-pressure air. Therefore, the average droplet SMD increases with increasing SMFR under an invariable air-inlet-pressure.

3.2.2. The Hydraulic Characteristics of the Liquid Film

1. The LFH and the LFV
Figure 9 shows the thermodynamic characteristic curves of the liquid film. Several points along the radius of the X-axis on the circular surface were taken to quantitatively analyze the LFH and LFV under different SMFRs.
Figure 9a shows the characteristic curves of the LFH. All four curves representing four different SMFRs display the similar tendency. It can be clearly observed that all LFH generated by the four different SMFRs are all under the value of 35 μ m . Within a radius of about 3 mm, the SMFR has little effect on the LFH; while the liquid film gradually expands outward along the radius (>3 mm), the average LFH will increase with the growth of SMFR. Additionally, the peak value of LFH also increases with increasing SMFR.
Figure 9b displays the characteristic curves of the LFV. The tendency of the curves is similar to that generated by the different SHs in Section 3.1.2. However, the four curves representing the four SMFRs also have similar values, meaning that the LFV with different SMFRs have similar velocity distribution, indicating that the SMFR has little effect on the LFV.
2. The LFGR
Figure 10 illustrates the curves of the LFGR. All four curves display the same tendency as time goes on under the different SMFRs, where all the four curves show an upward trend first then gradually flatten. This indicates that the formation of liquid film gradually completes to an invariable shape meaning that the LFH and the LFV are fixed in a range. The green curve stays at the highest position, indicating that the LFGR generated by the maximum SMFR of 0.035 kg/s has the highest ability. A greater SMFR means more coolant spray onto the target surface, thus the LFGR will increase with increasing SMFR. In addition, the LFGRs obtained by the four different SMFR are all lower than 0.5.

4. Conclusions

In the present paper, a three-dimensional model was developed to investigate both the thermodynamic characteristics of a liquid film with an air-blast nozzle for spray cooling and the spray quality. The model focuses on simulating the factors including the SMD of spray droplet size, the LFH, the LFV, and the LFGR, which all influence the heat transfer performance of spray cooling. The variation of factors under different SHs and SMFRs are discussed in detail. The primary conclusions can be formulated as follows:
(1)
The varied SH has little effect on the SMD with an invariable inlet condition; while the different SMFRs can influence the SMD of the spray droplet. With increasing SMFR, the SMD becomes larger. The average SMD of droplet size shows a sharp drop when sprayed from the nozzle in a short period of time (<1.5 ms), then approaching smoothness below a value of 40 μ m , the spray status tends to be stable;
(2)
Under all working conditions, the variations of LFH and LFV have a similar tendency. In the central area below the nozzle on the surface, the LFH is relatively low as well as the LFV while the LFV has a large speed fluctuation. Along the radius direction, the LFH increases first and then decreases, while the LFV increases gradually. When it approaches the circumferential edge, the velocity reaches its maximum;
(3)
A larger LFGR can be obtained by a lower SH. With the increase of SH, the average LFH decreases as well as the peak height but is higher at the edge-area. The LFV expands from the central region to the surrounding areas. A lower SH can be obtained a a higher average LFV and a lower speed fluctuation in the central area.
(4)
With the growth of SMFR, the LFGR and the average LFH become larger. In addition, the peak value of LFH also increases with increasing SMFR. The LFVs with different SMFRs have similar velocity distribution, indicating that the SMFR has little effect on the LFV.
Overall, the SH and the SMFR have major effects on the thermodynamic characteristics of the liquid film as well as the spray quality, which would have a large influence on the heat transfer performance of air-blast spray cooling technology. From the above results and conclusions, we find a promising approach to achieve a more efficient heat dissipation capacity with spray cooling technology.

Author Contributions

Conceptualization, Y.-Z.L.; methodology, J.-X.L.; formal analysis, J.-X.L.; writing—original draft preparation, J.-X.L. and T.L.; writing—review and editing, E.-H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The work uncovered in the present paper was conducted at the Spacecraft Thermal Control Laboratory, Beihang University.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

C constant τ RT breakup time
DDiameter (m) i ε Empirical (i = 1,2 and 3)
d l ligament diameterSubscripts
d p Particle diameter b buoyant
E specific energy including kinetic energy, intrinsic energy and potential energy D drag
F D drag force on droplet of unit mass (m/s2) g gas-phase
G b turbulence energy due to the buoyant force L ligament
G k turbulence energy due to the laminar velocity gradient l laminar
gGravitational acceleration (m/s2) p liquid droplet particle
J j component diffusion fluxRTRayleigh–Taylor model
K s wave number corresponding to the maximum growth rate t turbulence
K R T wave numberxx-axis
QHeat (W)yy-axis
ReReynold numberzz-axis
r c liquid radiusAcronyms
S h volume heat source due to the heat transfer between droplet and the continuous phaseCFDComputational Fluid Dynamics
TTemperature (K)DPMDiscrete Phase Model
u velocity of liquid-phase (m/s)LFHLiquid Film Height
Greek symbolsLFVLiquid Film Velocity
λ Thermal conductivity ( W m 1 K 1 )LFGRLiquid Film Generation Rate
ζ Uncertainty (%)LD-MOSFETLateral Diffused Metal Oxide Semiconductor Field Effect Transistor
ρ Density ( k g m 3 )SHSpray Height
σ Surface tension ( N m 1 )SMDSauter Mean Diameter
Ω unstable wave frequencySMFRSpray Mass Flow Rate
μ Dynamic viscosity ( P a s )3-DThree dimensional

References

  1. Zhang, H.; Li, Y.; Mao, Y.; Tian, S.; Liu, Y.; Yang, L. Optimization design of a novel closed-loop spray cooling system. In 2015 IEEE International Conference on Advanced Intelligent Mechatronics (AIM); IEEE: Piscataway, NJ, USA, 2015; pp. 324–329. [Google Scholar]
  2. Liping, P.; Kun, L.; Qi, G.; Shuxin, L.; Chao, Y. Numerical study of high-overload effect on liquid film of spray cooling. Appl. Therm. Eng. 2017, 127, 1015–1024. [Google Scholar] [CrossRef]
  3. Cotler, A.C.; Brown, E.R.; Dhir, V.; Shaw, M.C. Chip-level spray cooling of an LD-MOSFET RF power amplifier. IEEE Trans. Conmpon. Pack. Technol. 2004, 27, 411–416. [Google Scholar] [CrossRef]
  4. Silk, E.A.; Golliher, E.L.; Selvam, R.P. Spray cooling heat transfer: Technology overview and assessment of future challenges for micro-gravity application. Energy Conv. Manag. 2008, 49, 453–468. [Google Scholar] [CrossRef]
  5. Wang, J.X.; Li, Y.Z.; Zhang, Y.; Li, J.X.; Mao, Y.F.; Ning, X.W. A hybrid cooling system combining self-adaptive single-phase mechanically pumped fluid loop and gravity-immune two-phase spray module. Energy Convers. Manag. 2018, 176, 194–208. [Google Scholar] [CrossRef]
  6. Murshed, S.M.S.; De Castro, C.A.N. A critical review of traditional and emerging techniques and fluids for electronics cooling. Renew. Sustain. Energy Rev. 2017, 78, 821–833. [Google Scholar] [CrossRef]
  7. Pais, M.R.; Chow, L.C.; Mahefkey, E.T. Surface roughness and its effects on the heat transfer mechanism in spray cooling. ASME J. Heat Transf. 1992, 114, 211–219. [Google Scholar] [CrossRef]
  8. Bancoff, S.G. Stability and dynamics of thin heated liquid films. Int. J. Fluid Mech. Res. 1998, 25, 189–201. [Google Scholar] [CrossRef]
  9. Yang, B.H.; Wang, H.; Zhu, X.; Liao, Q.; Ding, Y.D.; Chen, R. Heat transfer enhancement of spray cooling with ammonia by microcavity surfaces. Appl. Therm. Eng. 2013, 50, 245–250. [Google Scholar] [CrossRef]
  10. Hsieh, S.S.; Fan, T.C.; Tsai, H.H. Spray cooling characteristics of water and R-134a, part II: Nucleate boiling. Int. J. Heat Mass Transf. 2004, V47, 5703–5712. [Google Scholar] [CrossRef]
  11. Zhou, N.; Chen, F.; Cao, Y.; Chen, M.; Wang, Y. Experimental investigation on the performance of a water spray cooling system. Appl. Therm. Eng. 2017, 112, 1117–1128. [Google Scholar] [CrossRef]
  12. Jedelsky, J.; Jicha, M. Energy conversion in effervescent atomization. In Proceedings of the 12th ICLASS, Heidelberg, Germany, 2–6 September 2012; pp. 2–9. [Google Scholar]
  13. Li, J.X.; Li, Y.Z.; Cai, B.Y.; Li, E.H. Experimental investigation of heat transfer mechanism of air-blast spray cooling system with a two-phase ejector loop for aeronautical application. Energies 2019, 12, 3963. [Google Scholar] [CrossRef] [Green Version]
  14. Beck, J.E.; Lefebvre, A.H.; Koblish, T.R. Liquid sheet disintegration by impinging air streams. Atom. Sprays 1991, 1, 155–170. [Google Scholar] [CrossRef]
  15. Gambaryan-Roisman, T.; Kyriopoulos, O.; Roisman, I.; Stephan, P.; Tropea, C. Gravity effect on spray impact and spray cooling. Microgravity Sci. Technol. 2007, 19, 3–4. [Google Scholar] [CrossRef]
  16. Sarkar, S.; Selvam, R.P. Direct numerical simulation of heat transfer in spray cooling through 3D multiphase flow modeling using parallel computing. J. Heat Transf. 2009, 131, 121007-1–121007-8. [Google Scholar] [CrossRef]
  17. Bussmann, M.; Mostaghimi, J.; Chandra, S. On a three-dimensional volume tracking model of droplet impact. Phys. Fluids 1999, 11, 1406–1417. [Google Scholar] [CrossRef]
  18. Issa, R.J.; Yao, S.C. Numerical model for spray-wall impaction and heat transfer at atmospheric conditions. J. Thermophys. Heat Transf. 2005, 19, 441–1117. [Google Scholar]
  19. Selvam, R.P.; Lin, L.; Ponnappan, R. Direct simulation of spray cooling: Effect of vapor bubble growth and liquid droplet impact on heat transfer. Int. J. Heat Mass Transf. 2006, 49, 4265–4278. [Google Scholar] [CrossRef]
  20. Selvam, R.P.; Hamilton, M.; Silk, E.A. Silk Spray cooling modeling: Liquid film thickness effect on heat transfer. AIP Conf. Proc. 2007, 880, 110–117. [Google Scholar]
  21. Gunjal, P.R.; Ranade, V.V.; Chaudhari, R.V. Dynamics of drop impact on solid surface: Experiments and VOF simulations. AICHE J. 2005, 51, 59–78. [Google Scholar] [CrossRef]
  22. Ganesan, S.; Rajasekaran, S.; Tobiska, L. Numerical modeling of the non- isothermal liquid droplet impact on a hot solid substrate. Int. J. Heat Mass Transf. 2014, 78, 670–687. [Google Scholar] [CrossRef]
  23. Šikalo, Š.; Wilhelm, H.D.; Roisman, I.V.; Jakirlić, S.; Tropea, C. Dynamic contact angle of spreading droplets: Experiments and simulations. Phys. Fluids 2005, 17, 62103. [Google Scholar]
  24. Sinha-Ray, S.; Yarin, A.L. Drop impact cooling enhancement on nano-textured surfaces. Part I: Theory and results of the ground (1 g) experiments. Int. J. Heat Mass Transf. 2014, 70, 1095–1106. [Google Scholar] [CrossRef]
  25. Wang, Z.; Xing, Y.; Liu, X.; Zhao, L.; Ji, Y. Computer modeling of droplets impact on heat transfer during spray cooling under vibration environment. Appl. Therm. Eng. 2016, 107, 453–462. [Google Scholar] [CrossRef]
Figure 1. Schematic view, (a) 3-D physical model and (b) mesh.
Figure 1. Schematic view, (a) 3-D physical model and (b) mesh.
Entropy 22 00308 g001
Figure 2. Flowchart of the calculation procedure.
Figure 2. Flowchart of the calculation procedure.
Entropy 22 00308 g002
Figure 3. Cloud pictures of spray process.
Figure 3. Cloud pictures of spray process.
Entropy 22 00308 g003
Figure 4. Variation of the SMD with different SH.
Figure 4. Variation of the SMD with different SH.
Entropy 22 00308 g004
Figure 5. Cloud pictures of the liquid film: (a) the LFH and (b) the LFV.
Figure 5. Cloud pictures of the liquid film: (a) the LFH and (b) the LFV.
Entropy 22 00308 g005
Figure 6. The thermodynamic characteristic curves of the liquid film: (a) the LFH and (b) the LFV.
Figure 6. The thermodynamic characteristic curves of the liquid film: (a) the LFH and (b) the LFV.
Entropy 22 00308 g006
Figure 7. The film generation rate.
Figure 7. The film generation rate.
Entropy 22 00308 g007
Figure 8. Variation of the average droplets SMD with different SMFRs.
Figure 8. Variation of the average droplets SMD with different SMFRs.
Entropy 22 00308 g008
Figure 9. The thermodynamic characteristic curves of the liquid film: (a) the LFH and (b) the LFV.
Figure 9. The thermodynamic characteristic curves of the liquid film: (a) the LFH and (b) the LFV.
Entropy 22 00308 g009
Figure 10. The film generation rate.
Figure 10. The film generation rate.
Entropy 22 00308 g010
Table 1. Description of simulating conditions.
Table 1. Description of simulating conditions.
CaseVarious Parameter Descriptions
SH (mm)SMFR (kg/s)
1150.025
2200.025
3250.025
4300.025
5200.020
6200.030
7200.035

Share and Cite

MDPI and ACS Style

Li, J.-X.; Li, Y.-Z.; Li, E.-H.; Li, T. Numerical Investigation on the Thermodynamic Characteristics of a Liquid Film upon Spray Cooling Using an Air-Blast Atomization Nozzle. Entropy 2020, 22, 308. https://doi.org/10.3390/e22030308

AMA Style

Li J-X, Li Y-Z, Li E-H, Li T. Numerical Investigation on the Thermodynamic Characteristics of a Liquid Film upon Spray Cooling Using an Air-Blast Atomization Nozzle. Entropy. 2020; 22(3):308. https://doi.org/10.3390/e22030308

Chicago/Turabian Style

Li, Jia-Xin, Yun-Ze Li, En-Hui Li, and Tong Li. 2020. "Numerical Investigation on the Thermodynamic Characteristics of a Liquid Film upon Spray Cooling Using an Air-Blast Atomization Nozzle" Entropy 22, no. 3: 308. https://doi.org/10.3390/e22030308

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop