• Atmosphere
• /
• v.17 no.4
• /
• pp.315-326
• /
• 2007

The effects of windbreaks on the reduction of suspended particles are investigated using a computational fluid dynamics (CFD) model with the ${\kappa}-{\varepsilon}$ turbulence closure scheme based on the renormalization group (RNG) theory. In the control experiment, the recirculation zones behind the storage piles are generated and, as a whole, relatively monotonous flow patterns appear. When the windbreaks with the 0% porosity are constructed, the recirculation zones are generated by the windbreaks and very complicated flow patterns appear due to the interference between the windbreaks and storage piles. The porosity of the windbreaks suppresses the generation of the recirculation zone and decreases the wind velocity in the windbreaks as well as that outside the windbreaks. As the emission of suspended particles from the storage piles are closely related with the friction velocity at the surfaces of the storage piles, variation of the friction velocity and total amount of the emission of the suspended particles with the height and porosity of the windbreaks are investigated. The results show that higher and more porous windbreaks emit less suspended particles and that the reduction effect of the porosity is still more effective than that of the height. In the case of the windbreak with 30 m height and 50% porosity, friction velocities above the storage piles are smaller than the critical friction velocity above which particles would be suspended. As a result, total amount of suspended particles are much fewer than those in other cases.

• Journal of Korean Society for Atmospheric Environment
• /
• v.22 no.1
• /
• pp.73-84
• /
• 2006

Hazardous Air Pollutants (HAPs) are characterized by being relatively heavier and denser than that of ambient air due to the various reasons such as higher molecular weight, low temperature and other complicated chemical transformations (Witlox, 1994). In an effort to investigate transport and diffusion from instantaneous emission of heavy gas, Lagrangian Particle Dispersion Model (LPDM) coupled with the RAMS output was employed. Both deposition process and buoyancy term were added on the atmospheric diffusion equations of LPDM, and the locations and concentrations of dense gas particle released from instantaneous single point source (emitting initially for 10 minutes only) were analyzed. The result overall shows that adding deposition process and buoyancy terms on the diffusion equation of LPDM has very small but detectable effect on the vertical and horizontal distribution of Lagrangian particles that especially transported for a fairly long traveling time. Also the slumping of dense gas can be found to be ignored horizontally compared to the advection by the horizontal wind suggesting that it was essential to couple the Lagrangian particle dispersion model coupled with the RAMS model in order to explain the dispersion of HAPs more accurately. However, during the initial time of instantaneous emission, buoyancy term play an important role on the vertical locations of dense particles for near surface atmosphere and around source area, indicating the importance of densities of HAPs in the beginning stage or short duration for the risk assessment of HAPs or management of heavy vapors during the explosive accidents.

• Wind and Structures
• /
• v.26 no.5
• /
• pp.267-277
• /
• 2018

Numerical simulations are performed of a long flexible cylinder undergoing vortex-induced vibration at a Reynolds number of 500. The cylinder is pinned at both ends, having an aspect ratio of 100 (cylinder length to cylinder diameter) and a mass ratio of 4.2 (structural mass to displaced fluid mass). Temporal and spatial information on the cross-flow (CF) and in-line (IL) vibrations is extracted. High modal vibrations up to the $6^{th}$ in the CF direction and the $11^{th}$ in the IL direction are observed. Both the CF and IL vibrations feature a multi-mode mixed pattern. Mode competition is observed. The $2^{nd}$ mode with a low frequency dominates the IL vibration and its existence is attributed to a wave group propagating back and forth along the span. Distributions of fluid force coefficients are correlated to those of the CF and IL vibrations along the span. Histograms of the x'-y motion phase difference are evaluated from the total simulation time and a complete vibration cycle representing the standing or travelling wave pattern. Correlations between the phase difference and the vibrations are discussed. Vortex structures behind the cylinder show an interwoven near-wake pattern when the standing wave pattern dominates, but an oblique near-wake pattern when the travelling wave pattern prevails.

• Wind and Structures
• /
• v.26 no.5
• /
• pp.293-304
• /
• 2018

