• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2002.04a
• /
• pp.35-40
• /
• 2002

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc compressor cascades and the results are compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE algorithm adopt pressure weighted interpolation method for non-staggered grid and hybrid scheme for the convertive terms. Turbulence modeling is very important for prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. In this paper k-$\varepsilon$ turbulence model with wall function is used to increase efficiency of computation times.

• Korea-Australia Rheology Journal
• /
• v.16 no.2
• /
• pp.101-108
• /
• 2004

The hemodynamics behavior of the blood flow is influenced by the presence of the arterial stenosis. If the stenosis is present in an artery, normal blood flow is disturbed. In the present study, the characteristics of pulsatile flow in the blood vessel with stenosis are investigated by the finite volume method. For the validation of numerical model, the computation results are compared with the experimental ones of Ojha et al. in the case of 45％ stenosis with a trapezoidal profile. Comparisons between the measured and the computed velocity profiles are favorable to our solutions. Finally, the effects of stenosis severity and wall shear stress are discussed in the present computational analysis. It can be seen, where the non-dimensional peak velocity is displayed for all the stenosis models at a given severity of stenosis, that it is exponentially increased. Although the stenosis and the boundary conditions are all symmetric, the asymmetric flow can be detected in the more than 57％ stenosis. The instability by a three-dimensional symmetry-breaking leads to the asymmetric separation and the intense swirling motion downstream of the stenosis.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.3
• /
• pp.29-36
• /
• 2002

A numerical analysis based on two-dimensional and three-dimensional incompressible Wavier-Stokes equations has been carried out for double-circular-arc compressor cascades and the results are compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE algorithm adopt pressure weighted interpolation method for non-staggered grid and hybrid scheme for the convective terms. Turbulence modeling is very important for prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. Considering computation times, $\kappa$-$\varepsilon$ turbulence model with wall function is used.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.3
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• Journal of Korean Society for Atmospheric Environment
• /
• v.14 no.3
• /
• pp.219-228
• /
• 1998

In the present investigation, a numerical model developed for the prediction of the wind flow over complex terrain is validated by comparing with the field experiments. For the solution of the Reynolds - Averaged Clavier- stokes equations which are the governing equations of the microscale atmospheric flow, the model is constructed based on the finite-volume formulation and the SIMPLEC pressure-correction algorithm for the hydrodynamic computation. The boundary- fitted coordinate system is employed for the detailed depiction of topography. The boundary conditions and the modified turbulence constants suitable for an atmospheric boundary- layer are applied together with the k- s turbulence model. The full- scale experiments of Cooper's Ridge, Kettles Hill and Askervein Hill are chosen as the validation cases . Comparisons of the mean flow field between the field measurements and the predicted results show good agreement. In the simulation of the wind flow over Askervein Hill , the numerical model predicts the three dimensional flow separation in the downslope of the hill including the blockage effect due to neighboring hills . Such a flow behavior has not been simulated by the theoretical predictions. Therefore, the present model may offer the most accurate prediction of flow behavior in the leeside of the hill among the existing theoretical and numerical predictions.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2010.10a
• /
• pp.27-27
• /
• 2010

There have been several studies about empirical relation between seismic source parameters(e.g., focal mechanisms, depths, magnitudes, etc.) and T-wave observation. In order to delineate the relation, numerical and theoretical approaches to figure out T-wave excitation mechanism are required. In an attempt to investigate source radiation and wave scattering effects in the oceanic crust on T-wave envelopes, we perform three-dimensional numerical modeling to synthesize T-wave envelopes. We first calculate seismic P- and SV-wave energy on the seafloor using the Direct Simulation Monte Carlo based on the Radiative Transfer Theory, which enables us to take into account both realistic seismic source parameters and wave scattering in heterogeneous media, and then estimate excited T-wave energy by normal mode computation. The numerical simulation has been carried out considering the following different conditions: source types (strike and normal faults), source depths (shallow and deep), and wave propagation through homogeneous and heterogeneous Earth media. From the results of numerical modeling, we confirmed that T-wave envelopes vary according to spatial seismic energy distributions on the seafloor for the various input parameters. Furthermore, the synthesized T-wave envelopes show directional patterns due to anisotropic source radiation, and the slope change of T-wave envelopes caused by focal depth. Seismic wave scattering in the oceanic crust is likely to control the shape of envelopes.

