• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.5
• /
• pp.731-740
• /
• 2001

A numerical study on the quantitative analogy between the fully developed turbulent flow in a straight square duct rotating about an axis perpendicular to that of the duct and that in a stationary curved duct of square cross-section is carried out. In order to clarify the similarity of the two flows, dimensionless parameters K(sub)TR=Re(sup)1/4/√Ro and Rossby number, Ro, in a rotating straight duct flow were used as a set corresponding to K(sub)TC=Re(sup)1/4/√λ and curvature ratio, λ, in a stationary curved duct flow so that they have the same dynamical meaning as those of the fully developed laminar flows. For the large values of Ro or λ, it is shown that the flow field satisfies the asymptotic invariance property, that is, there are strong quantitative similarities between the two flows such as flow patterns and friction factors for the same values of K(sub)TR and K(sub)TC.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.3
• /
• pp.201-209
• /
• 2013

Much numerical and experimental research has been done for the flow around an oscillating airfoil. The main research topics are vortex shedding, dynamic stall phenomenon, MAV's lift and thrust generation. Until now, researches mainly have been concentrated on analyzing the wake flow for the variation of frequency and amplitude at a low angle of attack. In this study, wake structures and acoustic wave propagation characteristics were studied for a plunging airfoil at high angle of attack. The governing equations are the Navier-Stokes equation with LES turbulence model. OHOC (Optimized High-Order Compact) scheme and 4th order Runge-Kutta method were used. The Mach number is 0.3, the Reynolds number is, and the angle of attack is from $20^{\circ}$ to $50^{\circ}$. The plunging frequency and the amplitude are from 0.05 to 0.15, and from 0.1 to 0.2, respectively. Due to the high resolution numerical method, wake vortex shedding and pressure wave propagation process, as well as the propagation characteristics of acoustic waves can be simulated. The results of frequency analysis show that the flow has the mixed characteristics of the forced plunging frequency and the vortex shedding frequency at high angle of attack.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.3
• /
• pp.31-39
• /
• 1999

In this paper, numerical experiment is attempted to analyze and compare the combustion efficiency of the burning sprays due to OFO, FOF triplet / FOOF split doublet injectors. Preconditioned Wavier-Stokes equation system with low Reynolds number $\kappa$-$\varepsilon$ model for turbulence closure, is LU-SGS time-integrated. Spray processes are modeled by DSF analysis with experimentally determined injection characteristics. n-heptane/air global reaction model approximates the combustion for simplicity, and the influence of turbulence on the chemical reaction is included using eddy dissipation model. The results showed the FOF triplet injector of highest combustion efficiency, whereas the OFO type of poet performance. It was also observed that the droplet mean diameter and the average gas temperature due to the mixing efficiency, are the representative parameters for the performance design of combustion.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1254-1264
• /
• 1999

A three-dimensional coupled turbulent fluid flow and solidification process were analyzed in a continuous casting process of a steel slab with Electromagnetic Brake(EMBR). A revised low-Reynolds number $k-{\varepsilon}$ turbulence model was used to consider the turbulent effects. The enthalpy-porosity relation was employed to suppress the velocity within a mushy region. The electromagnetic field was described by Maxwell equations. Tile application of EMBR to the mold region results in the decrease of the transfer of superheat to the narrow face, the increase of temperature in free surface region and most liquid of submold region, and the higher temperature gradient near the solidifying shell. The increasing magnetic flux density effects mainly to the surface temperature of the solidifying shell of narrow face, hardly to the one of wide face. It is seen that in the presence of EMBR a thicker solidifying shell is obtained at the narrow face of the slab.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.285-295
• /
• 2001

