• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 2002년도 추계학술발표회요약집
• /
• pp.167-167
• /
• 2002

• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• 대한기계학회논문집B
• /
• 제27권9호
• /
• pp.1262-1272
• /
• 2003

A two-dimensional direct numerical simulation is performed to investigate the dynamic behaviors of a single vortex in counter reacting and non-reacting flow field. A predictor-corrector-type numerical scheme with a low Mach number approximation is used in this simulation. A 16-step augmented reduced mechanism is adopted to treat the chemical reaction. The budget of the vorticity transport equation is examined to reveal a mechanism leading to the formation, destruction and transport of a single vortex according to the direction of vortex generation in reacting and non-reacting flows. The results show that air-side vortex has more larger strength than that of fuel-side vortex in both non-reacting and reacting flows. In reacting flow, the vortex is more dissipated than that in non-reacting flow as the vortex approach the flame. The total circulation in reacting flow, however, is larger than that in non-reacting flow because the convection transport of vorticity becomes much large by the increased velocity near the flame region. It is also found that the stretching and the convection terms mainly generate vorticity in non-reacting and reacting flows. The baroclinic torque term generates vorticity, while the viscous and the volumetric expansion terms attenuate vorticity in reacting flow. Furthermore, the contribution of volumetric expansion term on total circulation for air-side vortex is much larger than that of fuel-side vortex. It is also estimated that the difference of total circulation near stagnation plane according to the direction of vortex generation mainly attributes to the convection term.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 춘계 학술대회논문집
• /
• pp.33-38
• /
• 2000

A Numerical study of laminar vortex-shedding past a circular cylinder has been performed widely by many researchers. Many factors, such as numerical technique and domain size, number and shape of grid, affected predicting vortex shedding and Strouhal number. In the present study, the effect of convection scheme, time discretization methods and grid dependence were investigated. The present paper presents the finite volume solution of unsteady flow past circular cylinder at Re=200, 400. The Strouhal number was predicted using UDS, CDS, Hybrid, Power-law, LUDS, QUICK scheme for convection term, implicit and crank-nicolson methods for time discretization. The grid dependence was investigated using H-type mesh and O-type mesh. It also studied that the effect of mesh size of the nearest adjacent grid of circular cylinder. The effect of convection scheme is greater than the effect of time discretization on predicting Strouhal. It has been found that the predicted Strouhal number changed with mesh size and shape.

• 동력기계공학회지
• /
• 제13권5호
• /
• pp.11-17
• /
• 2009

A numerical study of unsteady mixed convection in a cavity with high viscous fluid is presented. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Rayleigh number ($10^3\;{\leq}\;Ra\;{\leq}\;10^5$), the Reynolds number (0 < Re $\leq$ 1), and the aspect ratio (0.5 $\leq$ AR $\leq$ 2). The fluid used is silicon oil, a high prandtl number fluid, Pr = 909.1. The results show velocity vectors and temperature distributions. It is found that the periodic flows in a cavity are observed at very low Reynolds numbers, and the period of periodic flow decreases with increasing Reynolds and Rayleigh numbers, and increases with increasing aspect ratio. Also, the Reynolds number range of periodic flow increases with increasing Rayleigh numbers and aspect ratio.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1997년도 추계 학술대회논문집
• /
• pp.49-54
• /
• 1997

Rice's monotone streamline upwind finite element method, which was proposed to treat convection-dominated flows, is applied to the linear triangular element. An alignment technique of unstructured grids with given velocity fields is used to prevent the interpolation error produced in evaluating the convection term in the upwind method. The alignment of grids is accomplished by optimizing a target function defined with the inner-product of a properly chosen side vector in the element with the velocity field. Two pure advection problems are considered to demonstrate the superiorities of the present approach in solving the convection-dominated flow on the unstructured grid. Solutions obtained with aligned grids are much closer to the exact solutions than those with initial regular grids. The capability of the present approach in predicting the appearance of the secondary vortex in the laminar confined jet impingement is shown by comparing streamlines to those produced by SIMPLE on a highly stretched grid toward the impingement plate.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제16권3호
• /
• pp.193-204
• /
• 2012

This paper is comparing numerical schemes for a differential equation with convection and fourth-order diffusion. Our model equation is $h_t+(h^2-h^3)_x=-(h^3h_{xxx})_x$, which arises in the context of thin film flow driven the competing effects of an induced surface tension gradient and gravity. These films arise in thin coating flows and are of great technical and scientific interest. Here we focus on the several numerical methods to apply the model equation and the comparison and analysis of the numerical results. The convection terms are treated with well known WENO methods and the diffusion term is treated implicitly. The diffusion and convection schemes are combined using a fractional step-splitting method.

• Journal of applied mathematics & informatics
• /
• 제26권5_6호
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• 대한기계학회논문집
• /
• 제11권3호
• /
• pp.528-536
• /
• 1987

본 연구에서는 간헐적인 유동에 의거한 체적 대류 모델을 설정하고, 구배 확 산 모델에 유동장의 중심부와 외부에서 서로 다른 무게값을 부여하는 혼성 확산 모델 (hybrid diffusion model)을 제안하며, 이 모델을 검증하기 위한 첫 단계로서 평면 제 트 유동에 대하여 수치 계산을 수행하다. 여기에는 간헐도에 대한 난류 전달 방정식 이 필요한데 이 방정식의 생성항(production term)은 곧 외부의 비회전 유동이 난류성 유체로 유입되는 정도를 의미하게 된다.

• 한국전산유체공학회지
• /
• 제19권1호
• /
• pp.15-20
• /
• 2014

Accurate simulation of free-surface wave flows around a ship is very important for better hull-form design. In this paper, a computational fluid dynamics (CFD) code, termed SNUFOAM, which is based on the open source libraries, OpenFOAM, was developed to predict the wave patterns around a ship. Additional anti-diffusion source term for minimizing a numerical diffusion, which was caused by convection differencing scheme, was considered in the volume-fraction transport equation. The influence of the anti-diffusion source term was tested by applying it to free-surface wave flow around the Wigley model ship. In results, the band width of the volume fraction contours between 0.1 to 0.9 at the hull surface was narrowed by considering the anti-diffusion term.