• Wind and Structures
• /
• v.17 no.1
• /
• pp.1-19
• /
• 2013

A multiscale finite element method is applied to the Spalart-Allmaras turbulence model based detached-eddy simulation (DES). The multiscale arises from a decomposition of the scalar field into coarse (resolved) and fine (unresolved) scales. It corrects the lack of stability of the standard Galerkin formulation by modeling the scales that cannot be resolved by a given spatial discretization. The stabilization terms appear naturally and the resulting formulation provides effective stabilization in turbulent computations, where reaction-dominated effects strongly influence near-wall predictions. The multiscale DES is applied in the context of high-Reynolds flow over the Commonwealth Advisory Aeronautical Council (CAARC) standard tall building model, for both uniform and turbulent inflows. Time-averaged pressure coefficients on the exterior walls are compared with experiments and it is demonstrated that DES is able to resolve the turbulent features of the flow and accurately predict the surface pressure distributions under atmospheric boundary layer flows.

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

A study of the shock wave turbulent boundary layer interaction is presented. The focus of the study is the interactions of the shock waves with the turbulent boundary layer on the falt plate. Three examples are investigated. The computations are performed, using mixed explicit-implicit generalized Galerkin finite element method. The linear equations at each time step are solved by a preconditioned GMRES algorithm. Numerical results indicate that the implicit scheme converges to the asymptotic steady state much faster than the explicit counterpart. The computed surface pressures and skin friction coefficients display good agreement with experimental data. The flowfield manifests a complex shock wave system and a pair of counter-rotating vortices.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.5 no.3
• /
• pp.25-32
• /
• 2001

A passive control method of the interaction between a weak normal shock-wave and a turbulent boundary-layer was simulated using two-dimensional Navier-Stokes computations. The inflow Mach number just upstream of the normal shock wave was 1.33. A porous plate wall having a cavity underneath was used to control the shock-wave/turbulent boundary-layer interaction. The flows through the porous holes and inside the cavity were investigated to get a better understanding of the flow physics involved in this kind of passive control method. The present computations were validated by some recent wind tunnel tests. The results showed that downstream of the rear leg of the $\lambda$-shock wave the main stream inflows into the cavity, but upstream of the rear leg of the $\lambda$-shock wave the flow proceeds from the cavity toward to the main stream. The flow through the porous holes did not choke fur the present shock/boundary layer interaction.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.2
• /
• pp.54-63
• /
• 2009

Large eddy simulation was conducted for flow development in a chamber with wall injection which simulates the cold flow in an idealized hybrid rocket motor. It was found that a peculiar timescale, roughly corresponding to St~0.5, resides in the flowfield resulting from the interaction between the main oxidizer and wall injected flows. However, the fact that this time characteristics is absent in the temperature field in the vicinity of the wall indicates that even a small regression rate renders the passive scalar, such as temperature, dissimilar to the velocity field. This implies that a classical approach, which assumes that constant turbulent Prandtl number, should be replaced by a more sophisticated turbulence models to accurately predict the temperature field in the hybrid motor.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.31-39
• /
• 2000

In this study, the numerical analysis for laminar and turbulent flows around a circular cylinder is carried out from Reynolds number $10^4\;to\;10^6$, where the dramatic change of drag coefficients exists around Reynolds number $4{\times}10^5$. The numerical results are compared with those of the experiments and they show reasonable agreements. The good results can be obtained without turbulent model by minimizing mesh size in the body normal direction and increasing mesh numbers in the body tangential direction. The future study for 2-dimensional flat plates with no pressure gradients and 3-dimensional bodies must be made to verify whether the same conclusions can be drawn.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.137-149
• /
• 2022

In this paper, the fluctuating lift and drag forces on 5:1 rectangular cylinders with two different geometric scales in three turbulent flow-fields are investigated. The study is particularly focused on understanding the influence of the ratio of turbulence integral length scale to structure characteristic dimension (the length scale ratio). The results show that both fluctuating lift and drag forces are influenced by the length scale ratio. For the model with the larger length scale ratio, the corresponding fluctuating force coefficient is larger, while the spanwise correlation is weaker. However, the degree of influence of the length scale ratio on the two fluctuating forces are different. Compared to the fluctuating drag, the fluctuating lift is more sensitive to the variation of the length scale ratio. It is also found through spectral analysis that for the fluctuating lift, the change of length scale ratio mainly leads to the variation in the low frequency part of the loading, while the fluctuating drag generally follows the quasi-steady theory in the low frequency, and the slope of the drag spectrum at high frequencies changes with the length scale ratio. Then based on the experimental data, two empirical formulas considering the influence of length scale ratio are proposed for determining the lift and drag aerodynamic admittances of a 5:1 rectangular cylinder. Furthermore, a simple relationship is established to correlate the turbulence parameter with the fluctuating force coefficient, which could be used to predict the fluctuating force on a 5:1 rectangular cylinder under different parameter conditions.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.2
• /
• pp.1-9
• /
• 1999

