• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1435-1442
• /
• 2014

This study presented numerical simulations of smooth-bed flows in the rectangular open-channel using the source code by OpenFOAM. For the analysis of the turbulent flow, Large Eddy Simulations were carried out and the dynamic sub-grid scale model proposed by Germano et al. (1991) is used to model the residual stress term. In order to analyze the coherent structure, the uw quadrant method proposed by Lu and Willmarth (1973) is used and the contribution rate and the fraction time of the instantaneous Reynolds stress are obtained in the Reynolds stress. The results by the present study are analyzed and compared with data from previous laboratory studies and direct numerical simulations. It is found that the contribution rate of the ejection events is larger than that of sweep events over the buffer layer in the open-channel flow over the smooth bed, however, the frequency of the sweep event is higher than that of the ejection events.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.291-300
• /
• 2002

An efficient large eddy simulation algorithm is used to compute surface pressure distributions on an eleven story (target) building on the NIST campus. Local meteorology, neighboring buildings, topography and large vegetation (trees) all play an important part in determining the flows and therefore the pressures experienced by the target. The wind profile imposed at the upstream surface of the computational domain follows a power law with an exponent representing a suburban terrain. This profile accounts for the flow retardation due to friction from the surface of the earth, but does not include fluctuations that would naturally occur in this flow. The effect of neighboring buildings on the time dependent surface pressures experienced by the target is examined. Comparison of the pressure fluctuations on the single target building alone with those on the target building in situ show that, owing to vortices shed by the upstream buildings, fluctuations are larger when such buildings are present. Even when buildings are lateral to or behind the target, the pressure disturbances generate significantly different flows around this building. A simple grid-free mathematical model of a tree is presented in which the trunk and the branches are each represented by a collection of spherical particles strung together like beads on a string. The drag from the tree, determined as the sum of the drags of the component particles, produces an oscillatory, spreading wake of slower fluid, suggesting that the behavior of trees as wind breakers can be modeled usefully.

• Wind and Structures
• /
• v.13 no.6
• /
• pp.557-580
• /
• 2010

A new inflow turbulence generation method and a combined dynamic SGS model recently developed by the authors were applied to evaluate the wind effects on 508 m high Taipei 101 Tower. Unlike the majority of the past studies on large eddy simulation (LES) of wind effects on tall buildings, the present numerical simulations were conducted for the full-scale tall building with Reynolds number greater than $10^8$. The inflow turbulent flow field was generated based on the new method called discretizing and synthesizing of random flow generation technique (DSRFG) with a prominent feature that the generated wind velocity fluctuations satisfy any target spectrum and target profiles of turbulence intensity and turbulence integral length scale. The new dynamic SGS model takes both advantages of one-equation SGS model and a dynamic production term without test-filtering operation, which is particular suitable to relative coarse grid situations and high Reynolds number flows. The results of comparative investigations with and without generation of inflow turbulence show that: (1) proper simulation of an inflow turbulent field is essential in accurate evaluation of dynamic wind loads on a tall building and the prescribed inflow turbulence characteristics can be adequately imposed on the inflow boundary by the DSRFG method; (2) the DSRFG can generate a large number of random vortex-like patterns in oncoming flow, leading to good agreements of both mean and dynamic forces with wind tunnel test results; (3) The dynamic mechanism of the adopted SGS model behaves adequately in the present LES and its integration with the DSRFG technique can provide satisfactory predictions of the wind effects on the super-tall building.

• Structural Engineering and Mechanics
• /
• v.46 no.2
• /
• pp.197-212
• /
• 2013

With more and more high-rise building being constructed in recent decades, bluff body flow with high Reynolds number and large scale dimensions has become an important topic in theoretical researches and engineering applications. In view of mechanics, the key problems in such flow are high Reynolds number turbulence and fluid-solid interaction. Aiming at such problems, a parallel fluid-structure interaction method based on socket parallel architecture was established and combined with the methods and models of large eddy simulation developed by authors recently. The new method is validated by the full two-way FSI simulations of 1:375 CAARC building model with Re = 70000 and a full scale Taipei101 high-rise building with Re = 1e8, The results obtained show that the proposed method and models is potential to perform high-Reynolds number LES and high-efficiency two-way coupling between detailed fluid dynamics computing and solid structure dynamics computing so that the detailed wind induced responses for high-rise buildings can be resolved practically.

• Wind and Structures
• /
• v.26 no.2
• /
• pp.89-98
• /
• 2018

