• Journal of Environmental Science International
• /
• v.23 no.7
• /
• pp.1223-1232
• /
• 2014

Density currents have been easily observed in environmental flows, for instance turbidity currents and pollutant plumes in the oceans and rivers. In this study, we explored the propagation dynamics of density currents using the FLOW-3D computational fluid dynamics code. The renormalization group (RNG) $k-{\varepsilon}$ scheme, a turbulence numerical technique, is employed in a Reynold-averaged Navier-Stokes framework (RANS). The numerical simulations focused on two different types of intrusive density flows: (1) propagating into a two-layer ambient fluid; (2) propagating into a linearly stratified fluid. In the study of intrusive density flows into a two-layer ambient fluid, intrusive speeds were compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting propagation speed of the density currents. We also numerically reproduced the effect of the ratio of current depth to the overall depth of fluid. The numerical model provided excellent agreement with the analytical values. It was also clearly demonstrated that RNG $k-{\varepsilon}$ scheme within RANS framework is able to accurately simulate the dynamics of density currents. Simulations intruding into a continuously stratified fluid with the various buoyancy frequencies are carried out. These simulations demonstrate that three different propagation patterns can be developed according to the value of $h_n/H$ : (1) underflows developed with $h_n/H=0$ ; (2) overflows developed when $h_n/H=1$ ; (3) intrusive interflow occurred with the condition of 0 < $h_n/H$ < 1.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.187-193
• /
• 2007

The flight vehicles have cavities such as wheel wells and bomb bays. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves. Resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, numerical analysis was performed for cavity flows by the unsteady compressible three dimensional Reynolds-Averaged Navier-Stokes (RANS) equations with Wilcox's ${\kappa}\;-\;{\omega}$ turbulence model. The cavity has the aspect ratios of 2.5, 3.5 and 4.5 for two-dimensional case, same aspect ratios with the W/D ratio of 2 for three-dimensional case. The Mach and Reynolds numbers are 0.53 and 1,600,000 respectively. The flow field is observed to oscillate in the "shear layer mode" with a feedback mechanism. Based on the SPL(Sound Pressure Level) analysis of the pressure variation at the cavity trailing edge, the dominant frequency was analyzed and compared with the results of Rossiter's formula. The MPI(Message Passing Interface) parallelized code was used for calculations by PC-cluster.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.4
• /
• pp.560-570
• /
• 2015

We present a simulation of a hybrid Reynolds-averaged Navier Stokes / Large Eddy Simulation (RANS/LES) based on detached eddy simulation (DES) for a Burrows and Kurkov supersonic planar mixing experiment. The preliminary simulation results are checked in order to validate the numerical computing capability of the current code. Mesh refinement studies are performed to identify the minimum grid size required to accurately capture the flow physics. A detailed investigation of the turbulence/chemistry interaction is carried out for a nine species 19-step hydrogen-air reaction mechanism. In contrast to the instantaneous value, the simulated time-averaged result inside the reactive shear layer underpredicts the maximum rise in $H_2O$ concentration and total temperature relative to the experimental data. The reason for the discrepancy is described in detail. Combustion parameters such as OH mass fraction, flame index, scalar dissipation rate, and mixture fraction are analyzed in order to study the flame structure.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.6
• /
• pp.1-7
• /
• 2005

Incompressible Reynolds-averaged Navier-Stokes equations with $\varepsilon{-SST}$ turbulence model are adopted for the investigation of the flow fields between the square cylinder and the ground. When the grounds moves, the diminish of the shear layer intensity on the ground promotes the interaction between the lower and the upper separated shear layer of the cylinder. Hence vortex shedding occurs at the lower gap height than stationary ground. In the moving ground, the secondary shedding frequency disappears due to the absence of the separation bubble on the ground which exists in the stationary ground. In addition, the shedding frequency and aerodynamic coefficients in the moving ground become higher than those of the stationary ground.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.30-37
• /
• 2004

