Benchmark solutions for twodimensional. transfer problems in irregular regions using multigrid method. Jingfa Li, Bo Yu and Min Wang.
|
|
- Dominic McKenzie
- 5 years ago
- Views:
Transcription
1 Special Issue Article Benchmark solutions for twodimensional fluid flow and heat transfer problems in irregular regions using multigrid method Advances in Mechanical Engineering 2015, Vol. 7(11) 1 17 Ó The Author(s) 2015 DOI: / aime.sagepub.com Jingfa Li, Bo Yu and Min Wang Abstract In the previous researches about the benchmark solutions in the irregular regions, there are some obvious shortcomings, such as small number of computational grids, simple computational domain, small characteristic number, and lack of mixed convection benchmark solution. In order to improve this situation, this article deeply studied the benchmark solutions for two-dimensional fluid flow and heat transfer problems in the irregular regions under the body-fitted coordinate system. Taking the lid-driven flow, natural convection, and mixed convection problems as the research objects and considering the influence of boundary configurations of the computational domains, Reynolds number, Prandtl number, and Grashof number, this article designed three groups of test cases. The program code is first verified through employing the nonlinear multigrid method based on the collocated finite volume method and then the numerical solutions of three test cases with grids and convergence criterion of are obtained to get the benchmark solutions estimated by the Richardson extrapolation method. Keywords Benchmark solutions, irregular regions, body-fitted coordinate system, multigrid, finite volume method Date received: 30 March 2015; accepted: 15 October 2015 Academic Editor: Chin-Lung Chen Introduction The fluid flow and heat transfer phenomena in practice always happen in the irregular regions, so how to achieve the accurate and efficient calculations of fluid flow and heat transfer problems in the irregular regions is no doubt of great importance. Currently, many articles published in the international academic journals are related to the numerical simulations of the fluid flow and heat transfer problems in the irregular regions. Thus, in order to verify the calculation results and provide references for different algorithms, it is quite necessary to study the benchmark solutions of fluid flow and heat transfer problems defined in the irregular regions. As one of the popular methods to discretize the irregular regions, the body-fitted coordinate grids have been widely used due to its good applicability and desirable ability to follow the boundaries, 1 which is mainly applied to internal flows over smoothly varying irregular boundaries. Therefore, this article aims to research National Engineering Laboratory for Pipeline Safety, Beijing Key Laboratory of Urban Oil and Gas Distribution Technology, China University of Petroleum Beijing, Beijing, China Corresponding author: Bo Yu, National Engineering Laboratory for Pipeline Safety, Beijing Key Laboratory of Urban Oil and Gas Distribution Technology, China University of Petroleum Beijing, Beijing , China. yubobox@vip.163.com Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).
2 2 Advances in Mechanical Engineering on the benchmark solutions of fluid flow and heat transfer problems in the two-dimensional (2D) irregular regions with the body-fitted coordinate grids. As a basic research topic, the benchmark solutions have been studied and reported systematically by researchers all over the world, especially the lid-driven flow 2,3 and natural convection 4,5 in a square cavity, and the calculated benchmark solutions have already been verified and proved. For the benchmark solutions in the irregular regions, the cases in the literature usually can be solved by utilizing the cylindrical coordinate system 6 or polar coordinate system, 7 but the method using the body-fitted coordinate system is rarely proposed. The limited literatures involving these issues are as follows: in the literature, Demirdzic et al. 8 showed the benchmark solutions of the lid-driven flow in an inclined square cavity with the angles of 30 and 45, natural convection in a 45 inclined square cavity, and a cylinder which is enclosed by a square duct. The Reynolds number is set to be 100 and 1000, respectively, in the case of lid-driven flow; the Prandtl number is set to be 0.1 and 10, respectively, while Rayleigh number is 10 6 in the case of natural convection. In that article, the multigrid (MG) method based on the collocated finite volume method (FVM) is used and coupling between pressure and velocity is treated by SIMPLE algorithm. And the benchmark solutions with grids are given. Besides, Oosterlee et al. 9 tested the correctness of the algorithm through comparing the results in the literature 8 and presented the benchmark solutions of L-shaped lid-driven cavity flow with grids and Reynolds number of 100 and Different from the literature, 8 Oosterlee et al. 9 employed the symmetric coupled alternating lines (SCAL) to deal with the coupling between pressure and velocity and solved the discretized equations by the means of nonlinear MG method based on the staggered FVM. Although the literature 8,9 has studied the benchmark solutions of 2D lid-driven cavity flow and natural convection in the irregular regions, there is still some room for improvement. For example, the maximum number of grids in the literature 8 is only in which case the numerical results may still be affected by the discretization error. And the maximum Reynolds number in the literature 8,9 is only 1000 and there is no research regarding the problems with high Reynolds number. In addition, the configurations in previous researches are relatively simple to indicate the influence of boundary shape on the accuracy and robustness of the numerical simulation. To the best knowledge of the authors, there are also no such literatures studying the benchmark solutions of mixed convection problems under the body-fitted coordinate system. Therefore, in order to improve this situation, this article aims to research on the benchmark solutions for 2D fluid flow and heat transfer problems in the irregular regions under the body-fitted coordinate system using MG method. The layout of this article is organized as follows. In section Test cases, the reasons and purposes to design the test cases are explained in detail, and the dimensionless general governing equations as well as the calculation conditions are given in this section. In section Numerical methods, we mainly focus on the numerical methods for solving the test cases. In section Program testing, the correctness of algorithms is tested through comparing the calculation results of liddriven flow and natural convection with those in the literature. 8 On this basis, section Benchmark solutions shows the benchmark solutions of different test cases under various calculation conditions. In the final section, the most important findings of this study are summarized. Test cases Description of test cases This article takes three typical cases of fluid flow and heat transfer problems, lid-driven flow, natural convection, and mixed convection, to be the research objects, noting them as test cases 1 3, respectively. The computational domains of the three test cases all contain the irregular configurations shown in Figure 1 where l m. The reasons and purposes of designing the three test cases are as follows: 1. Because previous literature did not discuss the benchmark solutions of the mixed convection problem under the body-fitted coordinate system, this article designs this kind of case to expand the application scope of the benchmark solutions. 2. The boundary shapes in the literature 8,9 are relatively simple. To investigate the influence of the boundary shape on the numerical accuracy, this article designs the regions with non-parallel east west boundaries as shown in Figure 1(b) and (c). The function of the west boundary shape in Figure 1(c) is x(y)=0:1 sin½p 2p(y=l)Š. To guarantee the reliability of the benchmark solutions, this article first compared the numerical results of lid-driven flow on the domain shown in Figure 1(a) with the reference results in the literature 8 to verify the reliability and accuracy of the present algorithms and code and then utilizes the solid program to carry out the subsequent numerical calculations.
3 Li et al. 3 Table 1. Parameters of the continuity equation and energy equation. Test cases Physical problems Continuity equation Energy equation F G S(j, h) F G S(j, h) 1 Lid-driven flow Natural convection Y 1/Pr 0 3 Mixed convection Y 1/Pr/Re 0 calculations in the governing equations. For the sake of concision, the details of dimensionless method are not provided. In the physical plane x y, the 2D dimensionless unsteady convection diffusion equation of the incompressible fluid is given as F t + (UF) X + Y + (VF) = Y X F G + R(X, Y ) Y G F X ð1þ where F, G, and R(X, Y) take different values in different test cases. Since the computational domains in Figure 1 are all irregular regions, this article transforms the dimensionless governing equation (1) in the physical plane x y into the computational plane j h for solution and the following governing equation can be obtained (Appendix 2) F t + 1 (U F) J + 1 (V F) j J h 1 G J j J (a F j b F h ) G J h J (g F h b F j ) + S(j, h) ð2þ Figure 1. Sketch map of the irregular computational domains: (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity. As shown in Tables 1 and 2, when F, G,andS(j, h) in equation (2) take different values, the dimensionless governing equations of lid-driven flow, natural convection, and mixed convection can be represented, respectively. Table 3 provides the dimensionless boundary conditions and calculation conditions of test cases 1 3. Dimensionless general governing equations and calculation conditions Different from the literature, 8,9 this article first applies the dimensionless method to the governing equations of the physical problems in test cases 1 3. The dimensionless treatment can not only make the numerical results more universal but also cover the magnitude gap between different physical variables to enhance the numerical accuracy and reduce the constant Numerical methods In order to guarantee the conservation of the dimensionless governing equation (2), this article applies the FVM to solve the three physical problems. The secondorder central difference scheme is employed to discretize the convection term and diffusion term. The deferred correction technology is used for the discretization of the convection term to better meet the diagonal dominant condition of coefficient matrix of the
4 4 Advances in Mechanical Engineering Table 2. Parameters of the momentum equations. Test cases Physical problems x momentum equation y momentum equation F G S(j, h) F G S(j, h) 1 Lid-driven flow U 1/Re (P h Y j P j Y h )=J V 1/Re (P j X h P h X j )=J 2 Natural convection U 1 (P h Y j P j Y h )=J V 1 (P j X h P h X j )=J + GrY 3 Mixed convection U 1/Re (P h Y j P j Y h )=J V 1/Re (P j X h P h X j )=J + RiY Table 3. Dimensionless boundary conditions and calculation conditions of test cases 1 3. Test cases Physical problems Dimensionless boundary conditions Dimensionless calculation conditions 1 Lid-driven flow North: U =1,V = 0; else: U = V =0 Re 00, 1000, Natural convection West: Y h ; east: Y c = 0; else: adiabatic Pr = 0.1, 10; Ra = Mixed convection North: U =1,V = 0; else: U = V =0 West: Y h ; east: Y c = 0; else: adiabatic Gr =10 6,10 7 ; Pr = 0.71; Re 000 discretized equations. To further satisfy the conservation of flux, the coupling between pressure and velocity is treated by the SIMPLE algorithm on collocated grids, where the dimensionless contravariant U* and V* are selected as the velocity components of the continuity equation in the computational plane and dimensionless U and V are selected as the velocity components of the momentum equations. The crossderivative terms in the continuity equation and momentum equations caused by the non-orthogonality of grids are treated as the source terms. To make the programming easier, this article utilizes the additional source term method to deal with the boundary conditions of the governing equation. To obtain the high-resolution numerical solutions and reduce the influence of discretization errors on the numerical results, the computational domain in computational plane is divided into uniform grids. And the convergence criterion is strictly set as the absolute residuals of the continuity equation less than The governing equations hold high nonlinearity because the lid-driven flow, natural convection, and mixed convection are all the multi-variable coupling problems. Besides, the number of computational grids is also large for grids adopted. As a result, the above factors would lead to the slow convergence rate of the computation process. To speed up computation, this article adopts the nonlinear MG method based on collocated FVM for the solution where the full approximation scheme (FAS) is chosen as the nonlinear version of the MG method. In FAS, V-cycle is adopted and strong implicit procedure (SIP) method is employed to solve the discretized equations. It is not necessary to obtain the precise solutions on coarse grid when MG method is adopted. Thus, the convergence criterion is only set on the finest mesh, while on other grid levels the number of smoothing is given. On the coarse grid, the numbers of pre-smoothing and postsmoothing for solving the momentum equations and energy equation are both set as 4, and set as 12 for solving the pressure correction equation. In the calculation, grid levels are determined by the coarsest grid number, which is set as in lid-driven flow and natural convection and in mixed convection, respectively. The under-relaxation factors for different test cases are as follows: the velocity and pressure under-relaxation factors in the case of lid-driven flow are 0.7 and 0.2, respectively; the velocity, pressure, and temperature under-relaxation factors in the cases of natural convection and mixed convection are 0.4, 0.1, and 0.5, respectively. To prove the computational efficiency of MG method, a set of numerical results such as computational cost, convergence rate, and convergence ratio in the computation of lid-driven flow in parallelogram cavity with Re 00 on grids are provided in Tables 4 and 5. Table 4 shows the computational cost for solving U-momentum equation, V-momentum equation, and P-correction equation, respectively, in a V-cycle of MG method. In Table 5, the convergence rate and convergence ratio are provided. The convergence rate is the maximum residual after each V-cycle, and the convergence ratio gives the ratio of residual of successive V-cycles. From the results, it can be seen that MG method is a wise choice in the present computation. Program testing In this section, the numerical results of lid-driven flow and natural convection are first given, and then the comparisons between the present numerical results and
5 Li et al. 5 Table 4. Computational cost of a V-cycle (y 1, y 2 ). Equation y 1 (Pre-smoothing) y 2 (Post-smoothing) WU U-momentum equation /3 V-momentum equation /3 P-correction equation WU: work unit. Table 5. Convergence rate and ratio. V-cycle Convergence rate Ratio V-cycle Convergence rate Ratio Table 6. Comparison between the numerical results of the lid-driven flow in this article with the reference results in the literature. 8 b = 45 8 Re 00 Re 000 Minimum Maximum Minimum Maximum c Present Demirdzic et al x Present Demirdzic et al y Present Demirdzic et al e Present Demirdzic et al the reference results in the literature 8 are conducted to verify the accuracy and reliability of the algorithms in this article, as shown in Tables 6 and 7 and Figure 2. From Table 6, Figure 2, and Table 7, it can be clearly seen that the numerical results in this article meet very well with the reference results in the literature, 8 and the numerical accuracy of the maximum Nusselt number and stream function values is higher than that in the literature. 8 In addition, the curves of U and V velocity profiles all go through the 15 selected points of the previous benchmark solutions in the literature. 8 Through the above comparisons, it can be reasonably concluded that the proposed algorithms are reliable and accurate to ensure the subsequent numerical calculations. Benchmark solutions In section Program testing, the algorithm in this article has been verified by the lid-driven flow and natural
6 6 Advances in Mechanical Engineering Table 7. Comparison between the numerical results of the natural convection in this article with the reference results in the literature. 8 Pr Pr Nu Present c max Present Demirdzic et al Demirdzic et al e Present x Present Demirdzic et al Demirdzic et al Nu h, max Present y Present Demirdzic et al Demirdzic et al y Present e Present Demirdzic et al Demirdzic et al Nu c, min Present c min Present Demirdzic et al. 8 Demirdzic et al y Present x Present Demirdzic et al. 8 Demirdzic et al Nu c Present y Present Demirdzic et al Demirdzic et al e Present Demirdzic et al Figure 2. Comparison between the velocity profiles along centerlines CL1 and CL2 for the lid-driven flow in this article with the reference results in the literature: 8 (a) U and (b) V. convection in the parallelogram cavity. In the following text, the benchmark solutions of the lid-driven flow, natural convection, and mixed convection will be calculated on the computational domains with grids. The detailed calculation conditions are presented in Table 3 of section Test cases. Lid-driven flow The numerical results of lid-driven flow with Reynolds number of 100 and 5000 are given below. Figure 3 is the streamlines for the lid-driven flow. It can be clearly seen that there is substantial difference between the contour map of Re = 5000 with that of Re =100. To be specific, the first vortex of Re = 5000 is much smaller than that of Re 00 and it is concentrated on the top right corner of the computational domain. However, the second vortex of Re = 5000 is bigger than that of Re 00 and it fills about three-fourths of the whole cavity. The reason for this result is that the driving force increases along with the Reynolds number and the fluid movement also becomes more intense. Thus, when the first vortex gets smaller, the second vortex has larger development space and the more intense fluid flow further enhances the full development of the second vortex. It can also be observed that the number of vortex increases with Reynolds number: when Re 00, there is only one main vortex and it fills almost the whole cavity; but when Re = 5000, there are three to five vortexes which means the fluid flow has been fully developed. Figure 3(a) (c) shows different development trends of the vortexes and the reason is the different computational domains and boundary shapes, leading to the vortex development in
7 Li et al. 7 Figure 3. Streamlines for the lid-driven flow when Re 00 and Re = 5000: (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity. Table 8. Minimum and maximum stream function values in vortex center and their position for the lid-driven flow on grids. Re Computational domain c min (x, y) c max (x, y) 100 (a) Parallelogram cavity (1.1136, ) ( , ) (b) Trapezoid cavity (1.4521, ) (2.0880, ) (c) Sine-type cavity ( , ) ( , ) 5000 (a) Parallelogram cavity (1.5075, ) ( , ) (b) Parallelogram cavity (1.5477, ) ( , ) (c) Sine-type cavity (1.0282, ) ( , ) different spaces. Figure 3 presents the fourth and fifth vortexes of the streamline which also suggests the high accuracy of this algorithm. Table 8 presents the minimum and maximum stream function values in vortex center and their position for the lid-driven flow on grids. And Figure 4 shows the V velocity profile at the centerline CL1 for the sake of qualitative comparison. Table 9 lists the selected values estimated by the Richardson extrapolation method for quantitative description of the V velocity profiles in Figure 4. Natural convection This section shows the numerical results of natural convection with Pr =10, Pr = 0.1, and Ra =10 6. Figure 5
8 8 Advances in Mechanical Engineering Figure 4. V velocity profile for the lid-driven flow along the centerline CL1: (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity. Figure 5. Streamlines for the natural convection when Ra =10 6 : (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity.
9 Li et al. 9 Table 9. Selected values of the velocity profiles in Figure 4. No. (a) Parallelogram cavity (b) Trapezoid cavity (c) Sine-type cavity Re 00 Re = 5000 Re 00 Re = 5000 Re 00 Re =5000 X V X V X V X V X V X V is the streamlines for the natural convection. We can see that there is only one main vortex when the Prandtl number is set to be 10 and 0.1, filling almost the whole computational domain with slight difference for the same area. The more intense flow when Pr = 0.1 than that when Pr 0 leads to the full development of the main vortex. It is obvious that Prandtl number has great influence on the shape of main vortex. Table 10 shows the minimum and maximum stream function values in vortex center and their position on grids in Figure 5. Figure 6 is the isotherms for the natural convection when Ra =10 6. It is shown that the isotherms are always vertical to the adiabatic boundaries, but for the same computational domain the curvature of isotherms differs from each other. When Pr = 0.1, the isotherms in the center of cavity have larger curvature than those when Pr 0, since when the Prandtl number is getting smaller, the heat diffusion ability of fluid flow becomes stronger than the momentum diffusion ability which results in the great buoyancy force. The convective movement between the cool and warm fluids becomes more intense, leading to the better mixing between fluids with different temperature. Figure 7 shows the Nusselt number along the hot wall when Ra =10 6 with Pr 0 and Pr = 0.1. Table 11 presents the selected values of the Nusselt number profiles in Figure 7. It is worth pointing out that the coordinate value Y of the selected points is the actual coordinate rather than the length of the hot wall. Mixed convection In order to expand the application scope of benchmark solutions, this article also studied the mixed convection problem. The mixed convection can be regarded as the combination between the lid-driven flow and natural convection; the driving forces of the fluid come from the force of the moving boundary as well as the buoyancy force. In the following text, the numerical results of mixed convection when Gr =10 6, Gr =10 7, and Pr = 0.71 are provided. Figure 8 is the streamlines for the mixed convection. It can be observed from Figure 8 that along with the increase in Grashof number, the number of vortex does not change but the shape changes a lot. When Gr =10 6, there is only one main vortex and the shape of the main vortex is concentric deformed circles which fill almost two-thirds of the cavity. When Gr =10 7, there is still one main vortex, but the shape is bigger than that of Gr =10 6 and presents as long strips which almost fill the whole region. The reason is that along with the increase in Grashof number, the intensity of natural convection becomes larger and the natural convection is dominant. In addition, the shape of
10 10 Advances in Mechanical Engineering Table 10. Minimum and maximum stream function values in vortex center and their position for the natural convection on grids. Pr Computational domain c min (x, y) c max (x, y) 10 (a) Parallelogram cavity ( , ) (1.7064, ) (b) Trapezoid cavity ( , ) (1.9771, ) (c) Sine-type cavity ( , ) (1.1201, ) 0.1 (a) Parallelogram cavity ( , ) (1.7064, ) (b) Trapezoid cavity (1.5217, ) ( , ) (c) Sine-type cavity ( , ) (1.1483, ) Figure 6. Isotherms for the natural convection when Ra =10 6 : (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity. streamlines when Gr =10 6 is similar to that in the liddriven flow when Re 000. Because when Gr =10 6 and Re 000, the Ri, which signifies that the influences of lid-driven flow and natural convection on the fluid flow are almost the same. The minimum and maximum stream function values in vortex center and their position are given in Table 12. Figure 9 is the isotherms for the mixed convection when Re 000 and Pr = It is shown that for different domains with various boundary shapes, the isotherms differ from each other greatly when Gr =10 6 and Gr =10 7. When Gr =10 6, the isotherms show the shape of vortex, while Gr =10 7,the shape of isotherms is similar to that of natural
11 Li et al. 11 Figure 7. Nusselt number along the hot wall for the natural convection when Ra =10 6 : (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity. Figure 8. Streamlines for the mixed convection when Re 000 and Pr = 0.71: (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity.
12 12 Advances in Mechanical Engineering Table 11. Selected values of the Nusselt number profile in Figure 7. No. (a) Parallelogram cavity (b) Trapezoid cavity (c) Sine-type cavity Pr =10 Pr = 0.1 Pr =10 Pr = 0.1 Pr =10 Pr = 0.1 Y Nuh Y Nuh Y Nuh Y Nuh Y Nuh Y Nuh
13 Li et al. 13 Figure 9. Isotherms for the mixed convection when Re 000 and Pr = 0.71: (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity. Table 12. Minimum and maximum stream function values in vortex center and their position for the mixed convection on grids. Gr Computational domain c min (x, y) c max (x, y) 10 6 (a) Parallelogram cavity ( , ) ( , ) (b) Trapezoid cavity (1.6385, ) (1.0374, ) (c) Sine-type cavity ( , ) (1.1280, ) 10 7 (a) Parallelogram cavity ( , ) (b) Trapezoid cavity (1.4432, ) (1.7962, ) (c) Sine-type cavity ( , ) (1.1426, ) convection when Gr =10 6. The reason for this phenomenon is the dominance of natural convection, and lid-driven flow on the mixed convection varies along with the different Grashof numbers. Although the three computational domains are quite different from each other, the isotherms are always vertical to the adiabatic boundaries. For the same calculation region, the curvature of isotherms when Gr =10 6 is larger than that when Gr =10 7, since the contour shape of the former situation is affected by lid-driven flow while the latter one is mainly influenced by the buoyancy force. Figure 10 shows the Nusselt number along the hot wall for the mixed convection when Re 000 and Pr = Table 13 gives the selected values of the Nusselt number profiles in Figure 10.
14 14 Advances in Mechanical Engineering Table 13. Selected values of the Nusselt number profiles in Figure 10. No. (a) Parallelogram cavity (b) Trapezoid cavity (c) Sine-type cavity Gr =10 6 Gr =10 7 Gr =10 6 Gr =10 7 Gr =10 6 Gr =10 7 Y Nuh Y Nuh Y Nuh Y Nuh Y Nuh Y Nuh
15 Li et al. 15 Figure 10. Nusselt number along the hot wall for the mixed convection when Re 000 and Pr = 0.71: (a) parallelogram cavity, (b) trapezoid cavity, and (c) sine-type cavity. Conclusion This article studied the benchmark solutions of the liddriven flow, natural convection, and mixed convection in 2D irregular regions under the body-fitted coordinate system. On the basis of verifying the reliability and correctness of the algorithms, this article employed the nonlinear MG method based on the collocated FVM to solve the above physical problems on grids under different calculation conditions. Based on the numerical results, the following conclusions can be obtained: 1. The results of lid-driven flow in parallelogram cavity verified the reliability of the present method and proved it to be more accurate than the benchmark solutions in the literature The benchmark solutions of lid-driven flow and natural convection on grids are given and this article also studied the influence of boundary shape, Reynolds number, Prandtl number, and Grashof number on the numerical results. 3. The benchmark solutions of the mixed convection are presented to expand the application scope of the benchmark solutions. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by the National Science Foundation of China (No , No and No ). References 1. Tao WQ. Numerical heat transfer. 2nd ed.xi an, China: Xi an Jiaotong University Press, Bruneau CH and Saad M. The 2D lid-driven cavity problem revisited. Comput Fluids 2006; 35: Botella O and Peyret R. Benchmark spectral results on the lid-driven cavity flow. Comput Fluids 1998; 27:
16 16 Advances in Mechanical Engineering 4. Hortmann M, Peric M and Scheuerer G. Finite volume multigrid prediction of laminar natural convection: benchmark solutions. Int J Numer Meth Fl 1990; 11: Barakos G, Mitsoulis E and Assimacopoulos D. Natural convection flow in a square cavity revisited: laminar and turbulent models with wall functions. Int J Numer Meth Fl 1994; 18: Mohamed AG, Valentine DT and Hassel RE. Numerical study of laminar separation over an annular backstep. Comput Fluids 1991; 20: Fuchs L and Tillmark N. Numerical and experimental study of driven flow in a polar cavity. Int J Numer Meth Fl 1985; 5: Demirdzic I, Lilek Z and Peric M. Fluid flow and heat transfer test problems for non-orthogonal grids: benchmark solutions. Int J Numer Meth Fl 1992; 15: Oosterlee CW, Wesseling P, Segal A, et al. Benchmark solutions for the incompressible Navier Stokes equations in general co-ordinates on staggered grids. Int J Numer Meth Fl 1993; 17: Appendix 1 Notation Jacobi factor, J = X j Y h X h Y j l length of computational domain (m) R(X, Y) source term in physical plane S(j, h) source term in computational plane U, V dimensionless velocities in the X- and Y-directions in physical plane U, V dimensionless contravariant velocities in the j- and h-directions in computational plane, U = UY h VX h, V = VX j UY j x, y coordinates in physical plane X, Y dimensionless coordinates in physical plane, X = x=l, Y = y=l J Greek symbols a, g metric coefficients in the j- and h-directions in computational plane, a = Xh 2 + Y h 2, g = Xj 2 + Y j 2 b angle of parallelogram computational domain b local orthogonality of grid in physical plane, b = X j X h + Y j Y h e relative error of numerical solution, e =(f j k + 1 f k )=f k + 1 j3 100 j, h coordinates in computational plane f k, f k + 1 numerical solutions on kth and k + 1th grid levels c c min, c max G Subscripts c cold wall h hot wall k grid level Characteristic number stream function minimum and maximum stream function values general dimensionless diffusion coefficient Gr Grashof number, Gr = gbðt h T c Þl 3 =n 2 Nu Nusselt number, Nu = hl=l Nu average Nusselt number, Nu = Ð b a Nu da= Ð b a da Nu c characteristic Nusselt number Nu c, min minimum Nusselt number along cold wall Nu h Nusselt number along hot wall Nu h, max maximum Nusselt number along hot wall Pr Prandtl number, Pr = rc p n=l Ra Rayleigh number, Ra = Pr Gr Re Reynolds number, Re = ul=n Ri Richardson number, Ri = Gr=Re 2 Appendix 2 Derivation of equation (2) The detailed derivation of equation (2) is as follows. In the physical plane x y, the two-dimensional dimensionless unsteady convection diffusion equation of the incompressible fluid is given below F t + (UF) X + Y + (VF) = Y X F G + R(X, Y ) Y G F X ð3þ To clearly show the derivation of equation (2), detailed transformation of every term in equation (3) is provided as follows 1. Unsteady term F t = F ð4þ t The unsteady term in computational plane is same as that in physical plane because it has nothing to do with spatial parameters.
17 Li et al Convection term (UF) X =(U jj X + U h h X )F + U(F j j X + F h h X ) = U jy h U h Y j J F + F jy h F h Y j J U J (U jf + F j U)Y h (U h F + F h U)Y j (UY h F) J (UY jf) ð5þ j h (VF) Y =(V jj Y + V h h Y )F + V(F j j Y + F h h Y ) = V hx j V j X h J F + F hx j F j X h J V J (V hf + VF h )X j (V j F + F j V)X h (VX j F) J (VX hf) ð6þ h j where J = X j Y h X h Y j, j X = Y h =J, h X = Y j =J, j Y = X h =J, and h Y = X j =J. Combining equations (5) and (6), we can obtain (UF) X + (VF) (UY h F) Y J (UY jf) j h + 1 (VX j F) J (VX hf) h j (UY h F) J (VX hf) j j + 1 (VX j F) J (UY jf) h h (UY h VX h )F J j + 1 (VX j UY j )F J h (U F) J + 1 (V F) j J h where U = UY h VX h and V = VX j UY j. 3. Diffusion term ð7þ F X G = X X G (F j j X + F h h X ) = G (F jj j X + F jh h X )j X + F j j XX + G (F hj j X + F hh h X )h X + F h h XX = G (F jj j 2 X + F hhh 2 X + 2F jhh X j X + F j j XX + F h h XX ) = G F jj Yh 2 + F hhyj 2 2F! jhy h Y j + F J 2 j j XX + F h h XX ð8þ F Y G = Y Y G (F j j Y + F h h Y ) = G (F jj j Y + F jh h Y )j Y + f j j YY + G (F hj j Y + F hh h Y )h Y + F h h YY = G (F jj j 2 Y + F hhh 2 Y + 2F jhh Y j Y + F j j YY + F h h YY ) = G F jj Xh 2 + F hhxj 2 2F! hjx h X j + F J 2 j j YY + F h h YY X Combining equations (8) and (9), we can obtain G F + Y F G Y ð9þ X h i = G J F jj(x 2 2 h + Y h 2 )+F hh(xj 2 + Y j 2 ) 2F jh(x j X h + Y j Y h ) + G F j ðj XX + j YY Þ+ F h ðh XX + h YY Þ = G J 2 (a F jj 2b F jh + g F hh )+0 = G J 2 (a F jj b F jh )+ G J 2 (g F hh b F jh ) 1 G J j J (a F j b F h ) G J h J (g F h b F j ) ð10þ where a = X 2 h + Y 2 h, b = X j X h + Y j Y h, and g = X 2 j + Y 2 j. 4. Source term R(x, y)=s(j, h) ð11þ Substituting equations (4), (7), (10), and (11) into equation (3), the governing equation in the computational plane j h can be obtained as follows F t + 1 (U F) J j 1 G J j + 1 J 1 h + 1 (V F) J h J (a F j b F h ) G J (g F h b F j ) + S(j, h) ð12þ
Natural 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 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 informationAbstract. Introduction
Combined forced and natural convection in a square cavity - numerical solution and scale analysis A.T. Franco/ M.M. Ganzarolli'' "DAMEC, CEFET, PR 80230-901, Curitiba, PR Brasil >>DE, FEM, UNICAMP 13081-970,
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 informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
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 informationNumerical study of heat transfer enhancement by rectangular winglet vortex generator pair in a channel
Special Issue Article Numerical study of heat transfer enhancement by rectangular winglet vortex generator pair in a channel Advances in Mechanical Engineering 2016, Vol. 8(5) 1 11 Ó The Author(s) 2016
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 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 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 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 information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
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 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 informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
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 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 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 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 informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationEFFECT OF THE INLET OPENING ON MIXED CONVECTION INSIDE A 3-D VENTILATED CAVITY
THERMAL SCIENCE: Year 2018, Vol. 22, No. 6A, pp. 2413-2424 2413 EFFECT OF THE INLET OPENING ON MIXED CONVECTION INSIDE A 3-D VENTILATED CAVITY by Hicham DOGHMI *, Btissam ABOURIDA, Lahoucin BELARCHE, Mohamed
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 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 informationUNSTEADY MIXED CONVECTION IN A POROUS MEDIA FILLED LID-DRIVEN CAVITY HEATED BY A SEMI-CIRCULAR HEATERS
THERMAL SCIENCE, Year 2015, Vol. 19, No. 5, pp. 1761-1768 1761 UNSTEADY MIXED CONVECTION IN A POROUS MEDIA FILLED LID-DRIVEN CAVITY HEATED BY A SEMI-CIRCULAR HEATERS by Md Mustafizur RAHMAN a,c, Hakan
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 informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 90 (014 ) 517 53 10th International Conference on Mechanical Engineering, ICME 013 Numerical simulation on buoyancy-driven heat
More informationAvailable online at ScienceDirect. Procedia Engineering 105 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 15 (215 ) 176 183 6th BSME International Conference on Thermal Engineering (ICTE 214) Natural Convection Flow and Heat Transfer
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
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 informationA STUDY ON NATURAL CONVECTION HEAT TRANSFER IN COMPLEX BOUNDARIES
http://doi.org/10.4867/jpe-017-01-11 JPE (017) Vol.0 (1) Mohapatra, C. R. Preliminary Note A STUDY ON NATURAL CONVECTION HEAT TRANSFER IN COMPLEX BOUNDARIES Received: 3 February 017 / Accepted: 01 April
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 informationAdvanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell
Laminar external natural convection on vertical and horizontal flat plates, over horizontal and vertical cylinders and sphere, as well as plumes, wakes and other types of free flow will be discussed in
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 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 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 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 informationNOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7 Number Pages 3 39 c Institute for Scientific Computing and Information NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL
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 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 informationAli J. Chamkha* Ahmed Ali Shaker
Progress in Computational Fluid Dynamics, Vol. 12, No. 5, 2012 309 Conduction-combined forced and natural convection in a lid-driven parallelogram-shaped enclosure divided by a solid partition Ali J. Chamkha*
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES
Journal of Mathematics and Statistics 9 (4): 339-348, 2013 ISSN: 1549-3644 2013 doi:10.3844/jmssp.2013.339.348 Published Online 9 (4) 2013 (http://www.thescipub.com/jmss.toc) ENERGY PERFORMANCE IMPROVEMENT,
More informationConvection Workshop. Academic Resource Center
Convection Workshop Academic Resource Center Presentation Outline Understanding the concepts Correlations External Convection (Chapter 7) Internal Convection (Chapter 8) Free Convection (Chapter 9) Solving
More informationAdaptability analysis of radiative transport diffusion approximation in planar-graded-index media
Research Article Adaptability analysis of radiative transport diffusion approximation in planar-graded-index media Advances in Mechanical Engineering 2018, ol. 10(11) 1 6 Ó The Author(s) 2018 DOI: 10.1177/1687814018809613
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationStudy on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow
Research Article Study on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow Advances in Mechanical Engineering 1 8 Ó The Author(s) 2015 DOI: 10.1177/1687814015573829
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 informationFREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER
Free Convective Heat Transfer From an Object at Low Rayleigh Number FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER Md. Golam Kader and Khandkar Aftab Hossain * Department of Mechanical
More informationNUMERICAL STUDY OF MIXED CONVECTION AND THERMAL RADIATION IN A SQUARE CAVITY WITH AN INSIDE INCLINED HEATER
NUMERICAL STUDY OF MIXED CONVECTION AND THERMAL RADIATION IN A SQUARE CAVITY WITH AN INSIDE INCLINED HEATER N. HAMICI a, D. SADAOUI a a. Laboratory of Mechanic, Materials and Energy (L2ME), University
More informationMIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER
Suranaree J. Sci. Technol. Vol. 17 No. 2; April - June 2010 139 MIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER Md. Mustafizur Rahman 1 *, M. A. Alim 1 and Sumon Saha
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 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 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 informationThe Natural Convective Heat Transfer in Rectangular Enclosure Containing Two Inclined Partitions
AASCIT Journal of Energy 2015; 2(3): 36-43 Published online June 10, 2015 (http://www.aascit.org/journal/energy) The Natural Convective Heat Transfer in Rectangular Enclosure Containing Two Inclined Partitions
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 informationMIXED CONVECTION HEAT TRANSFER OF NANOFLUIDS IN A LID DRIVEN SQUARE CAVITY: A PARAMETRIC STUDY
International Journal of Mechanical and Materials Engineering (IJMME), Vol. 8 (2013), No. 1, Pages: 48-57. MIXED CONVECTION HEAT TRANSFER OF NANOFLUIDS IN A LID DRIVEN SQUARE CAVITY: A PARAMETRIC STUDY
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 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 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 informationNumerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles
Numerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles O. A. SHAHRUL Faculty of Mechanical & Manufacturing Engineering Universiti
More informationThe Influence of the Angle of Inclination on Laminar Natural Convection in a Square Enclosure
The Influence of the Angle of Inclination on Laminar Natural Convection in a Square Enclosure Kumar Dookhitram 1* Roddy Lollchund 2 1. Department of Applied Mathematical Sciences, School of Innovative
More informationIntroduction to Heat and Mass Transfer. Week 14
Introduction to Heat and Mass Transfer Week 14 HW # 7 prob. 2 Hot water at 50C flows through a steel pipe (thermal conductivity 14 W/m-K) of 100 mm outside diameter and 8 mm wall thickness. During winter,
More informationFree Convective Heat Transfer From A Vertical Surface For The Case Of Linearly Varying Thermal Potential
American Journal of Engineering Research (AJER) e-issn : 232-847 p-issn : 232-936 Volume-2, Issue-9, pp-71-75 www.ajer.org Research Paper Open Access Free Convective Heat Transfer From A Vertical Surface
More informationNUMERICAL ANALYSIS OF NATURAL CONVECTION IN A RIGHT- ANGLED TRIANGULAR ENCLOSURE
Frontiers in Heat and Mass Transfer Available at www.thermalfluidscentral.org NUMERICAL ANALYSIS OF NATURAL CONVECTION IN A RIGHT- ANGLED TRIANGULAR ENCLOSURE Manoj Kr. Triveni *, Dipak Sen, RajSekhar
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 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 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 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 informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
More informationFinite Element Analysis of MHD Natural Convection in a Rectangular Cavity and Partially Heated Wall
Engineering and Applied Sciences 2017; 2(3): 53-58 http://www.sciencepublishinggroup.com/j/eas doi: 10.11648/j.eas.20170203.12 Finite Element Analysis of MHD Natural Convection in a Rectangular Cavity
More informationLaminar Mixed Convection in the Entrance Region of Horizontal Quarter Circle Ducts
Proceedings of the 5th IASME/WSEAS Int. Conference on Heat Transfer Thermal Engineering and Environment Athens Greece August 5-7 007 49 Laminar Mixed Convection in the Entrance Region of Horizontal Quarter
More informationMeysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI
THERMAL SCIENCE: Year 014, Vol. 18, No., pp. 479-49 479 NUMERICAL INVESTIGATION OF ENTROPY GENERATION IN LAMINAR FORCED CONVECTION FLOW OVER INCLINED BACKWARD AND FORWARD FACING STEPS IN A DUCT UNDER BLEEDING
More informationNUMERICAL STUDY OF THE THERMOSOLUTAL CONVECTION IN A THREE- DIMENSIONAL CAVITY SUBMITTED TO CROSS GRADIENTS OF TEMPERATURE AND CONCENTRATION
NUMERICAL STUDY OF THE THERMOSOLUTAL CONVECTION IN A THREE- DIMENSIONAL CAVITY SUBMITTED TO CROSS GRADIENTS OF TEMPERATURE AND CONCENTRATION Meriem OUZAOUIT 1 *, Btissam ABOURIDA 1, Lahoucine BELARCHE
More informationTHREE-DIMENSIONAL MIXED CONVECTION HEAT TRANSFER IN A PARTIALLY HEATED VENTILATED CAVITY. Corresponding author;
THREE-DIMENSIONAL MIXED CONVECTION HEAT TRANSFER IN A PARTIALLY HEATED VENTILATED CAVITY Hicham DOGHMI 1 *, Btissam ABOURIDA 1, Lahoucin BELARCHE 1, Mohamed SANNAD 1, Meriem OUZAOUIT 1 1 National School
More informationA new computational method for threaded connection stiffness
Research Article A new computational method for threaded connection stiffness Advances in Mechanical Engineering 2016, Vol. 8(12) 1 9 Ó The Author(s) 2016 DOI: 10.1177/1687814016682653 aime.sagepub.com
More informationHeat transfer and flow resistance performance of shutter baffle heat exchanger with triangle tube layout in shell side
Research Article Heat transfer and flow resistance performance of shutter baffle heat exchanger with triangle tube layout in shell side Advances in Mechanical Engineering 2016, Vol. 8(3) 1 8 Ó The Author(s)
More informationConvection in Three-Dimensional Separated and Attached Flow
Convection in Three-Dimensional Separated and Attached Flow B. F. Armaly Convection Heat Transfer Laboratory Department of Mechanical and Aerospace Engineering, and Engineering Mechanics University of
More informationTHE EFFECT OF VARIATION OF BASE FLUID ON NATURAL CONVECTION IN CAVITY FILLED WITH NAOFLUID IN THE PRESENCE OF MAGNETIC FIELD
THE EFFECT OF VARIATION OF BASE FLUID ON NATURAL CONVECTION IN CAVITY FILLED WITH NAOFLUID IN THE PRESENCE OF MAGNETIC FIELD Salma H. RADWAN 1*, Mohamed TEAMAH 2, Mohamed M. ABO ELAZM 3, Wael El-MAGHLANY
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 information3D Numerical Study on Laminar Forced Convection in V-Baffled Square Channel
American Journal of Applied Sciences 10 (10): 1287-1297, 2013 ISSN: 1546-9239 2013 Boonloi and Jedsadaratanachai, This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0
More informationTHERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW
THERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW Luigi De Giorgi, Carlo Cima, Emilio Cafaro Dipartimento di Energetica, Politecnico di Torino, Torino, Italy Volfango Bertola School
More informationA Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer
Copyright c 2008 ICCES ICCES, vol.6, no.1, pp.13-18 A Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer K agamani Devi 1,D.W.Pepper 2 Summary Over the past few years, efforts have
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 informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More informationFREE CONVECTION LIMIT TO HEAT TRANSFER FROM HEAT SINKS WITH SQUARE FINS ON A HORIZONTAL BASE Comparison with Cylindrical Fins
THERMAL SCIENCE, Year 2015, Vol. 19, No. 5, pp. 1565-1574 1565 FREE CONVECTION LIMIT TO HEAT TRANSFER FROM HEAT SINKS WITH SQUARE FINS ON A HORIZONTAL BASE Comparison with Cylindrical Fins by Subhas Chandra
More informationº Õߪí À À Ë ËÕπ â«ß Õ µ «â«µ Àπàß ÿ æ Ë Ëß øíß å π µ «Õ ªìπøíß å π π ß»
««æ π. ªï Ë 31 Ë 4 µÿ - π«2551 643 º Õߪí À À Ë ËÕπ â«ß Õ µ «â«µ Àπàß ÿ æ Ë Ëß øíß å π µ «Õ ªìπøíß å π π ß» µ π» «1 À µ å»ÿπ å ß µ ÕßÀ «ß ª ÿ π 12121 ËÕ 1 惻 π 2550 µõ ËÕ 20 ÿπ π 2551 àõ «µ Àπàß ÿ æ Ë
More informationEffects of Conduction on Magneto-Hydrodynamics Convection Flow in Triangular Enclosures
Effects of Conduction on Magneto-Hydrodynamics Mixed Convection Flow in Triangular Enclosures 1 Muhammad Sajjad Hossain*, M. A. Hoque and 3 Md. Kamruzzaman 1 Assistant Professor of Mathematics, Lecturer
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 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 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 informationEffects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
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 informationSTUDY OF CONVECTION HEAT TRANSFER ENHANCEMENT INSIDE LID DRIVEN CAVITY UTILIZING FINS AND NANOFLUID
STUDY OF CONVECTION HEAT TRANSFER ENHANCEMENT INSIDE LID DRIVEN CAVITY UTILIZING FINS AND NANOFLUID Arash Mirabdolah Lavasani a*, Mousa Farhadi b, A.A Rabienataj Darzi c a Department of Mechanical Engineering,
More informationHEAT TRANSFER BY CONVECTION. Dr. Şaziye Balku 1
HEAT TRANSFER BY CONVECTION Dr. Şaziye Balku 1 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 the
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 informationAn underrelaxation control method
The current issue and full text archive of this journal is available at www.emeraldinsight.com/064-4401.htm An under-relaxation for accelerating the iteration convergence of flow field simulation C.H.
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 informationThe Study of Natural Convection Heat Transfer in a Partially Porous Cavity Based on LBM
Send Orders for Reprints to reprints@benthamscience.ae 88 The Open Fuels & Energy Science Journal, 014, 7, 88-93 Open Access The Study of Natural Convection Heat Transfer in a Partially Porous Cavity Based
More information¼ outer radial distance of the outer. r o,oc
The current issue and full text archive of this journal is available at http://www.emeraldinsight.com/0264-4401.htm Numerical and heat transfer in multi-channel, narrow-gap fuel element Huang Jun, Wang
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 information