• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.911-918
• /
• 1997

The turbulent viscous wake flows behind a single airfoil, two-dimensional stationary blade row and three-dimensional rotating blade row were calculated, and the numerical results were compared with experimental ones. The numerical technique was based on the SIMPLE algorithm using three turbulent closure models, standard k-.epsilon. model(WFM), low Reynolds number k-.epsilon. model(LRN) and Reynolds stress model (RSM). In the case of a single airfoil, WFM, LRN and RSM presented fairly good velocity distributions in the wake compared with experimental data. In the case of the stationary blade row, LRN and RSM presented better results than WFM for wake velocity distribution, and especially LRN showed best results among these three turbulent models. In the case of the rotating blade row, WFM and LRN showed fairly good agreement with experimental data of the three-dimensional velocity component distributions in the range from hub to mid span region. LRN was also superior to WFM in accuracy of prediction for the wake velocity distribution as same with the cases of a airfoil and the stationary blade row.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.88-92
• /
• 2009

When exposed to uniform magnetic fields externally applied, paramagnetic particles acquire dipole moments and the induced moments interacting with each other lead to the formation of chainlike structures or clusters of particles aligned with the field direction. A direct simulation method, based on the Maxwell stress tensor and a fictitious domain method, is applied to solve flows with magnetic chains in simple shear flow. We assumed that the particles constituting the chains are paramagnetic, and inertia of both flow and magnetic particles is negligible. The numerical scheme enables us to take into account both hydrodynamic and magnetic interactions between particles in a fully coupled manner, enabling us to numerically visualize breakup and reformation of the chains by the combined effect of the external field and the shear flow. Simple shear flow with suspended magnetic chains is solved in a periodic domain for a given magnetic field. Dynamics of interacting magnetic chains is found to be significantly affected by a dimensionless parameter called the Mason number, the ratio of the viscous force to the magnetic force in the shear flow. The effect of particle area fraction on the chain dynamics is investigated as well.

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.22-31
• /
• 2001

The two-dimensional unsteady incompressible viscous flow behind rectangular bluff bodies between two parallel plates was analyzed numerically. The steady state flow and the vortex flow behind rectangular bluff bodies in the channel were investigated for two regimes i.e., the laminar(Re = 100, 300, 500) and the turbulent flows(Re = 10⁴∼10/sup 6/). The vortex shedding was generated by a physical disturbance(6%) numerically imposed at the rear of the bluff bodies for a short time. It was observed that the perturbed flow became periodic after a transient period. And in the case of unsteady inflow, the sinusoidal pulsatile flow was applied as the inlet condition in the turbulent flow of Reynolds number of 1.0×10/sup 5/. FLUENT code was employed to solve the problems. The power-law scheme was used to get stable linearized equations and the PISO algorithm was applied to finding the solution of them.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.1
• /
• pp.84-94
• /
• 2003

It is known that tip clearance flows reduce the pressure rise, flow range and efficiency of the turbomachinery. So, the clear understanding about flow fields in the tip region is needed to efficiently design the turbomachinery. The Navier-Stokes code with the proper treatment of the boundary conditions has been developed to analyze the three-dimensional steady viscous flow fields in the transonic rotating blades and a numerical study has been conducted to investigate the detail flow physics in the tip region of transonic rotor, NASA Rotor 67. The computational results in the tip region of transonic rotors show the leakage vortices, leakage flow from pressure side to suction side and their interaction with a shock. Depen ding on the operating conditions, toad distributions and the position of shock-wave on the blade surface are very different close to the blade tip of the transonic compressor rotor. The load distribution and the shock-wave position close to the blade tip had the close relationship with the starting position of leakage vortex and the direction of leakage flow.

• 유체기계공업학회:학술대회논문집
• /
• 1998.12a
• /
• pp.42-48
• /
• 1998

Numerical study is presented for the analysis of three-dimensional incompressible turbulent flows in multiblade centrifugal fan. Reynolds-averaged Navier-Stokes equations with standard k - $\epsilon$ turbulence model are transformed to non-orthogonal curvilinear coordinates, and are discretized with finite volume approximations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. The computational area is divided into three blocks; core, impeller and scroll, which are linked by multi-block method. The flow inside of the fan is regarded as steady flow, and mathematical formula established from the cascade theory and empirical coefficient are employed to simulate tile flow through the impeller. From comparisons between the computational results and the experimental data, the validity of the mathematical formula for the blade forces was examined and good results were obtained qualitatively. Hence, we can get the flow characteristics of multi-blade centrifugal fan and it will be a corner stone of the development of the multiblade centrifugal fan.

