• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.614-617
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulations of three dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence free flow field at each physical time step. The one equation Spalart-Allmaras turbulence model has been adopted to solve the high-Reynolds number flow fields. This method has been applied to calculate the steady flow fields around submarine configurations and unsteady flow fields around a 3-D infinite cylinder.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.517-518
• /
• 2006

Transient and asymmetric density distributions have been investigated by three-dimensional digital speckle tomography with a novel integration method. Multiple CCD images captured movements of speckles in three angles of view simultaneously because the flows were asymmetric and unsteady. The speckle movements which have been formed by a ground glass between no flow and downward butane flow from a circular half opening have been calculated by a cross-correlation tracking method so that those distances can be transferred to deflection angles of laser rays for density gradients. A novel integration method has been developed to obtain projection data from the deflection angles for the speckle tomography.'The three-dimensional density fields have been reconstructed from the accurate projection values by a real-time multiplicative algebraic reconstruction technique (MART) with the developed integration method.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.352-357
• /
• 2010

A parametric study has been carried out to elucidate the characteristics of channel flow with a streamwise-periodic array of cylinders. This flow configuration is relevant to heat exchanger applications. The presence of cylinders in channel flow causes the attached wall boundary layer to separate, leading to significant change in flow instabilities. There exist two kinds of instabilities; flow undergoes a primary instability (Hopf bifurcaiton) at a lower Reynolds number, and the unsteady two-dimensional flow becomes unstable to three-dimensional disturbances at a higher Reynolds number. We report here the dependencies of the primary instability as well as the flow characteristics of the subsequent unsteady flow including flow-induced forces and Strouhal number of vortex shedding, on the distance between the cylinder and the channel wall.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.6
• /
• pp.767-776
• /
• 1985

Unsteady thermally stratified flow caused by two-dimensional surface discharge of warm water into a oblong channel was investigated. Experimental study was focused on the rapidly developing thermal diffusion at small Richardson number. The basic objective were to study the interfacial mixing between a flowing layer of warm water and an underlying body of cold water and to accumulate experimental data to test computational turbulence models. Mean velocity field measurements were carried out by using NMR-CT (Nuclear Magnetic Resonance-Computerized Tomography). It detects quantitative flow image of any desired section in any direction of flow in short time. Results show that at small Richardson number warm layer rapidly penetrates into the cold layer because of strong turbulent mixing and instability between the two layers. It is found that the transfer ofheat across the interface is more vigorous than that of momentum. It is also proved that the NMR-CT technique is a very valuable tool to measure unsteady three dimensional flow field.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.12
• /
• pp.1774-1783
• /
• 1998

The transition from steady laminar to chaotic convection in a glass melting furnace specified by upper surface temperature distribution has been studied by the direct numerical analysis of the two and three-dimensional time dependent Navier-Stokes equations. The thermal instability of convection roll may take place when modified Rayleigh number($Ra_m$) is larger than $9.71{\times}10^4$. It is shown that the basic flows in a glass melting furnace are steady laminar, unsteady periodic, quasi-periodic or chaotic flow. The dimensionless time scale of unsteady period is about the viscous diffusion time, ${\tau}_d=H^2/{\nu}_0$. Through primary and secondary instability analyses the fundamental unsteady feature in a glass melting furnace is well defined as the unsteady periodic or weak chaotic flow.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.54-59
• /
• 2010

A parametric study has been carried out to elucidate the characteristics of channel flow with a streamwise-periodic array of cylinders. This flow configuration is relevant to heat exchanger applications. The presence of cylinders in channel flow causes the attached wall boundary layer to separate, leading to significant change in flow instabilities. There exist two kinds of instabilities; flow undergoes a primary instability (Hopf bifurcation) at a lower Reynolds number, and the unsteady two-dimensional flow becomes unstable to three-dimensional disturbances at a higher Reynolds number. We report here the dependencies of the primary instability as well as the flow characteristics of the subsequent unsteady flow, including flow-induced forces and Strouhal number of vortex shedding, on the distance between the cylinder and the channel wall.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.34-37
• /
• 2009

Deforming mesh should be used when bodies are deforming or moving relative to each other due to the presence of aerodynamic forces and moments. Also, the flow solver for such a flow problem should satisfy the geometric conservation law to ensure the accuracy of the solutions. In this paper, a RANS(Reynolds Averaged Navier-Stokes) solver including automatic mesh capability using TFI(Transfinite Interpolation) method and GCL is developed and applied to flows induced by oscillating wings with given frequencies. The computations are performed both on deforming meshes and on rigid meshes. The computational results are compared with experimental data, which shows a good agreement.

• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.37-45
• /
• 2010

Deforming mesh should be used when bodies are deforming or moving relative to each other due to the presence of aerodynamic forces and moments. Also, the flow solver for such a flow problem should satisfy the geometric conservation law to ensure the accuracy of the solutions. In this paper, a RANS(Reynolds Averaged Navier-Stokes) solver including automatic mesh capability using TFI(Transfinite Interpolation) method and GCL is developed and applied to flows induced by oscillating wings with given frequencies. The computations are performed both on deforming meshes and on rigid meshes. The computational results are compared with experimental data, which shows a good agreement.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.550-555
• /
• 2003

When a shock wave arrives at an open end of tube, an impulse wave is discharged from the tube exit and complicated flow is formed near tube exit. The flow field is influenced by the cross-sectional geometry of tube exit, such as circular, square, rectangular, trapezoid and etc. In the current study, three-dimensional propagation characteristics of impulse wave discharged from the tube exit with non-circular cross section are numerically investigated using a CFD method. Total variation diminishing (TVD) scheme is used to solve the three-dimensional, unsteady, compressible Euler equations. Computations are performed for the Mach numbers of the incident shock wave $M_{s}$ below 1.5. The results obtained show that the peak pressure of the impulse wave and propagation directivity depends on the cross-sectional geometry of tube exit and the Mach number of incident shock wave.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.