• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.39-51
• /
• 1998

This paper describes the analysis of flow field around an automobile. The governing equations of the 3-D unsteady incompressible Navier-Stokes equations are solved by the iterative time marching scheme. The Chimera grid technique has been applied to efficiently simulate the flow around the side-view mirror. To validate the capability of simulating the flow around a ground vehicle, the flows around the Ahmed body with 12.5$^{\circ}$ and 30$^{\circ}$ of slant angles are simulated and good agreements with experiment and other numerical results are achieved. To validate Chimera grid technique, the flow field around a cylinder was also calculated. The computed results are also well agreed with other numerical results and experiment. After code validations, the flow phenomena around the ground vehicle are evidently shown. The flow around the side-view mirror is also well simulated using the Chimera grid technique.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.1
• /
• pp.11-25
• /
• 1991

The steady, incompressible developing 3-dimensional turblent flow in a square sectioned curved duct has been investigated by using partially-parabolic equation and Finite Analytic Method. The calculation of turbulent flow field is performed using 2-equation K-$\epsilon$ turbulence model, modified wall function, simpler algorithm and numerically generated body fitted coordinates. Iso-mean velocity contours at the various sections are compared with the existing experimental data and elliptic solutions by other authors. In the region of $0^{\circ}<{\theta}<71^{\circ}$, present results agree with the experimental data much better than the elliptic solution for the similar number of grid points. Furthermore, for the same tolerance, the present solution converges four times faster than the elliptic solution.

• The Bulletin of The Korean Astronomical Society
• /
• v.38 no.1
• /
• pp.47.1-47.1
• /
• 2013

Supersonic turbulence is believed to decay rapidly within a flow crossing time irrespective of the degree of magnetization. However, this consensus of decaying magnetohydrodynamic (MHD) turbulence relies on local isothermal simulations, which are unable to investigate the role of global magnetic fields and structures. Utilizing three-dimensional MHD simulations including interstellar cooling and heating, we investigate decaying MHD turbulence within cold neutral medium sheets embedded in warm neutral medium. Early evolution is consistent with previous studies characterized rapid decay of turbulence with the decaying time shorter than a flow crossing time and power-law temporal decay of turbulent kinetic energy with slope of -1. If initial magnetic fields are strong and perpendicular to the sheet, however long term evolutions of kinetic energy shows that a significant amount of turbulent energy still remains even after ten flow crossing times, and decaying rate is reduced as field strengths increase. We analyse power spectra of remaining turbulence to show that incompressible, in-plane motions dominate.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.754-759
• /
• 2001

This paper presents unsteady computational investigations and wind tunnel tests on the flow field around a square cylinder with a gap between the body and the ground plane. Two-dimensional unsteady, incompressible Navier-Stokes codes are developed for the computation of the viscous turbulent flows. By computing the flow around a square cylinder without ground effect, three two-equation turbulence models are evaluated and the developed code is validated. The results show a good agreement with experimental values and other computational results. Critical gap height at which the formation of Karman vortex streets is interrupted, is demonstrated and another transition regime is pointed out

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.3
• /
• pp.375-385
• /
• 1996

Turbulent flow and heat transfer characteristics around staggered tube banks were studied using the 3-D Navier-Stokes equations and energy equation governing a steady incompressible flow, which were reformulated in a non-orthogonal coordinate system with cartesian velocity components and discretized by the finite volume method with a non-staggered variable arrangement. The predicted turbulent kinetic energy using RNG $k-{\varepsilon}$ model was lower than that of standard $k-{\varepsilon}$ model but showed same result for mean flow field quantities. The prediction of the skin friction coefficient using RNG $k-{\varepsilon}$ model showed better trend with experimental data than standard $k-{\varepsilon}$ model result. The inclined flow showed higher velocity and skin friction coefficient than transverse flow because of extra strain rate ($\frac{{\partial}w}{{\partial}y}$). Also, this was why the inclined flow showed higher local heat transfer coefficient than the transverse flow.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.14 no.12
• /
• pp.1004-1013
• /
• 2002

Turbulent carbon dioxide flows and cooling heat transfers under supercritical state in a straight duct with a square cross-section are numerically analyzed employing low Reynolds number $\kappa-\varepsilon$ model and algebraic stress model. The flow is assumed to be (quasi-incompressible. Predicted Nusselt number and friction factor are compared with the experimental data, Blasius correlation for friction factor and Dittus-Boelter correlation for Nusselt number. Computational results for the Fanning's friction factor agree well with the all Rohsenow and Choi's correlation, Liou and Hwang's experimental data and Blasius correlation. The results obtained by algebraic stress model agree more with the Liou and Hwang's experimental data, while the results obtained by low Reynolds number $\kappa-\varepsilon$ model agree more with Blasius correlation. In the computation of Nusselt number, Dittus-Boelter correlation can not exactly fit the computational results. Therefore we propose the new correlation$Nu=0.053Re^{0.73}Pr^{0.4}$for the turbulent cooling heat transfer of carbon dioxide under supercritical state.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.267-277
• /
• 1993

A finite element method(FEM) is presented and applied to the two-dimensional bottom turbulent boundary layer. The time-dependent incompressible motion of a viscous fluid is formulated by using the well-known Navier-Stokes equations and vorticity equation in terms of the velocity and pressure fields. The general numerical formulation is based on Galerkin method and solved by introducing the mixing length theory of Prandtl for eddy kinematic viscosity of a turbulent flow field. Numerical computations of the transport of sediment on an arbitrary sea-bed due to wave motion in the turbulent boundary layer are carried out. The results obtained by the FEM made clear the difference in characteristic features between the boundary layer due to oscillatory flow and the boundary layer due to wave motion.

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.22-31
• /
• 2001

The two-dimensional unsteady incompressible viscous flow behind rectangular bluff bodies between two parallel plates was analyzed numerically. The steady state flow and the vortex flow behind rectangular bluff bodies in the channel were investigated for two regimes i.e., the laminar(Re = 100, 300, 500) and the turbulent flows(Re = 10⁴∼10/sup 6/). The vortex shedding was generated by a physical disturbance(6%) numerically imposed at the rear of the bluff bodies for a short time. It was observed that the perturbed flow became periodic after a transient period. And in the case of unsteady inflow, the sinusoidal pulsatile flow was applied as the inlet condition in the turbulent flow of Reynolds number of 1.0×10/sup 5/. FLUENT code was employed to solve the problems. The power-law scheme was used to get stable linearized equations and the PISO algorithm was applied to finding the solution of them.