• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.380-388
• /
• 2010

A thermal-hydraulic code, named CUPID, has been developed for the analysis of transient two-phase flows in nuclear reactor components. A two-fluid three-field model was used for steam-water two-phase flows. To obtain numerical solutions, the finite volume method was applied over unstructured cell-centered meshes. In steam-water two-phase flows, a phase change, i.e., evaporation of condensation, results in a great change in the flow field because of substantial density difference between liquid and vapor phases. Thus, two-phase flows are very sensitive to the local pressure that determines the phase change. This in turn puts emphasis on the accurate evaluation of local pressure gradient. This paper presents a new numerical scheme to evaluate the pressure gradient at cell centers on unstructured meshes. The results of the new scheme for a simple test function a gravity-driven cavity, and a wall boiling two-phase flow are compared with those of the previous schemes in the cupid code.

• Journal of Advanced Marine Engineering and Technology
• /
• 제30권1호
• /
• pp.58-64
• /
• 2006

This work focuses on reducing the noise and vibration levels of an LPi fuel pump, which are generated from the dynamic motions of pump elements and non-uniform flow of fuel. The noise and vibration levels increase as the revolution speed of the cam goes up. The fuel pump consists of five cavity cells, plungers and diaphragms, which are driven by the cam. The optimal design of the cam profile is performed to decrease the accelerations of moving Parts and to obtain a smooth hydraulic force through a dynamic analysis of a cam-plunger mechanism. The cam-Plunger with a cavity is modeled as a 2 degrees of freedom system having non-linear contacts, the cam profile being represented in terms of Fourier series in order to determine the optimal shape of the cam. From the optimized cam Profile, the acceleration of the diaphragm is reduced in $78\%$, the hydraulic force becoming smoother in case that the hydraulic force is rapidly dropped.

• 해양환경안전학회:학술대회논문집
• /
• 해양환경안전학회 2005년도 추계학술대회지
• /
• pp.159-163
• /
• 2005

채널 내부의 비정상 흐름을 규명하기 위하여 실제 사용 중인 판형 열교환기 모델을 대상으로 가시화 실험과 입자영상유속계(PIV)를 사용한 계측을 수행하였다. 실험은 기복을 가지고 있는 채널의 높이와 내부유동의 부분 평균속도벡터에 따른 7가지 종류의 레이놀즈수를 적용하였고, 순간속도벡터분포와 유동특성을 고찰하였다. 실험에서 삼각형상 그루브 채널은 하부 채널과 채널 흐름에 의해서 받고 있는 전단응력의 그루브 흐름과의 관계에서 복합적인 흐름의 형태로 나타났다. 삼각형상 그루브와 채널 사이의 경계인 전단혼합층은 채널에서 난류강도가 상승하는 주흐름에 영향을 미쳤다.

• KIEE International Transactions on Electrophysics and Applications
• /
• 제4C권4호
• /
• pp.137-141
• /
• 2004

In this paper, a thermopneumatic PMDS (polydimethlysiloxane) micropump with nozzle/diffuser elements is presented. The micropump is composed of nozzle/diffuser elements as dynamic valves, an actuator consisting of a circular PDMS diaphragm and a Cr/Au heater on a glass substrate. Four PDMS layers are used for fabrication of an actuator chamber, actuator diaphragm by a spin coating process, spacer layer, and nozzle/diffuser by the SU-8 molding process. The radius and thickness of the actuator diaphragm is 2 mm and 30 ${\mu}{\textrm}{m}$, respectively. The length and the conical angle of the nozzle/diffuser elements are 3.5 mm and 20$^{\circ}$, respectively. The actuator diaphragm is driven by the air cavity pressure variation caused by ohmic heating and natural cooling. The flow rate of the micropump in the frequency domain is measured for various duty cycles of the square wave input voltage. When the square wave input voltage of 5 V DC is applied to the heater, the maximum flow rate of the micropump is 44.6 ${mu}ell$/min at 100 Hz with a duty ratio of 80% under the zero pressure difference.

• 한국전산유체공학회지
• /
• 제13권4호
• /
• pp.13-23
• /
• 2008

This paper is an extension of previous study[1] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite (serendipity) cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires four degrees-of-freedom at each element corners. Those degrees-of-freedom are the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational basis functions from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[2].

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.33-41
• /
• 2008

This paper is an extension of previous study[9] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite serendipity cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires at each element corners four degrees-of-freedom such as the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational bases from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[7].

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.33-41
• /
• 2008

This paper is an extension of previous study[9] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite serendipity cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires at each element corners four degrees-of-freedom such as the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational bases from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[7].

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2004년도 학술대회지
• /
• pp.216-221
• /
• 2004

In two dimensional incompressible flaws, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint and hence achieve an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flaw analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flaw showed the present procedure to be stable, very efficient and useful in flaw analyses in conjunction with hierarchical elements.

• 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.

• 대한수학회지
• /
• 제50권5호
• /
• pp.1129-1163
• /
• 2013

Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.