• Particle and aerosol research
• /
• v.13 no.4
• /
• pp.173-182
• /
• 2017

An Eulerian-Lagrangin approach is used to compute particle dispersion from a power plant chimney. For air flow, three-dimensional incompressible filtered Navier-Stokes equations are solved with a subgrid-scale model by integrating the Newton's equation, while the dispersed phase is solved in a Lagrangian framework. The velocity ratios between crossflow and a jet of 0.455 and 0.727 are considered. Flow fields and particle distribution of both cases are evaluated and compared. When the velocity ratio is 0.455, it demonstrates a Kelvin-Helmholtz vortex structure above the chimney caused by the interaction between crossflow and a jet, whereas the other case shows flow structures at the top of the chimney collapsed by fast crossflow. Also, complex wake structures cause different particle distributions behind the chimney. The case with the velocity ratio of 0.727 demonstrates strong particle concentration at the vortical region, whereas the case with the velocity ratio of 0.455 shows more dispersive particle distribution. The simulation result shows similar tendency to the experimental result.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.10
• /
• pp.955-966
• /
• 2009

DES method is applied to an axisymmetric base flow at supersonic mainstream. The model is based on the Spalart-Allmaras (S-A) turbulence model in the RANS mode, and is based on the subgrid scale model in the Large-eddy simulation (LES) mode. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology which is less expensive than LES. Flow properties at the edge of base, such as boundary layer thickness, momentum thickness and skin fraction are compared with Dutton et al [experimental data to proper prediction of base flowfiled. From the present results, The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region and small eddy motions inside the recirculating region. Moreover, The present results of using an empirical constant $C_{DES}$ of 1.2 shows good agreement with experimental data than conventional empirical constant $C_{DES}$ of 0.65.

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

• Journal of computational fluids engineering
• /
• v.18 no.3
• /
• pp.72-78
• /
• 2013

The turbulent flow behavior of air supply and exhaustion in the Shin-gum-ho subway station is analyzed for ordinary and emergency state. The depth of Shin-gum-ho station is 43.6m which consists of the island-type platform(8th floor in underground) and a two-story lobby (first & second floor in underground). An emergency stairway connects between the platform and the lobby. Ventilation operation mode for ordinary state is set up as a combination of air supply and exhaustion in the lobby and platform, while for emergency state it is set up as a full air supply in the lobby and a full exhaustion in the platform. The entire station is covered for simulation. The ventilation diffusers are modeled as 95 square shapes of $0.6m{\times}0.6m$ in the lobby and as 222 square shapes of $0.6m{\times}0.6m$ and 4 rectangular shapes of $1.2m{\times}0.8m$ in the platform. The total of 7.5million grids are generated and whole domain is divided to 22 blocks for MPI efficiency of calculation. Large eddy simulation(LES) is applied to solve the momentum equation and Smagorinsky model($C_s$=0.2) is used as SGS(subgrid scale) model. The time-averaged velocity fields are compared to experimental data and show a good agreement with it.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.1
• /
• pp.9-15
• /
• 2015

We investigated the Reynolds number effects on the flow over a twisted offshore structure in the range of 3×10³≤ Re ≤ 1 × 10⁴. To analyze the effect of the twisted surface treatment, a large eddy simulation (LES) with a dynamic subgrid model was employed. A simulation of the cylindrical structure was also carried out to compare the results with those of the twisted offshore structure. As Re increased, the mean drag and lift coefficient of the twisted offshore structure increased with the same tendency as those of the cylindrical structure. However, the increases in the mean drag and lift coefficient of the twisted offshore structure were much smaller than those of the cylindrical structure. Furthermore, elongated shear layer and suppressed vortex shedding from the twisted offshore structure occurred compared to those of the cylindrical cylinder, resulting in a drag reduction and suppression of the vortex-induced vibration (VIV). In particular, the twisted offshore structure achieved a significant reduction of over 96% in VIV compared with that of the cylindrical structure, regardless of increasing Re. As a result, we concluded that the twisted offshore structure effectively controlled the flow structures with reductions in the drag and VIV compared with the cylindrical structure, irrespective of increasing Re.

• Atmosphere
• /
• v.15 no.2
• /
• pp.91-99
• /
• 2005

This study compares the summer monsoon circulations during a heavy rainfall period over the Korean peninsular from 11 to 18 July 2004, simulated by three widely used regional models; WRF, MM5, and RSM. An identical model setup is carried out for all the experiments, except for the physical option differences in the RSM. The three models with a nominal resolution of about 50 km over Korea are nested by NCEP-DOE reanalysis data. Another RSM experiment with the same cumulus parameterization scheme as in the WRF and MM5 is designed to investigate the importance of the representation of subgrid-scale parameterized convection in reproducing monsoonal circulations in East Asia. All thee models are found to be capable of reproducing the general distribution of monsoonal precipitation, extending northeastward from south China across the Korean peninsula, to northern Japan. The results from the WRF and MM5 are similar in terms of accumulated precipitation, but a slightly better performance in the WRF than in the MM5. The RSM improves the bias for precipitation as compared to those from the WRF and MM5, but the pattern correlation is degraded due to overestimation of precipitation in northern China. In the comparison of simulated synoptic scale features, the RSM is found to reproduce the large-scale features well compared to the results from the MM5 and WRF. On the other hand, the simulated precipitation from the RSM with the convection scheme used in the MM5 and WRF is closer to that from the WRF and MM5 simulations, indicating the significant dependency of simulated precipitation in East Asia on the cumulus parameterization scheme.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.3
• /
• pp.529-561
• /
• 2014

The two-phase turbulent flow past an interface-piercing circular cylinder is studied using a high-fidelity orthogonal curvilinear grid solver with a Lagrangian dynamic subgrid-scale model for large-eddy simulation and a coupled level set and volume of fluid method for air-water interface tracking. The simulations cover the sub-critical and critical and post critical regimes of the Reynolds and sub and super-critical Froude numbers in order to investigate the effect of both dimensionless parameters on the flow. Significant changes in flow features near the air-water interface were observed as the Reynolds number was increased from the sub-critical to the critical regime. The interface makes the separation point near the interface much delayed for all Reynolds numbers. The separation region at intermediate depths is remarkably reduced for the critical Reynolds number regime. The deep flow resembles the single-phase turbulent flow past a circular cylinder, but includes the effect of the free-surface and the limited span length for sub-critical Reynolds numbers. At different Froude numbers, the air-water interface exhibits significantly changed structures, including breaking bow waves with splashes and bubbles at high Froude numbers. Instantaneous and mean flow features such as interface structures, vortex shedding, Reynolds stresses, and vorticity transport are also analyzed. The results are compared with reference experimental data available in the literature. The deep flow is also compared with the single-phase turbulent flow past a circular cylinder in the similar ranges of Reynolds numbers. Discussion is provided concerning the limitations of the current simulations and available experimental data along with future research.

• Wind and Structures
• /
• v.28 no.5
• /
• pp.285-298
• /
• 2019

For super long-span bridges, the aerodynamic forces induced by the flow passing the box girder should be considered carefully. And the Reynolds number sensitively of aerodynamic characteristics is one of considerable issue. In the study, a numerical study on the Reynolds number sensitivity of aerodynamic characteristic (flow pattern, pressure distribution and aerodynamic forces) of a twin-box girder were carried out using large eddy simulation (LES) with the dynamic Smagorinsky-Lilly subgrid model. The results show that the aerodynamic characteristics have strong correlation with the Reynolds number. At the leading edge, the flow experiences attachment, departure, and reattachment stages accompanying by the laminar transition into turbulence, causing pressure plateaus to form on the surface, and the pressure plateaus gradually shrinks. Around the gap, attributing that the flow experiences stages of laminar cavity flow, the wake with alternate shedding vortices, and turbulent cavity flow in sequence with an increase in the Reynolds number, the pressures around the gap vary greatly with the Reynold number. At the trailing edge, the pressure gradually recovers as the flow transits to turbulence (the flow undergoes wake instability, shear layer transition-reattachment station), In addition, at relative high Reynolds numbers, the drag force almost does not change, however, the lift force coefficient gradually decreases with an increase in Reynolds number.