• Wind and Structures
• /
• v.11 no.4
• /
• pp.275-289
• /
• 2008

Large eddy simulations have been performed within and over different types of urban building arrays. This paper adopted three dimensionless parameters, building frontal area density (${\lambda}_f$) the variation degree of building height (${\sigma}_h$), and the staggered degree of building range ($r_s$), to study the systematic influence of building spacing, height and layout on wind and turbulent characteristics. The following results have been achieved: (1) As ${\lambda}_f$ decrease from 0.25 to 0.18, the mean flow patterns transfer from "skimming" flow to "wake interference" flow, and as ${\lambda}_f$ decrease from 0.06 to 0.04, the mean flow patterns transfer from "wake interference" flow to "isolated roughness" flow. With increasing ${\lambda}_f$, wind velocity within arrays increases, and the vortexes in front of low buildings would break, even disappear, whereas the vortexes in front of tall buildings would strengthen and expand. Tall buildings have greater disturbance on wind than low buildings do. (2) All the wind velocity profiles and the upstream profile converge at the height of 2.5H approximately. The decay of wind velocity within the building canopy was in positive correlation with ${\lambda}_f$ and $r_s$. If the height of building arrays is variable, Macdonald's wind velocity model should be modified through introducing ${\sigma}_h$, because wind velocity decreases at the upper layers of the canopy and increases at the lower layers of the canopy. (3) The maximum of turbulence kinetic energy (TKE) always locates at 1.2 times as high as the buildings. TKE within the canopy decreases with increasing ${\lambda}_f$ and $r_s$ but the maximum of TKE are very close though ${\sigma}_h$ varies. (4) Wind velocity profile follows the logarithmic law approximately above the building canopy. The Zero-plane displacement $z_d$ heighten with increasing ${\lambda}_f$, whereas the maximum of and Roughness length $z_0$ occurs when ${\lambda}_f$ is about 0.14. $z_d$ and $z_0$ heighten linearly with ${\sigma}_h$ and $r_s$, If ${\sigma}_h$ is large enough, $z_d$ may become higher than the average height of buildings.

• Journal of Korea Water Resources Association
• /
• v.47 no.8
• /
• pp.703-715
• /
• 2014

This study aims to investigating the effect of Schmidt number (${\sigma}_c$) on sediment suspension and hydrodynamics calculation. The range of ${\sigma}_c$ is also studied based on the flux Richardson number ($Ri_f$) and gradient Richardson number ($Ri_g$). Numerical experiments are carried out by 1 dimensional vertical model. Both cohesive and non-cohesive sediments are tested under the conditions of pure current and oscillatory flow. The turbulence damping effect due to sediment suspension is examined considering ${\sigma}_c$ as a constant for the damping effect. The results of this study show the consistent effect of ${\sigma}_c$ on sediment suspension regardless of hydrodynamic condition. It is also found that the model overestimates the flow velocity and turbulent kinetic energy when the damping effect is not considered. Under the conditions of $Ri_f$ and $Ri_g$ causing density stratification, it is known that the vertical mixing of sediment is reasonably calculated in the range of ${\sigma}_c$ from 0.3 to 0.5.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.11a
• /
• pp.417-423
• /
• 2007

Oxy-gasification or oxygen-blown gasification, enables a clean and efficient use of coal and opens a promising way to CO2 capture. The coal gasification process of a slurry feed type, entrained-flow coal gasifier was numerically predicted in this paper. The purposes of this study are to develop an evaluation technique for design and performance optimization of coal gasifiers using a numerical simulation technique, and to confirm the validity of the model. By dividing the complicated coal gasification process into several simplified stages such as slurry evaporation, coal devolatilization, mixture fraction model and two-phase reactions coupled with turbulent flow and two-phase heat transfer, a comprehensive numerical model was constructed to simulate the coal gasification process. The influence of turbulence on the gas properties was taken into account by the PDF (Probability Density Function) model. A numerical simulation with the coal gasification model is performed on the Conoco-Philips type gasifier for IGCC plant. Gas temperature distribution and product gas composition are also presented. Numerical computations were performed to assess the effect of variation in oxygen to coal ratio and steam to coal ratio on reactive flow field. The concentration of major products, CO and H2 were calculated with varying oxygen to coal ratio (0.2-1.5) and steam to coal ratio(0.3-0.7). To verify the validity of predictions, predicted values of CO and H2 concentrations at the exit of the gasifier were compared with previous work of the same geometry and operating points. Predictions showed that the CO and H2 concentration increased gradually to its maximum value with increasing oxygen-coal and hydrogen-coal ratio and decreased. When the oxygen-coal ratio was between 0.8 and 1.2, and the steam-coal ratio was between 0.4 and 0.5, high values of CO and H2 were obtained. This study also deals with the comparison of CFD (Computational Flow Dynamics) and STATNJAN results which consider the objective gasifier as chemical equilibrium to know the effect of flow on objective gasifier compared to equilibrium. This study makes objective gasifier divided into a few ranges to study the evolution of the gasification locally. By this method, we can find that there are characteristics in the each scope divided.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

