• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1624-1629
• /
• 2008

The governing equations in generalized curvilinear coordinates for a 3D pulsatile flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. The computational technique implements the implicit approximate factorization method of the Beam and Warming method (1978), which is the extension of the Alternate Direction Implicit (ADI) method. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Rindt & Steenhoven, 1991).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1392-1396
• /
• 2007

The governing equations in generalized curvilinear coordinates for 3D laminar flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms. and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. Multigrid methods are also applied because solving the equations on the coarse grids requires much less computational effort per iteration than on the fine grid. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Humphrey et al., 1977).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.98-103
• /
• 2000

The cut-off is a unique feature associated with duct acoustics due to the presence of duct walls. Duct geometry effect on sound radiation is another issue of duct acoustics. The radiation of duct internal noise to ambient from duct open ends with various geometries is studied via numerical methods. The linearized Euler's equations in generalized curvilinear coordinates are solved by the DRP finite difference scheme. A number of accurate boundary conditions are used at boundaries for the computational domain to minimize the non-physical reflections. The far field sound pressure levels are computed by the Kirchhoff integration method. We investigate the cut off phenomana and duct geometry effects on sound radiation with numerical results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.1010-1015
• /
• 2002

Over the last few decades, noise has played a major role in the development of aircraft engines. The dominant noise is generated by the wake interactions of fan and downstream stator. Engine inlet and exhaust ducts are being fitted with liner materials that aid in damping fan related noise. In this paper, the radiation of duct internal noise from duct open ends with liners is studies via numerical methods. The linearized Euler's equations in generalized curvilinear coordinates are solved by the DRP scheme. The far field sound pressure levels are computed by the Kirchhoff integration method. Through comparison of sound directivity from bell-mouth duct with and without liners, it is shown that radiation from engine inlet is affected by liner effects or a soft wall boundary condition.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.186-191
• /
• 1999

An Euler solver is developed to predict accurate aerodynamic data such as lift coefficient, drag coefficient, and moment coefficient. The conservation law form of the compressible Euler equations are used in the generalized curvilinear coordinates system. The Euler solver uses a finite volume method and the second order Roe's flux difference splitting scheme with min-mod flux limiter to calculate the fluxes accurately. An implicit scheme which includes the boundary conditions is implemented to accelerate the convergence rate. The multi-block grid is integrated into the flow solver for complex geometry. The flowfields are analyzed around NACA 0012 airfoil in the cases of $M_{\infty}=0.75,\;\alpha=2.0\;and\;M_{\infty}=0.80,\;\alpha=1.25$. The numerical results are compared with other numerical results from the literature. The final goal of this research is to prepare a robust and an efficient Navier-Stokes solver eventually.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.2
• /
• pp.9-17
• /
• 2002

For turbulent spray combustion flows a coupled numerical procedure was developed, This method was discretized by using generalized curvilinear coordinates to handle complex geometries. The preconditioning and eigenvalue rescaling techniques were employed to provide efficient convergences over a wide range of subsonic Mach numbers. The accuracy was validated by simulating the laminar cavity flow. The film cooling effect of a liquid rocket engine (KSR-III) were investigated by a spray combustion analysis. The film cooling showed a negative effect on the combustion efficiency. In the combustion chamber wall, the film cooling effect was revealed to be promoted by the production of fuel rich zone.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1996.06a
• /
• pp.59-62
• /
• 1996

A numerical study of a two-dimensional acoustic field is carride ort for a spinning vortex pair located neat a wall to investigate the effect of the wall from the spinning quadrupole source in unsteady vortical flows. Based on the known incompressible flow field, the perturbed compressible acoustic terms derived from the Euler equations are calculated. Non-reflecting boundary conditions on the free field and the solid boundary conditions are developed for a generalized curvilinear coordinates system to investigate the effect of a curced wall. It is concluded that the sound generated by the quadrupole sources of unsteady vortical flows in the presence of a flat wall or a circular cylinder can be calculated by using the source terms of hydrodynamic flow fluctuations in both near and far acoustic fields simultaneously.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.414-426
• /
• 1997

Numerical computations were conducted to characterize the three-dimensional laminar flow through an injector orifice having an inclined angle of 30 .deg.. For this study, the incompressible Navier-Stokes equations in generalized curvilinear coordinates, using a pseudocompressibility approach for continuity equation, were solved. The computations were performed using the finite difference implicit, approximately factored scheme of Beam and Warming and multi-block grids of complete continuity at block interfaces. The multi-block computations were validated for the steady state using direct comparison of multi-block solutions with equivalent single-block ones, including 2-D 180.deg. TAD and 3-D 90.deg. pipe bend. The comparisons between the numerical solutions and the flow field measurements for a tube with sudden contraction were presented in this work for solution validation. Computational results showed the nature of complex flow fields within the inclined injector orifice, including strong pressure-driven secondary flows in the cross stream induced by the effect of streamline curvature. In addition, asymmetric secondary flows were induced in the Reynolds number range above assumed laminar flow regime considered. However, turbulence calculations and grid dependency studies are needed for more accurate computations.

• 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.