• 대한조선학회논문집
• /
• 제41권2호
• /
• pp.79-89
• /
• 2004

The viscous free surface flow around KRISO container ship (KCS) is computed using the finite volume based multi block RANS code, WAVIS developed at KRISO. The free surface is captured with the Level-Set method and the realizable k-$\varepsilon$ model is employed for turbulence closure. The computations are carried out at model scale. For accurate free surface solution and its stable convergence the computations are performed with a suitable grid refinement around the free surface by applying an implicit discretization method based on a finite volume method to the Level-Set formulation. In all computational cases the numerical results agree well with experimental measurements.

• 대한조선학회논문집
• /
• 제42권6호
• /
• pp.583-592
• /
• 2005

The turbulent free surface flow around a self-propelled KRISO 138K LNG Carrier is numerically simulated using the finite volume based multi-block RANS code, WAVIS developed at HRISO. The realizable k-$\varepsilon$ turbulence model with a wail function is employed for the turbulence closure. The free surface is captured with the Level-Set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. In order to obtain an accurate free surface solution and stable convergence, the computations are executed with a proper fine grid refinement around the free surface and with an adoption of implicit discretization scheme for the Level-Set formulation. The computed velocity vectors at the several stations and wave patterns show a good agreement with the experimental results measured at the KRISO towing tank.

• 한국추진공학회지
• /
• 제6권2호
• /
• pp.75-84
• /
• 2002

고체 추진 로켓 내부 연소실의 비정상 유동을 수치적으로 해석하였다. 완전 보존식을 이용하여 2 차원 축-대칭 연소실 안의 연소 불안정을 해석하기 위한 수치 기법을 구성하였는데 비정상 유동을 해석하기 위한 수정된 $\kappa$-$\varepsilon$ 난류 모델이 사용되었다. 이산화한 지배 방정식은 연관된 경계 조건을 포함하여 dual time-stepping 방법으로 시적분 하였다. 정상 상태의 계산을 기반으로 연소실 내의 천이 압력파의 비정상 상태를 수치적으로 모사하기 위하여 압력 펄스 및 압력 변동을 연소실 상단에 부과하였다. 로켓 모터 연소실 내의 다양한 정상 상태 및 비정상 상태의 특성을 계산 및 해석하였다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.62-68
• /
• 2004

A comprehensive numerical analysis has been carried out for both non-reacting and reacting flows in a scramjet engine combustor with and without a cavity. The theoretical formulation treats the complete conservation equations of chemically reacting flows with finite-rate chemistry of hydrogen-air. Turbulence closure is achieved by means of a k-$\omega$ two-equation model. The governing equations are discretized using a MUSCL-type TVD scheme, and temporally integrated by a second-order accurate implicit scheme. Transverse injection of hydrogen is considered over a broad range of injection pressure. The corresponding equivalence ratio of the overall fuel/air mixture ranges from 0.167 to 0.50. The work features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous studies. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the .underlying physical mechanisms. Much of the flow unsteadiness is related not only to the cavity, but also to the intrinsic unsteadiness in the flow-field. The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The roles of the cavity, injection pressure, and heat release in determining the flow dynamics are examined systematically.

• 한국정밀공학회지
• /
• 제15권9호
• /
• pp.117-126
• /
• 1998

For investigating rubber pad sheet metal forming process, the rubber pad deformation characteristics as well as the contact problem of rubber pad-sheet metal has been analyzed. In this paper, the behavior of the rubber deformation is represented by hyper-elastic constitutive relations based on a generalized Mooney-Rivlin model. Finite element procedures for the two-dimensional responses, employing total Lagrangian formulations are implemented in an implicit form. The volumetric incompressibility condition of the rubber deformation is included in the formulation by using penalty method. The sheet metal is characterized by elasto-plastic material with strain hardening effect and analyzed by a commercial code. The contact procedure and interface program between rubber pad and sheet metal are implemented. Inflation experiment of circular rubber pad identifies the behaviour of the rubber pad deformation during the process. The various form dies and scaled down apparatus of the rubber-pad forming process are fabricated for simulating realistic forming process. The obtaining experimental data and FEM solutions were compared. The numerical solutions illustrate fair agreement with experimental results. The forming pressure distribution according to the dimensions of sheet metal and rubber pads, various rubber models and rubber material are also compared and discussed.

• 한국컴퓨터그래픽스학회논문지
• /
• 제4권1호
• /
• pp.79-88
• /
• 1998

현재까지 불과 같은 유체를 렌더링하기 위한 많은 모델들이 개발되었다. 특히 메타볼은 단순성과 유동성의 특정으로 인해 눈이나 물과 같은 유체의 곡면을 모델링하기 위한 기법으로 많이 사용되어졌다. 본 논문은 중력장내에서 변형되는 물방울을 모델링하기 위한 새로운 방법을 제안한다. 이전의 연구에서 물방울은 단순한 메타볼에 의해 표현되었다. 물방울에 중력이 고려되지 않았기 때문에 렌더링 결과는 실제적인 물방울의 형태를 표현할 수 없었다. 우리는 기존의 메타볼에 실세계의 중력과 마찰력을 고려함으로서 평탄한 지면위에 놓여진 물방울의 형태를 생성하는 새로운 방법을 보여준다. 본 연구의 새로운 메타볼 모델은 이전의 메타볼이 가지는 스칼라 공간의 isosurface값을 중력을 고려한 벡터 공간으로 이동하여 얻어진다. 결과로서 광선 추적 기법과 물방울의 그림자 형성을 통해 실제적인 물방울의 형태를 렌더링한다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2003년도 제20회 춘계학술대회 논문집
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.