• The KSFM Journal of Fluid Machinery
• /
• v.18 no.6
• /
• pp.63-70
• /
• 2015

This paper proposes a novel airfoil named "KA2" for the blade of the wind turbine systems. Dynamic loads characteristics are analyzed and compared using aerodynamic data of ten airfoils including the proposed airfoil. The blade is divided into the sixteen elements in the longitudinal direction of the blade for applying the Blade Element Method Theory (BEMT) method, and in each element, torque, thrust, and pitching moment are calculated using turbulent time varying wind speed and aerodynamic data of each wing. Additionally, each force and torque is accumulated in the whole region of the blade for the estimation of representative values. The magnitude of such forces is comparatively analyzed for different airfoils. The angle of attack is constant below the rated wind speed due to the fact that the tip speed ratio is kept at the constant value, and it increases in the region of over rated wind speed as the tip speed ratio decreasing with constant rated rpm and increasing wind speed. Such increase in the angle of attack causes the changes of the force acting on the airfoil with different characteristics of lift and drag in the stall region of each different airfoil. Even though the mean wind speed is in the rated speed in a given time, because of the turbulence, it has either the over rated or under rated speed most of the time. Furthermore, the dynamic properties of each force are analyzed in this rated wind speed in order to objectively understand the dynamic properties of the blades which are designed based on the different airfoils. These dynamic properties are also compared by the standard deviation of time varying characteristics. Moreover, the output characteristics of the wind turbine are investigated with different airfoils and wind speeds. Based on these investigations, it was revealed that the proposed airfoil (KA2) is well applicable to the blade with passive pitch control system.

• Journal of Korean Society of Environmental Engineers
• /
• v.31 no.1
• /
• pp.29-34
• /
• 2009

This study was performed under four operational conditions for nitrogen removal in metal finishing wastewater. The conditions include electrode gap, reducing agent, the recycling of treated wastewater in 1st step and the simultaneous treatment of nitrate and other materials. Result showed that the removal efficiency of $NO_3{^-}-N$ was highest at the electrode gap of 10 mm. As the electrode gap was shorter than 10 mm, the removal efficiency of $NO_3{^-}-N$ decreased due to increasing in concentration polarization on electrode. And, in case that the electrode gap was longer than 10 mm, the removal efficiency of $NO_3{^-}-N$ increased with an increase in energy consumption. Because hydrogen ions are consumed when nitrate is reduced, reducing reaction of nitrate was effected more in acid solution. As 1.2 excess amount of zinc was injected, the removal efficiency of $NO_3{^-}-N$ increased due to increasing in amount of reaction with nitrate. As the effluent from 1st step in the reactor was recycled into the 1st step, the removal efficiency of $NO_3{^-}-N$ increased. Because the zinc were detached from the cathode and concentration-polarization was decreased due to formation of turbulence in the reactor. The presence of $NH_4{^+}-N$ did not affect the removal efficiency of $NO_3{^-}-N$ but the addition of heavy metal decreased the removal efficiency of $NO_3{^-}-N$. As chlorine is enough in wastewater, the simultaneous treatment of nitrate and ammonia nitrogen may be possible. The problem that heavy metal decrease the removal efficiency of $NO_3{^-}-N$ may be solved by increasing current density or using front step of electrochemical process for heavy metal removal.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.4
• /
• pp.648-656
• /
• 2022

As environmental regulations such as the International Maritime Organization (IMO)'s strategy to reduce greenhouse gases(GHG) are strengthened, technology development such as eco-friendly ships and alternative fuels is expanding. As part of this, ship propulsion technology using energy reduction and wind propulsion technology is emerging, especially in shipping companies and shipbuilders. By securing wind propulsion technology and introducing empirical research into shipbuilding and shipping, a high value-added market using eco-friendly technology can be created. Moreover, by reducing the fuel consumption rate of operating ships, GHG can be reduced by 6-8%. Rotor Sail (RS) technology is to generate a hydrodynamic lift in the vertical direction of the cylinder when the circular cylinder rotates at a constant speed and passes through the fluid. This is called the Magnus effect, and this study attempted to propose a plan to increase propulsion efficiency through a numerical analysis study on turbulence flow characteristics around RS, a wind power assistance propulsion system installed on a ship. Therefore, C_L and C_D values according to SR and AR changes were derived as parameters that affect the aerodynamic force of the RS, and the flow characteristics around the rotor sail were compared according to EP application.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.2
• /
• pp.93-103
• /
• 2023

