• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.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.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.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.

• Journal of computational fluids engineering
• /
• v.17 no.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.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.17 no.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.

• Water for future
• /
• v.6 no.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.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• 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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.3
• /
• pp.247-255
• /
• 2016

The purpose of this study is to investigate unsteady-state characteristics of a pintle thruster with various pintle speed. Based on steady state experimental results, non-linear pintle stroke equation is obtained and applied to the unsteady state experimental system. For the unsteady state experiments, three different pintle speeds are used: 3.10 mm/s, 5.65 mm/s, 10.83 mm/s, respectively. Results show that backward pintle stroke results in faster convergence time because of high chamber pressure during backward pintle stroke sequence. During the forward and backward process, thrust curve shows singular points. These phenomenons is caused by variation of mass flow rate, which is mainly due to changes of both chamber pressures and nozzle throat area. This behavior becomes distinctive for a faster pintle speed case.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• 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.