• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.316-323
• /
• 2000

A nonlinear low-Reynolds-number heat transfer model is developed to predict turbulent flow and heat transfer in separated and reattaching flows. The $k-{\varepsilon}-f_{\mu}$ model of Park and Sung (1997) is extended to a nonlinear formulation, based on the nonlinear model of Gatski and Speziale (1993). The limiting near-wall behavior is resolved by solving the $f_{\mu}$ elliptic relaxation equation. An improved explicit algebraic heat transfer model is proposed, which is achieved by applying a matrix inversion. The scalar heat fluxes are not aligned with the mean temperature gradients in separated and reattaching flows; a full diffusivity tensor model is required. The near-wall asymptotic behavior is incorporated into the $f_{\lambda}$ function in conjunction with the $f_{\mu}$ elliptic relaxation equation. Predictions of the present model are cross-checked with existing measurements and DNS data. The model preformance is shown to be satisfactory.

• Ocean Science Journal
• /
• v.43 no.2
• /
• pp.81-90
• /
• 2008

A study of estuarine flows during a neap tide was performed using 13-hour roving acoustic Doppler current profiles (ADCP) and conductivity-temperature-depth (CTD) profiles in the Altamaha River estuary, Georgia, U.S.A. The least-squared harmonic analysis method was used to fit the tidal ($M_2$) component and separate the flow into two components: the tidal and residual ($M_2$-removed) flows. We applied this method to depth-averaged data. Results show that the $M_2$ component demonstrates over 95% of the variability of observation data. As the flow was dominated by the $M_2$ tidal component in a narrow channel, the tidal ellipse distribution was essentially a back-and-forth motion. The amplitude of $M_2$ velocity component increased slightly from the river mouth (0.45 m/sec) to land (0.6 m/sec) and the phase showed fairly constant values in the center of the channel and rapidly decreasing values near the northern and southern shoaling areas. The residual flow and transport calculated from depth-averaged flow shows temporal variability over the tidal time scale. Strong landward flows appeared during slack waters which may be attributed to increased baroclinic forcing when turbulent mixing decreases.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.39-45
• /
• 2010

The turbulent flows in a tunnel mock-up($10L{\times}0.5W{\times}0.25H$ m3 : scale reduction 1/20) with rectangular cross section were investigated. The instantaneous velocity fields of Re = 49,029, 89,571 were measured by the 2-D PIV system which is consisted of double pulsed Nd:Yag laser and the tracer particles in the straight-duct mock-up where the flows were fully developed. The mean velocity profiles were taken from the ensemble averages of 1,000 instantaneous velocity fields. Simultaneously, numerical simulations(RANS) were performed to compare with experimental data using STREAM code. Non-linear eddy viscosity model (NLEVM : Abe-Jang-Leschziner Eddy Viscosity Model) was employed to resolve the turbulent flows in the duct. The calculated mean velocity profiles were well compared with PIV results. In the log-law profiles, the experimental data were in good agreement with numerical simulations all the way to the wake region except the viscous sub-layer (near wall region).

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.815-818
• /
• 2002

The elliptic relaxation model(ERM) with the inhomogeneous correction intermediate between near wall with and far from the wall. The source of the ERM usually was appled quasi-homogeneous pressure-strain correlation in homogeneous situations. This formulation was easily applied to the linear model or non-linear pressure-strain model. It is observed that the boundary conditions of the relaxation operator dominate the homogeneous pressure-strain model in the near wall region. While looking at high-Reynolds number flows, it was found necessary to modify the effect of the relaxation operator throughout the log region by accounting for gradients of the flatness variable and turbulent length scales. These effects are kinematic blocking of the wall normal velocity fluctuation and pressure reflections from the surface. This model is wall distances and unit vectors which make the model applicable to flows boundary by a complex geometry. Inhomogeneous correction model is computed inertial and non-inertial channel flow These are compared DNS(Kim et at., Kristofffrsen & Andersson) for channel flow. The present model could be predicted well for rotating flows.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.2
• /
• pp.22-31
• /
• 1997

