• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.842-845
• /
• 2005

A method for the numerical simulation of two-dimensional free-surface flow resulting from the propagation of regular gravity waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid finite-differences and spectral-method scheme. The formulation includes a boundary-fitted transformation, and is suitable for extension to incorporate large-eddy simulation (LES) and large-wave simulation (LWS) terms for turbulence and breaking wave modeling, respectively. Results are presented for the simulation of the free-surface flow over two different bottom topographies, with constant slope values of 1:10 and 1:20, two different inflow wave lengths and two different inflow wave heights. An absorption outflow zone is utilized and the results indicate minimum wave reflection from the outflow boundary. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion, while wave heights remain more or less unchanged. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape.

• Journal of computational fluids engineering
• /
• v.4 no.3
• /
• pp.52-62
• /
• 1999

The development of a parallelized aerodynamic simulation process involving moving bodies is presented. The implementation of this process is demonstrated using a fully systemized Chimera methodology for steady and unsteady problems. This methodology consists of a Chimera hole-cutting, a new cut-paste algorithm for optimal mesh interface generation and a two-step search method for donor cell identification. It is fully automated and requires minimal user input. All procedures of the Chimera technique are parallelized on the Cray T3E using the MPI library. Two and three-dimensional examples are chosen to demonstrate the effectiveness and parallel performance of this procedure.

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

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.13-20
• /
• 1997

A fundamental study for the numerical scheme to simulate unsteady nonlinear waves by solving Euler equations is presented. First a conservation form and a non-conservation form of the Euler equations with a free surface fitted coordinate system are compared. Next, a time splitting fractional step method and an alternating direction implicit(ADI) method for the time integration are compared. For the comparative study, flow calculations around a bottom bump in a channel and a NACA 0012 hydrofoil in a flume are performed. The results show that the ADI method with a third order upwind differencing scheme is very efficient in reducing the computing time with keeping the accuracy, And, there is no distinct difference between two expression forms except that the non-conservative form shows faster wave propagating velocity than the conservation form. Some results are compared with experiments and show good agreement.

• Journal of the Korean Mathematical Society
• /
• v.45 no.2
• /
• pp.559-573
• /
• 2008

In this paper, two classes of functions, involving a parameter and the classical Euler gamma function, and two functions, involving the classical Euler gamma function, are verified to be logarithmically completely monotonic in $(-\frac{1}{2},\infty)$ or $(0,\infty)$; some inequalities involving the classical Euler gamma function are deduced and compared with those originating from certain problems of traffic flow, due to J. Wendel and A. Laforgia, and relating to the well known Stirling's formula.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.99-106
• /
• 2005

Performance of first, second and third order accurate methods for calculation of in viscid fluxes in fluid flow governing equations are investigated here. For the purpose, an upwind method based on Roe's scheme is used to solve 2-dimensional Euler equations. To increase the accuracy of the method two different schemes are applied. The first one is a second and third order upwind-based algorithm with the MUSCL extrapolation Van Leer (1979), based on primitive variables. The other one is an upwind-based algorithm with the Chakravarthy extrapolation to the fluxes of mass, momentum and energy. The results show that the thickness of shock layer in the third order accuracy is less than its value in second order. Moreover, applying limiter eliminates the oscillations near the shock while increases the thickness of shock layer especially in MUSCL method using Van Albada limiter.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.253-256
• /
• 2006

A numerical study for spray flow of fuel and oxidizer droplets in the combustion chamber has been conducted prior to the analysis of spray combustion of the liquid rocket engine. As the spray combustion model, DSF model and Euler-Lagrange scheme have been used. While the coupling effects of the droplets between gas phase and evaporated vapor have been calculated using PSIC model, SIMPLER algorithm and QUICK scheme have been used as numerical schemes. As the results, the calculations have shown velocity and temperature distribution in combustion chamber as well as mole fraction of fuel and oxidizer.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.26-34
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geon unstructured triangular meshes. The flux terms are discretized based on a cell-centered formulation with the Roe's flux-difference splitting. The solution is advanced in time us backward-Euler time-stepping scheme. At each time step, the linear system of equation approximately solved wi th the Gauss-Seidel relaxation scheme. The effect of turbulence is with a standard k-ε two-equation model which is solved separately from the mean flow equation the same backward-Euler time integration scheme. The triangular meshes are generated advancing-front/layer technique. Validations are made for flows over the NACA 0012 airfoil. Douglas 3-element airfoil. Good agreements are obtained between the numerical result experiment.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.51-58
• /
• 1997

A reliable computational solver has been developed for the analysis of three-dimensional inviscid compressible flows around a nacelle of a high bypass ratio turbofan engine, The numerical algorithm is based on the modified Godunov scheme to allow the second order accuracy for space variables, while keeping the monotone features. Two step time integration is used not only to remove time step limitation but also to provide the second order accuracy in a time variable. The multi-block approach is employed to calculate the complex flow field, using an algebraic, conformal, and elliptic method. The exact solution of Riemann problem is used to define boundary conditions. The accuracy of the developed solver is validated by comparing its results around the isolated nacelle in the cruise flight regime with the solution obtained using a commercial code "RAMPANT. "