Flow characteristics of the inflow ahead of a rotating propeller attached to a container ship model were investigated using a two-frame PIV (Particle Image Velocimetry) technique. Ensemble-averaged mean velocity fields were measured at four different blade phases. The mean velocity fields show the acceleration of inflow due to the rotating propeller and the velocity deficit in the near-wake region. The axial velocity distribution of inflow in the upper plane of propeller is quite different from that in the lower plane due to the thick hull boundary layer. The propeller inflow also shows asymmetric axial velocity distribution in the port and starboard side. As the inflow moves toward the propeller, the effect of phase angle variation of propeller blade on the inflow becomes dominant. In the upper plane above the propeller axis the inflow has very low axial velocity and large turbulent kinetic energy, compared with the lower plane. The boundary layer developed along the bottom surface of stern hull forms a strong shear layer affecting vortex structure of the propeller near-wake.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.356-368
• /
• 1992

The strongly interactive flow field near a missile afterbody containing a centered exhaust jet is numerically investigated. The thin shear layer and full formulation of compressible, Reynolds I averaged Navier-Stokes equations are solved. A time-dependent implicit numericals algorithm is used to obtain solution for a variety of flow conditions. Turbulence closure is implemented by the Baldwin-Lomax algebraic eddy viscosity model. An adaptive grid technique is adopted to resolve flow regimes with large gradients and to improve the accuracy and efficiency of the computation, Numerical results show good agreemement with experimental data in all regimes.

• Journal of Ship and Ocean Technology
• /
• v.6 no.2
• /
• pp.37-50
• /
• 2002

In the present numerical simulation of viscous free surface flow around a ship, two-fluids in-compressible Reynolds-averaged Navier-Stokes equations with the standard $\textsc{k}-\varepsilon$turbulence model are discretized on a regular grid by using a finite volume method. A local level-set method is introduced for capturing the free surface movement and the influence of the viscous layer and dynamic boundary condition of the free surface are implicitly considered. Partial differential equations in the level-set method are discretized with second order ENO scheme and explicit Euler scheme in the space and time integration, respectively. The computational results for the Series-60 model with $C_B=0.6$ show a good agreement with the experimental data, but more validation studies for commercial complicated hull forms are necessary.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Journal of Korea Water Resources Association
• /
• v.50 no.1
• /
• pp.55-64
• /
• 2017

This study presents a numerical model for simulating turbidity currents using the ULTIMATE scheme. For this, the layer-averaged model is used. The model is applied to laboratory experiments, where the flume is composed of sloping and flat parts, and the characteristics of propagating turbidity currents are investigated. Due to the universal limiter of the ULTIMATE scheme, the frontal part of the turbidity currents at a sharp gradient without numerical oscillations is computed. Simulated turbidity currents propagate super-critically to the end of the flume, and internal hydraulic jumps occur at the break-in-slope after being affected by the downstream boundary. It is found that the hydraulic jumps are computed without numerical oscillations if Courant number is less than 1. In addition, factors that affect propagation velocity of turbidity currents is studied. The particle size less than $9{\mu}m$ does not affect propagation velocity but the buoyancy flux affects clearly. Finally, it is found that the numerical model computes the bed elevation change due to turbidity currents properly. Specifically, a discontinuity in the bed elevation, arisen from the hydraulic jumps and resulting difference in sediment entrainment, is observed.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.2
• /
• pp.73-83
• /
• 2021

In this study, the numerical investigation of the flow around the SUBOFF submarine model is performed by using the Detached Eddy Simulation (DES) method which is developed based on the SST k-ω turbulence model. At the DES analysis level, complex vortical flows around the submarine model are caused mainly by the vortices due to the appendages and their interactions with the flows from the hull boundary layer and other appendages. The complexity and scale of the vortical flow obtained from the numerical simulations are highly dependent on the grid. The computed local flow properties of the submarine model are compared with the available experimental data showing a good agreement. The DES analysis more reasonably estimates the physical phenomena inherent in the experimental result in a low radius of the propeller plane where vortical flows smaller than the RANS scale are dominant.