A numerical study is conducted on the flow characteristics of a rectangular cylinder (chord-to-width ratio C/W = 2 - 10) mounted close to a rigid wall at gap-to-width ratios G/W = 0.25 - 6.25. The effects of G/W and C/W on the Strouhal number, vortex structure, and time-mean drag and lift forces are examined. The results reveal that both G/W and C/W have strong influences on vortex structure, which significantly affects the forces on the cylinder. An increase in G/W leads to four different flow regimes, namely no vortex street flow (G/W < 0.75), single-row vortex street flow ($0.75{\leq}G/W{\leq}1.25$), inverted two-row vortex street flow ($1.25<G/W{\leq}2.5$), and two-row vortex street flow (G/W > 2.5). Both Strouhal number and time-mean drag are more sensitive to C/W than to G/W. For a given G/W, Strouhal number grows with C/W while time-mean drag decays with C/W, the growth and decay being large between C/W = 2 and 4. The time-mean drag is largest in the single-row vortex street regime, contributed by a large pressure on the front surface, regardless of C/W. A higher C/W, in general, leads to a higher time-mean lift. The maximum time-mean lift occurs for C/W = 10 at G/W = 0.75, while the minimum time-mean lift appears for C/W = 2 at the same G/W. The impact of C/W on the time-mean lift is more substantial in single-row vortex regime. The effect of G/W on the time-mean lift is larger at a larger C/W.

• Wind and Structures
• /
• v.30 no.5
• /
• pp.465-474
• /
• 2020

Vortex-induced vibrations of a yawed flexible cylinder near a plane boundary are numerically investigated at a Reynolds number Re_n= 500 based on normal component of freestream velocity. Free to oscillate in the in-line and cross-flow directions, the cylinder with an aspect ratio of 25 is pinned-pinned at both ends at a fixed wall-cylinder gap ratio G/D = 0.8, where D is the cylinder diameter. The cylinder yaw angle (α) is varied from 0° to 60° with an increment of 15°. The main focus is given on the influence of α on structural vibrations, flow patterns, hydrodynamic forces, and IP (Independence Principle) validity. The vortex shedding pattern, contingent on α, is parallel at α=0°, negatively-yawed at α ≤ 15° and positively-yawed at α ≥ 30°. In the negatively- and positively-yawed vortex shedding patterns, the inclination direction of the spanwise vortex rows is in the opposite and same directions of α, respectively. Both in-line and cross-flow vibration amplitudes are symmetric to the midspan, regardless of α. The RMS lift coefficient C_L,rms exhibits asymmetry along the span when α ≠ 0°, maximum C_L,rms occurring on the lower and upper halves of the cylinder for negatively- and positively-yawed vortex shedding patterns, respectively. The IP is well followed in predicting the vibration amplitudes and drag forces for α ≤ 45° while invalid in predicting lift forces for α ≥ 30°. The vortex-shedding frequency and the vibration frequency are well predicted for α = 0° - 60° examined.

• Wind and Structures
• /
• v.30 no.5
• /
• pp.533-546
• /
• 2020

Large Eddy Simulation (LES) is used to study the effects of steady slot suction on the aerodynamic forces of and flow around a wall-mounted finite-length square cylinder. The aspect ratio H/d of the tested cylinder is 5, where H and d are the cylinder height and width, respectively. The Reynolds number based on free-stream oncoming flow velocity U_∞ and d is 2.78×10⁴. The suction slot locates near the leading edge of the free end, with a width of 0.025d and a length of 0.9d. The suction coefficient Q (= U_s/U_∞) is varied as Q = 0, 1 and 3, where U_s is the velocity at the entrance of the suction slot. It is found that the free-end steady slot suction can effectively suppress the aerodynamic forces of the model. The maximum reduction of aerodynamic forces occurs at Q = 1, with the time-mean drag, fluctuating drag, and fluctuating lift reduced by 3.75%, 19.08%, 40.91%, respectively. For Q = 3, all aerodynamic forces are still smaller than those for Q = 0 (uncontrolled case), but obviously higher than those for Q = 1. The involved control mechanism is successfully revealed, based on the comparison of the flow around cylinder free end and the near wake for the three tested Q values.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.2 no.3
• /
• pp.80-86
• /
• 2000

