• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.434-440
• /
• 2011

Oscillating airfoil haw been challenged for the dynamic stalls of airfoil am wind turbines at high angle of attach. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance am safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed for the oscillating airfoil at high angle of attack around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.2 and Reynolds number of $1.2{\times}10^4$. The lift, drag, pressure distribution, etc. are analyzed according to the pitching oscillation. Unsteady velocity field, periodic vortex shedding, the unsteady pressure distribution, and the acoustic fields are analyzed. The effects of these unsteady characteristics in the aerodynamic coefficients are analyzed.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.13-16
• /
• 2008

After the collapse of the Tacoma bay bridge at Tacoma Washington, the accurate prediction of aerodynamics became crucial to the sound design of bridges. CFD(Computational Fluid Dynamics) becomes important tool for the prediction on wind effects on the bridge due to the recent development of CFD. The usage of CFD is further prompted by the advantages in using CFD, such as low-cost and fast feed-back of design. In this paper, an unsteady compressible Reynolds averaged Navier-Stokes code is used for the computation of the flow over bridges. Coakley's ��q-${\omega}$ �� two-equation turbulence model is used for the turbulent eddy viscosity. For accurate and stable computations, the local preconditioning method is adapted to the code. Aerodynamic characteristics of a couple bridges are presented to show the validity and the accuracy of the method.

• Ocean and Polar Research
• /
• v.24 no.2
• /
• pp.135-146
• /
• 2002

The upper ocean in the western equatorial Pacific warm pool during TOGA-COARE IMET IOP was simulated using a one-dimensional turbulence closure ocean mixed-layer model, which considered recent observations, such as the remarkable enhancement of turbulent kinetic energy near the ocean surface. The shoaling/deepening of the mixed layer and warming/cooling subsurface water in the model were in reasonable agreement with the observations. There was a significant improvement in simulating the cooling trend of the sea surface temperature under a westerly wind burst with heavy rainfall over previous simulations using bulk mixed-layer models. By contrast the simulated sea surface salinity (SSS) departed significantly from the observed SSS, especially during a westerly burst and the subsequent restratification period, which might be due to 3-D control processes, such as downwelling/upwelling or advection.

• Journal of the Korean Society of Visualization
• /
• v.2 no.1
• /
• pp.32-38
• /
• 2004