The transport phenomena in a wedge-shaped pool of twin-roll continuous caster are affected by the various operating parameters such as the melt-feed pattern, roll-gap thickness, melt-superheat, and casting speed. A computer program has been developed for analyzing the two-dimensional, steady conservation equations for transport phenomena during twin-roll continuous casting process in order to estimate the turbulent melt-flow, temperature fields, and solidification in the wedge-shaped pool. The turbulent characteristics of the melt-flow were considered using a low-Reynolds-number K-$\xi$ turbulence model. Based on the computer program, the effects of the different melt-feed patterns, roll-gap thicknesses, and superheats of melt on the variations of the velocity and temperature distributions, and the mushy solidification were examined. The results show that the liquidus line is located considerably at the upstream region, and in the lower region appear the well-mixed melt-flow and most widely developed mushy zone. Besides, the variation of melt-flow due to varying melt-feed patterns, affects mainly the liquidus line, and scarcely has effects on the solidus line in the outlet region.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.1 no.1
• /
• pp.1-12
• /
• 2009

Wing-in-Ground vehicles and aerodynamically assisted boats take advantage of increased lift and reduced drag of wing sections in the ground proximity. At relatively low speeds or heavy payloads of these craft, a flap at the wing trailing-edge can be applied to boost the aerodynamic lift. The influence of a flap on the two-dimensional NACA 4412 airfoil in viscous ground-effect flow is numerically investigated in this study. The computational method consists of a steady-state, incompressible, finite volume method utilizing the Spalart-Allmaras turbulence model. Grid generation and solution of the Navier-Stokes equations are completed using computer program Fluent. The code is validated against published experimental and numerical results of unbounded flow with a flap, as well as ground-effect motion without a flap. Aerodynamic forces are calculated, and the effects of angle of attack, Reynolds number, ground height, and flap deflection are presented for a split and plain flap. Changes in the flow introduced with the flap addition are also discussed. Overall, the use of a flap on wings with small attack angles is found to be beneficial for small flap deflections up to 5% of the chord, where the contribution of lift augmentation exceeds the drag increase, yielding an augmented lift-to-drag ratio.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.8
• /
• pp.799-807
• /
• 2010

Direct numerical simulations (DNS) of turbulent channel flow for which the drag is reduced by using polymer additives have been performed by a pseudo-spectral method. The Reynolds number based on the friction velocity and half-channel height is 395, and the polymeric stresses due to the polymer additives are evaluated using the FENE-P (finitely extensible nonlinear elastic-Peterlin) model. The numerical results show that the drag reduction rate is significantly affected by the parameters used in the FENE-P model, such as the maximum extensibility and relaxation time of the polymer molecules. The turbulence data for both low- and high-drag reduction regimes are analyzed. In addition, the effects of FENE-P model parameters on the flow characteristics have been investigated for the same drag reduction rate due to the polymer additives. Finally, the present DNS results have been used to verify the correlation between rheological parameters and the extent of drag reduction, which was suggested by Li et al. (2006).

• Journal of Mechanical Science and Technology
• /
• v.16 no.10
• /
• pp.1303-1313
• /
• 2002

Effects of the reduced frequency of upstream wake on downstream unsteady boundary layer flow were simulated by using a Wavier-Stokes code. The Wavier-Stokes code is based on an unstructured finite volume method and uses a low Reynolds number turbulence model to close the momentum equations. The geometry used in this paper is the MIT flapping foil experimental set-up and the reduced frequency of the upstream wake is varied in the range of 0.91 to 10.86 to study its effect on the unsteady boundary layer flow. Numerical solutions show that they can be divided into two categories. One is so called the low frequency solution, and behaves quite similar to a Stokes layer. Its characteristics is found to be quite similar to those due to either a temporal or spatial wave. The low frequency solutions are observed clearly when the reduced frequency is smaller than 3.26. The other one is the high frequency solution. It is observed for the reduced frequency larger than 7.24. It shows a sudden shift of the phase angle of the unsteady velocity around the edge of the boundary layer. The shift of phase angle is about 180 degree, and leads to separation of the boundary layer flow from corresponding outer flow. The high frequency solution shows the characteristics of a temporal wave whose wave length is half of the upstream frequency. This characteristics of the high frequency solution is found to be caused by the strong interaction between unsteady vortices. This strong interaction also leads to destroy of the upstream wake strips inside the viscous sublayer as well as the buffer layer.