• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.984-994
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a spectral, static error analysis. To investigate the effect of numerical dissipation on LES solutions, power spectra of discretization errors are evaluated for isotropic turbulence models in both continuous and discrete wavevector spaces. Contrary to the common belief, the aliasing errors from upwind-biased schemes are larger than those from comparable non-dissipative schemes. However, this result is the direct consequence of the definition of the power spectral density of the aliasing error, which poses the limitation of the static error analysis for upwind schemes.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.995-1006
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a dynamic analysis. Large eddy simulation of isotropic turbulence is performed with various dissipative and non-dissipative schemes to investigate the effect of numerical dissipation on the resolved solutions. It is shown by the present dynamic analysis that upwind schemes reduce the aliasing error and increase the finite differencing error. The existence of optimal upwind scheme that minimizes total numerical error is verified. It is also shown that the finite differencing error from numerical dissipation is the leading source of numerical errors by upwind schemes. Simulations of a turbulent channel flow are conducted to show the existence of the optimal upwind scheme.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.3
• /
• pp.22-30
• /
• 1996

In-cylinder flow of a purpose-built small HSDI Hydra Diesel engine was investigated by laser Doppler velocimetry(LDV) during induction and compression processes. The flow was quantified in terms of ensemble-averaged axial and swirl velocities, normalized by the mean piston speed, at a plane located 12mm from the cylinder head and corresponding to the mid-plane of the diametrically-opposed quartz windows at an enigne speed of 1000rpm. The formation of toroidal vortices during the intake process and the evolution and decay of swirl motion during the compression process were observed. Turbulence at around TDC of compression became homogeneous and isotropic.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.13 no.12
• /
• pp.1184-1195
• /
• 2001

Analysis of fluid flow in rectangular ducts has been conducted since it has a wide application. The purpose is to provide experimental data for the comparison with computational results. Velocity distributions inside a rectangular duct with $90^{\circ}$ mitered elbow are measured by 5W laser doppler velocity meter for Reynolds numbers of 4,049, 8,104, and 12,186. Flow rates obtained by the integration of measured velocity profile at three cross-sections, which are inlet, middle section after the elbow, and outlet, have errors less than 0.9% among them. Turbulent fluctuation components in two directions are found to have almost similar magnitude each other at a certain location due to the isotropic characteristic of turbulence.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.1
• /
• pp.42-53
• /
• 2008

The two dimensional PIV (particle image velocimetry) measurement technique is applied to water flow in a narrow cavitation tunnel. The nearly homogeneous and isotropic turbulent flows are generated by the honeycomb installed in the tunnel and visualized with a PIV technique. The velocities in the measurement plane at the tunnel centerline 184cm downward from the honeycomb were measured and calculated by an image correlation technique. The turbulent properties are evaluated and each term in the turbulent kinetic energy equation is calculated for the conditions with different internal pressures. Lowering the internal pressure gives an effect on the turbulent flow due to growing bubbles which are resolved in the water. The turbulent kinetic energy in the measurement plane is decayed much slower than those of other research results carried out with wind tunnels. With decreasing the tunnel internal pressures the turbulent intensities are increased about 1.5 times and the anisotropic tendency is also increased.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.4
• /
• pp.519-526
• /
• 1996

A 2D-LDV system was employed to investigate the flow field characteristics in fully developed drag reducing turbulent channel flows. The additive used in this study was Habon-G which showed splendid drag reduction effect and minimum mechanical degradation trend in the closed flow circulation loop. In order to have better understanding of the drag reduction mechanism, the instantaneous velocities were carefully measured under various experimental conditions and the flow characteristics including time-averaged velocity, turbulent intensity and Reynolds shear stresses were carefully assessed. The time-averaged velocity profiles of surfactant flows showed more parabolic shape(typically shown in a laminar flow) together with significant suppression of turbulent production, yielding the shear induced micelle structure orienting in the flow direction due to its isotropic characteristics. Especially it was observed that the maximum intensity for drag reducing flows was shifted away from the wall and that the streamwise and normal turbulent intensities were strongly altered. This phenomenon strongly suggests that the viscous sublayer becomes thicker with addition of surfactant. Turbulent momentum transport was drastically suppressed across the whole drag reducing channel flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.8 s.227
• /
• pp.902-909
• /
• 2004

A steady-state Reynolds averaged Navier-Stokes simulation was conducted to investigate the distribution of the Reynolds stress tensor inside tip leakage vortex of a linear compressor cascade. Two different inlet flow angles ${\beta}=29.3^{\circ}$(design condition) and $36.5^{\circ}$(off-design condition) at a constant tip clearance size of $1\%$ blade span were considered. Classical methods of solid mechanics, applied to view the Reynolds stress tensor in the principal direction system, clearly showed that the high anisotropic feature of turbulent flow field was dominant at the outer part of tip leakage vortex near the suction side of the blade and endwall flow separation region, whereas a nearly isotropic turbulence was found at the center of tip leakage vortex. There was no significant difference in the anisotropy of the Reynolds normal stresses inside tip leakage vortex between the design and off-design condition.

• Journal of the Korean Society of Visualization
• /
• v.17 no.2
• /
• pp.48-57
• /
• 2019

A straightforward strategy for particle image velocimetry (PIV) interrogation and post-processing has been proposed, aiming at reducing errors and clarifying vortex structures. The interrogation window size should be kept small to reduce bias error and improve spatial resolution. A spatial filter is then applied to the velocity field to reduce random error and clarify flow structure. The performance of three popular spatial filters were assessed: box filter, median filter, and local quadratic polynomial regression filter. In order to quantify random uncertainty, the image matching (IM) method is applied to an experimental dataset of homogeneous and isotropic turbulence (HIT) obtained by 2D-PIV. We statistically analyze the uncertainty propagation through the spatial filters, and verify the reduction in random uncertainty. Moreover, we illustrate that the spatial filters help clarify vortex structures using vortex identification criteria. As a result, PIV random uncertainty was reduced and the vortex structures became clearer by spatial filtering.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.24 no.5
• /
• pp.219-232
• /
• 2020

The slip-weakening model developed by Ohnaka and Yamashita is extended over the breakdown zone by equating the scaling relationships for the breakdown zone and the whole rupture area. For the extension, the study uses the relationship between rupture velocity and radiation efficiency, which was derived in the theory of linear elastic fracture mechanics, and the definition of f_max given in the specific barrier model proposed by Papageorgiou and Aki. The results clearly show that the extended scaling relationship is governed by the ratio of rupture velocity to S wave velocity, and the velocity ratio can be determined by the ratio of characteristic frequencies of a Fourier amplitude spectrum, which are corner frequency, f_c, and source-controlled cut-off frequency, f_max, or vice versa. The derived relationship is tested by using the characteristic frequencies extracted from previous studies of more than 130 shallow crustal events (focal depth less than 25 km, M_W 3.0~7.5) that occurred in Japan. Under the assumption of a dynamic similarity, the rupture velocity estimated from f_max/f_c and the modified integral timescale give quite similar scale-dependence of the rupture area to that given by Kanamori and Anderson. Also, the results for large earthquakes show good agreement to the values from a kinematic inversion in previous studies. The test results also indicate the unavailability of the spectral self-similarity proposed by Aki because of the scale-dependent rupture velocity and the rupture velocity-dependent f_max/f_c; however, the results do support the local similarity asserted by Ohnaka. It is also remarkable that the relationship between the rupture velocity and f_max/f_c is quite similar to Kolmogorov's hypothesis on a similarity in the theory of isotropic turbulence.