• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.541-550
• /
• 1997

Computational and experimental investigations on the three-dimensional flowfield through an automotive cooling fan are carried out in this work. Steady, incompressible, three-dimensional, turbulent flow through a rotating axial-flow fan is analyzed with Reynolds averaged Navier-Stokes equations and standard k-.epsilon. turbulence model. The governing equations are discretized with finite-volume approximations in non-orthogonal curvilinear coordinates. Computational static pressures on the casing wall agree well with the experimental data which are measured in this work. And, they are sensitive to the change of tip clearance. The flowfield is not significantly affected by the thickness of the blade. The k-.omega. model gives the static pressure rise on the casing wall which is similar to that with the k-.epsilon. model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.7
• /
• pp.907-915
• /
• 2000

The preconditioned Krylov subspace methods were applied to the incompressible Navier-Stoke's equations for convergence acceleration. Three of the Krylov subspace methods combined with the five of the preconditioners were tested to solve the lid-driven cavity flow problem. The MILU preconditioned CG method showed very fast and stable convergency. The combination of GMRES/MILU-CG solver for momentum and pressure correction equations was found less dependency on the number of the grid points among them. A guide line for stopping inner iterations for each equation is offered.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.11 no.4
• /
• pp.55-68
• /
• 2007

First-order least-squares method of a distributed optimal control problem for the incompressible Navier-Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing velocity-flux variables and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $L^2$ norm of the residual equations of the system. The optimal error estimates for least-squares finite element approximations are obtained.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• Journal of Power System Engineering
• /
• v.2 no.3
• /
• pp.15-20
• /
• 1998

An integrated computing system has been developed for a Navier-Stokes design procedure of an axial-flow compressor blades. The process is done on the four separate steps, i.e., determination of the basic profiles, generation of computational grids, cascade flow simulation and analysis of the computed results in design sense. Applications are made to the blade design of the LP compressor. Computational results are analyzed with respect to the flow-field characteristics and are compared with the expected design requirements. The present system are coupled with the design procedure of the turbomachinery blades using the Navier-Stokes technique.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.397-403
• /
• 2010

The Immersed boundary method(IBM) is one of CFD techniques which can simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In these cases, pressure buildup near IB was found to occur when linear interpolation and stadard mass conservation is used and the interpolation scheme became complicated when higher order of interpolation is adopted. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. Bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Proceedings of the Korean Society for History of Mathematics Conference
• /
• 2003.11a
• /
• pp.4.2-4.2
• /
• 2003

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.1
• /
• pp.30-45
• /
• 2013

The present work focuses on the unsteady aerodynamics and aeroelastic properties of a small-medium sized wind-turbine blade operating under ideal conditions. A tapered/twisted blade representative of commercial blades used in an experiment setup at the National Renewable Energy Laboratory is considered. The aerodynamic loads are computed using Computational Fluid Dynamics (CFD) techniques. For this purpose, FLUENT$^{(R)}$, a commercial finite-volume code that solves the Navier-Stokes and the Reynolds-Averaged Navier-Stokes (RANS) equations, is used. Turbulence effects in the 2D simulations are modeled using the Wilcox k-w model for validation of the CFD approach. For the 3D aerodynamic simulations, in a first approximation, and considering that the intent is to present a methodology and workflow philosophy more than highly accurate turbulent simulations, the unsteady laminar Navier-Stokes equations were used to determine the unsteady loads acting on the blades. Five different blade pitch angles were considered and their aerodynamic performance compared. The structural dynamics of the flexible wind-turbine blade undergoing significant elastic displacements has been described by a nonlinear flap-lag-torsion slender-beam differential model. The aerodynamic quasi-steady forcing terms needed for the aeroelastic governing equations have been predicted through a strip-theory based on a simple 2D model, and the pertinent aerodynamic coefficients and the distribution over the blade span of the induced velocity derived using CFD. The resulting unsteady hub loads are achieved by a first space integration of the aeroelastic equations by applying the Galerkin's approach and by a time integration using a harmonic balance scheme. Comparison among two- and three- dimensional computations for the unsteady aerodynamic load, the flap, lag and torsional deflections, forces and moments are presented in the paper. Results, discussions and pertinent conclusions are outlined.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.