We developed a general purpose program for the analysis of flows in the combustor with recirculating flow regime and simulated the flows. The program uses non-staggered grids based on finite volume method and the primitive variables are cartesian velocities. The combustion model is irreversible one step reaction with infinite chemistry The Favre averaged governing equations are considered and the clipped gaussian distribution is considered as a probability density function of the conserved scalar. We calculated turbulent diffusion flame with recirculating flow regime. Simulation shows two recirculating regions like experimental results. Velocity, turbulent kinetic energy, temperature and concentration distribution in simulation agree well with experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.7
• /
• pp.1030-1040
• /
• 1998

One of the main unresolved issues in large-eddy simulation(LES) of wall-bounded turbulent flows is the requirement of high spatial resolution in the near-wall region, especially in the spanwise direction. Such high resolution required in the near-wall region is generally used throughout the computational domain, making simulations of high Reynolds number, complex-geometry flows prohibitive. A grid-embedding strategy using a nonconforming spectral domain-decomposition method is proposed to address this limitation. This method provides an efficient way of clustering grid points in the near-wall region with spectral accuracy. LES of transitional and turbulent channel flow has been performed to evaluate the proposed grid-embedding technique. The computational domain is divided into three subdomains to resolve the near-wall regions in the spanwise direction. Spectral patching collocation methods are used for the grid-embedding and appropriate conditions are suggested for the interface matching. Results of LES using the grid-embedding strategy are promising compared to LES of global spectral method and direct numerical simulation. Overall, the results show that the spectral domain-decomposition grid-embedding technique provides an efficient method for resolving the near-wall region in LES of complex flows of engineering interest, allowing significant savings in the computational CPU and memory.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.1
• /
• pp.11-23
• /
• 1989

A computer code has been developed to simulate three-dimensional viscous flows over a ship-stern. Semi-elliptic forms of Reynolds equations are adopted and numerically generated body-fitted coordinate systems are used to resolve complex geometries of the ship-hull. A standard form of $k-\varepsilon$ turbulence model is adopted for evaluation of the Reynolds stresses. Turbulent flows on a model with 3:1 elliptic sections and the SSPA-720 container ship model are predicted by using the code. Calculated pressure distributions of hull-surfaces and mean velocity distributions are generally in good agreements with measured values in wind-tunnels. But turbulent kinetic energies tend to be over-estimated near the stern in comparison with measured data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.1
• /
• pp.230-240
• /
• 1990

An exprimental research has been carried out to find the intermittent flow pattern in the transition region of a turbulent round jet in order to elucidate detailed turbulence structure and to accumulate basic data necessary for computational turbulence modelling. Turbulent signals were processed digitally to obtain conventional or conditional velocity components. The high-order conditional correlations obtained in this study showed similar trends as those of other free shear flows. It was found that the non-turbulent fluid contributes negligibly to the turbulent kinetic energy production and its diffusive transport and that the diffusion by bulk convection has the same order of magnitude as the gradient diffusion in the free boundary region. The statistical analyses such as flatness factor, skewness factor and probability density functions of turbulent and non-turbulent zone durations have also been performed.

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

Two-phase(air-solid, air-liquid droplet) turbulent round jet has been analyzed numerically using two equation turbulence model. The mean motion of suspending particles in air has been treated as the secondary fluid with virtual density and eddy viscosity. In this paper, the local mean velocity of secondary fluid is not assumed to be the same as that of the primary one. Dissipation rate of turbulent kinetic energy which arises because the particles can not catch up with the turbulent fluctuations of the primary fluid has been modelled by using the concept of Kolmogorov's spectral energy transfer. Numerical computations were performed for flows with different volume fraction of the dispersed phase and the diameter of particle. Results show that the total rate of turbulent energy dissipation, turbulent intensities and spreading rate of jets are reduced by the increase of volume fraction of dispersed phase. However it does not show consistent tendency with increasing the particle diameter. This investigation also shows that presence of particles in the fluid modifies the structure of the primary fluid flow significantly. Predicted velocity profiles and turbulence properties qualitatively agree with available data.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2002.04a
• /
• pp.50-54
• /
• 2002

With a view to estimating the effect of flows on evolving sprays and combustion in ramjet combustor and corresponding extent of combustion, ram air flows in combustion chamber is numerically experimented. Preconditioned three dimensional Navier-Stokes system of equations per transient, compressible, turbulent flows in IRR(Integral Rocket Ramjet) combustor is numerically integrated. Flow properties in the side-dump ramjet combustor, rectangular duct with two 60-deg curved inlets located radially at an angle of 180-deg, are addressed in terms of mixing quality and extent of combustion efficiency.