• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2000년도 제20회 KOSCO SYMPOSIUM 논문집
• /
• pp.117-121
• /
• 2000

We developed a general purpose program for the analysis of flows in a gas turbine combustor. The program uses non-staggered grids based on finite volume method and the cartesian velocities as primitive variables. We calculated a flow inside the C-type diffuser to check the boundary fitted coordinate. The velocity profiles at cross section agree well with experimental results. We calculated turbulent diffusion flame behind a bluff body for the combustion simulation. Simulation shows two recirculating region like experimental results. Simulated velocity, turbulent kinetic energy, temperature and concentration distribution agree well with experimental data. Finally, simulation of axisymmetric flows with swirl shows two recirculating region like experimental results.

• 한국연소학회지
• /
• 제5권1호
• /
• pp.55-66
• /
• 2000

A general purpose program for the analysis of flows in a gas turbine combustor is developed. The program uses non-staggered grids based on finite volume method and the cartesian velocities as primitive variables. A flow inside the C-type diffuser is simulated to check the boundary fitted coordinate. The velocity profiles at cross section agree well with experimental results. A turbulent diffusion flame behind a bluff body is simulated for the combustion simulation. Simulated results show good agreement with experimental data. Finally, a turbulent flow with swirl in a gas turbine combustor was simulated. The results show two recirculating region and simulated velocity fields agree well with experimental data. The distance between two recirculating regions becomes shorter as swirl angle increases. Swirl angle changes angular momentum and streamlines in flow fields.

• 한국전산유체공학회지
• /
• 제12권4호
• /
• pp.44-53
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed for a 3D component of a nuclear system code and a component-scale analysis tool. A two-fluid three-field model is used for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. An unstructured grid is adopted for realistic simulations of the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been applied to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the modified numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing very well.

• 한국안전학회지
• /
• 제15권3호
• /
• pp.14-22
• /
• 2000

A numerical analysis has been performed on the two-dimensional rectangular gallium melting problem using the enthalpy method. The major advantage of this method is that the physical domain is discretized with fixed grids without transforming variables and the interface conditions of phase change are accounted for the definition of suitable source terms in the governing equations. But in the fixed method, there is some ambiguity in defining the porosity constant which has no physical interpretation. If the velocity correction is included in the momentum equation, for the appropriate range of porosity constant, the realistic predictions are obtained. The object of the present work is to predict the phase change within the molten steel with thin riser slab using the modified enthalpy-porosity method. The computational procedures for predicting velocity and temperature are based on the finite volume method and the non-staggered grid system. The influence of natural convection on the melting process is considered. A comparison with the experimental results shows that the modified method is better than the previous one.

• Journal of Ship and Ocean Technology
• /
• 제5권1호
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• 한국수자원학회논문집
• /
• 제44권6호
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 한국추진공학회지
• /
• 제1권2호
• /
• pp.22-31
• /
• 1997

비직교 일반좌표계를 사용하는 범용 연소유동해석 프로그램을 개발하여 재순환영역을 가지는 연소유동 등을 해석하였다. 본 프로그램은 유한체적법에 근간을 둔 비엇갈림격자계를 사용하며 직교좌표 속도성분을 종속변수로 하였다. 연소모델은 무한반응속도의 1단계 비가역 반응을 고려하였으며, 밀도가중평균된 지배방정식을 고려하였다. Conserved Scalar의 확률밀도함수를 clipped gauss 분포로 가정하였다 본 프로그램을 사용하여 재순환영역을 가지는 난류확산화염을 계산하였다. 계산결과는 실험결과와 마찬가지로 두개의 재순환영역을 보였으며, 속도, 난류운동에너지, 온도 및 농도 등이 실험결과와 양호하게 일치하는 결과를 얻을 수 있었다.