A new PIV technique was developed for decreasing optical error which was created during simultaneous measurements of velocity fields at a wall-normal plane and wall-parallel plane by a plane PIV and a Stereo-PIV systems. For experimental study, two different Reynolds number based on momentum thickness, Re$_{ =514 and 934 were generated in a blowing type wind tunnel under the condition of zero pressure gradient. The two Nd:Yag laser systems and three CCD cameras were synchronized. to obtain instantaneous velocity fields at the same time. To avoid optical noise at the crossing line by the two laser light sheets, a new optical arrangement using polarization was applied. The obtained velocity fields show the existence of hairpin packet structure vividly and the idealized hairpin vortex signature is confirmed by experiment.

• Proceedings of KOSOMES biannual meeting
• /
• 2001.10a
• /
• pp.137-144
• /
• 2001

A Plankton ecosystem model was developed to investigate effects of hydrodynamic processes including advection and diffusion on size-structured phytoplankton dynamics in the mesohaline zone of the York River estuarine system, Virginia, USA. The model included 12 state variables representing the distribution of carbon and nutrients in the surface mixed layer. Groupings of autotrophs and heterotrophs were based on cell site and ecological hierarchy Forcing functions included incident radiation, temperature, wind stress, mean How and tide which includes advective transport and turbulent mixing. The ecosystem model was developed in FORTRAN using differential equations that were solved using the 4th order Runge-Kutta technique. The model showed that microphytoplankton blooms during winter-spring resulted from a combination of vertical advection and diffusion of phytoplankton cells rather than in-situ production in the lower York River estuary.

• Proceedings of the KSR Conference
• /
• 2001.05a
• /
• pp.123-128
• /
• 2001

Numerical simulation is conducted to clarify the heat transfer and fluid flow characteristics of brake shoe for freight car. High order up-wind scheme for governing equations, k-epsilon turbulent model and SIMPLEC algorithm based on finite volume method are used to solve the physical shoe model. The governing equations are solved by TDMA(Tri-Diagonal Matrix Algorithm) with line-by-line method and block correction. From the results of simulation, the characteristics of cooling pattern is strongly affected by the velocity of train and the material of shoe. The face lift of shoe affects on the temperature distribution of rear surface of shoe as well as the front surface of that. Due to the grooves in shoe, it will be expected to cool the frictional heat and result in the reduction of maintenance efforts.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.313-324
• /
• 1992

Two-dimensional Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system, then mapped into a generalized body-fitted coordinate system. All direction of viscous terms are incorporated and turbulent effects are modeled using the Baldwin-Lomax algebraic model. Equation are discretized using finite difference method on the C-type grids and solved using implicit LU-ADI decomposition scheme. Calculations are made at a VKI turbine cascade flow in a transonic wind-tunnel and compared to experimental data. Present numerical scheme is shown to be in good agreement with the previous experimental results and simulates the two-dimensional viscous flow phenomena.

• Wind and Structures
• /
• v.1 no.1
• /
• pp.59-66
• /
• 1998

The vibrations of bodies subjected to fluid flow can cause modifications in the flow conditions, giving rise to interaction forces that depend primarily on displacements and velocities of the body in question. In this paper the linearized equations of motion for bodies of arbitrary prismatic or cylindrical cross-section in two-dimensional cross-flow are presented, considering the three degrees of freedom of the body cross-section. By restraining the rotational motion, equations applicable to circular tubes, pipes or cables are obtained. These equations can be used to determine stability limits for such structural systems when subjected to non uniform cross-flow, or to evaluate, under the quasi static assumption, their response to vortex or turbulent excitation. As a simple illustration, the stability of a pipe subjected to a bidimensional flow in the direction normal to the pipe axis is examined. It is shown that the approach is extremely powerful, allowing the evaluation of fluid-structure interaction in unidimensional structural systems, such as straight or curved pipes, cables, etc, by means of either a combined experimental-numerical scheme or through purely numerical methods.

• Journal of the Korean earth science society
• /
• v.37 no.7
• /
• pp.448-464
• /
• 2016

The purpose of this study is to optimize the parameters related to consensus coherency within the PCL 1300, the operating program of wind profiler, from a validation of wind data between rawinsonde and wind profiler at Chupungryeong ($36^{\circ}13^{\prime}$, $127^{\circ}59^{\prime}$) site in Korea. It is then to analyze the diurnal and seasonal characteristics of the turbulence energy dissipation rate (${\varepsilon}$) in clear and rainy days from March 2009 to February 2010. In comparison of the wind data between wind profiler and rawinsonde during April 22-23, 2010, it was shown in a big error more than $10ms^{-1}$ over the height of 3,000 meters in the zonal (u) and meridional (v) wind components. When removing more than $10ms^{-1}$ in each wind speed difference of u an v components between the two instruments, the correlation coefficients of these wind components were 0.92 and 0.88, respectively, and the root mean square errors were 3.07 and $1.06ms^{-1}$. Based on these results, when the data processing time and the minimum available data within the PCL 1300 program were adjusted as 30 minutes and 60%, respectively, the bias errors were small. In addition, as a result of an analysis of sensitivity to consensus coherency of u and v components within the PCL1300 program, u components were underestimated in radial coherency, instantaneous and winbarbs coherency, whereas v components were overestimated. Finally by optimizing parameters of the PCL1300 program, the diurnal and seasonal means of ${\varepsilon}$ at each height were higher in rainy days than those in clear days because of increasing in the vertical wind speed due to upward and downward motions. The mean ${\varepsilon}$ for clear and rainy days in winter was lower than those of other seasons, due to stronger horizontal wind speed in winter than those in other seasons. Consequently, when the turbulence energy dissipation rates in the vertical wind speed of more than ${\pm}10cm\;s^{-1}$ were excluded for clear and rainy days, the mean ${\varepsilon}$ in rainy days was 6-7 times higher than that in clear days, but when considering them, it was 4-5 times higher.