Eddy covariance method is widely used in measuring vertical fluxes of mass and energy between the atmosphere and the biosphere. In this method, scalar concentration is measured with either open-path or closed-path sensors. For the latter, fluctuations of scalar concentration are attenuated as the sample travels through a long tube, resulting in flux loss. To quantify this tube attenuation, water vapor concentrations measured with both closed-path and open-path sensors were analyzed. Our statistical analysis showed that the power spectral density obtained from the closed-path sensor was different from that from the open-path sensor in the frequency range of > 0.5 Hz. The loss of water vapor flux due to tube attenuation was < 5% during midday. At nighttime, however, the flux loss increased significantly because of the low wind speeds and the weak turbulence sources. Theoretical calculation for the tube attenuation showed a small bias in high frequency range probably because of the interaction of sticky water vapor with a tube wall.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.5
• /
• pp.405-412
• /
• 2015

This paper represents a filter design to estimate the airspeed of a spin-stabilized, trajectory-correctible artillery ammunition. Due to the limited power and space in operational point of view, the airspeed sensor is not installed, and thus the airspeed need to be estimated using limited sensor measurements. The only IMU measurements(three-axis specific forces and angular rates) are used in this application. The extended Kalman filter algorithm is applied since a linear filter can not cover the its wide operational range in airspeed and altitude. In the implementation of the EKF, the state and measurement equations are transformed into the no-roll frame for simple form of Jacobian matrix. The simulation study is conducted to evaluate the performance of the filter under various environment conditions of sensor noise and wind turbulence. In addition, the effect of the choice in filter design parameters, i.e. process error covariance matrices is analyzed on the performance of the estimation of airspeed and angular rates.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.259-267
• /
• 2017

Flow fields around a KARI-11-180 airfoil, SDM and transonic body are numerically simulated by using an unstructured meshes based compressible flow solver developed at KAIST. RANS equations are solved to analyse the flow fields and Roe's FDS method is adopted to evaluate convective fluxes. Turbulence effect of the flow fields is modeled by a SA model, SST model and ${\gamma}-{\widetilde{Re}}_{{\theta}t}$ model. It is found that smaller drag coefficients are predicted for the KARI-11-180 airfoil when a transition phenomenon is considered and small deviations exist between CFD and EFD results. For the SDM, flow separation is observed at a leading edge and calculated aerodynamic properties show similar tendencies to experimental results. A shock wave on main wings of the transonic body is successfully captured by the present flow solver at a Mach number 0.9. Estimated pressure profiles by means of the present CFD method also agree well with those of wind tunnel results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.8 no.5
• /
• pp.442-449
• /
• 1984

Turbulent jet flow has been studied in many ways; a plane jet, a rectangular jet, an annular jet, a round jet, a wall jet, a parallel jet, a valve jet, a cross jet, a slit jet and etc. In this report, a 45.deg. cross jet flow was tried by using two same dimensioned nozzels(dia..phi.20)which were set up at the exit of the subsonic wind tunnel. Each jet flows to the direction of 22.5.deg. to the axis of downstream of the mixed flow. The centerline of each jet meets at the distance of 217.3mm and their mixing flow could be imagined to develop beyond that distance, so the measurement was effectuated at X/X$_{0}$=1.2-1.5. The section of the mixed flow a elliptic circle which is formed by the 22.5.deg. inclined flows to the X direction. This experimental study aimed at the investigation of the turbulent mixing process of two jets; the mean velocities, the turbulent shear stresses, the correlation coefficients, and the momentum were respectively measured. The mean velocity distribution profiles of the down-stream component measured in the Y direction coincide well with the empirical equation of Gortler and those measured in the Z direction agree with the equation of H. Schlichting. Other mean velocities V over bar and W over bar components were randomly distributed. The higher values with same order of the intensity of turbulence were largely distributed at the central part of the flow. The momentum was decreased up to 70% by the shock losses and the development of intense turbulences, but it kept its value constantly beyond X/d=14. Two-channel hot-wire anemometer systems (model 1050 series), X-type hot-wire made of tungsten (dia. .phi.e.mu.m, long 3mm, model 0252 T5), a computer(model HP 9845B0, and a plotter (model HP 9872C) were used for the experiments and the analyses.s.