• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.4
• /
• pp.459-472
• /
• 1996

An experimental investigation on the wind flow over smooth bell-shaped two-dimensional hills with hill slopes (the ratio of height to half width) of 0.3 and 0.5 is performed in an atmospheric boundary-layer wind tunnel. Two categories of the models are used in the present investigation; six two-dimensional single-hills, and four continuous double-hills. The measurements of the flow field and surface static-pressure distribution are carried out over the Reynolds number (based on the hill height) of 1.9 $\times 10^4, 3.3 \times 10^4, and 5.6 \times 10^4$. The velocity profiles and turbulence characteristics are measured by the pitot-tube and X-type hot-wire anemometer, respectively. The undisturbed boundary-layer profile on the bottom surface of the wind tunnel is reasonably consistent with the power-law profile with $\alpha = 7.0 (1/\alpha$ is the power-law exponent) and shows good spanwise uniformities. The profiles of turbulent intensity are found to be consistent along the centerline of the wind tunnel. The measured non-dimensional speed-up profiles at the hill crest show good agreements with the predictions of Jackson and Hunt's linear theory. The flow separation occurs in the hill slope of 0.5, and the oil-ink dot method is used to find the reattachment points in the leeside of the hill. The measured reattachment points are compared with the numerical predictions. Comparisons of the mean velocity profiles and surface pressure distributions between the numerical predictions and the experimental results show good agreements.

• Wind and Structures
• /
• v.30 no.6
• /
• pp.547-558
• /
• 2020

With the development of economy and construction technology, more and more bridges are built in complex mountainous areas. Accurate assessment of wind parameters is important in bridge construction at complex terrain. In order to investigate the wind characteristics in the high-altitude difference area, a complex mountain terrain model with the scale of 1:2000 was built. By using the method of wind tunnel tests, the study of wind characteristics including mean wind characteristics and turbulence characteristics was carried out. The results show: The wind direction is affected significant by the topography, the dominant wind direction is usually parallel to the river. Due to the sheltering effect of the mountain near the bridge, the wind speed and wind attack angle along the bridge are both uneven which is different from that at flat terrain. In addition, different from flat terrain, the wind attack angle is mostly negative. The wind profiles obey exponential law and logarithmic law. And the fitting coefficient is consistent with the code which means that it is feasible to use the method of wind tunnel test to simulate complex terrain. As for turbulence characteristics, the turbulence intensity is also related to the topography. Increases sheltering effect of mountain increases the degree of breaking up the large-scale vortices, thereby increasing the turbulence intensity. Also, the value of turbulence intensity ratio is different from the recommended values in the code. The conclusions of this study can provide basis for further wind resistance design of the bridge.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1417-1426
• /
• 2001

Turbulent flow characteristics in the near wake of a square cylinder have been studied experimentally by using a Digital PIV method. Experiments are performed at the Reynolds numbers of 1600 and 3900 based on the free-stream velocity and the square height. The ensemble averaged turbulence statistics are acquired from 2030 realizations of instantaneous fluctuating velocity field after the conventional Reynolds decomposition. The differences in turbulent intensity and Reynolds shear stress profiles fur both oases indicate that the effect of Reynolds number seems to be descernible mainly due to the occurrence of transition in the separated shear layer. Because of the periodic nature of vortex shedding process, transverse velocity fluctuations contribute dominantly , to turbulent kinetic energy distribution. A comparison with previous LDV data obtained at much higher Reynolds number shows a fairly good agreement each other. It turns out that the effect of Reynolds number diminishes as increasing Reynolds number, which is a well-known feature of a sharp-edged bluff body wake. The streamwise variation of turbulence intensities are compared with those from a circular cylinder along the centerline at the same Reynolds number. The overall magnitudes and the decay rates of turbulence intensities are quite similar, but some differences are noticeble especially in the transverse intensity variation.