• The KSFM Journal of Fluid Machinery
• /
• v.19 no.1
• /
• pp.11-17
• /
• 2016

In this paper, a numerical analysis was performed to investigate the effect of the pitch angle of a helical nozzle on the performance characteristics of a vortex tube. Three-dimensional numerical simulation has been performed with standard $k-{\varepsilon}$ turbulence model by using FLUENT 13.0. The effect of the pitch angle of helical nozzle was described in term of ${\beta}$. A CFD analysis was performed on ${\beta}=0^{\circ}$, $5^{\circ}$, $10^{\circ}$, $15^{\circ}$. In order to realize the influence of ${\beta}$ on performances of the vortex tube. Computation results were expressed by the ${\beta}-{\Delta}T_{h,c}$ graph and radial profiles of axial velocity and swirl velocity. The results showed that ${\beta}$ which improves energy separation capacity of vortex tube was $5^{\circ}$ at ${\alpha}=0.33$, 0.5 and $10^{\circ}$ at ${\alpha}=0.33$. Besides, It was confirmed that the results were closely related to axial velocity and swirl velocity.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.73-79
• /
• 2005

Micro sensor is very useful for flow measurements in a number of engineering applications. Especially, it is necessary for the development of MEMS. This paper presents the 3D numerical simulation of flows around a micro flow sensor, which is mounted on a flat plate. The effects of the sensor configuration (i.e. bottom gap) and the Reynolds number on the flow field are numerically investigated. The numerical results indicate that the bottom gap clearly affects the flow fields over the top surface of the sensor. The Reynolds numbers also show a significant influence on the flow nature, especially on the recirculation zone at downstream of the sensor. The present results illustrate a certain improvement on the flow field for the sensor installed at O.5mm above the wall with four pillars, comparing with that directly mounted on the wall.

• Journal of Advanced Marine Engineering and Technology
• /
• v.30 no.2
• /
• pp.275-283
• /
• 2006

Experiment and three dimensional numerical investigations of incompressible turbulent flow through square channels with one- and two-sided ribbed walls are performed to determine pressure drop and heat transfer. The CFX(version 5.7) software package is used for the computation. The ribbed walls have a $45^{\circ}$ inclined square rib. Uniform heat flux is maintained on whole inner heat transfer channel area. The numerical results coincide with experimental data that obtained for $7,600{\le}Re{\le}24.900$, the pitch-to-rib height ratio (p/e) of 8.0. and the rib height-to-channel hydraulic diameter ratio ($e/D_h$) of 0.0667. The results show that values of local heat transfer coefficient and friction factor in the channel with two-sided ribbed wall are higher than those in the channel with one-sided ribbed walls.

• Journal of computational fluids engineering
• /
• v.14 no.1
• /
• pp.70-77
• /
• 2009

In this study, time-dependent numerical analysis was carried out to investigate the plasma jet impingement on a flat plate, and a compressible form of two-dimensional inviscid gas dynamics equations were solved using the flux corrected transport algorithm. The mathematical modeling of Joule heating in the polycarbonate capillary bore and the mass ablation from the bore wall was incorporated in the numerical analysis and the series of computation was performed for three cases depending on the distance of the opposing plate from the capillary exit. The computational results reveal that the presence of the opposing plate does not affect the flow conditions inside the capillary when compared to the case of open-air plasma discharge. In the exterior region, the flow structure shows the typical supersonic underexpanded jet which consists of the strong Mach disk in front of the opposing plate and the barrel shock at the side of the jet. It is found that the shock evolution becomes more quasi-steady when the plate distance decreases. Also, the effects of the distance between the capillary bore exit and the opposing plate on the flow conditions along the opposing plate are investigated and the pressure variation on the plate shows more complicated interaction between the plasma discharge and the opposing plate when the location of plate becomes closer to the capillary exit.