Warships widely spread numerous chaffs using a blast, which form chaff clouds that create false radar cross-sections to deceive enemy radars. In this study, we established a numerical framework based on a one-way coupling of computational fluid dynamics and discrete element method to simulate the spatiotemporal distribution of chaff clouds for warships in the air. Using the framework, we investigated the effects of wind, initial chaff cartridge angle, and blast pressure on the distribution of chaff clouds. We observed three phases for the chaff cloud diffusion: radial diffusion by the explosion, omnidirectional diffusion by turbulence and collision, and gravity-induced diffusion by the difference in the fall speed. The wind moved the average position of the chaff clouds, and the diffusion due to drag force did not occur. The direction of radial diffusion by the explosion depended on the initial angle of the cartridge, and a more vertical angle led to a wider distribution of the chaffs. As the blast pressure increased, the chaff clouds spread out more widely, but the distribution difference in the direction of gravity was not significant.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.9 no.2
• /
• pp.109-119
• /
• 2006

This study investigated the effect of artificially enhanced mesozooplankton on the phytoplankton dynamics during fall blooming period using a mesocosm in Jangmok bay located in the Southern Sea of Korea in 2001. The four bags with 2,500 liter seawater containment were directly filled with the ambient water. And then, abundances of mesozooplankton in two experimental bags were treated 6 times higher than those in control bags by towing with net($300{\mu}m$) through the ambient water. Phytoplankton community between control and experimental bags were not significantly different in terms of chlorophyll-a(chl-a) concentration and standing crop (one-way ANOVA, p>0.05) during the study period. Initial high standing crop and chl-a concentration of phytoplankton drastically decreased and remained low until the end of the experiment in all bags. Diatoms, accounting for most of the phytoplankton community, consisted of Skeletonema costatum, Pseudo-nitzschia seriata, Chaetoceros curvisetus, Ch. debilis, Cerataulina pelagica, Thalassiosira pacifica, Cylindrotheca closterium, and Leptocylindrus danicus. Noctiluca scintillans dominated the temporal variation of mesozooplankton abundances, which peaked on Day 10 in the control and experimental bags, while the next dominant copepods showed their peak on Day 7. Shortly after mesozooplankton addition, copepod abundance in the experimental bags was obviously higher than that in the control bags on Day 1, however, it became similar to that in the control bags during the remnant period. It was supported by the higher abundance and length of both ctenophores and hydromedusae in experimental bags relative to the control bags. However, the cascading trophic effect, commonly leading to re-increase of phytoplankton abundance, was not found in the experimental bags, indicating that copepods were not able to control the phytoplankton in the bags based on the low grazing rate of Acartia erythraea. Besides that, rapidly sunken diatoms in the absence of natural turbulence as well as N-limited condition likely contributed the no occurrence of re-increased phytoplankton in the experimental bags.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.10 no.3
• /
• pp.82-93
• /
• 2008

In order to ensure a standardized data analysis of the eddy covariance measurements, Hong and Kim's quality control program has been updated and used to process eddy covariance data measured at two levels on the main flux tower at Gwangneung site from January to May in 2005. The updated program was allowed to remove outliers automatically for $CO_2$ and latent heat fluxes. The flag system consists of four quality groups(G, D, B and M). During the study period, the missing data were about 25% of the total records. About 60% of the good quality data were obtained after the quality control. The number of record in G group was larger at 40m than at 20m. It is due that the level of 20m was within the roughness sublayer where the presence of the canopy influences directly on the character of the turbulence. About 60% of the bad data were due to low wind speed. Energy balance closure at this site was about 40% during the study period. Large imbalance is attributed partly to the combined effects of the neglected heat storage terms, inaccuracy of ground heat flux and advection due to local wind system near the surface. The analysis of wind direction indicates that the frequent occurrence of positive momentum flux was closely associated with mountain valley wind system at this site. The negative $CO_2$ flux at night was examined in terms of averaging time. The results show that when averaging time is larger than 10min, the magnitude of calculated $CO_2$ fluxes increases rapidly, suggesting that the 30min $CO_2$ flux is influenced severely by the mesoscale motion or nonstationarity. A proper choice of averaging time needs to be considered to get accurate turbulent fluxes during nighttime.