• 한국해양공학회지
• /
• 제7권2호
• /
• pp.150-161
• /
• 1993

원초변수를 이용한 Navier-Stokes 방정식의 수치계산기법을 개발하고, 이를 응용하여 backward facing step의 층류 유동을 계산하였다. 직교좌표계에서의 비압축성 Navier-Stokes방정식을 풀기위해 시간과 공간항을 2차 정도의 유한 차분을 사용하여 이산화하였고 비교차격자계를 사용하여 양해법으로 수치 계산하였다. 운동량방정식과 연속방정식으로 부터 유도된 압력방정식(pressure-poisson equation)을 이용하여 무발산 조건을 만족시켰ㄲ다. Backward facing step의 층류 유동을 100.$\leq$R$_e$$\leq$1000 범위에 대해서 수치 계산하였으며 실험결과와 잘 일치하는 결과를 구할 수 있었다. 특히 step뒤에서 생기는 박리구간의 길이는 다른 계산결과들보다 실험치에 가까운 값을 얻을 수 있었으며, Re가 600보다 클때는 위쪽 벽에 또 다른 박리 유동이 발생되는 현상이 예측되었다.

• 한국정밀공학회지
• /
• 제25권4호
• /
• pp.100-107
• /
• 2008

Bellows is an important component in Joule-Thomson cryocooler, which minimize the excessive flow of the cryogenic gas. The bellows is made of Monel 400 and its geometry is an axial symmetric shell. During cool-down process, the pressure and volume within bellows must be satisfied with Benedict-Webb-Rubin state equation. Moreover, Poisson's ratio of Monel 400 is nearly constants, but its Young's modulus varies for a drop in temperature. Under these conditions, bellows contracts in the axial direction like a spring. To evaluate deformation of bellows at cryogenic temperature, the numerical calculation of the volume within bellows and finite element analysis are iteratively used in this research. the numerical results show that deformation of the bellows is approximately linear for change of temperature.

• Journal of Advanced Marine Engineering and Technology
• /
• 제21권5호
• /
• pp.557-563
• /
• 1997

The unsteady flow in three-dimensional cubic cavity with narrow channel at upper region is investigated experimentally for three kinds of Reynolds number, 1*10/sup 4/, 3*10/sup 4/ and 5*10/sup 4/ based on the cavity width and cavity inlet mean flow velocity. Instant velocity vectors are obtained simultaneously at whole field by PIV(Particle Image Velocimetry). Wall pressure distributions are estimated using Poisson equation from the velocity data. Results of PIV reveal that severe unsteady flow fluctuation within the cavity are remarkable at all Reynolds numbers and sheared mixing layer phenomena are also found at the region where inlet driving flow is collided with the clock-wise rotating main primary vortex. Instant velocity profiles reveal that deformed forced vortex formation is observed throughout the entire region and spanwise kinetic energy migration is conspicuous.

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

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.8-12
• /
• 1997

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• Journal of Advanced Marine Engineering and Technology
• /
• 제22권6호
• /
• pp.935-941
• /
• 1998

Experiments were carried out for a cubic cavity flow. Contrinuous shear stress is supplied by driven flow for high Reynolds number and three kinds of aspect ratios. Velocity vectors are obtained by PIV and they are used as velocity components for Poisson equation for pressure, Related boundary conditions and no-slip condition at solid wall and the linear velocity extrapolation on the upper side of cavity are well examined for the present study. For calculation of pressure resolution of grid is basically $40{\times}40$ and 2-dimensional uniform mesh using MSC staggered grid is adopted. The flow field within the cavity maintains a forced-vortex formation and almost of the shear stress from the driving inflow is transformed into rotating flow energy and the size of the distorted forced-vortex increases with increment of Reynolds number

• Journal of Advanced Marine Engineering and Technology
• /
• 제20권5호
• /
• pp.30-38
• /
• 1996

A PIV system consisting of an image board slit into personal computer and 2-D sheet light projector is developed and related techniques to improve its performance are discussed. A grey-level cross correlation method capable of overcoming particle seeding limitation is suggested. And a sub-pixed interpolation method in determining the vector terminal is preposed by considering 8-neighbours correlation distributing patterns. Furthermore, pressure estimation from the acquired velocity vectors by applying the Poisson equation is persented with reasonable feasibility. As a practical application of the present system, evaporator flows are analysed and attained instantaneous velocity vectors reveal that the flow phenomena maintain turbulent fluctuation.

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

The development of unsteady three-dimensional incompressible viscous solver based on unsteady physical curvilinear coordinate system is presented. A 12-point finite analytic scheme based on local uniform grid spacing is extended for nonuniform grid spacing. The formulation of a condition is suggested to avoid the oscillation of the series summations produced by the application of the method of separation of variables. SIMPLER and pressure Poisson equation techniques are used for solving a velocity-pressure coupled problem. The matrix is solved using the Generalized Minimal RESidual (GMRES) method to enhance the convergence rate of unsteady flow solver and the Kinematic boundary condition of a free surface flow. It is demonstrated that the numerical solutions of these equations are less mesh sensitive.

• Journal of Electrical Engineering and Technology
• /
• 제7권1호
• /
• pp.97-102
• /
• 2012

To predict accurately the breakdown voltage in air at atmospheric pressure, a fully coupled finite element analysis combining the hydrodynamic diffusion-drift equations with Poisson's equation is proposed in the current paper. As three kinds of charged transport particles are nonlinearly coupled with spatial electric fields, the equations should be solved by an iterative numerical scheme, in which secondary effects, such as photoemission and photoionization, are considered. The proposed method has been successfully applied to evaluate the breakdown voltage in circular parallel-plane electrodes. Its validity has been proved through the comparison of the predicted and experimental results. The effects of numerical conditions of the initial charge, photoemission, and background ionization on the discharge phenomena are quantitatively assessed through Taguchi's design of experiment method.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.410-415
• /
• 2001

An analytical and numerical examination of second-order fractional-step methods and boundary condition for the incompressible Navier-Stokes equations is presented. In this study, the compatibility condition for pressure Poisson equation and its boundary conditions, stability, and numerical accuracy of canonical fractional-step methods has been investigated. It has been found that satisfaction of compatibility condition depends on tentative velocity and pressure boundary condition, and that the compatible boundary conditions for type D method and approximately compatible boundary conditions for type P method are proper for divergence-free velocity for type D and approximately divergence-free for type P method. Instability of canonical fractional-step methods is induced by approximation of implicit viscous term with explicit terms, and the stability criteria have been founded with simple model problems and numerical experiments of cavity flow and Taylor vortex flow. The numerical accuracy of canonical fractional-step methods with its consistent boundary conditions shows second-order accuracy except $D_{MM}$ condition, which make approximately first-order accuracy due to weak coupling of boundary conditions.