• 대한조선학회논문집
• /
• 제54권6호
• /
• pp.470-475
• /
• 2017

Since the intake air of gas turbine engine of marine purpose contains water particles, inertial separator for separating the air and water particles are provided. Saw type and wave type separator are now used to separate inflow water particle from the gas. In this paper, the design parameters of saw type separator are studied by numerical simulations. Using the commercial CFD program, Star-CCM+, Lagrangian-Eulerian method was used to perform the analysis of two phase flow of the mist in the air. This method solves Reynolds-Averaged Navier-Stokes equations in Eulerian framework for the continuous phase, while solves equation of motion for individual particles in Lagrangian framework. Lagrangian multiphase method was applied to monitor the particles of different sizes and shapes and to verify collision between particles by chasing particles. Water particles were injected through injectors located at the inlet of the separator and escape mode was used which assumes that the particles attached on the surface of inertial separator were removed from the simulation, effectively escaping the solution domain. Through the numerical computations with the inlet condition of constant water particle size in the wetness fraction of 85%, efficiency of eliminating the water particle and the pressure drop between the inlet and outlet were examined.

• 대한기계학회논문집B
• /
• 제38권5호
• /
• pp.431-436
• /
• 2014

컨덕턴스는 유동저항의 반대되는 개념으로써, 광범위하게 유동지표로써 현재 사용되고 있다. 하지만, 유동 컨덕턴스 연구들은 지금까지 매우 드물며, 압축 공기의 표준화 장치에 대한 체계적인 연구가 필요하다. 본 연구에서는 2차 방정식 난류모델을 이용한 압축성 N-S 방정식을 적용하여 수치해석적(CFD) 연구를 수행하였다. 본 CFD 결과는 기존의 실험 데이터를 통해 검증되었으며 공압 배관의 입구 및 출구에서의 컨덕턴스와 마찰 계수의 값은 유량을 평가하기 위해 사용되었다. 본 결과는 컨덕턴스가 공압 배관의 입구 및 출구에서의 압력 비율에 의존한다는 것을 보여준다.

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

본 연구에서는 2 차원 캐비티 유동에서 경계조건 부여 방법의 차이에 따른 유체유동의 이력이 최종적인 정상상태에 미치는 효과를 수치실험을 통하여 관찰하였다. 그 결과 유동의 Reynolds 수가 작은 경우는 유동장이 가지는 이력이 시간의 경과와 함께 소멸하여 하나의 정상유동 상태를 가지는데 반하여, 큰 Reynolds 수에서는 유동장 이력이 소멸하지 않고 유지됨으로 인하여 다수의 정상유동 형태가 나타남이 확인이 되었다. 엄밀한 이론적 해석은 현재 이루워지지 않았으나 유체 운동량의 증가에 따라 Navier-Stokes 방정식의 비선형성이 영향을 미치고 있다고 여겨진다.

• 대한기계학회논문집B
• /
• 제20권12호
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

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

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 대한기계학회논문집B
• /
• 제26권12호
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2010년도 제39회 하계학술대회 초록집
• /
• pp.74-74
• /
• 2010

The aqua-plasma is the non-thermal plasma in electrical conductive electrolyte by generates the vapor film layer on the immersed metal electrode surface. This plasma can generate the hydroxyl radical by dissociate the water molecule with the plasma electron. To develop the plasma discharge device for high efficiency in the hydroxyl radical generation, proper model for estimation of plasma power is necessary. In this work, the 1-D spherical model was developed, considering temperature dependence material constants. The relation between the plasma power and hydroxyl generation was also studied by the comparison between the optical emission intensity from the hydroxyl radical using monochromator and estimated plasma power. First, the thickness of vapor layer thickness was estimated using the Navier-Stokes fluid equation in order to calculate the discharge E-field inside vapor layer. Using the E-field magnitude and power balance on the plasma generation, it was possible to estimate the plasma power. The plasma power was assumed to uniformly fill the vapor layer and the temperature of vapor layer was fixed in the boiling temperature of electrolyte, 375K. In the experiment, the aqua-plasma was discharged in the saline by applied the voltage on the bipolar electrode. The range of applied voltage was 234 to 280V-rms in the frequency of 380 kHz. Two type electrodes were produced with two ${\Phi}0.2$ tungsten. The plasma power was estimated from the V-I signal from the two high voltage probes and current probe. The estimated plasma power agreed with the profile of emission intensity when the plasma discharged between the metal electrode and vapor layer surface. However, when the plasma discharged between the metal electrodes, the increasing rate of emission intensity was lower than the increase of plasma power. It implies that the surface reaction is more sufficient rather than the volume reaction in the radical generation, due to the high density of water molecule in the liquid.

• 대한기계학회논문집B
• /
• 제31권11호
• /
• pp.895-903
• /
• 2007

CMBT(Curved Moving Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the curved solid wall of moving obstacle in a flow field. In our research CMBT was used to perform LBM simulation of a flow over a moving circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of CMBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a fixed circular cylinder in a channel and the results were compared against the solution of Navier-Stokes equation with deforming mesh technique. The simulations were performed in a moderate range of reynolds number at each moving cylinder to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is ar Re=250 and the result is the same as the case of fixed cylinder. As the cylinder approaching to one wall, the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. As the velocity ratio increase the third vortex are generated by interacting with the 2nd vortexes developed on the upper and lower wall boundary layer. The resultant $C_d$ decrease as reynolds number increasing and the Cd approached to a value when Re>1000.

• 대한조선학회논문집
• /
• 제48권3호
• /
• pp.222-232
• /
• 2011

This article examines the scale effect of the flow characteristics, resistance and propulsion performance on a 317k VLCC. The turbulent flows around a ship in both towing and self-propulsion conditions are analyzed by solving the Reynolds-averaged Navier-Stokes equation together with the application of Reynolds stress turbulence model. The computations are carried out in both model- and full-scale. A double-body model is applied for the treatment of free surface. An asymmetric body-force propeller is used. The speed performances including resistance and propulsion factors are obtained from two kinds of methods. One is to analyze the computational results in model scale through the revised ITTC' 78 method. The other is directly to analyze the computational results in full scale. Based on the computational predictions, scale effects of the resistance and the self-propulsion factors including form factor, thrust deduction fraction, effective wake fraction and various efficiencies are investigated. Scale effects of the streamline pattern, hull pressure and local flow characteristics including x-constant sections, propeller and center plane, and transom region are also investigated. This study presents a useful tool to hull-form and propeller designers, and towing-tank experimenters to take the scale effect into consideration.