• International Journal of Aeronautical and Space Sciences
• /
• 제3권1호
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

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

The propagation of the impulse wave discharged from the Inclined exit of a pipe is investigated through shock tube experiment and numerical computations. The pressure histories and directivities of the impulse wave propagating outside from the exit of pipe with several different configurations are analyzed for the range of the incident shock wave Mach number between 1.1 and 1.4. In the shock tube experiments, the impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. Computations using the two-dimensional. unsteady, compressible, Euler equations are carried out to represent the experimented impulse waves. Computed Schlieren images predict the experimented impulse waves with a good accuracy. The results obtained show that for the radial direction the peak pressure of the impulse wave discharged depends upon the Inclined angle of the exit of the pipe. but for the axial direction it is almost constant regardless of the inclined angle of the pipe exit.

• 한국전산유체공학회지
• /
• 제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권4호
• /
• pp.382-394
• /
• 1988

Dynamic characteristics of thermally-forced stratification process in a square enclosure with a linear temperature profile at the side walls have been investigated through flow visualization experiment and numerical analysis. The experiment was performed on air with the Rayleigh numbers of order $10^5$. A particle tracer method is used for the flow visualization and to obtain a sudden linear temperature profile at the side walls copper blocks which already have a linear temperature profile are come into contact with the thin copper plates of the test section. Immediately a meridional circulation is developed and heat transfer takes place from the wall to the interior region by circulation of fluid and finally a thermal stratification is achieved. In the numerical study, QUICK scheme for convective terms, SIMPLE algorithm for pressure correction, and the implicit method for the time marching are adopted for the integration of conservation equations. Comparison of flow visualization and numerical results shows that the developing flow patterns are very similar in dynamic nature even though there is a time lag due to the inevitable time delay in setting up a linear temperature profile. For high Rayleigh numbers, the oscillatory motion is likely to take place and stratified region is extended. However, initial temperature adjustment process is much slower than that for low Rayleigh numbers.

• 물과 미래
• /
• 제6권1호
• /
• pp.29-35
• /
• 1973

For the simplicity in the analytical solution, the simple surge tank has been chosen for the test where an unsteady flow is porduced by suddenly closing the valve controlling the discharge. The valve is loated immediately downstream from the surge tank. Momentumn equations in the penstock and in the surge column are measured recored on the oscillograph and then the calibration of surge column heights and scale readings on the oscillograph trace are made. The diameter of the penstock are determined by the trial and error method. The water levels in the surge column are computed by numerical integration of the differential equation of the surge tank. The relationships between the results from the experiment and numerical computation are figured, compared and discussed.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집B
• /
• pp.535-540
• /
• 2000

The objective of the present study is to visualize the pulsatile flow fields by using three-dimensional computer simulation and the PIV system. A closed flow loop system was built for the steady and unsteady experiments. The Harvard pulsatile pump was used to generate the pulsatile pressure and velocity waveforms. Conifer powder as the tracing particles was added to water to visualize the flow field. Two consecutive particle images were captured by a CCD camera for the image processing. The cross-correlation method in combination with the moving searching area algorithm was applied for the image processing of the flow visualization. The pulsatile flow fields were visualized effectively by the PIV system in conjunction with the applied algorithm. The range validation and the area interpolation methods were used to obtain the final velocity vectors with high accuracy. The finite volume predictions were used to analyze three-dimensional flow patterns in the bifurcation model. The results of the PIV experiment and the computer simulation are in good agreement and the results show the recirculation zones and formation of the paired secondary flow distal to the apex of the bifurcated model. The results also show that the branch flow is pushed strongly to the inner wall due to the inertial force effect and helical motions are generated as the flow proceeds toward the outer wall.

• 한국항공우주학회지
• /
• 제44권3호
• /
• pp.247-255
• /
• 2016

본 논문에서는 핀틀 추력기의 비정상상태 특성을 파악하기 위하여 구동기의 속도를 측정하였으며, 정상상태 실험을 통해 비정상상태 실험 시스템을 구성하였다. 비정상상태 실험은 총 3 가지의 구동 속도(3.01 mm/s, 5.65 mm/s, 10.83 mm/s)를 이용하여 진행하였다. 그 결과 핀틀이 후진하는 경우가 전진하는 경우보다 더 빨리 명령 압력 값에 수렴하였으며, 이는 핀틀이 후진하는 경우가 연소실의 압력이 높은 상태로 형성되어 있기 때문이다. 핀틀이 전 후진하는 경우에 추력 곡선에 특이점들이 나타났으며, 이는 연소실 압력과 노즐 목 면적 변화에 기인하여 나타난 질유량 변화로 인한 것이다. 구동 속도가 빠를수록 이러한 현상이 뚜렷하게 나타났다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.748-753
• /
• 2001

The objective of the present study is to visualize the pulsatile flow fields by using three-dimensional computer simulation and the PIV system. A closed flow loop system was built for the steady and unsteady experiments. The Harvard pulsatile pump was used to generate the pulsatile pressure and velocity waveforms. Conifer powder as the tracing particles was added to water to visualize the flow field. Two consecutive particle images were captured by a CCO camera for the image processing at several cross section. The range validation and the area interpolation methods were used to obtain the final velocity vectors with high accuracy. The finite volume predictions were used to analyze three-dimensional flow patterns in the bifurcation model. The results of the PIV experiment and the computer simulation are in good agreement and the results show the recirculation zones and formation of the paired secondary flow distal to the apex of the bifurcated model. The results also show that the branch flow is pushed strongly to the inner wall due to the inertial force effect and helical motions are generated as the flow proceeds toward the outer wall.

• 대한기계학회논문집
• /
• 제18권4호
• /
• pp.1019-1030
• /
• 1994

Natural convection heat transfer from a horizontal ice cylinder immersed in quiescent cold pure water was studied experimentally. The experiment was conducted for the ambient water temperatures ranging from $2.0^{\cric}C$ to $10.0^{\circ}C$. The flow fields around an ice cylinder and its melting shapes were visualized and local Nusselt numbers obtained. Especially, its attention was focused on the density maximum effects and stagnation point Nusselt number. From the visualized photographs of flow fields, three distinct flow patterns were observed with the ambient water temperature variation. The melting shapes of ice cylinder are various in shape with flow patterns. Steady state upflow was occured at the range of $2.0^{\circ}C \leq T_{\infty} \leq 4.6^{\circ}C$ and steady state downflow was occured at $T_{\infty} \geq 6.0^{\circ}C$. In the range of $4.7^{\circ}C < T_{\infty} < 6.0^{\circ}C$, three-dimensional unsteady state flow was observed. Especially, the melting shapes of ice cylinder have formed the several spiral flutes for the temperatures ranging from $5.5^{\circ}C$ to $5.8^{\circ}C$. For upflow regime, the maximum stagnation point Nusselt number exists at $T_{\infty} = 2.5^{\circ}C$ and as the ambient water temperature increases the Nusselt number decreases. At ambient water temperature of about $5.7^{\circ}C$, Nusselt number shows its minimum value.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.729-731
• /
• 2008

Three-dimensional unsteady simulation using RNG $\kappa-\varepsilon$ turbulence model is used in full flow passage of model Kaplan turbine. Then the pressure pulsation is obtained. Monitoring data of pressure pulsation in the turbine is obtained through experiment and is compared with the numerical simulation. And a good coherence is shown between the simulation and the experiment. Then the regularity of the pressure pulsation s distribution and transmission in the turbine is discussed in detail.