• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제19권2호
• /
• pp.103-121
• /
• 2015

In this paper, we briefly review and describe a projection algorithm for numerically computing the two-dimensional time-dependent incompressible Navier-Stokes equation. The projection method, which was originally introduced by Alexandre Chorin [A.J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comput., 22 (1968), pp. 745-762], is an effective numerical method for solving time-dependent incompressible fluid flow problems. The key advantage of the projection method is that we do not compute the momentum and the continuity equations at the same time, which is computationally difficult and costly. In the projection method, we compute an intermediate velocity vector field that is then projected onto divergence-free fields to recover the divergence-free velocity. Numerical solutions for flows inside a driven cavity are presented. We also provide the source code for the programs so that interested readers can modify the programs and adapt them for their own purposes.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.1287-1292
• /
• 2007

Recently, the frequency of unexpecting heavy rains has been increased due to abnormal climate and extreme rainfall. There was a limit to analyze 1D or 2D stream flow that was applied simple momentum equation and fixed energy conservation. Therefore, hydrodynamics flow analysis in rivers has been needed 3D numerical analysis for correct stream flow interpretation. In this study, CFD model on FLOW-3D was applied to stream flow analysis, which solves three dimenson RANS(Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behavior and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as ${\kappa}-{\varepsilon}$, RNG ${\kappa}-{\varepsilon}$ and LES. Those numerical analysis results have been illustrated by the turbulence energy effects, velocity of flow distributions, water level pressure distributions and eddy flows around the piers at Jangwall bridge in urbarn stream.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회B
• /
• pp.2460-2464
• /
• 2008

The present paper suggests new method to know the effects of molecular diffusion and the helicity of microchannel flows on mixing in passive micromixers, which are essential components of a microfluidic chip. In this study, 'Helix Index' is newly defined as the magnitude of chaotic advection. Relationship between Helix Index and Mixing Index is analyzed numerically such as the wide range of Peclet and Reynolds numbers in three dimensional serpentine microchannel when using soluble solutions (water/glycerol). As a result, a simple algebraic equation is derived by this relationship based on a regression analysis. The algebraic equation is found to be able to accurately predict the mixing performance without solving the coupled, complex momentum and mass transfer equations.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제4권3호
• /
• pp.256-266
• /
• 2012

Cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, torpedoes, etc. The present work has developed the base code to solve the cavitating flows past the axisymmetric bodies with several forebody shapes. The governing equation is the Navier-Stokes equation based on homogeneous mixture model. The momentum is in the mixture phase while the continuity equation is solved in liquid and vapor phase, separately. The solver employs an implicit preconditioning algorithm in curvilinear coordinates. The computations have been carried out for the cylinders with hemispherical, 1-caliber, and 0-caliber forebody and, then, compared with experiments and other numerical results. Fairly good agreements with experiments and numerical results have been achieved. It has been concluded that the present numerical code has successfully accounted for the cavitating flows past axisymmetric bodies. The present code has also shown the capability to simulate ventilated cavitation.

• Coupled systems mechanics
• /
• 제8권2호
• /
• pp.147-168
• /
• 2019

This work aims to simulate three-dimensional heavy oil flow in a reservoir with heater-wells. Mass, momentum and energy balances, as well as correlations for rock and fluid properties, are used to obtain non-linear partial differential equations for the fluid pressure and temperature, and for the rock temperature. Heat transfer is simulated using a two-equation model that is more appropriate when fluid and rock have very different thermal properties, and we also perform comparisons between one- and two-equation models. The governing equations are discretized using the Finite Volume Method. For the numerical solution, we apply a linearization and an operator splitting. As a consequence, three algebraic subsystems of linearized equations are solved using the Conjugate Gradient Method. The results obtained show the suitability of the numerical method and the technical feasibility of heating the reservoir with static equipment.

• 대한설비공학회지:설비저널
• /
• 제14권4호
• /
• pp.293-299
• /
• 1985

This paper describes an economical prediction procedure for heat transfer phenomenon through a channel containing an abrupt asymmetric expansion in flow cross-seetional area. Numerical solutions for the flow field are obtained by the finite difference numerical method applied to the modified boundary layer equations. Modified boundary energy equation is used to analyze heat transfer as modified boundary momentum equation. Predictions of the method compare very favorable with exprimental data. Results of this study by modified boundary layer equation are as follows : 1. The computation time required for the scheme is at least an order of magnitude less than for the numerical solution of the full Navier-stokes and Energy eguations. 2. In laminar flow, the maximum heat transfer occurs downstream of the reattachment point.

• 대한기계학회논문집B
• /
• 제23권2호
• /
• pp.243-253
• /
• 1999

The dual reciprocity boundary element method (DRBEM) is used to solve the Graetz problem of laminar flow inside circular duct. In this method the domain integral tenn of boundary integral equation resulting from source term of governing equation is transformed into equivalent boundary-only integrals by using the radial basis interpolation function, and therefore complicate domain discretization procedure Is completely removed. Velocity profile is obtained by solving the momentum equation first and then, using this velocities as Input data, energy equation Is solved to get the temperature profile by advancing from duct entrance through the axial direction marching scheme. DRBEM solution is tested for the uniform temperature and heat flux boundary condition cases. Local Nusselt number, mixed mean temperature and temperature profile inside duct at each dimensionless axial location are obtained and compared with exact solutions for the accuracy test Solutions arc in good agreement at the entry region as well as fully developed region of circular duct, and their accuracy are verified from error analysis.

• 한국전기전자재료학회:학술대회논문집
• /
• 한국전기전자재료학회 2001년도 추계학술대회 논문집 Vol.14 No.1
• /
• pp.141-144
• /
• 2001

An accurate cross sections set are necessary for the quantitatively understanding and modeling of plasma phenomena. By using the electron swarm method. we determine an accurate electron cross sections set for objective atoms or molecule at low electron energy range. In previous paper, we calculated the electron transport coefficients in pure $CF_4$ molecular gas by using two-term approximation of the Boltzmann equation. And by using this simulation method. we confirmed erroneous calculated results of transport coefficients for $CF_{4}$ molecule treated in this paper having 'C2v symmetry' as $C_{3}H_{8}$ and $C_{3}F_{8}$ which have large vibrational excitation cross sections which may exceed elastic momentum transfer cross section. Therefore, in this paper, we calculated the electron transport coefficients(W and $ND_L$) in pure $CF_4$ gas by using multi-term approximation of the Boltzmann equation by Robson and Ness which was developed at lames-Cook university, and discussed an application and/or validity of the calculation method by comparing the calculated results by two-term and multi-term approximation code.

• 한국전기전자재료학회:학술대회논문집
• /
• 한국전기전자재료학회 2001년도 추계학술대회 논문집
• /
• pp.141-144
• /
• 2001

An accurate cross sections set are necessary for the quantitatively understanding and modeling of plasma phenomena. By using the electron swarm method, we determine an accurate electron cross sections set for objective atoms or molecule at low electron energy range. In previous paper, we calculated the electron transport coefficients in pure CF$_4$ molecular gas by using two-term approximation of the Boltzmann equation. And by using this simulation method, we confirmed erroneous calculated results of transport coefficients for CF$_4$ molecule treated in this paper having 'C2v symmetry'as C$_3$H$_{8}$ and C$_3$F$_{8}$ which have large vibrational excitation cross sections which may exceed elastic momentum transfer cross section. Therefore, in this paper, we calculated the electron transport coefficients(W and ND$_{L}$) in pure CF$_4$ gas by using multi-term approximation of the Boltzmann equation by Robson and Ness which was developed at James-Cook university, and discussed an application and/or validity of the calculation method by comparing the calculated results by two-term and multi-term approximation code.e.

• Tribology and Lubricants
• /
• 제20권5호
• /
• pp.284-291
• /
• 2004

We present a three-dimensional average flow model considering elastic deformation of pad asperities for chemical mechanical planarization. To consider the contact deformation of pad asperities in the calculation of the flow factor, three-dimensional contact analysis of a semi-infinite solid based on the use of influence functions is conducted from computer generated three dimensional roughness data. The average Reynolds equation and the boundary condition of both force and momentum balance are used to investigate the effect of pad roughness and external pressure conditions on film thickness and wafer position angle.