• Journal of Ocean Engineering and Technology
• /
• v.28 no.4
• /
• pp.268-273
• /
• 2014

To investigate the turbulence characteristics within the boundary layer over a flat plate, an experimental study was performed using a PIV technique in a circular water channel. For two water velocities, 0.92 and 1.99 m/s, the water velocity profiles were taken and analyzed to determine turbulent characteristics such as the Reynolds stress, Taylor micro-length scale, and Kolmogorov length scale within the defect law region of the boundary layer. These analysis methods may be applied to research on the friction drag reduction technology using micro-bubbles or an air sheet over the surface of a ship's hull, because the physical reason for the friction drag reduction could be found by understanding the variation of the turbulence characteristics and structures in the boundary layer.

• 한국신재생에너지학회:학술대회논문집
• /
• 2011.05a
• /
• pp.59.2-59.2
• /
• 2011

Although there are many activities on the construction of wind farm to produce amount of power from the wind, in practice power productions are not as much as its expected capabilities. This is because a lack of both the prediction of wind resources and the aerodynamic analysis on turbines with far wake effects. In far wake region, there are velocity deficits and increases of the turbulence intensity which lead to the power losses of the next turbine and the increases of dynamic loadings which could reduce system's life. The analysis on power losses and the increases of fatigue loadings in the wind farm is needed to prevent these unwanted consequences. Therefore, in this study velocity deficits have been predicted and aerodynamic analysis on turbines in the far wake is carried out from these velocity profiles. Ainslie's eddy viscosity wake model is adopted to determine a wake velocity and aerodynamic analysis on wind turbines is predicted by the numerical methods such as blade element momentum theory(BEMT) and vortex lattice method(VLM). The results show that velocity recovery is more rapid in the wake region with higher turbulence intensity. Since the velocity deficit is larger when the turbine has higher thrust coefficient, there is a huge aerodynamic power loss at the downstream turbine.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Journal of the korean society of oceanography
• /
• v.32 no.2
• /
• pp.75-84
• /
• 1997

Relation between the length scale and the wall proximity function in the Mellor-Yamada level 2.5 turbulence closure model has been investigated through various experiments using a range of wall proximity functions. The model performance has been evaluated quantitatively by comparing with laboratory data for wind-driven flow (Baines and Knapp, 1965) and for open-channel flows without and with adverse wind action (Tsuruya, 1985). Comparison shows that a symmetric wall proximity function used by Blumberg and Mellor(1987) gives rise to current profiles with better accuracy than asymmetric wall proximity functions considered. It is noted that in modelling homogeneous flows the length scale 1= 0.31${\|}$z${\|}$(1+z/h) can be used with tolerable accuracy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.891-899
• /
• 1997

This paper presents the results of a detailed experimental examination of fully developed asymmetric flows between annular tubes with square-ribbed surface roughness. The main emphasis of the research has been on establishing the turbulence structure, particularly in the central region of the channel where the two dissimilar wall flows interact. Measurements have included profiles of time mean velocities, turbulence intensities, turbulent shear stresses, triple velocity correlations, skewness, and flatness. The region of greatest interaction is characterized by strong diffusional transport of turbulent shear stress and kinetic energy from rough toward the smooth wall region, giving rise to an appreciable separation between the planes of zero shear stresses depending on positions of roughness on the walls.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.1-7
• /
• 2001

Wind flow perturbations, recirculations and turbulence generated by buildings often dominate air pollutant distributions around buildings. This paper describes dispersion of contaminants in the vicinity of a building by solving the concentration equation based on previously simulated wind flow field. Turbulence closure is achieved by using the standard k-ε two-equation model. The paper shows application of the CIP method for solving a species concentration equation of contaminant gas around a rectangular building for two different sources under conditions of neutral atmospheric stratification. Results have been compared to the experimental data and the previous numerical results by hybrid scheme. The computational results of concentration profiles by the CIP method agree well with experimental data.