• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.27-34
• /
• 2015

In this study, preconditioned adjoint equations for the quasi one-dimensional Euler equations are developed, and their computational benefit at all speed is assessed numerically. The preconditioned adjoint equations are derived without any assumptions on the preconditioning matrix. The dissipation for Roe type numerical flux is also suggested to scale the dissipation term properly at low Mach numbers as well as at high Mach numbers. The new preconditioned method is validated against analytical solutions. The convergence characteristics over wide range of Mach numbers is evaluated. Finally, several inverse designs for the nozzle are conducted and the applicability of the method is demonstrated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2347-2355
• /
• 1996

Purpose of the present study is to develop a computational design program for shape optimization, combining the numerical optimization technique with the flow analysis code. The present methodology is then validated in three cases of aerodynamic shape optimization. In the numerical optimization, a feasible direction optimization algorithm and shape functions are considered. In the flow analysis, the Navier-Stokes equations are discretized by a cell-centered finite volume method, and Roe's flux difference splitting TVD scheme and ADI method are used. The developed design code is applied to a transonic channel flow over a bump, and an external flow over a NACA0012 airfoil to minimize the wave drag induced by shock waves. Also a separated subsonic flow over a NACA0024 airfoil is considered to determine a maximum allowable thickness of the airfoil without separation.

• Journal of the Korean Society of Combustion
• /
• v.4 no.1
• /
• pp.117-129
• /
• 1999

A method for the mixing augmentation of the transverse injection of fuel in a Scramjet combustor is suggested. Based on the fact that the main factor controlling the mixing characteristics in transverse injection is the effective back-pressure (the pressure around the injection hole) it is tried to make a flow expansion near the injection port in order to reduce effective back-pressure. A three dimensional Navier-Stokes code adopting the upwind method of Roe#s flux difference splitting scheme is used. The k-w SST turbulence model turbulence model is used to calculate the turbulent viscosity. It is shown that the reducing of the effective back-pressure make it possible to increase the penetration distance and to increase the mixing rate without excessive losses of stagnation pressure. Also, it is shown that the streamwise vorticity generated by baroclinic torque has great influences on the mixing process.

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

To investigate the flow inside the centrifugal impeller, computer program which can solve Three-dimensional compressible turbulent flow has been developed. The Navier-Stokes equations were chosen as the governing equation for viscous flow while Euler equations for inviscid case. Time marching method was incorporated with the Flux Difference Splitting method suggested by Roe to capture the steep gradients such as a shock. For high order of accuracy, MUSCL approach was adopted while differentiable limiter to ensure TVD property. For turbulence closure, Baldwin- Lomax model was applied due to its simplicity. To demonstrate the capabilities of present program, several validation problems have been solved and compared with experiments and other available data. From the above calculations generally good agreements were obtained. Finally, the developed code was applied to Eckardt's impeller and the performance prediction was carried out. Some important aspects on boundary condition for successful simulation were discussed and the remedy was also introduced.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.93-98
• /
• 2001

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.

• International Journal of Aeronautical and Space Sciences
• /
• v.8 no.2
• /
• pp.44-53
• /
• 2007

A three-dimensional viscous flow solver has been developed for the prediction of the aerodynamic performance of hovering helicopter rotor blades using unstructured hybrid meshes. The flow solver utilized a vertex-centered finite-volume scheme that is based on the Roe's flux-difference splitting with an implicit Jacobi/Gauss-Seidel time integration. The eddy viscosity are estimated by the Spalart- Allmaras one-equation turbulence model. Calculations were performed at three operating conditions with varying tip Mach number and collective pitch setting for the Caradonna-Tung rotor in hover. Additional computations are made for the UH-60A rotor in hover. Reasonable agreements were obtained between the present results and the experiment in both blade loading and overall rotor performance. It was demonstrated that the present vertex-centered flow solver is an efficient and accurate tool for the assessment of rotor performance in hover.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.48-54
• /
• 1998

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.25-28
• /
• 2008

An aerodynamic calculation in hover of KUH main rotor blade is performed using a three-dimensional unstructured hybrid mesh viscous flow solver. The flow solver utilizes a vertex-centered finite-volume scheme that is based on the Roe's flux-difference splitting with an implicit Jacobi/Gauss-Seidel time integration. The eddy viscosity are estimated by the Spalart-Allmaras one-equation turbulence model. A solution-adaptive mesh refinement technique is used for efficient capturing of the tip vortex. Calculations are performed at several operating conditions with varying collective pitch setting for KUH main rotor blade in hover. Good agreements are obtained between the present and other results using HOST and CAMRAD II in overall rotor performance. It is demonstrated that the present vertex-centered flow solver is an efficient and accurate tool for the assessment of rotor performance in hover.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.60-68
• /
• 2005

For study on the unsteady wall interference effect, flows around a forced oscillating airfoil in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The Spalart-Allmaras one-equation model is employed for the turbulence effect. The computed results of the oscillating airfoil having a thin wake showed that the lift curve slope is increased and the magnitude of hysteresis loop is reduced by the interference effects. Since the vortex around the airfoil is generated and convected downstream faster than the free-air condition, the phase of lift, drag and pitching moment coefficients was shifted. The pressure on the test section wall shows harmonic terms having the oscillating frequency contained in the wail effect.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.