• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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.

• 한국소음진동공학회논문집
• /
• 제12권5호
• /
• pp.346-354
• /
• 2002

The present study is to investigate the propagation characteristics of the impulse waves discharged from the exit of the convergent and divergent pipes. An experiment is carried out using a shock tube with an open end and is compared to the computation of the axisymmetric, compressible, unsteady Euler equations, which are solved by the second-order total variation diminishing (TVD) scheme. For the computational work, several initial compression waves are assumed inside the pipe so that those are the same to the experimental ones of the shock tube. The results show that the peak pressures of the impulse waves discharged from the exit of convergent and divergent pipes decrease with an increase in the wavelength of the initial compression wave. All of the impulse waves have a strong directivity toward the pipe axis, regardless of the exit type of the pipe employed. The impulse waves discharged from the divergent pipe are stronger than those from the straight pipe, while the impulse waves of the convergent pipe are weaker than those from the straight pipe. It is found that the convergent pipe can play a role of a passive control to reduce the peak pressure of the impulse wave. The present computations represent the experimented impulse waves with a good accuracy.

• Water Engineering Research
• /
• 제3권1호
• /
• pp.57-62
• /
• 2002

This research is intended to integrate long-term operation rules and real time operation policy for conservation & flood control in a reservoir. The familiar Yield model has been modified and used to provide long-term rule curves. The model employs linear programming technique under given physical conditions, i.e., total capacity, dead storage, spillways, outlet capacity and their respective elevations to find required and desired minimum storage fur different demands. To investigate the system behavior resulting from the above-mentioned operating policy, i.e., the rule curves, the simulation model was used. Results of the simulation model show that the results of the optimization model are indeed valid. After confirmation of the above mentioned rule curves by the simulation models, gate operation procedure was merged with the long term operation rules to determine the optimum reservoir operating policy. In the gate operation procedure, operating policy in downstream flood plain, i.e., determination of damaging and non-damaging discharges in flood plain, peak floods, which could be routed by reservoir, are determined. Also outflow hydrograph and variations of water surface levels for two known hydrographs are determined. To examine efficiency of the above-mentioned models and their ability in determining the optimum operation policy, Esteghlal reservoir in Iran was analyzed as a case study. A numerical model fur the solution of two-dimensional dam break problems using fractional step method is developed on unstructured grid. The model is based on second-order Weighted Averaged Flux(WAF) scheme with HLLC approximate Riemann solver. To control the nonphysical oscillations associated with second-order accuracy, TVD scheme with SUPERBEE limiter is used. The developed model is verified by comparing the computational solutions with analytic solutions in idealized test cases. Very good agreements have been achieved in the verifications.

• 대한조선학회논문집
• /
• 제49권6호
• /
• pp.461-468
• /
• 2012

In this study, a Cartesian-grid method based on finite volume approach is applied to simulate the ship motions in large amplitude waves. Fractional step method is applied for pressure-velocity coupling and TVD limiter is used to interpolate the cell face value for the discretization of convective term. Water, air, and solid phases are identified by using the concept of volume-fraction function for each phase. In order to capture the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with weighed line interface calculation (WLIC) method which considers multidimensional information. The volume fraction of solid body embedded in the Cartesian grid system is calculated using a level-set based algorithm, and the body boundary condition is imposed by a volume weighted formula. Numerical simulations for the two-dimensional barge type model and Wigley hull in linear waves have been carried out to validate the newly developed code. To demonstrate the applicability for highly nonlinear wave-body interactions such as green water on the deck, numerical analysis on the large-amplitude motion of S175 containership is conducted and all computational results are compared with experimental data.

• 천문학회보
• /
• 제39권1호
• /
• pp.70.2-70.2
• /
• 2014

We developed a three-dimensional (3D) magnetohydrodynamic (MHD) simulation code to reproduce the structure of a solar wind and the propagation of a coronal mass ejection (CME) through it. This code is constructed by a finite volume method based on a total variation diminishing (TVD) scheme using an unstructured grid system (Tanaka 1994). The grid system can avoid the singularity arising in the spherical coordinate system. In this study, we made an improvement of the code focused on the propagation of a CME through a solar wind, which extends a previous work done by Nakamizo et al. (2009). We first reconstructed a solar wind in a steady state from physical values obtained at 50 solar radii away from the Sun via an MHD tomography applied to interplanetary scintillation (IPS) data (Hayashi et al. 2003). We selected CR2057 and inserted a spheromak-type CME (Kataoka et al. 2009) into a reconstructed solar wind. As a result, we found that our simulation well captures the velocity, temperature and density profiles of an observed solar wind. Furthermore, we successfully reproduce the general characteristics of an interplanetary coronal mass ejection (ICME) obtained by the Helios 1/2 spacecraft (R. J. FORSYTH et al. 2006).

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

추진기관은 노즐을 통해 추력을 발생하며 축소-확대 형상의 초음속 노즐에서는 노즐의 설계 팽창비가 내부 유동의 전압력과 배압의 압력비보다 매우 클 때 충격파의 발생과 함께 경계층 박리를 유발한다. 노즐 내부에서의 충격파 발생과 유동의 박리는 주어진 유동의 압력비에 가장 적합한 노즐형상을 구현하는 것으로 실제의 구조적 노즐의 형상보다 짧은 노즐에서 나타나는 유동과 같은 현상을 보인다. 수치 해석적 방법으로 고정된 형상의 2차원 노즐 내부의 충격파와 경계층 박리 현상에 관한 연구를 수행하였고 Hunter가 행한 실험적 연구와 비교하였다. 수치해석은 TVD 기법을 이용한 압축성 유체 해석 코드와 SST 2방정식 난류 모델을 이용하여 수행되었다. 낮은 압력비에서의 충격파와 경계층과의 상호작용에 의한 $\lambda$형태의 충격파 시스템을 잘 보여주고 있고 추력 값의 비교를 통해 고정된 형상의 노즐을 이용하여 필요한 운용범위를 충족할 수 있음을 알 수 있었다.

• Journal of Mechanical Science and Technology
• /
• 제20권8호
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence 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 with 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 ratios 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 first 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.

• 한국수자원학회논문집
• /
• 제43권10호
• /
• pp.911-920
• /
• 2010

저수지와 같은 갇혀진 수체는 상류에서 유입되는 오염물질 뿐만 아니라 성층현상에 의해서도 오염될 수 있다. 갇혀진 수체에서의 연직순환은 이러한 오염을 줄이는데 중요한 역할을 하는데, 연직순환을 일으키는 인자로는 빛의 입사, 바람, 물의 온도 및 열의 확산 등이 있으며, 그중에서도 가장 중요한 것은 바람의 영향이다. 그러므로 성층화된 흐름에서 바람에 의해 발생하는 연직순환에 대한 수치모형을 개발하고 적용하는 것이 필요하다. 본 연구는 수온성층흐름을 해석할 수 있는 3차원 수치모형을 제시하였다. 유속성분은x-축과 y-축 방향에서의 운동량방정식으로부터 3단계에 걸쳐 계산되고, 자유수면 변위와 온도변화 등의 스칼라양은 각각 자유수면방정식과 이송-확산 방정식으로부터 계산된다. 본 연구에서 제시한 모형의 정확도를 검증하기 위하여 정사각형수조에서 진동하는 자유수면의 해석해와 비교하였고, 성층화된 흐름에서 발생하는 연직순환에 대하여 수치모의를 실시하였다. 그 결과, 본 연구에서 개발된 수치모형이 흐름 내부의 현상을 잘 묘사함을 알 수 있었다.