• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.5
• /
• pp.733-743
• /
• 2000

We solve the integral representation of the Navier-Stokes equations in a lagrangian view by tracking the particles, which have vortex strengths. We simulate the unsteady viscous flow around an impulsively started cylinder using the vortex particle method. Particles are advanced via the Biot-Savart law for a lagrangian evolution of particles. The particle strength is modified based on the scheme of particle strength exchange. The solid boundary satisfies the no-slip boundary condition by the vorticity generation algorithm. We newly modify the diffusion scheme and the boundary condition for simulating an unsteady flow efficiently. To save the computation time, we propose the mixed scheme of particle strength exchange and core expansion. We also use a lot of panels to ignore the curvature of the cylinder, and not to solve the evaluation of the surface density. Results are compared to those from other theoretical and experimental works.

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

A program based on a 2-D CFD code has been developed to simulate a gas turbine engine. 2-D Navier-Stokes implicit code with $k-\omega$ turbulent model is used in compressor and turbine. Lumped method chemical equilibrium code with 10 species of molecular is applied to combustor with assuming perfect mixture and 100% combustion efficiency at constant pressure state. Fluid properties are shared on interfaces between engine components. Compressor supplies outlet temperature and pressure to combustor. At the same time, combustor also carries temperature and pressure to turbine. The back pressure of compressor outlet is transferred by inlet pressure of turbine. Unsteady phenomena in rotor-stator are covered by mixing-plane method. The running condition of engine can be determined only by given the inlet condition of compressor, the outlet condition of turbine, equivalence ratio and rotating speed.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2384-2389
• /
• 2007

The 2nd order response surface method (RSM) has been carried out to get optimum thermal design for enhanced heat transfer on square channel with bleed holes. The RSM was used as an optimization technique with Reynolds-averaged navier-stokes equation. Turbulence model for heat transfer analysis used RNG k-epsilon model. The wall function used enhanced wall function. Numerical local heat transfer coefficients were similar to the experimental tendency. Two design variables such as attack angle of rib (${\alpha}$), rib pitch-to-rib height ratio (p/e) were chosen. Operation condition considered bleeding ratio per bleed hole ($BR_{hole}$). A response surface were constructed by the design variables and operation condition. As a result, adjusted $R^2$ was more than 0.9. Optimization results of various objective function were similar to heat transfer in channel with and without bleed flow. But friction factor was lower than channel without bleed flow.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.157-164
• /
• 2015

The present work performs three-dimensional flow calculations based on Reynolds Averaged Navier-Stokes (RANS) and Delayed Detached Eddy Simulation (DDES) to investigate the flow field of a transonic rotor (NASA Rotor 37) at near-stall condition. It is found that the DES approach is likely to predict well the complex flow characteristics such as secondary vortex or turbulent flow phenomenon than RANS approach, which is useful to describe the flow mechanism of a transonic compressor. Especially, the DES results show improvement of predicting the flow field in the wake region and the model captures reasonably well separated regions compared to the RANS model. Besides, it is discovered that the three-dimensional vortical flows after the vortex breakdown from the rotor tip region are widely distributed and its vortex structures are clearly present. Near the rotor leading edge, a part of the tip leakage flows in DES solution spill over into next passage of the blade owing to the separation vortex flow and the backflow is clearly seen around the trailing edge of rotor tip. Furthermore, the DES solution shows strong turbulent eddies especially in the rotor hub, rotor tip section and the downstream of rotor trailing edge compared to the RANS solution.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.104-111
• /
• 2009

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, the boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Using the newly modified code, numerical simulations were performed and compared with other computational results as well as the experimental data for the supersonic flows over an oblique shock with a bleed region.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• Journal of Mechanical Science and Technology
• /
• v.18 no.6
• /
• pp.1026-1035
• /
• 2004

An immersed boundary method for solving the Navier-Stokes and thermal energy equations is developed to compute the heat transfer over or inside the complex geometries in the Cartesian or cylindrical coordinates by introducing the momentum forcing, mass source/sink, and heat source/sink. The present method is based on the finite volume approach on a staggered mesh together with a fractional step method. The method of applying the momentum forcing and mass source/sink to satisfy the no-slip condition on the body surface is explained in detail in Kim, Kim and Choi (2001, Journal of Computational Physics). In this paper, the heat source/sink is introduced on the body surface or inside the body to satisfy the iso-thermal or iso-heat-flux condition on the immersed boundary. The present method is applied to three different problems : forced convection around a circular cylinder, mixed convection around a pair of circular cylinders, and forced convection around a main cylinder with a secondary small cylinder. The results show good agreements with those obtained by previous experiments and numerical simulations, verifying the accuracy of the present method.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.1
• /
• pp.91-98
• /
• 2017

Performance of a low speed axial compressor is enhanced through a proper configuration of blade row tip injection and casing treatment of groove type. Air injectors were mounted evenly spaced upstream of the blade row within the casing groove and were all aligned parallel to the compressor axis. The groove, which covers all the blade tip chord length, extends all-round the casing circumference. Method of investigation is based on solution of the unsteady form of the Navier-Stokes equations utilizing $k-{\omega}$ SST turbulence model. Extensive parametric studies have been carried out to explore effects of injectors' flow momentums and yaw angles on compressor performance, while being run at different throttle valve setting. Emphasis has been focused on situations near to stall condition. Unsteady numerical analyses for untreated casing and no-injection case for near stall condition provided to discover two well-known criteria for spike stall inception, i.e., blade leading edge spillage and trailing edge back-flow. Final results showed that with only 6 injectors mounted axially in the casing groove and at yaw angle of 15 degrees opposite the direction of the blade row rotation, with a total mass flow rate of only 0.5% of the compressor main flow, surprisingly, the stall margin improves by 15.5%.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.5
• /
• pp.523-529
• /
• 2008

The course stability performance for tankers is evaluated by computational fluid dynamics. In the present work, a Reynolds averaged Navier-Stokes (RANS) code is applied to a maneuvering problem covering the pure drift and yaw motions. The purposes of this study are to evaluate the hydrodynamic force in the bare hull (AFRAMAX) in pure drift and yaw motion and to provide information about the trends in the forces and moments when the rudder angles are varied. The flow simulation is performed by FLUENT. The CFD code is examined to find the optimistic computational condition such as size of grid, turbulence model and initial condition. The hydrodynamic derivatives in drift and pure yaw motion are estimated by the numerical simulation, and then the stability levers are calculated. It is confirmed that the computations show the superiority and inferiority of course stability performance according to the hull forms. Finally, the CFD code is applied to the estimation of the rudder forces when the rudder angles are varied. The propeller effect expressed by the body force distribution is also included.