• Korea-Australia Rheology Journal
• /
• v.16 no.3
• /
• pp.117-128
• /
• 2004

The effect of molecular parameters on the steady vortex behaviors in the inertial viscoelastic flow past a cylinder has been investigated. FENE-CR model was considered as a constitutive equation. A recently developed iterative solution method (Kim et al., (in press)) was found to be successfully applicable to the computation of inertial viscoelastic flows. The high-resolution computations were carried out to understand the detailed flow behaviors based on the efficient iterative solution method armed with ILU(0) type pre-conditioner and BiCGSTAB method. The discrete elastic viscous split stress-G/streamline upwind Petrov Galerkin (DEVSS-G/SUPG) formulation was adopted as a stabilization method. The vortex size decreased as elasticity increases. However, the vortex enhancement was also observed in the case of large extensibility, which means that the vortex behavior is strongly dependent upon the material parameters. The longitudinal gradient of normal stress was found to retard the formation of vortex, whereas the extensional viscosity played a role in the vortex enhancement. The present results are expected to be helpful for understanding the inertial vortex dynamics of viscoelastic fluids in the flow past a confined cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.11
• /
• pp.1371-1378
• /
• 1999

Flow through turbomachinery has a very complex structure and Is Intrinsically unsteady. In addition, trend to highly loaded turbomachinery makes the flow extremely complex due to the interaction between rotor and stator. In this study, flows through UTRC LSRR turbine are numerically analyzed using 2 dimensional Navier-Stokes equations. The convective terms of the governing equations are discretized using the Van-Leer's FVS(Flux vector splitting) with an upwind TVD scheme. The conventional central differencing is used to discretize the diffusion terms on the finite volume. The accurate unsteady motion is achieved by using a 2nd order accurate, 3-point Euler implicit scheme. The quasi-conservative zonal scheme is used for calculating the flow variables on the zonal interface between the rotor and stator. The axial gap between stator and rotor has been configured in two variations, 15% and 65% of average chord length. The analysis program is validated using experimental results and the effect of axial gap is examined. The numerical analysis results are presented by time averaged pressure coefficient and pressure magnitude coefficient and compared with experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.1
• /
• pp.9-16
• /
• 2014

This paper presents the theoretical analysis of a flow driven by surface tension and gravity in an inclined circular tube. A governing equation is developed for describing the displacement of a non-Newtonian fluid(Power-law model) that continuously flows into a circular tube owing to surface tension, which represents a second-order, nonlinear, non-homogeneous, and ordinary differential form. It was found that quantitatively, the theoretical predictions of the governing equation were in excellent agreement with the solutions of the equation for horizontal tubes and the past experimental data. In addition, the predictions compared very well with the results of the force balance equation for steady.

• Journal of Ship and Ocean Technology
• /
• v.10 no.2
• /
• pp.1-10
• /
• 2006

Numerical analysis of viscous flow around twin-skeg hull forms was conducted according to the variations of distance between skegs and vertical skeg inclinations by using a hydrodynamic analysis system, WAVIS. Six twin-skeg hull forms were derived by combining three distances between skegs (16m, 20m, 24m) and four vertical skeg angles ($0^{\circ},\;10^{\circ},\;15^{\circ},\;20^{\circ}$). It is found that the better resistance performance can be obtained with larger vertical skeg angle and smaller skeg distance for the present test cases. It also can be seen that the same trend is true for the nominal wake distributions in the propeller plane. Those tendencies were confirmed by the experimental results of MOERI. It is shown that numerical analysis can be a useful and practical tool for the evaluation and improvement of hydrodynamic performances for the complex stern hull forms with twin skegs.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.328-333
• /
• 2014

A hybrid particle-mesh method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations is a combination of the Vortex-In-Cell(VIC) method for convection and the penalization method for diffusion. The key feature of the numerical methods is to determine velocity and vorticity fields around a solid body on a temporary grid, and then the time evolution of the flow is computed by tracing the convection of each vortex element using the Lagrangian approach. Assuming that the vorticity and velocity fields are to be computed in time domain analysis, pressure fields are estimated through a complete set of solutions at present time step. It is possible to obtain vorticity and velocity fields prior to any pressure calculation since the pressure term is eliminated in the vorticity-velocity formulation. Therefore, pressure field is explicitly treated by solving a suitable Poisson equation. In this paper, we propose a simple way to numerically implement the vorticity-velocity-pressure formulation including a penalty term. For validation of the proposed numerical scheme, we illustrate the early development of viscous flows around an impulsive started circular cylinder for Reynolds number of 9500.