• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• Journal of the Korean Society of Visualization
• /
• v.9 no.4
• /
• pp.16-21
• /
• 2011

The effects of stroke change on turbulent kinetic energy for the in-cylinder flow of a four-valve SI engine were studied. For this study, the same intake manifold, head, cylinder, and the piston were used to examine turbulence characteristics in two different strokes. In-cylinder flow measurements were conducted using three dimensional LDV system. The measurement method, which simultaneously collects 3-D velocity data, allowed a evaluation of turbulent kinetic energy inside a cylinder. High levels of turbulent kinetic energy were found in regions of high shear flow, attributed to the collisions of intake flows. These specific results support the more general conclusion that the inlet conditions play the dominant role in the generation of the turbulence fields during the intake stroke. However, in the absence of two counter rotating vortices, this intake generated turbulent kinetic energy continues to decrease but at a much faster rate.

• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.52-61
• /
• 2004

In the present study, a solution algorithm for the computation of unsteady flows on unstructured meshes is presented. Dual time stepping is incorporated to achieve the second-order temporal accuracy while reducing errors associated with linearization and factorization. This allows any time step size, which is suitable for considering physical phenomena of interest. The Gauss-Seidel scheme is used to solve the linear system of equations. A special treatment based on spring analogy is made to handle meshes with high aspect-ratio cells. The present method was validated by comparing the results with experimental data and those obtained from rigid motion.

• Journal of Advanced Marine Engineering and Technology
• /
• v.24 no.1
• /
• pp.31-40
• /
• 2000

This paper describes the numerical computations of the interaction between the streamwise vortex and a flat plate 3-D turbulent boundary layer. In the present study, the main interest is in the behavior of the streamwise vortices introduced in turbulent boundary layers. The flow behind a vortex generator is modeled by the information that is avilable from studies on the dalta winglet. An algorithm of the solution of the incompressible Navier-Stokes equations for three-dimensional turbulent flows, together with a two layer turbulent model to resolve the near-wall flow, is based on the method of artificial compressibility. The present results show boundary layer distortion due to vortices, such as strong spanwise flow divergence and boundary thinning, and have a good agreement with the experimental data.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.42 no.1
• /
• pp.38-43
• /
• 2006

The 3-D modeling of ocean finite element method(OFEM) using $k-{\varepsilon}$ turbulent model and tetrahedron grids has been used to investigate the internal wave generation during the expansion of the deep water from the open sea to the shelf with a simple shape, which can be widely used in the fields of submarine development, ocean environment and meteorology, etc. In this paper, the detailed configuration of internal wave with its length and height and also the distribution of salinity and turbulent kinematic energy, etc. were derived. It is hoped that this OFEM method can be successfully applied to the numerical calculation of internal wave for and the oceanographic problems (tidal flows around underwater hill, plateau, Georges Bank, etc.) and ocean engineering problems(flow past artificial sea reefs) in future.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.12 no.3
• /
• pp.25-39
• /
• 2004

An experimental study was quantitatively carried out in order to investigate the influence of flow conditions on a boundary layer in the near-wake of a flat plate. Tripping wires attached at various positions were selected to change flow conditions of a boundary layer in the vicinity of trailing edge. The flows such as laminar, transitional, and turbulent boundary layer at 0.98C from the leading edge are imposed to investigate the evolution of symmetric and asymmetric wake. Measurements were made at freestream velocity of 6.0m/s, and the corresponding Reynolds number is $2.8{\times}10^5$. An x-type hot-wire probe(55P61) was employed to measure at 8 stations in the near-wake region. Test results show that the near-wake of the flat plate for the case of a laminar and transitional boundary layer is sensitive to mean flow shear generated after separation but for the case of turbulent boundary layer is insensitive.

• Journal of Power System Engineering
• /
• v.4 no.2
• /
• pp.11-16
• /
• 2000

Simultaneous quantitative measurements are made of the velocity fields for two-dimensional natural convection in a rectangular enclosure using PIV(Particle Image Velocimetry). The experiments are performed at a Prandtl number of 6.62, an aspect ratio of 1.0, Rayleigh numbers from $1.294{\times}10^6\;to\;3.8841{\times}10^6$, and angles of inclination of $0^{\circ},\;30^{\circ}\;and\;60^{\circ}$ inside a $30mm{\times}30mm{\times}8mm$ cavity made of an acrylic glass 10mm, with two isothermal copper walls kept at a prescribed temperature. The experimental results agreed very well with the numerical results. It was found that the flow consisted of a large double convection cell at angle of inclination of $60^{\circ}$.

• Journal of Advanced Marine Engineering and Technology
• /
• v.15 no.1
• /
• pp.71-74
• /
• 1991

Uncertainty analyses accompanied by the measurement of the velocity vectors in 3-D cavity flows are carried out. Twenty-one elemental errors are esimated or calculated according to the ANIS/ASME uncertainty analysis manual. Error components associater with the PIV(Particle Imaging Velocimetry) are reasonably small and the errors caused by the flow characteristics are fairly large, which confirm the reliability of the PIV measurement and also give good information to the planning phase of the experiment by discriminating the most critical parameter. The present study reveals that vector length expressed by pixels is the most influential. Calculated relative uncertainty for the all experimental conditions is ranging about 5-10% in terms of the representative velocity 0.5U. U is here the belt velocity on the cavity apparatus. Approximating equations to show the relative rss uncertainties are given and graphic representations are followed for the easier understanding of the uncertainty, existing in the velocity profiles of the cavity flow.

• Wind and Structures
• /
• v.4 no.1
• /
• pp.31-44
• /
• 2001

Mean concentrations of ammonia gas released as a tracer from an isolated low-rise building have been measured and predicted. Predictions were calculated using computational fluid dynamics (CFD) and two dispersion models: a diffusion model and a Lagrangian particle tracking technique. Explicit account was taken of the natural variation of wind direction by a technique based on the weighted summation of individual steady state wind direction results according to the probability density function of the wind direction. The results indicated that at distances >3 building heights downstream the weighted predictions from either model are satisfactory but that in the near wake the diffusion model is less successful. Weighted solutions give significantly improved predictions over unweighted results. Lack of plume spread is identified as the main cause of inaccuracies in predictions and this is linked to inadequate resolution of flow features and mixing in the CFD model. Further work on non-steady state simulation of wake flows for dispersion studies is recommended.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.7
• /
• pp.1866-1872
• /
• 1994

Turbulent flow in a channel with a square rib periodically mounted on one wall is studied by large-eddy simulation(LES). An efficient 3D Navier-Stokes solver has been written for this geometry using a fractional step method and a multi-grid technique. The Reynolds number considered is 82, 000 based on the mean velocity above the obstacle height. Near-wall turbulence is approximated by a wall-layer model based on the turbulence intensity at the grid point nearest a solid wall. The results show a good qualitative agreement with experiments currently available for a single rib, indicating that LES can be a useful tool in simulating complex turbulent flows.