• Journal of applied mathematics & informatics
• /
• 제13권1_2호
• /
• pp.53-65
• /
• 2003

Here we present a study of small-amplitude, shallow water waves on sloping beds. The beds considered in this analysis are linear and quadratic in nature. First we start with stating the relevant governing equations and boundary conditions for the theory of water waves. Once the complete prescription of the water-wave problem is available based on some assumptions (like inviscid, irrotational flow), we normalize it by introducing a suitable set of non-dimensional variables and then we scale the variables with respect to the amplitude parameter. This helps us to characterize the various types of approximation. In the process, a summary of equations that represent different approximations of the water-wave problem is stated. All the relevant equations are presented in rectangular Cartesian coordinates. Then we derive the equations and boundary conditions for small-amplitude and shallow water waves. Two specific types of bed are considered for our calculations. One is a bed with constant slope and the other bed has a quadratic form of surface. These are solved by using separation of variables method.

• 한국전산유체공학회지
• /
• 제5권2호
• /
• pp.38-46
• /
• 2000

It has been indicated that the rotor/stator interaction has distinct causes of unsteadiness, such as the viscous vortex shedding, wake/stator interaction and potential rotor/stator interaction. In this paper, the mechanism of unsteady potential interaction in one stage axial compressor is numerically investigated for blade row ratio 1:1 and 2:3 at design point and for blade row ratio 2:3 at off-design point in two-dimensional view point. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting(FVS) and Cubic spline interpolation is applied on zonal interface. In this study the flow unsteadiness due to potential interaction are found to be larger in blade row ratio 2:3 than in 1:1. The total pressure rise in blade row ratio 2:3 is closer to the real value in design point than that in 1:1. The change of unsteady pressure amplitude according to the variation of stator exit pressure is very small.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.42-49
• /
• 2008

A wavelet method is presented in order to improve the computational efficiency of two dimensional unsteady flow problems while maintaining the order of accuracy of conventional CFD schemes. First, by using the interpolating wavelet transformation including decomposition and thresholding, an adaptive dataset to a solution is constructed. Then, inviscid and viscous fluxes are calculated only at the points within an adaptive dataset, which enhances the computational efficiency. Second, thresholding step is modified to maintain the spatial and temporal accuracy of conventional CFD schemes automatically by selecting the threshold value between user-defined value and the magnitude of spatial or temporal truncation error. The wavelet method suggested in this study is successfully applied to various unsteady flow problems and it is shown that the computational efficiency is enhanced with maintaining the computational accuracy of CFD schemes.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.8-15
• /
• 2012

A code developed using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method is applied to simulate three-dimensional internal waves. The material interface is regarded as a moving contact discontinuity and is captured on the basis of mass conservation without any additional treatment across the interface. Inviscid fluxes are estimated using the flux-difference splitting scheme for incompressible fluids of different density. The hybrid Cartesian/immersed boundary method is used to enforce the boundary condition for a moving three-dimensional body. Immersed boundary nodes are identified within an instantaneous fluid domain on the basis of edges crossing a boundary. The dependent variables are reconstructed at the immersed boundary nodes along local normal lines to provide the boundary condition for a discretized flow problem. The internal waves are simulated, which are generated by an pitching ellipsoid near an material interface. The effects of density ratio and location of the ellipsoid on internal waves are compared.

• 대한기계학회논문집B
• /
• 제21권11호
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• Journal of Mechanical Science and Technology
• /
• 제18권2호
• /
• pp.313-324
• /
• 2004

A rich phenomenon in the dynamics of azimuthal vortices in a circular cylinder caused by the inertial oscillation is investigated numerically at high Reynolds numbers and moderate Rossby numbers. In the actual spin-up flow where both the Ekman circulation and the bottom friction effects are included, the first appearance of a seed vortex is generated by the Ekman boundary-layer on the bottom wall and the subsequent roll-up near the corner bounded by the side wall. The existence of the small vortex then rapidly propagates toward the inviscid region and induces a complicated pattern in the distribution of azimuthal vorticity, i.e. inertial oscillation. The inertial oscillation however does not deteriorate the classical Ekman-pumping model in the time scale larger than that of the oscillatory motion. Motions of single vortex and a pair of vortices are further investigated under a slip boundary-condition on the solid walls. For the case of single vortex, repeated change of the vorticity sign is observed together with typical propagation of inertial waves. For the case of a pair of vortices with a two-step profile in the initial azimuthal velocity, the vortices' movement toward the outer region is resisted by the crescent-shape vortices surrounding the pair. After touching the border between the core and outer regions, the pair vortices weaken very fast.

• 대한기계학회논문집B
• /
• 제24권3호
• /
• pp.382-389
• /
• 2000

The instantaneous and ensemble averaged flow characteristics of a round jet issuing normally into a crossflow was studied using a flow visualization technique and Particle Image Velocimetry measurements. Experiments were performed at a jet-to-crossflow velocity ratio, 3.3, and two Reynolds numbers, 1050 and 2100, based on crossflow velocity and jet diameter. Instantaneous laser tomographic images of the vertical center plane of the crossflow jet showed that there exist very different natures in the flow structures of the near field jet even though the velocity ratio is the same. It was found that the shear layer becomes much thicker when the Reynolds number is 2100 due to the strong entrainment of the inviscid fluid by turbulent interaction between the jet and crossflow. The mean and second order statistics were calculated by ensemble averaging over 1000 realizations of instantaneous velocity fields. The detail characteristics of mean flow field, stream wise and vertical r.m.s. velocity fluctuations, and Reynolds shear stress distributions were presented. The new PlV results were compared with those from previous experimental and LES studies.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2007년도 제29회 추계학술대회논문집
• /
• pp.35-38
• /
• 2007

액체로켓엔진(Liquid Rocket Engine) 의 천이성능 예측을 위해 선행연구 되었던 LRE 시스템 모듈화 프로그램의 결과를 살펴보고, 일본의 로켓엔진 동적 해석 프로그램(Rocket Engine Dynamic Simulator)의 엔진 시스템 동적 해석 방법과 모델링에 대해 고찰하였다. LRE 시스템 모듈화 프로그램에서는 각 구성품에 대한 설계 인자를 수학적으로 모델링하였고 구성품 간의 유량과 압력을 매칭시켜 통합하여, 로켓엔진 시스템의 요구조건을 만족하는 각 구성품에 대한 주요 설계 파라미터를 도출하는 과정에 관하여 논의하였다. 로켓엔진 시스템을 유한한 배관요소들의 연결로 모델링하고, 시간의 함수로 표현되는 보존방정식을 적용하여 터보펌프, 밸브, 오리피스，추력실 등 유체기기의 작동 특성을 모사하는 동특성 설계 과정에 관하여 고찰한다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2007년도 제29회 추계학술대회논문집
• /
• pp.49-54
• /
• 2007

오일러나 Navier-Stokes방정식을 통한 터빈 캐스케이드 유동 해석은 비교적 정확한 해를 구할 수 있으나 많은 계산 시간을 필요로 한다. 비점성, 비압축성 유동에 적용이 가능한 패널법은 빠르고 합리적인 유동 정보를 얻을수 있지만 고속 유동의 경우 압축성 보정이 반드시 이뤄져야한다. 본 논문에서는 압축성이 보정된 패널법을 이용하여 터빈 블레이드 표면의 속도 분포를 계산하였다. 그 결과, 압축성이 보정된 패널법의 결과는 실험이나 유한 체적법에 의해 계산된 결과와 잘 일치하였다.