• 대한기계학회논문집B
• /
• 제28권7호
• /
• pp.809-817
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• Journal of applied mathematics & informatics
• /
• 제33권5_6호
• /
• pp.635-648
• /
• 2015

We consider a weakly coupled system of singularly perturbed convection-diffusion equations with discontinuous source term. The diffusion term of each equation is associated with a small positive parameter of different magnitude. Presence of discontinuity and different parameters creates boundary and weak interior layers that overlap and interact. A numerical method is constructed for this problem which involves an appropriate piecewise uniform Shishkin mesh. The numerical approximations are proved to converge to the continuous solutions uniformly with respect to the singular perturbation parameters. Numerical results are presented which illustrates the theoretical results.

• 천문학회지
• /
• 제18권2호
• /
• pp.86-99
• /
• 1985

A generalization of the original ${\ddot{O}}pik's$ cellular convection theory has been made to accomodate a rotating convective medium. With the use of the formulation, a set of rotating model envelopes of the sun and late type main sequence stars have been constructed under three different rotation periods. Their thermal structures are presented and characteristics of their convection are discussed in the context of stellar dynamo. In the present study it is noted that the rotational angular velocity increases in wards with depth, and its increase turns out to be about 6% at the bottom of the solar convection zone.

• 한국마린엔지니어링학회:학술대회논문집
• /
• 한국마린엔지니어링학회 2005년도 전기학술대회논문집
• /
• pp.210-215
• /
• 2005

A three-dimensional numerical simulation is performed to investigate on a low Reynolds number mixed convection in a horizontal rectangular channel with the upper part cooled and the lower part heated uniformly. The three-dimensional governing equations are solved using a finite volume method. For convective term, the central differencing scheme is used and for the pressure correction, the PISO algorithm is used. Solutions are obtained for A=4, Pr=0.72, 10, 909, the Reynolds number ranging from $2.1{\times}10^{-2}$ to $1.2{\times}10^{-1}$, the Rayleigh number is $3.5{\times}10^4$. It is found that vortex roll structures of mixed convection in horizontal rectangular channel can be classified into three roll structures which affected by Prandtl number and Reynolds number.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집B
• /
• pp.51-56
• /
• 2000

The Finite Element Solutions Is reported on solid-liquid phase change in porous media with natural convection including freezing. The model is based on volume averaged transport equations, while phase change is assumed to occur over a small temperature range. The FEM (Finite Element Method) algorithm used in this study is 3-step time-splitting method which requires much less execution time and computer storage the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the energy equation. For natural convection including melting and solidification the numerical results show reasonable agreement with FDM (Finite Difference Method) results.

• 동력기계공학회지
• /
• 제11권4호
• /
• pp.56-64
• /
• 2007

A numerical simulation is carried out mixed convection in horizontal channel with a heat source from below of rectangular cavity. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Reynolds number ($10^{-2}{\leq}Re{\leq}50$), the Rayleigh number ($10^3{\leq}Ra{\leq}2.06{\times}10^5$), the Prandtl number ($0.72{\leq}Pr{\leq}909$), the aspect ratio ($0.5{\leq}AR=W/H{\leq}2$) and the angle of inclination ($0^{\circ}{\theta}60^{\circ}$). Mean Nusselt number distributions were obtained and effect of Reynolds number, Rayleigh number and Prandtl number on mixed convection in the horizontal channel with rectangular cavity were investigated.

• 대한기계학회논문집B
• /
• 제25권5호
• /
• pp.619-626
• /
• 2001

The present experimental and numerical study investigate flow and natural convection heat transfer characteristics of a composed rectangular-parallelogrammic enclosure with a guide vane. The governing equations for the two-dimensional, laminar, natural convection process in an enclosure are discretized by the control volume approach which insures the conservative characteristics to be satisfied in the calculation domain, and solved by a modified SIMPLE algorithm. The momentum and energy equations are coupled through the buoyancy term. In this results of experimental study, the natural convection heat transfer characteristics was well coincided with conclusions of other earlier experimental researches and numerical analysis.

• 설비공학논문집
• /
• 제12권1호
• /
• pp.33-39
• /
• 2000

There are many under going researchs for the natural convection and fluid flow in rectangular enclosure. In this paper, the optimal model that is the most frequently used for the analysis of a turbulent natural convection in rectangular enclosure is suggested by comparing with the result of Cheesewright's experiment. As We can see the distribution of the velocity, temperature, and turbulent kinitic energy, ST model tends to exaggerate the result of the experiment. The LS model generates better experimental result than the ST and DA's. Therefore, it is resonable to adopt the LS model that contains explicit physical meanings of each term in eouation of turbulent kinitic energy.

• 한국전산유체공학회지
• /
• 제7권1호
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2005년도 제36회 하계학술대회 논문집 B
• /
• pp.939-941
• /
• 2005

This research presents a 3D multi-objective approach regarding both magnetic and thermal characteristics associated with design of C-core actuator. The adjoint variable topology sensitivity equations are derived using the continuum method for three dimension. The sensitivity is verified using the Finite Difference Method(FDM). Convection interpolation function is proposed for density method of topologies such that convection term can be taken into consideration for practical design in the process of the optimization.