• Journal of Ocean Engineering and Technology
• /
• v.10 no.3
• /
• pp.50-61
• /
• 1996

The breaking phenomenon of regular periodic waves generated by a numerical wave maker is simulated by finite-difference method which can cope with strong interface motions. The air and water flows are simultaneously solved in the time-marching solution procedure for the Navier-Stokes equation. A density-function technique is devised for the implemenation of the interface conditions. The accuracy is examined and applied to the simulation of two-dimensional breaking phenomena of periodic gravity waves.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.46-53
• /
• 2016

Direct numerical simulation is conducted to observe the behavior of microbubbles in isotropic turbulence. Navier-Stokes equation and the motion of equation for microbubbles are solved with periodic boundary condition in a cube domain. Vorticity contour, enstrophy ratio, relative reduction of bubble rise velocity, and the closest distance of particles are investigated for various Stokes numbers and gravity factors to understand clustering of microbubbles. Also, clustering due to the effect of the lift force is investigated.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.441-444
• /
• 2006

Helicity in isotropic turbulence was well known to have intermittent fashion in their statistics. But its exact explanation about the onset of intermittency of helicity in turbulence did not give clearly yet. Most probable causes comes from the vortical motion of the fluids. Distribution of the angle between fluid velocity and vorticity have alignment tendency. This may be a clue to investigate intermittency of helicity. In this study, we aim to review and establish approaches to reveal the mechanism and the origin of intermittency of helicity in the isotropic turbulence. To do those work, we look for some quantities like helicity, enstrophy, acceleration and its flatness. And also correlations among them are sought.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.11a
• /
• pp.110-113
• /
• 2007

LES code was developed and verified to investigate the turbulence evolution in the gas generator equipped with turbulence ring. We compared mixing enhancement of various turbulence rings having different shape and size. To qualify the mixedness, we used two perturbation values of energy and enstrophy.

• Journal of Korean Society for Atmospheric Environment
• /
• v.19 no.5
• /
• pp.569-581
• /
• 2003

A method of synthetic time series generation was developed and applied to the simulation of homogeneous turbulence in a periodic 3 - D box and the hourly wind data simulation. The method can simulate almost exact sample auto and cross correlations of multiple time series and control non-Gaussian distribution. Using the turbulence simulation, influence of correlations, non-Gaussian distribution, and one-direction anisotropy on homogeneous structure were studied by investigating the spatial distribution of turbulence kinetic energy and enstrophy. An hourly wind data of Typhoon Robin was used to illustrate a capability of the method to simulate sample cross correlations of multiple time series. The simulated typhoon data shows a similar shape of fluctuations and almost exactly the same sample auto and cross correlations of the Robin.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.9
• /
• pp.1032-1041
• /
• 2004

In this study we decompose the 3-dimensional velocity field of isotropic turbulent flow into the coherent and the incoherent structure using the discrete wavelet. It is shown that the coherent structure, 3% wavelet modes, has 98% energy and 88% enstrophy and its statistical characteristics are almost same as the original turbulence structure. And it is confirmed that the role of the coherent structure is that it produces the turbulent kinetic energy at the inertia range then transfers energy to the dissipation range. The incoherent structure, with residual wavelet modes, is uncorrelated and has the Gaussian probability density function but it dissipates the kinetic energy in dissipation range. On the procedure, we propose a new but easy way to get the threshold by applying the energy partition percentage concept about coherent structure. The vorticity field extracted from the wavelet-decomposed velocity field has the same structure as the result of the precedent studies which decomposed vorticity field directly using wavelet. Therefore it has been shown that velocity and vorticity field are on the interactive condition.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.819-822
• /
• 2002

It is necessary for the numerical simulation of 3-dimensional incompressible isotropic decaying turbulence to construct 3-dimensional initial velocity field which resembles the fully developed turbulence. Although the previous velocity field generation method proposed by Rogallo(1981) satisfies continuity equation and 3-dimensional energy spectrum, it has limitation, as indicated in his paper, that it does not produce the higher velocity moments(e. g. velocity derivative skewness) characteristic of real turbulence. In this study, a new velocity field generation method which is able to control velocity derivative skewness of initial velocity field is proposed. Brief descriptions of the new method and a few parameters which is used to control velocity derivative skewness are given. A large eddy simulation(LES) of isotropic decaying turbulence using dynamic subgrid-scale model is carried out to evaluate the performance of the initial velocity field generated by the new method. It was shown that the resolved turbulent kinetic energy decay curve and the resolved enstrophy decay curve from the initial field of new method were more realistic than those from the initial field of Rogallo's method. It was found that the dynamic model coefficient from the former was initially half the stationary value and experienced relatively short transition period, though that from the latter was initially zero and experienced relatively longer transition period.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.391-398
• /
• 2008

LES analysis was conducted with in-house CFD code to investigate the turbulence evolution and interaction due to turbulence ring and splash plate in the gas generator. Though chemical reaction was not accounted for, the results can be useful in determining the turbulence characteristics generated by ring and plate. The calculation results show that the installation of turbulence ring can introduce additional turbulences and improve turbulent mixing in the downstream flow. However, the addition of splash plate in the downstream of TR brings totally different shape of perturbation energy and enstrophy distribution for turbulent mixing. This enhancement can be done by the formation of the intensively strong vorticity production and mixing behind the plate. Pressure drop was found to be a reasonable level of about 1% or less of initial pressure in all calculation cases. Also, calculation results revealed that the variation of TR shape and intrusion length did not change the characteristics of turbulent mixing in the chamber. Even though the effect of installation location of splash plate on the turbulent mixing is not investigated yet, calculation results conclude the addition of splash plate leads to the increase in turbulent mixing with an acceptable pressure drop.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.12
• /
• pp.1171-1179
• /
• 2008

LES analysis was conducted with in-house CFD code to investigate the turbulence evolution and interaction due to turbulence ring and splash plate in the gas generator. The calculation results show that the installation of turbulence ring can introduce additional turbulences and significantly improve turbulent mixing in the downstream flow. However, the addition of splash plate in the downstream of TR(Turbulence Ring) brings totally different shape of perturbation energy and enstrophy distribution into turbulent mixing. This enhancement can be done by the formation of the intensively strong vorticity and mixing behind the plate. Pressure drop was found to be a reasonable level of about 1% or less of initial pressure in all calculation cases. Also, calculation results revealed that the variation of shape and intrusion length of TR did not greatly affect the characteristics of turbulent mixing in the chamber. Even though the effect of installation location of splash plate on the turbulent mixing was not investigated yet, calculation results conclude the addition of splash plate leads to the increase in turbulent mixing with an acceptable pressure drop.