Numerical Analysis of Laminar Forced Convection in a Parallel Plate Channel Using Artificial Compression and Finite Volume Method

• Main Authors: Wang, Shu-Lung, 王述龍
• Other Authors: Lei, Hsien-Yu
• Format: Others
• Language: zh-TW
• Published: 2013
• Online Access: http://ndltd.ncl.edu.tw/handle/23791996515160399514

碩士 === 國立臺灣海洋大學 === 機械與機電工程學系 === 102 === Finite-volume, finite-difference, and finite-element are three numerical methods generally used in the simulation of engineering problems, in which the finite volume method (Finite Volume Method, FVM) is the most common approach in the field of thermal and fluid flow analysis. The area of solution is divided into many non-overlapping control volumes, each grid node is surrounded by a control volume, in the FVM. An integral process is performed for each control volume so that conservation laws (such as mass, momentum, and energy) could be satified within each control volume specified. Therefore, the FVM is applied in the discretization of the governing equations in this thesis. The artificial compressibility method proposed by Chorin is used in the FVM approach to obtain velocity and pressure in the forced flow field, and the temperature field is then solved by FVM by plugging in its velocity. The staggered grid is used in the numerical solutions written in FORTRAN. A steady two-dimensional laminar forced convection between a parallel plate channel with constant wall temperature is analyzed. Local friction coefficient and Nusselt Number are computed numerically. Hydrodynamic and thermal entry length are also obtained numerically. The results of fully developed flow and thermal entry cases are compared with those of analytic solutions. For the developing flow case, results are compared with those of empirical correlations and the numerical solution of the commercial software package Fluent. The applicability of the FORTRAN program developed in this thesis is justified through the comparisons for solving this type of forced convection problem. The optimization of grid size is studied and used to analyze the forced flow and convection of with higher Reynolds numbers up to 800. Velocity, pressure, temperature, entry length, and local friction coefficient are compared well with those of Fluent results. Different thermal boundary conditions (such as constant wall heat flux), three-dimensional rectangular or circular pipe flow could be studied in the future.