AN INTERFACIAL TRACKING MODEL FOR CONVECTION-CONTROLLED MELTING PROBLEMS
|
|
- Charla Mosley
- 5 years ago
- Views:
Transcription
1 Numerical Heat Transfer, Part B, 59: , 2011 Copyright # Taylor & Francis Group, LLC ISSN: print= online DOI: / AN INTERFACIAL TRACKING MODEL FOR CONVECTION-CONTROLLED MELTING PROBLEMS Qicheng Chen 1, Yuwen Zhang 2, and Mo Yang 1 1 College of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai, People s Republic of China 2 Department of Mechanical and Aerospace Engineering, University of Missouri, Columbia, Missouri, USA An interfacial tracking model is used to simulate melting in an enclosure in the presence of natural convection. It obtains the melting-front location by calculating the energy balance at the solid liquid interface and is a simple and convenient method to solve the solid liquid phase-change problem. It combines the advantages of the both deforming and fixed grids method and can handle natural convection controlled melting and solidification problems. It is demonstrated through comparisons with various experiments and numerical results that the interfacial tracking model can be used to solve melting problems controlled by natural convection at different Rayleigh numbers ranging from 10 5 to INTRODUCTION Modeling and numerical simulation for solid liquid phase-change problems has become an active area in the last several decades due to its wide application in energy systems [1] as well as thermal manufacturing, such as laser drilling [2], laser welding [3], and selective laser sintering [4], building systems [5], such as encapsulated materials used in plasterboard or packed beds [6], and in thermal energy storage systems [7 9]. Numerical simulation of convection-dominated phase-change problems is a challenging task. These problems are not only strongly nonlinear, i.e., the fluid flow in the melt dictates the heat transfer, they also involve moving boundaries. A variety of conventional numerical techniques have been developed for solving these problems: enthalpy [10, 11], apparent heat capacity [12], isotherm migration [13], and coordinate transformation methods [14 19]. These methods have been introduced by researchers to overcome the difficulties in handling moving boundaries. Previous works on multidimensional moving boundary problems include those of Received 28 July 2010; accepted 24 November Support for this work by the U.S. National Science Foundation under Grant CBET and the Chinese National Natural Science Foundation under Grants and is gratefully acknowledged. Address correspondence to Yuwen Zhang, Department of Mechanical and Aerospace Engineering, University of Missouri, E3411 Lafferre Hall, Columbia, MO 65211, USA. zhangyu@missouri.edu 209
2 210 Q. CHEN ET AL. NOMENCLATURE a nb coefficient in Eq. (22) a p coefficient in Eq. (22) A ratio of height and width (¼W=H) b term in Eq. (22) C p heat capacity, J=m 3 K d e coefficient in Eq. (23) d n coefficient in Eq. (24) f p liquid fraction Fo Fourier number (¼n t=h 2 ) g gravitational acceleration, m=s 2 Gr Grashof number (¼gbDTH 3 =n 2 ) h s latent heat of fusion, J=kg H height of wall, m k thermal conductivity in the liquid, W=mK K dimensionless thermal conductivity (¼k=k ) K s ratio of thermal conductivity (¼k s =k ) bk modified dimensionless thermal conductivity Nu Nusselt number at the heated wall (¼hH=k ) p pressure, Pa P dimensionless pressure [¼ðp þ q 1 gyþh 2 =qn 2 l ] Pr Prandtl number of the liquid PCM (¼n =a ) Ra Rayleigh number (¼Gr Pr) s location of solid liquid interface, m S dimensionless location of solid liquid interface (¼s=H) S 0 dimensionless location of solid liquid interface at the last time step Ste Stefan number [¼c (T h T m )=h s ] t time, s T temperature, K T c temperature of the cold wall, K T h temperature of the heated wall, K T i initial temperature, K melting point, K T m u velocity component in the x direction, m=s u I solid liquid interfacial velocity (¼qs=qt), m=s U dimensionless velocity in the x direction (¼uH=n ) U I dimensionless solid liquid interfacial velocity (¼u I H=n ) v velocity component in the y direction, m=s V dimensionless velocity in the y direction (¼vH=n ) Vol volume of the fusion Vol 0 volume of the enclosure W width of the enclosure, m x dimensional coordinate, m X dimensionless coordinate (¼x=H) y dimensional coordinate, m Y dimensionless coordinate (¼y=H) a thermal diffusivity, m 2 =s h dimensionless temperature [¼(T T m )= (T h T m )] m viscosity of the liquid PCM, kg=ms n kinematic viscosity, m 2 =s q density, kg=m 3 s dimensionless time (¼n t=h 2 ) Subscripts e east face of the control volume E east i initial I interface liquid m melting point n north face of the control volume N north s solid S south w west face of the control volume W west Duda et al. [20], Saitoh [21], Gong and Mujumdar [22], Cao et al. [23], Khillarkar et al. [24], Chatterjee and Prasad [25], and Beckett et al. [26]. These numerical models can be divided into two groups [27]: deforming-grid schemes (or strong numerical solutions) and fixed-grid schemes (or weak numerical solutions). Deforming-grid schemes transform the solid and liquid phases into fixed regions by using a coordinate transformation technique. The governing equations and boundary conditions are complicated due to the transformation. These schemes have been used successfully to solve multidimensional problems with or without natural convection. Ho and Viskanta [28] solved heat transfer during melting in a
3 INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS 211 rectangular cavity heated from the side isothermally using a deforming-grid scheme. Their results clearly revealed that the rates of melting and heat transfer were greatly affected by the buoyancy-driven convection in the liquid domain. Gadgil and Gobin [29] used a deforming-grid scheme to analyze the two-dimensional melting in a rectangular enclosure in the presence of natural convection. The disadvantage of deforming grid schemes is that they require significant amounts of computational time. On the other hand, the fixed-grid schemes use one set of governing equations for the whole computational domain, including both liquid and solid phases, and the solid liquid interface is later determined from the temperature distribution. This simplicity makes the computation much faster than deforming-grid schemes, while it still provides reasonably accurate results [30]. There are two major methods in the fixed-grid schemes: the enthalpy method [31] and the equivalent heat capacity method [32, 33]. The enthalpy method can handle solid liquid phase changes at a fixed melting point as well as in a range of temperature. When it is used to model phase change at a fixed temperature, it has difficulty with temperature and interfacial location oscillation. The equivalent heat capacity method must assume that phase change occurs in a range of temperature. In order to simulate phase change at fixed melting point, the temperature range must be very small. However, the equivalent heat capacity method encounters convergence problems when the range of phase-change temperature is small. Cao and Faghri combined the advantages of the enthalpy and equivalent heat capacity methods and proposed a temperature-transforming model [34]. Similar to the equivalent heat capacity method, it is also assumed that solid liquid phase change occurs in a range of temperatures. However, it works for phase changes occurring in any ranges of phase-change temperature. Therefore, the model can properly handle phase change occurring over a temperature range or at a single temperature. For conduction-controlled solid liquid phase-change problems, the converged solution can be obtained with any grid size and time step. For convection-controlled phase-change problems, Ma and Zhang [35] suggested that in order to use the temperature-transforming model efficiently and obtain convergent and reasonable results, the grid size must be selected with a suitable time step (which should not be too small). In this article, a fixed-grid interfacial tracking method [36] is used to solve the solid liquid phase-change problem in a rectangular enclosure. The solid liquid phase change is a moving-boundary problem due to the existence of the moving solid liquid interface. The interfacial tracking method also combines the advantages of the deforming-grid and fixed-grid methods. The location of the solid liquid interface is obtained by energy balance at the solid liquid interface at every time step. The results obtained using the interfacial tracking method will be compared with existing experimental and numerical results. PHYSICAL MODEL Figure 1 shows the physical model of the problem: melting inside a rectangular enclosure. The top and bottom walls are insulated, whereas the left and right walls are kept at constant temperatures of T h and T c, respectively. The initial temperature is equal to T i, which is below the melting point, T m.
4 212 Q. CHEN ET AL. The following assumptions are made: 1. The phase-change material (PCM) is pure and homogeneous. 2. The liquid phase of the PCM is a Newtonian and incompressible fluid. 3. The thermophysical properties of the solid and liquid PCM are independent of temperature, except for the temperature-dependent liquid density in buoyancy term (Boussinesq approximation). 4. The volume change due to solid liquid phase change is negligible. 5. The natural convection in the liquid phase is laminar and two-dimensional. The dimensional governing equations are as follows. Continuity equation: Momentum equations: qðquþ qt qðqvþ qt Figure 1. Physical model. qu qx þ qv qy ¼ 0 þ qðqu2 Þ qx þ qðquvþ ¼ qp qy qx þ q qx þ qðquvþ qx þ qðqvvþ ¼ qg qp qy qy þ q qx m qu þ q qx qy m qu qy m qv þ q qx qy m qv qy ð1þ ð2þ ð3þ Energy equation: qc p qt qt þ u qt qx þ v qt qy ¼ q qx k qt qx þ q qy k qt qy ð4þ
5 The boundary conditions of Eqs. (1) (4) are as follows. Left vertical wall: Right vertical wall: Bottom horizontal wall: INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS 213 x ¼ 0 T ¼ T h u ¼ 0 v ¼ 0 ð5þ x ¼ W T ¼ T c u ¼ 0 v ¼ 0 ð6þ y ¼ 0 qt qy ¼ 0 u ¼ 0 v ¼ 0 ð7þ Top horizontal wall: Melting front: y ¼ H " x ¼ s 1 þ qs # 2 qy qt qy ¼ 0 u ¼ 0 v ¼ 0 ð8þ x ¼ s T ¼ T m ð9þ k s qt s qx k qt qx ¼ q h s qs qt Equation (10) is the energy balance at the solid liquid interface, and the term 1 þðqs=qyþ 2 is to account for the effect of inclination of the solid liquid interface. Introducing the following nondimensional variables, X ¼ x H s ¼ t n H 2 Y ¼ y H h ¼ T T m T h T m S ¼ s H K s ¼ k s k ð10þ U ¼ u H n V ¼ v H n P ¼ ðp þ q 1gyÞH 2 Ste ¼ c ðt h T m Þ h s qn 2 Sc ¼ T c T m T h T c ð11þ the governing equations can be nondimensionalized as qu qx þ qv qy ¼ 0 qu qs þ qu 2 qx þ qðuvþ qy ¼ q2 U qx 2 þ q2 U qy 2 qv qs þ quv qx þ qv 2 qy ¼ q2 V qx 2 þ q2 V qy 2 þ Gr h ð12þ ð13þ ð14þ
6 214 Q. CHEN ET AL. qh qs þ U qh qx þ V qh qy ¼ 1 Pr! q 2 h qx 2 þ q2 h qy 2 ð15þ where Gr is the Grashof number based on height, Gr ¼ gbdth 3 =n 2 The boundary conditions of Eqs. (12) (15) are X ¼ 0 h ¼ 1 U ¼ 0 V ¼ 0 ð16þ X ¼ A h ¼ S c U ¼ 0 V ¼ 0 ð17þ Y ¼ 0 Y ¼ 1 Ste Pr qh qy ¼ 0 U ¼ 0 V ¼ 0 ð18þ qh qy ¼ 0 U ¼ 0 V ¼ 0 ð19þ X ¼ S h ¼ 0 ð20þ " 1 þ qs # 2 qy qh qx K s qh s ¼ U I K qx NUMERICAL SOLUTION AND INTERFACIAL TRACKING Discretization of Governing Equations The above two-dimensional governing equations are discretized by applying a finite-volume method [37], in which conservation laws are applied over finite-sized control volumes around grid points and the governing equations are then integrated over the control volume. A staggered grid arrangement is used in the discretization of the computational domain in the momentum equations. A power-law scheme is used to discretize convection=diffusion terms in the momentum and energy equations. The main algebraic equation resulting from this control-volume approach is in the form of ð21þ a P / P ¼ X a nb / nb þ b ð22þ where / P represents the value of general variable / (U, V,orh) at the grid point P, / nb are the values of the variable at P s neighbor grid points, and a P, a nb,andb are corresponding coefficients and terms derived from the original governing equations. The numerical simulation is accomplished using the SIMPLE algorithm [37]. The velocity-correction equations for corrected U and V in the algorithm are U e ¼ U e þ d eðp 0 P P0 E Þ V n ¼ V n þ d nðp 0 P P0 N Þ ð23þ ð24þ
7 INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS 215 where e and n represent the control-volume faces between grid P and its east neighbor E and grid P and its north neighbor N, respectively. In the present work, the governing equations are used for the entire computational domain. The velocity in the solid region is set to zero by letting a p ¼ and b ¼ 0 in Eq. (22) for the momentum equation. Interfacial Tracking Method Wang and Matthys [38] proposed an effective interface tracking method by introducing an additional node at the interface, which divides the control volume containing the interface into two small control volumes. In this work, an interfacial tracking method [36] that was developed for conduction-controlled melting of metal film under irradiation by a femtosecond laser will be extended to be able to handle convection-controlled solid liquid phase-change problems. This method is an alternative approach that does not require dividing the control volume containing the interface but can still accurately account for the energy balance at the interface. For a control volume that contains a solid liquid interface, the dimensionless temperature h P is numerically set as the interfacial temperature (h I ¼ 0) by letting a P ¼ and b ¼ 0 in Eq. (22) with / ¼ h. The above treatment yields an accurate result when the solid liquid interface is exactly at grid point P, as shown in Figure 2a. When the interfacial location within the control volume is not at grid point P, there are two scenarios as shown in Figure 2b: The interface is on the right side of the grid point; or (c) the interface is on the left side of the grid point. With scenario (b), a modified dimensionless thermal conductivity, bk w, at the face of the control volume w is introduced by equating the actual heat flux across the face of the control volume w, based on the position and temperature of the main grid point P [39], Considering h P ¼ h I, Eq. (25) becomes K w ðh W h I Þ K ¼ b w ðh W h P Þ þðf P 0:5ÞðDXÞ P ð25þ bk w ¼ þðf P 0:5ÞðDXÞ P K w ð26þ Similarly, a modified thermal conductivity at face e of the control volume can be obtained as bk e ¼ ðdxþ e ðdxþ e ðf P 0:5ÞðDXÞ P K e ð27þ With the scenario shown in Figure 2c, a modified dimensionless thermal conductivity, bk w, at the face of control volume w is obtained by a similar procedure: bk w ¼ ð0:5 f P ÞðDXÞ P K w ð28þ
8 216 Q. CHEN ET AL. Figure 2. Grid system near the liquid solid interface. A modified thermal conductivity at face e of the control volume is bk e ¼ ðdxþ e ðdxþ e þð0:5 f P ÞðDXÞ P K e ð29þ Combining the situations of Figures 2b and 2c, the universal modified dimensionless thermal conductivity can be expressed as bk w ¼ ð0:5 f P ÞðDXÞ P K w ð30þ
9 INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS 217 ðdxþ bk e ¼ e K e ðdxþ e þð0:5 f P ÞðDXÞ P ð31þ The modified thermal conductivities defined by Eqs. (30) and (31) are used to obtain the coefficients for grid points W and E, which allows the temperature at the main grid P to be used in the computation regardless of the location of the interface within the control volume. To determine the interfacial location, the energy balance at the interface, Eq. (21), can be discretized, and the solid liquid interfacial velocity can be obtained as ( U I ¼ Ste 1 þ S S ) 2 S Pr ðdyþ s K w ðh W h I Þ K e ðh I h E Þ ð0:5 f P ÞðDXÞ P ðdxþ e þð0:5 f P ÞðDXÞ P ð32þ where S S is the interfacial location at the grid at the south of grid P. The interfacial location is then determined using S ¼ S 0 þ U I Ds and the liquid fraction in the control volume that contains the interface is f ¼ S X P ðdxþ P =2 ðdxþ P Numerical Solution Procedure The numerical solution starts from time s ¼ 0. Once the temperature at the first control volume from the heated surface obtained exceeds the melting point, the temperature of the first control volume sets at the melting point by letting a P ¼ and b ¼ 0inEq.(22) with / ¼ h. After melting is initiated, the following iterative procedure is employed to solve for the interfacial velocity and the interfacial location at each time step. ð33þ ð34þ 1. Assume an interface velocity U I, using the velocity U I of the last time step as initial value. 2. Determine the new interface location, S, from Eq. (33). 3. Obtain the modified dimensionless thermal conductivities, bk w and bk e, at the faces of control volumes w and e from Eqs. (30) and (31). 4. Solve Eqs. (12) (15) to obtain the temperature distributions and obtain the new interface velocity U I by using Eq. (32). 5. Compare the newly obtained U I and the assumed value in step 1. If the difference is less than 10 5 and the maximum difference between the temperatures obtained from two consecutive iterative steps is less than 10 5, the interfacial location for the current step is obtained. If not, the process is repeated until the convergence criterion is satisfied.
10 218 Q. CHEN ET AL. RESULTS AND DISCUSSION The interfacial tracking method will be validated by comparing its results with experimental results as well as other numerical results. The first numerical simulation is performed under conditions that are the same as the experiments carried out by Okada [40]. The subcooling parameter is 0.01, which means that the initial temperature is very close to the melting point. The Rayleigh number is , the Stefan number is 0.045, and the Prandtl number is After a grid number and time step test, the grid number used in the simulation was and the time step was 0.1. Figures 3 and 4 show the positions of melting fronts obtained by the interfacial tracking method compared with experimental results by Okada [40] at different dimensionless times. At the early time of s ¼ 39.9, the melting interfacial velocity is slower than the experimental result at the bottom of the enclosure; however, the melting interfacial velocity obtained by the numerical solution is faster than the experimental results at the top of the enclosure. The volume faction of the liquid obtained by the interfacial tracking method is about 3.3% lower than that of the experimental results. At a longer time of s ¼ 78.68, the interfacial tracking method gives results very close to the experimental results. The difference between the predicted and measured liquid volume fraction is only 0.2%. The results in Figures 3 and 4 demonstrate that the interfacial tracking method can correctly solve melting problems at low Rayleigh number. In order to make sure that the interfacial tracking method is also valid at high Rayleigh numbers, an additional numerical simulation was performed based on the conditions specified by Ho and Viskanta [28]. The aspect ratio of the enclosure is A ¼ The subcooling parameter is equal to 0.01, the Rayleigh number is , the Stefan number is 0.09, and the Prandtl number is Since the top surface of the PCM is a free surface in [28], the velocity boundary condition at the top Figure 3. Comparison of the locations of the melting fronts (Ra ¼ , s ¼ 39.9).
11 INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS 219 Figure 4. Comparison of the locations of the melting fronts (Ra ¼ , s ¼ 78.68). specified in Eq. (19) is changed to qu Y ¼ 1 qy ¼ 0 ð35þ Since the boundary layer at the heated wall and the solid liquid interface decreases with increasing Rayleigh number, a finer grid must be used for the case of higher Rayleigh number to capture the behavior of the natural convection. After a grid number test, the grid number used in the simulation was 80 80, which was different from Ho and Viskata s Figure 5 shows the comparison between the locations of solid liquid interfaces obtained using the present model and that of Ho and Viskanta [28]. At early dimensionless time of s ¼ , the agreement is very good at the bottom and top of the rectangular cavity. As melting continues to the dimensionless time of s ¼ , the effect of the natural convection makes the top portion of the liquid region wider and the results obtained by the numerical and experimental results once again agree very well. When the dimensionless time s ¼ , the effect of natural convection on the melting becomes very strong and the interfacial tracking method is capable of simulating the strong natural convection during melting. Comparisons among all the numerical results obtained by the interfacial tracking method and experimental results indicate that the interfacial tracking method can obtain very good results even at higher Rayleigh number. The molten volume fraction of liquid is a very important parameter to assess the overall melting rate. The predicted instantaneous molten volume fraction, Vol=Vol 0, is evaluated from the solid liquid interface position by a numerical
12 220 Q. CHEN ET AL. Figure 5. Comparison of the locations of the melting fronts (Ra ¼ 10 8 ). integration of the instantaneous position and is defined as Vol=Vol 0 ¼ 1 Z H sdy¼ 1 Z A SdY ð36þ HW 0 A 0 Figure 6 shows a comparison of the molten volume fraction obtained from the interfacial tracking method and the experiments [28]. It can be seen that the liquid volume fraction obtained by numerical simulation agrees very well with the measured results. This comparison once again demonstrates that the interfacial tracking method can be used to solve melting problems at high Rayleigh number. Additional simulation was carried out under the same conditions as Gadgil and Gobin s [29] numerical solution based on a deformed-grid approach. Subcooling was not considered, and the Rayleigh number was The Stefan number was 0.2 and the Prandtl number was 50. The grid number was and the time step was A larger grid number was used because the Stefan number of 0.2 is higher than the 0.09 for Figure 5, and natural convection under higher Stefan number is stronger. Figure 7 shows a comparison of results obtained from the interfacial tracking method and Gadgil and Gobin s [29] numerical solution based on a deformed grid. The results obtained by the present numerical solution are represented by various lines, while the results obtained by Gadgil and Gobin s [29] are represented by symbols. At early dimensionless time of s ¼ 0.54, the agreement of the results is excellent. As time progresses to s ¼ 0.76, the natural convection plays a more important role
13 INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS 221 Figure 6. Predicted and experimental molten volume fractions as a function of dimensionless time. and the two results are still very close to each other. At later dimensionless time of s ¼ 1.10, the result obtained by the interfacial tracking method is slower than that of the Gadgil and Gobin s [29] numerical result in the lower portion of the rectangular. Gadgil and Gobin carried out their simulation by dividing the transient process into a large number of quasi-static steps. In each quasi-static step, steady-state natural convection in the liquid phase was calculated. On the other hand, the present interfacial tracking method does not require the quasi-steady-state assumption, Figure 7. Time evolution of the melting-front position.
14 222 Q. CHEN ET AL. Figure 8. Heat transfer from the heated wall as a function of time. and transient natural convection was considered. Therefore, the location of the melting front obtained by the interfacial tracking method is slower than that obtained by Gadgil and Gobin [29]. However, the overall melting rates obtained by the present method and that of [29] for all three different instances agreed very well. The Nusselt number at the heated wall was further calculated. The Nusselt number gives the heat transfer from the heated wall to the liquid domain, and it can be obtained as Nu ¼ hh k where the heat transfer coefficient h can be obtained from k qt qx ¼ hðt h T m Þ x¼0 Substituting Eq. (38) into Eq. (37), the dimensionless Nusselt number can be obtained as Nu ¼ qh qx X¼0 ð37þ ð38þ ð39þ The average Nusselt number on the heated wall is Z 1 qh Nu ¼ qx dy X¼0 0 ð40þ Figure 8 shows a quantitative comparison of the average Nusselt number at the heated wall obtained by different methods. During early time, two results are close to each other. After the dimensionless time reaches 1, the Nusselt number obtained by Gadgil and Gobin [29] becomes constant very quickly. On the contrary, the Nusselt number obtained by the present method decreases slowly after reaching
15 INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS 223 dimensionless time 1. This different trend may also be attributed to the different ways that the transient effect is treated. CONCLUSIONS An interfacial tracking method was proposed and tested for two dimensional convection-controlled melting problems at different Rayleigh numbers and Stefan numbers. The interfacial tracking method combines the advantages of both fixedand deforming-grid schemes. It can successfully solve multidimensional problems with or without natural convection. It retains the simplicity of grid construction of fixed-grid schemes while maintaining the accuracy of the energy balance at the solid interface of deforming grids. This simplicity makes the computation much faster than using the deforming-grid method. It also eliminates the need for of assuming a range of phase-change temperature as in other fixed grid methods. It does not product nonlinear oscillation of the temperature and interfacial location as enthalpy and equivalent heat capacity methods do. The interfacial tracking method can solve melting problems accurately at different Rayleigh and Stefan numbers. The results are much closer to the experimental results and other numerical results, and it also yields better results for cases with higher Rayleigh number. However, the interfacial tracking method requires more grids to deal with energy balance at the solid liquid interface. Although the method is presented for 2-D melting in a rectangular cavity, it can be extended to any other geometric configuration as long as the shape of the interface can be described by a continuous differentiable function. REFERENCES 1. Y. Zhang and A. Faghri, Heat Transfer Enhancement in Latent Heat Thermal Energy Storage System by Using an External Radial Finned Tube, J. Enhanced Heat Transfer, vol. 3, pp , Y. Zhang and A. Faghri, Vaporization, Melting and Heat Conduction in the Laser Drilling Process, Int. J. Heat Mass Transfer, vol. 42, pp , K. Mundra, T. DebRoy, and K. M. Kelkar, Numerical Prediction of Fluid Flow and Heat Transfer in Welding with Moving Heat Source, Numer. Heat Transfer A, vol. 29, pp , Y. Zhang, A. Faghri, C. W. Buckley, and T. L. Bergman, Three-Dimensional Sintering of Two-Component Metal Powders with Stationary and Moving Laser Beams, ASME J. Heat Transfer, vol. 122, pp , V. V. Tyagi and D. Buddhi, PCM Thermal Storage in Buildings: A State of Art, Renewable Sustainable Energy Rev., vol. 11, pp , A. F. Regin, S. C. Solanki, and J. S. Saini, Heat Transfer Characteristics of Thermal Energy Storage System Using PCM Capsules: A Review, Renewable Sustainable Energy Rev., vol. 12, pp , A. Sharma, V. V. Tyagi, C. R. Chen, and V. Buddhi, Review on Thermal Energy Storage with Phase Change Materials and Applications, Renewable Sustainable Energy Rev., vol. 13, pp , M. M. Farid, A. M. Khudhair, S. A. K. Razack, and S. Al- Hallaj, A Review on Phase Change Energy Storage: Materials and Applications, Energy Conversion Management, vol. 45, pp , 2004.
16 224 Q. CHEN ET AL. 9. M. M. Farid, Solar Energy Storage with Phase Change, J. Solar Energy Res., vol. 4, pp , N. Shamsundar and E. M. Sparrow, Effect of Density Change on Multidimensional Conduction Phase Change, J. Heat Transfer, vol. 97, pp. 550, A. B. Crowley, Numerical Solution of Stefan Problems, Int. J. Heat Mass Transfer, vol. 21, p. 215, C. Bonacina, G. Comini, A. Fasano, and M. Primicerio, Numerical Solution of Phase-Change Problems, Int. J. Heat Mass Transfer, vol. 16, p. 1825, J. Crank and R. S. Gupta, Isotherm Migration Method in Two Dimensions, Int. J. Heat Mass Transfer, vol. 18, p. 1101, C. F. Hsu, E. M. Sparrow, and S. V. Patankar, Numerical Solution of Moving Boundary Problems by Boundary Immobilization and a Control-Volume Based Finite Difference Scheme, Int. J. Heat Mass Transfer, vol. 24, p. 1335, E. M. Sparrow and W. Chuck, An Implicit=Explicit Numerical Solution Scheme for Phase-Change Problems, Numer. Heat Transfer, vol. 7, p. 1, E. M. Sparrow, S. Ramadhyani, and S. V. Patankar, Effect of Subcooling on Cylindrical Melting, J. Heat Transfer, vol. 100, p. 395, F. B. Cheung, T. C. Chawla, and D. R. Pedersen, The Effect of Heat Generation and Wall Interaction on Freezing and Melting in a Finite Slab, Int. J. Heat Mass Transfer, vol. 27, p. 29, P. Rattanadecho, Experimental and Numerical Study of Solidification Process in Unsaturated Granular Packed Bed, AIAA J. Thermophys. Heat Transfer, vol. 18, pp , P. Rattanadecho, The Theoretical and Experimental Investigation of Microwave Thawing of Frozen Layer Using Microwave Oven (Effects of Layered Configurations and Layered Thickness), Int. J. Heat Mass Transfer, vol. 47, pp , J. L. Duda, M. F. Malone, R. H. Notter, and J. S. Vrentas, Analysis of Two Dimensional Diffusion Controlled Moving Boundary Problems, Int. J. Heat Mass Transfer, vol. 18, p. 901, T. Saitoh, Numerical Method for Multi-Dimensional Freezing Problems in Arbitrary Domains, J. Heat Transfer, vol. 100, p. 294, Z.-X. Gong and A. S. Mujumdar, Flow and Heat Transfer in Convection Dominated Melting in a Rectangular Cavity Heated from Below, Int. J. Heat Mass Transfer, vol. 41, pp , W. Cao, W. Huang, and R. D. Russell, An r-adaptive Finite Element Method Based Upon Moving Mesh PDEs, J. Comput. Phys., vol. 149, p. 221, D. B. Khillarkar, Z. X. Gong, and A. S. Mujumdar, Melting of a Phase Change Material in Concentric Horizontal Annuli of Arbitrary Cross-Section, Appl. Thermal Eng., vol. 20, p , A. Chatterjee and V. Prasad, A Full 3-Dimensional Adaptive Finite Volume Scheme for Transport and Phase-Change Processes Part I: Formulation and Validation, Numer. Heat Transfer A, vol. 37, pp , G. Beckett, J. A. MacKenzie, and M. L. Robertson, A Moving Mesh Finite Element Method for the Solution of Two-Dimensional Stefan Problems, J. Comput. Phys., vol. 168, pp , V. R. Voller, An Overview of Numerical Methods for Solving Phase Change Problems, in W. J. Minkowycz and E. M. Sparrow (eds.), Advances in Numerical Heat Transfer, Vol. 1, Taylor & Francis, Basingestoke, UK, C. J. Ho and R. Viskanta, Heat Transfer during Melting from an Isothermal Vertical Wall, J. Heat Transfer, vol. 106, pp , A. Gadgil and D. Gobin, Analysis of Two-Dimensional Melting in Rectangular Enclosures in Presence of Convection, J. Heat Transfer, vol. 106, pp , 1984.
17 INTERFACIAL TRACKING MODEL FOR MELTING PROBLEMS K. Sasaguchi, A. Ishihara, and H. Zhang, Numerical Study on Utilization of Melting of Phase Change Material for Cooling of a Heated Surface at a Constant Rate, Num. Heat Transfer A, vol. 29, pp , B. Binet and M. Lacroix, Melting from Heat Sources Flush Mounted on a Conducting Vertical Wall, Int. J. Numer. Meth. Heat Fluid Flow, vol. 10, pp , K. Morgan, A Numerical Analysis of Freezing and Melting with Convection, Comput. Meth. Appl. Mech. Eng., vol. 28, pp , J. S. Hsiao, An Efficient Algorithm for Finite Difference Analysis of Heat Transfer with Melting and Solidification, ASME Paper 84-WA=HT-42, Y. Cao and A. Faghri, A Numerical Analysis of Phase Change Problem Including Natural Convection, J. Heat Transfer, vol. 112, pp , Z. Ma and Y. Zhang, Solid Velocity Correction Schemes for A Temperature Transforming Model for Convection Phase Change, Int. J. Numer. Meth. Heat Fluid Flow, vol. 16, pp , Y. Zhang and J. K. Chen, An Interfacial Tracking Method for Ultrashort Pulse Laser Melting and Resolidification of a Thin Metal Film, J. Heat Transfer, vol. 130, pp , S. V. Pantankar, Numerical Heat Transfer and Fluid Flow, McGraw-Hill, New York, G. X. Wang and E. F. Matthys, Numerical Melting of Phase Change and Heat Transfer During Rapid Solidification Processes: Use of Control Volume Integral with Element Subdivision, Int. J. Heat Mass Transfer, vol. 35, pp , E. D. Olsson and T. L. Bergman, Reduction of Numerical Fluctuations in Fixed Grid Solidification Simulations, AIAA=ASME Thermal Physics and Heat Transfer Conference, Seattle, WA, June 18 20, 1990, ASME HTD-vol. 130, pp M. Okada, Analysis of Heat Transfer during melting from a Vertical Wall, Int. J. Heat Mass Transfer, vol. 27, pp , 1984.
NUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS
NUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS Christiano Garcia da Silva Santim, chrisoff22@yahoo.com.br Luiz Fernando Milanez, milanez@fem.unicamp.br Universidade
More informationNUMERICAL SIMULATION OF UNSTEADY NATURAL CONVECTION FROM HEATED HORIZONTAL CIRCULAR CYLINDERS IN A SQUARE ENCLOSURE
Numerical Heat Transfer, Part A, 65: 715 731, 2014 Copyright # Taylor & Francis Group, LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407782.2013.846607 NUMERICAL SIMULATION OF UNSTEADY NATURAL
More informationNumerical Study of the Moving Boundary Problem During Melting Process in a Rectangular Cavity Heated from Below
American Journal of Applied Sciences 4 (4): 25-256, 2007 ISSN 546-9239 2007 Science Publications Corresponding Author: Numerical Study of the Moving Boundary Problem During Melting Process in a Rectangular
More informationNumerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall
Heat Transfer Research, 2011, Vol. 42, No. 3 Numerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall SAEID JANI, 1* MEYSAM AMINI, 2 and MOSTAFA MAHMOODI
More informationTransactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X
Numerical analysis of melting process in a rectangular enclosure J. Goscik Department of Mechanical Engineering, Technical University of Bialystok, 15-351 Bialystok, Wiejska 45C, Poland Abstract In the
More informationAnalysis of Forced Convection Heat Transfer in Microencapsulated Phase Change Material Suspensions
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 9, No. 4, October-December 1995 Analysis of Forced Convection Heat Transfer in Microencapsulated Phase Change Material Suspensions Yuwen Zhangh and Amir
More informationFarah SOUAYFANE, Pascal-Henry BIWOLE, Farouk FARDOUN Colloque des doctorants LJAD 2017 Barcelonnette - May 2017
Farah SOUAYFANE, Pascal-Henry BIWOLE, Farouk FARDOUN Colloque des doctorants LJAD 2017 Barcelonnette - May 2017 fsouayfane@unice.fr Outline I. Introduction II. Objectives III. Physical phenomena IV. Case
More informationCorresponding author.
EFFECTS OF BEAM SIZE AND PULSE DURATION ON THE LASER DRILLING PROCESS Nazia Afrin, Pengfei Ji, J. K. Chen and Yuwen Zhang 1 Department of Mechanical and Aerospace Engineering University of Missouri Columbia,
More informationSELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
More informationNumerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 541 552 Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus M. Sankar
More informationNatural Convection in Parabolic Enclosure Heated from Below
www.ccsenet.org/mas Modern Applied Science Vol. 5, No. 3; June 011 Natural Convection in Parabolic Enclosure Heated from Below Dr. Ahmed W. Mustafa (Corresponding auther) University of Tikrit, College
More informationANALYSIS OF NATURAL CONVECTION IN HORIZONTAL CONCENTRIC ANNULI OF VARYING INNER SHAPE
Numerical Heat Transfer, Part A, 68: 1155 1174, 2015 Copyright # Taylor & Francis Group, LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407782.2015.1032016 ANALYSIS OF NATURAL CONVECTION IN
More informationNATURAL CONVECTION HEAT TRANSFER IN PARTIALLY OPEN ENCLOSURES CONTAINING AN INTERNAL LOCAL HEAT SOURCE
Brazilian Journal of Chemical Engineering ISSN 0104-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 24, No. 03, pp. 375-388, July - September, 2007 NATURAL CONVECTION HEAT TRANSFER IN PARTIALLY OPEN
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationNATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A FLUSH MOUNTED HEATER ON A SIDE WALL
Journal of Naval Architecture and Marine Engineering December, 2010 DOI: 10.3329/jname.v7i2.3292 http://www.banglajol.info NATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A
More informationCombined Natural Convection and Thermal Radiation in an Inclined Cubical Cavity with a Rectangular Pins Attached to Its Active Wall
Periodicals of Engineering and Natural Sciences ISSN 2303-4521 Vol.5, No.3, November 2017, pp. 347~354 Available online at:http://pen.ius.edu.ba Combined Natural Convection and Thermal Radiation in an
More informationSigma J Eng & Nat Sci 36 (1), 2018, Sigma Journal of Engineering and Natural Sciences Sigma Mühendislik ve Fen Bilimleri Dergisi
Sigma J Eng & Nat Sci 36 (1), 2018, 49-62 Sigma Journal of Engineering and Natural Sciences Sigma Mühendislik ve Fen Bilimleri Dergisi Research Article THE ENERGY EFFICIENT CONFIGURATIONS OF NATURAL CONVECTION
More informationINSTRUCTOR: PM DR MAZLAN ABDUL WAHID
SMJ 4463: HEAT TRANSFER INSTRUCTOR: PM ABDUL WAHID http://www.fkm.utm.my/~mazlan TEXT: Introduction to Heat Transfer by Incropera, DeWitt, Bergman, Lavine 5 th Edition, John Wiley and Sons Chapter 9 Natural
More informationULTRAFAST PHASE CHANGE DURING FEMTOSECOND LASER INTERACTION WITH GOLD FILMS: EFFECT OF FILM THICKNESS
Numerical Heat Transfer, Part A, 57: 893 910, 2010 Copyright # Taylor & Francis Group, LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407782.2010.489491 ULTRAFAST PHASE CHANGE DURING FEMTOSECOND
More informationNatural convection adjacent to a sidewall with three fins in a differentially heated cavity
ANZIAM J. 48 (CTAC2006) pp.c806 C819, 2007 C806 Natural convection adjacent to a sidewall with three fins in a differentially heated cavity F. Xu 1 J. C. Patterson 2 C. Lei 3 (Received 31 August 2006;
More informationThis file contains the author's copy of:
This file contains the author's copy of: Tasnim, S.H., and Collins, M.R., "Suppressing Natural Convection in a Differentially Heated Square Cavity with an Arc Shaped Baffle", International Communications
More informationStudy on the flow and thermal characteristics of a heat storage system
THE ASIAN SYMPOSIUM ON COMPUTATIONAL HEAT TRANSFER AND FLUID FLOW - 2011, 22 26 SEPTEMBER 2011, KYOTO, JAPAN Study on the flow and thermal characteristics of a heat storage system Chung-Jen Tseng, Tzu-Yu
More informationNatural Convection Heat Transfer from a Rectangular Block Embedded in a Vertical Enclosure
Chapter 5 Natural Convection Heat Transfer from a Rectangular Block Embedded in a Vertical Enclosure Xiaohui Zhang Additional information is available at the end of the chapter http://dx.doi.org/10.5772/52666
More informationMELTING OF A PHASE CHANGE MATERIAL IN A HORIZONTAL ANNULUS WITH DISCRETE HEAT SOURCES
THERMAL SCIENCE, Year 2015, Vol. 19, No. 5, pp. 1733-1745 1733 MELTING OF A PHASE CHANGE MATERIAL IN A HORIZONTAL ANNULUS WITH DISCRETE HEAT SOURCES by Mohammad MASTIANI a*, Abdolrahman DADVAND b, Hooshyar
More informationEffectiveness NTU performance of finned PCM storage unit K.A.R. Ismail, M.M. Gongalves
Effectiveness NTU performance of finned PCM storage unit K.A.R. Ismail, M.M. Gongalves Campinas, SP, Brazil ABSTRACT This paper presents a mathematical model based upon two dimensional formulation of the
More informationNATURAL CONVECTIVE HEAT TRANSFER FROM A RECESSED NARROW VERTICAL FLAT PLATE WITH A UNIFORM HEAT FLUX AT THE SURFACE
HEFAT2007 5 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics Sun City, South Africa Paper number: OP2 NATURAL CONVECTIVE HEAT TRANSFER FROM A RECESSED NARROW VERTICAL FLAT
More informationHEFAT th International Conference on Heat Transfer, Fluid Mechanics, and Thermodynamics September 2005, Cairo, Egypt AA10
HEFAT5 4 th International Conference on Heat Transfer, Fluid Mechanics, and Thermodynamics 9- September 5, Cairo, Egypt AA Numerical Study of Natural Convection Heat Transfer in Enclosures with Conducting
More informationCharacteristics of forced convection heat transfer in vertical internally finned tube B
International Communications in Heat and Mass Transfer 32 (2005) 557 564 www.elsevier.com/locate/ichmt Characteristics of forced convection heat transfer in vertical internally finned tube B A. Al-Sarkhi*,
More informationMYcsvtu Notes HEAT TRANSFER BY CONVECTION
www.mycsvtunotes.in HEAT TRANSFER BY CONVECTION CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in
More informationNumerical Study for Free Convection Heat Transfer inside an Inclined Cavity with Cylindrical Obstacles
International ournal of Engineering and Technology Volume 3 No. 5, May, 13 Numerical Study for Free Convection eat Transfer inside an Inclined Cavity with Cylindrical Obstacles Khudheyer S. Mushatat College
More informationNatural Convection in Vertical Channels with Porous Media and Adiabatic Extensions
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions Assunta Andreozzi 1,a, Bernardo Buonomo 2,b, Oronzio Manca 2,c and Sergio Nardini 2,d 1 DETEC, Università degli Studi
More informationFINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
More informationNumerical Investigation of Conjugate Natural Convection Heat Transfer from Discrete Heat Sources in Rectangular Enclosure
Proceedings of the World Congress on Engineering 4 Vol II, WCE 4, July - 4, 4, London, U.K. Numerical Investigation of Conjugate Natural Convection Heat Transfer from Discrete Heat Sources in Rectangular
More informationPHYSICAL MECHANISM OF NATURAL CONVECTION
1 NATURAL CONVECTION In this chapter, we consider natural convection, where any fluid motion occurs by natural means such as buoyancy. The fluid motion in forced convection is quite noticeable, since a
More informationNumerical Analysis of Laminar Natural Convection in a Quadrantal Cavity with a Solid Adiabatic Fin Attached to the Hot Vertical Wall
Journal of Applied Fluid Mechanics, Vol., No., pp. 01-10, 2013. Available online at www.jafmonline.net, ISSN 13-32, EISSN 13-3. Numerical Analysis of Laminar Natural Convection in a Quadrantal Cavity with
More informationY. L. He and W. Q. Tao Xi an Jiaotong University, Xi an, China. T. S. Zhao Hong Kong University of Science and Technology, Kowloon, Hong Kong, China
Numerical Heat Transfer, Part A, 44: 399 431, 2003 Copyright # Taylor & Francis Inc. ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407780390206625 STEADY NATURAL CONVECTION IN A TILTED LONG CYLINDRICAL
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationEffect of an adiabatic fin on natural convection heat transfer in a triangular enclosure
American Journal of Applied Mathematics 2013; 1(4): 78-83 Published online November 10, 2013 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.20130104.16 Effect of an adiabatic fin on
More informationMaximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
More informationAnalysis of a Fluid Behavior in a Rectangular Enclosure under the Effect of Magnetic Field
Analysis of a Fluid Behavior in a Rectangular Enclosure under the Effect of Magnetic Field Y.Bakhshan and H.Ashoori Abstract In this research, a 2-D computational analysis of steady state free convection
More informationAn experimental investigation of the melting process in a rectangular enclosure
International Journal of Heat and Mass Transfer 42 (1999) 3659±3672 www.elsevier.com/locate/ijhmt An experimental investigation of the melting process in a rectangular enclosure Y. Wang, A. Amiri, K. Vafai*
More information6.9 Applications of SolidLiquid Phase Change
6.9 Applications of Solid-Liquid Phase Change 6.9.1 Latent Heat Thermal Energy Storage The theoretical model employed in this study is shown in Fig. 6.31. At the very beginning of the process (t =0), the
More informationSTUDY ON THE ROUTE TO CHAOS OF NATURAL CONVECTION IN CYLINDRICAL ENVELOPE WITH AN INTERNAL CONCENTRIC CYLINDER WITH SLOTS
Proceedings of the 14th International Heat Transfer Conference IHTC14 August 8-13, 010, Washington, DC, SA Proceedings of the International Heat Transfer Conference IHTC14 August 8-13, 010, Washington,
More informationEntropy 2011, 13, ; doi: /e OPEN ACCESS. Entropy Generation at Natural Convection in an Inclined Rectangular Cavity
Entropy 011, 13, 100-1033; doi:10.3390/e1305100 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Generation at Natural Convection in an Inclined Rectangular Cavity Mounir
More informationChapter 9 NATURAL CONVECTION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition in SI Units Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 9 NATURAL CONVECTION PM Dr Mazlan Abdul Wahid Universiti Teknologi
More informationNumerical simulation of transient forced convection in a square enclosure containing two heated circular cylinders
Numerical simulation of transient forced convection in a square enclosure containing two heated circular cylinders Fariborz Karimi a, Hongtao Xu a, Zhiyun Wang a, Mo Yang a, *, Yuwen Zhang b a School of
More informationNATURAL CONVECTION AND RADIATION IN CIRCULAR AND ARC CAVITY
Proceedings of the International Conference on Mechanical Engineering 9 (ICME9) - 8 December 9, Dhaka, Bangladesh ICME9-TH- NATURAL CONVECTION AND RADIATION IN CIRCULAR AND ARC CAVITY Naheed Ferdous, Md.
More informationand Training, Shuweikh, Kuwait
This article was downloaded by:[university of Minnesota] On: 26 July 2007 Access Details: [subscription number 731602659] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered
More informationA Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder
American Journal of Computational Mathematics, 2015, 5, 41-54 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/10.4236/ajcm.2015.51003 A Finite Element Analysis
More informationENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD
HEFAT 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 6 8 July Malta ENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD Nawaf
More informationHeat transfer increase with thin fins in three dimensional enclosures
157 Heat transfer increase with thin fins in three dimensional enclosures R. L. Frederick & S. Samper Universidad de Chile, Departamento de Ingeniería Mecánica, Santiago, Chile Abstract Heat transfer enhancement
More informationCOMPUTACIONAL SIMULATION OF TWO-DIMENSIONAL TRANSIENT NATURAL CONVECTION IN VOLUMETRICALLY HEATED SQUARE ENCLOSURE
COMPUTACIONAL SIMULATION OF TWO-DIMENSIONAL TRANSIENT NATURAL CONVECTION IN VOLUMETRICALLY HEATED SQUARE ENCLOSURE Camila Braga Vieira, camila@lasme.coppe.ufrj.br Jian Su, sujian@lasme.coppe.ufrj.br Universidade
More informationNatural Convection Inside a Rectangular Enclosure with Two Discrete Heat Sources Placed on Its Bottom Wall
Proceedings of the 2 nd World Congress on Momentum, Heat and Mass Transfer (MHMT 17) Barcelona, Spain April 6 8, 2017 Paper No. ENFHT 108 ISSN: 2371-5316 DOI: 10.11159/enfht17.108 Natural Convection Inside
More informationCOMPUTATIONAL ANALYSIS OF ENCAPSULATED THERMAL ENERGY PHASE CHANGE STORAGE SYSTEM: CYLINDRICAL AND SPHERICAL GEOMETRY
International Journal of Mechanical Engineering and Technology (IJMET) Volume 9, Issue 5, May 2018, pp. 662 668, Article ID: IJMET_09_05_073 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=9&itype=5
More informationMIXED CONVECTION HEAT TRANSFER FROM A PARTICLE IN SUPERCRITICAL WATER
THERMAL SCIENCE, Year 2016, Vol. 20, No. 2, pp. 483-492 483 MIXED CONVECTION HEAT TRANSFER FROM A PARTICLE IN SUPERCRITICAL WATER by Liping WEI, Youjun LU*, and Jinjia WEI State Key Laboratory of Multiphase
More informationA Mixed Finite Element Formulation for Solving Phase Change Problems with Convection
A Mixed Finite Element Formulation for Solving Phase Change Problems with Convection Youssef Belhamadia 1, Abdoulaye S. Kane 2, and André Fortin 3 1 University of Alberta, Campus Saint-Jean and Department
More informationNumerical Study of Laminar Free Convection About a Horizontal Cylinder with Longitudinal Fins of Finite Thickness
Published in International Journal of Thermal Sciences 007 This is author version post-print Archived in Dspace@nitr http://dspace.nitrkl.ac.in/dspace Numerical Study of Laminar Free Convection About a
More informationINFILTRATION AND SOLID-LIQUID PHASE CHANGE IN POROUS MEDIA. A Dissertation. Presented to the Faculty of the Graduate School of
INFILTRATION AND SOLID-LIQUID PHASE CHANGE IN POROUS MEDIA A Dissertation Presented to the Faculty of the Graduate School of University of Missouri-Columbia In Partial Fulfillment of the Requirements for
More informationRiyadh 11451, Saudi Arabia. ( a b,c Abstract
Effects of internal heat generation, thermal radiation, and buoyancy force on boundary layer over a vertical plate with a convective boundary condition a Olanrewaju, P. O., a Gbadeyan, J.A. and b,c Hayat
More informationEffect of Buoyancy Force on the Flow Field in a Square Cavity with Heated from Below
International Journal of Discrete Mathematics 017; (): 43-47 http://www.sciencepublishinggroup.com/j/dmath doi: 10.11648/j.dmath.01700.13 Effect of Buoyancy Force on the Flow Field in a Square Cavity with
More informationNUMERICAL STUDY OF HEAT TRANSFER IN A FLAT PLAT THERMAL SOLAR COLLECTOR WITH PARTITIONS ATTACHED TO ITS GLAZING. Adel LAARABA.
NUMERICAL STUDY OF HEAT TRANSFER IN A FLAT PLAT THERMAL SOLAR COLLECTOR WITH PARTITIONS ATTACHED TO ITS GLAZING Adel LAARABA. Department of physics. University of BATNA. (05000) Batna, Algeria Ccorresponding
More informationAnalysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
More informationAnalysis of Natural Convection Flow in a Trapezoidal Cavity Containing a Rectangular Heated Body in Presence of External Oriented Magnetic Field
Publications Available Online J. Sci. Res. 10 (1), 11-23 (2018) JOURNAL OF SCIENTIFIC RESEARCH www.banglajol.info/index.php/jsr Analysis of Natural Convection Flow in a Trapezoidal Cavity Containing a
More informationRAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8
HEFAT01 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 01 Malta RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8 Leong S.S. School of Mechanical
More informationNATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION
NATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION Éliton Fontana, eliton_fontana@hotmail.com Universidade Federal de Santa Catarina - UFSC Adriano
More informationNUMERICAL STUDIES OF TRANSITION FROM STEADY TO UNSTEADY COUPLED THERMAL BOUNDARY LAYERS
International Journal of Computational Methods Vol. 11, Suppl. 1 (214) 13442 (15 pages) c World Scientific Publishing Company DOI: 1.1142/S2198762134427 NUMERICAL STUDIES OF TRANSITION FROM STEADY TO UNSTEADY
More informationCOMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS
HEFAT214 1 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 214 Orlando, Florida COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES
More informationInternational Journal of Thermal Sciences
International Journal of Thermal Sciences 105 (2016) 137e158 Contents lists available at ScienceDirect International Journal of Thermal Sciences journal homepage: www. elsevier. com/ locate/ ijts Natural
More informationNumerical Study of Natural Convection in. an Inclined L-shaped Porous Enclosure
Adv. Theor. Appl. Mech., Vol. 5, 2012, no. 5, 237-245 Numerical Study of Natural Convection in an Inclined L-shaped Porous Enclosure S. M. Moghimi 1 *, G. Domairry 2, H. Bararnia 2, Soheil Soleimani 2
More informationNatural Convection from a Long Horizontal Cylinder
Natural Convection from a Long Horizontal Cylinder Hussein Awad Kurdi Saad Engineering Technical College of Al Najaf, Al-Furat Al-Awsat Technical University, Iraq ABSTRACT: Natural convection from a Long
More informationThermal conductivity measurement of two microencapsulated phase change slurries
Thermal conductivity measurement of two microencapsulated phase change slurries Xiaoli Ma (corresponding author), Siddig Omer, Wei Zhang and S. B. Riffat Institute of Sustainable Energy Technology, School
More informationNumerical simulation heat transfer by natural convection in liquid metal with a sinusoidal temperature
EPJ Web of Conferences 114, 02077 (2016) DOI: 10.1051/ epjconf/ 2016114 02077 C Owned by the authors, published by EDP Sciences, 2016 Numerical simulation heat transfer by natural convection in liquid
More informationNumerical Study of PCM Melting in Evacuated Solar Collector Storage System
Numerical Study of PCM Melting in Evacuated Collector Storage System MOHD KHAIRUL ANUAR SHARIF, SOHIF MAT, MOHD AFZANIZAM MOHD ROSLI, KAMARUZZAMAN SOPIAN, MOHD YUSOF SULAIMAN, A. A. Al-abidi. Energy Research
More informationNumerical solution of hyperbolic heat conduction in thin surface layers
International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,
More informationChapter 3 NATURAL CONVECTION
Fundamentals of Thermal-Fluid Sciences, 3rd Edition Yunus A. Cengel, Robert H. Turner, John M. Cimbala McGraw-Hill, 2008 Chapter 3 NATURAL CONVECTION Mehmet Kanoglu Copyright The McGraw-Hill Companies,
More informationMechanical Engineering. Postal Correspondence Course HEAT TRANSFER. GATE, IES & PSUs
Heat Transfer-ME GATE, IES, PSU 1 SAMPLE STUDY MATERIAL Mechanical Engineering ME Postal Correspondence Course HEAT TRANSFER GATE, IES & PSUs Heat Transfer-ME GATE, IES, PSU 2 C O N T E N T 1. INTRODUCTION
More informationInternational Journal of Multidisciplinary and Current Research
International Journal of Multidisciplinary and Current Research Research Article ISSN: 2321-3124 Available at: http://ijmcr.com in a Two-Phase Closed Cylindrical Thermosiphon in Conditions of Convective
More informationResearch Article Study of the Transient Natural Convection of a Newtonian Fluid inside an Enclosure Delimited by Portions of Cylinders
Research Journal of Applied Sciences, Engineering and Technology 7(5): 3069-3074, 04 DOI:0.906/rjaset.7.644 ISSN: 040-7459; e-issn: 040-7467 04 Maxwell Scientific Publication Corp. Submitted: August 7,
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationNatural Convection in Porous Triangular Enclosure with a Circular Obstacle in Presence of Heat Generation
American Journal of Applied Mathematics 2015; 3(2): 51-58 Published online March 20, 2015 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.20150302.14 ISSN: 2330-0043 (Print); ISSN:
More informationNatural Convection and Entropy Generation in a Porous Enclosure with Sinusoidal Temperature Variation on the Side Walls
Avestia Publishing Journal of Fluid Flow, Heat and Mass Transfer Volume 1, Year 14 Journal ISSN: 368-6111 DOI: 1.11159/jffhmt.14.4 Natural Convection and Entropy Generation in a Porous Enclosure with Sinusoidal
More informationFREE CONVECTION OF A NANOFLUID IN A SQUARE CAVITY WITH A HEAT SOURCE ON THE BOTTOM WALL AND PARTIALLY COOLED FROM SIDES
Abbasian Arani A. A., et al.: Free Convection of a Nanofluid in a Square Cavity S283 FREE CONVECTION OF A NANOFLUID IN A SQUARE CAVITY WITH A HEAT SOURCE ON THE BOTTOM WALL AND PARTIALLY COOLED FROM SIDES
More informationOn the influence of tube row number for mixed convection around micro tubes
Thessaloniki, Greece, 22-24 August 211 On the influence of tube row number for mixed convection around micro tubes Chuanshan DAI, Qiuxiang WANG, Biao LI * Corresponding author: Tel.: +86-22274183; Fax:
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (24 ) 55 556 th International Conference on Mechanical Engineering, ICME 23 Analysis of heat transfer and flow due to natural
More informationSpecific heat capacity. Convective heat transfer coefficient. Thermal diffusivity. Lc ft, m Characteristic length (r for cylinder or sphere; for slab)
Important Heat Transfer Parameters CBE 150A Midterm #3 Review Sheet General Parameters: q or or Heat transfer rate Heat flux (per unit area) Cp Specific heat capacity k Thermal conductivity h Convective
More informationEFFECT OF VARYING THE HEATED LOWER REGION ON FLOW WITHIN A HORIZONTAL CYLINDER
ISTP-1, 5, PRAGUE 1 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA EFFECT OF VARYING THE HEATED LOWER REGION ON FLOW WITHIN A HORIZONTAL CYLINDER S. S. Leong School of Mechanical and Manufacturing Engineering
More informationEFFECT OF HEATED WALL POSITION ON MIXED CONVECTION IN A CHANNEL WITH AN OPEN CAVITY
Numerical Heat Transfer, Part A, 43: 259 282, 2003 Copyright # 2003 Taylor & Francis 1040-7782/03 $12.00 +.00 DOI: 10.1080/10407780390122664 EFFECT OF HEATED WALL POSITION ON MIXED CONVECTION IN A CHANNEL
More informationC ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127
C ONTENTS Preface xviii Nomenclature xxvi CHAPTER ONE BASICS OF HEAT TRANSFER 1 1-1 Thermodynamics and Heat Transfer 2 Application Areas of Heat Transfer 3 Historical Background 3 1-2 Engineering Heat
More informationNATURAL CONVECTION OF AIR IN TILTED SQUARE CAVITIES WITH DIFFERENTIALLY HEATED OPPOSITE WALLS
Proceedings of the International onference on Mechanical Engineering 0 (IME0 8-0 December 0, Dhaka, Bangladesh IME- NATURAL ONVETION OF AIR IN TILTED SQUARE AVITIES WIT DIFFERENTIALLY EATED OPPOSITE WALLS
More informationA Numerical Study of Forced Convection Heat Transfer for Staggered Tube Banks in Cross-Flow
A Numerical Study of Forced Convection Heat Transfer for Staggered Tube Banks in Cross-Flow T. A. Tahseen 1, M. Ishak 1,2 and M. M. Rahman 1,2 1 Faculty of Mechanical Engineering, University Malaysia Pahang
More informationSIMULATION OF MIXED CONVECTIVE HEAT TRANSFER USING LATTICE BOLTZMANN METHOD
International Journal of Automotive and Mechanical Engineering (IJAME) ISSN: 2229-8649 (Print); ISSN: 2180-1606 (Online); Volume 2, pp. 130-143, July-December 2010 Universiti Malaysia Pahang DOI: http://dx.doi.org/10.15282/ijame.2.2010.3.0011
More informationOPTIMAL POSITIONING OF STRIPS FOR HEAT TRANSFER REDUCTION WITHIN AN ENCLOSURE
Numerical Heat Transfer, Part A, 66: 17 40, 2014 Copyright Taylor & Francis Group, LLC ISSN: 1040-7782 print/1521-0634 online DOI: 10.1080/10407782.2013.869081 OPTIMAL POSITIONING OF STRIPS FOR HEAT TRANSFER
More informationFALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
Proceedings of the International Conference on Mechanical Engineering 2 (ICME2) 8-2 December 2, Dhaka, Bangladesh ICME-TH-6 FALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by:[cdl Journals Account] On: 6 January 2008 Access Details: [subscription number 785022370] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered
More informationLAMINAR NATURAL CONVECTION IN VERTICAL 2D GLAZING CAVITIES
Mechanical and Industrial Engineering University of Massachusetts, Amherst AMINAR NATURA CONVECTION IN VERTICA 2D GAZING CAVITIES Bhaskar Adusumalli ABSTRACT Finite element predictions of natural convection
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationFlow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers ( )
Advances in Fluid Mechanics VII 391 Flow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers (10 6 10 9 ) R. L. Frederick & S. Courtin
More informationHeat Transfer Convection
Heat ransfer Convection Previous lectures conduction: heat transfer without fluid motion oday (textbook nearly 00 pages) Convection: heat transfer with fluid motion Research methods different Natural Convection
More informationOscillatory MHD Mixed Convection Boundary Layer Flow of Finite Dimension with Induced Pressure Gradient
Journal of Applied Fluid Mechanics, Vol. 9, No., pp. 75-75, 6. Available online at www.jafmonline.net, ISSN 75-57, EISSN 75-65. DOI:.8869/acadpub.jafm.68.5.876 Oscillatory MHD Mixed Convection Boundary
More informationCOMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
Suranaree J. Sci. Technol. Vol. 20 No. 4; October - December 2013 257 COMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
More information