For most full-scale tall buildings the Reynolds number of a flow field around a circular cylinder under strong wind is usually greater than $2{\times}10^7$, which is difficult to achieve in most wind tunnel tests. To explore the wind characteristics of tall cylindrical buildings with equidirectional grooves from subcritical to transcritical flow ($6.6{\times}10^4{\leq}Re{\leq}3.3{\times}10^5$ and $9.9{\times}10^6{\leq}Re{\leq}7.2{\times}10^7$), wind tunnel tests and full-scale large eddy simulations were carried out. The results showed that the rectangular-grooves narrow the wake width due to the downstream movement of the separation point and the deeper grooves cause smaller mean and fluctuating pressure while the peak pressure is little affected. Furthermore, the grooves lead to lower frequency of vortex shedding but the Strouhal number remains at the range from 0.15 to 0.35. The drag coefficient of the cylinders with grooves was found to be 2~3 times as large as that of smooth cylinders.

• International Journal of High-Rise Buildings
• /
• v.8 no.4
• /
• pp.255-264
• /
• 2019

This paper describes large eddy simulation of wind pressures on a square cylinder in a uniform flow and a high-rise building immersed in an atmospheric turbulent boundary layer. For the atmospheric boundary layer case, the inflow turbulence is generated by a numerical wind tunnel. In the numerical simulation, particular attention is devoted to the performance of an auto hexahedral non-structural mesh. Both simulations are performed for three grid systems: an auto hexahedral non-structured grid, a structured Cartesian grid and a non-structured triangular prism grid, and for three grid numbers. The present study shows that the auto hexahedral unstructured mesh achieves the best simulation results for wind pressures on the square cylinder and the high-rise building. When the grid number is sufficiently large, the differences among the results obtained from the three investigated grid systems are not significant. However, the advantage of the auto hexahedral unstructured mesh becomes clear when the grid number decreases, because it enables a balanced distribution of orthogonal grids. The results described in this paper demonstrate that the auto hexahedral non-structured mesh has good potential applicability to simulation of urban flows.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.3
• /
• pp.833-840
• /
• 2014

This study presents numerical simulation of turbulent flows in a rectangular open-channel that has a width-to-depth ratio of 2 using the source code provided by OpenFOAM. Large eddy simulations are carried out by solving the filtered continuity and momentum equations numerically. For the non-isotropic residual stress term, Smagorinsky's (1963) model is used. The flow in the open-channel whose width-to-depth ratio is 2, from experiment of Tominaga et al. (1989), is simulated numerically. Simulation results are compared with measured data by Tominga et al. (1989) and Nezu and Rodi (1985) and with LES data by Shi et al. (1999). Comparisons revealed that the model simulates the mean flow and turbulence statistics well. Specifically, the model reproduced the inner secondary currents located at the corner of sidewall and free surface successfully. In addition, the vortical component of turbulence intensity shows bulged contours towards the bottom edge.

• Wind and Structures
• /
• v.10 no.6
• /
• pp.511-532
• /
• 2007

The present study aims to generate turbulent inflow data to more accurately represent the turbulent flow around a square cylinder when the inflow turbulence level is significant. The modified random flow generation (RFG) technique in conjunction with a previously developed LES code is successfully adopted into a finite element based fluid flow solver to generate the required inflow turbulence boundary conditions for the three-dimensional (3-D) LES computations of transitional turbulent flow around a square cylinder at Reynolds number of 22,000. The near wall region is modelled without using wall approximate conditions and a wall damping coefficient is introduced into the calculation of sub-grid length scale in the boundary layer of the cylinder wall. The numerical results obtained from simulations are compared with each other and with the experimental data for different inflow turbulence boundary conditions in order to discuss the issues such as the synthetic inflow turbulence effects on the 3-D transitional flow behaviour in the near wake and the free shear layer, the basic mechanism by which stream turbulence interacts with the mean flow over the cylinder body and the prediction of integral flow parameters. The comparison among the LES results with and without inflow turbulence and the experimental data emphasizes that the turbulent inflow data generated by the present RFG technique for the LES computation can be a viable approach in accurately predicting the effects of inflow turbulence on the near wake turbulent flow characteristics around a bluff body.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.16 no.1
• /
• pp.55-63
• /
• 2012

The turbulent flow characteristics in a coaxial injector were investigated by the nonlinear $k-{\varepsilon}-f_{\mu}$ model of Park et al.[1] and large eddy simulation (LES). In order to analyze the geometric effects on the scalar mixing for nonreacting variable-density flows, several recessed lengths and momentum flux ratios are selected at a constant Reynolds number. The nonlinear $k-{\varepsilon}-f_{\mu}$�� model proposed the meaningful characteristics for various momentum flux ratios and recess lengths. The LES results showed the changes of small-scale structures by the recess. When the inner jet was recessed, the development of turbulent kinetic energy became faster than that of non-recessed case. Also, the mixing characteristics were mainly influenced by the variation of shear rates, but the local mixing was changed by the adoption of recess.

• 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.