Applicability of the pressure gradient method which is formulated based on pressure gradient is verified against turbulent flow analysis. In the pressure gradient method, pressure gradient instead of pressure itself is obtained using continuity constraint. Since correct pressure gradient is found only when mass conservation is satisfied, pressure gradient method can reflect physics of flow field properly The pressure gradient method is formulated with semi-staggered grid system which locates each primitive variables on the same grid point but evaluates pressure gradient in-between. This grid system ensures easy programming and reflection of correct physics in analysis. For verifying applicability of this method, the pressure gradient method is applied to turbulent flow analysis with low Reynolds number $\kappa$-$\varepsilon$ model. Turbulent flows include fully developed channel flow, backward-facing step flow, and conical diffuser flow. Prediction results show that the pressure gradient method can be applied to turbulent flow analysis. However, the pressure gradient method requires somewhat long computation time. Proper way to find optimum under-relaxation factor, $\gamma$, is also need to be developed.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• Journal of the Korean Society of Safety
• /
• v.9 no.1
• /
• pp.100-109
• /
• 1994

In this study, it is investigated that the possibility of the numerical simulation for the incineration of the hazardous material, crbon tetrachloride($CCl_4$). A 3-dimensional numerical technology is applied for turbulent reacting flows of the full-scale Dow Chemical incinerator. The calculations are made by a CRAY-2S, super computer. The major parameters considered in this study are kiln revolution rate (rpm), filling ratio of the solid waste(f), burner Injection velocity and angle, and turbulent air jets for swirl. And the employed turbulent reaction model is the eddy break-up model which is a kind of fast chemistry model assuming general equilibrium and used for a premixed flame. The calculated flow fields are presented and discussed. 1) The presence of turbulent air nozzles for swirl gives rise to visible increase of the convective motion over the region of the solid waste. This implies the possibility to enhance the mixing of the waste with the surrounding all and thereby to reduce thermal and species stratification, which were reported in a large rotary kiln operation. 2) Considering that the location of the recirculation region has a strong relation with the heating rate of the solid waste, the control of the recirculation region by the burner injection angle Is quite desirable in the sense of the flexible design of the rotary kiln incinerator for a carbon tetrachloride. 3) Finally, it is found that the eddy break-up model Is not suitable for carbon tetrachloride($CCl_4$) because this model is not incorporated the flame inhibition trend due to the presence $CCl_4$compound.

• International Journal of Ocean System Engineering
• /
• v.2 no.1
• /
• pp.50-56
• /
• 2012

This study investigates flow fields and energy dissipation due to regular wave interaction with a perforated vertical breakwater, through velocity data measurement in a two-dimensional wave tank. As the waves propagate through the perforated breakwater, the incoming wave energy is reflected back to the ocean, dissipated due to very turbulent flows near the perforations and inside the chamber, and transmitted through the perforations of the breakwater. This transmitted energy is further reduced due to the presence of the perforated back wall. Hence most of the energy is either reflected or dissipated in the vicinity of the structure, and only a small amount of the incoming wave energy is transmitted through the structure. In this study, particle image velocimetry (PIV) technique was employed to measure two-dimensional instantaneous velocity fields in the vicinity of the structure. Measured velocity data was treated statistically, and used to calculate mean flow fields, turbulence intensity and turbulent kinetic energy. For investigation of the flow pattern, time-averaged mean velocity fields were examined, and discussed using the cross-sections through slot and wall for comparison. Flow fields were obtained and compared for various cases with different regular wave conditions. In addition, turbulent kinetic energy was estimated as an approach to understand energy dissipation near the perforated breakwater. The turbulent kinetic energy was distributed against wave height and wave period to see the dependence on wave conditions.