• Journal of Mechanical Science and Technology
• /
• 제15권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.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회(2)
• /
• 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.

• 한국전산유체공학회지
• /
• 제4권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.

• 한국전산유체공학회지
• /
• 제20권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.

• 한국전산유체공학회지
• /
• 제2권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.

• 대한수학회지
• /
• 제45권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.

• 한국전산유체공학회지
• /
• 제10권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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2006년도 제26회 춘계학술대회논문집
• /
• pp.253-256
• /
• 2006

액체로켓의 연소실내 분무연소 현상을 해석하기 위해 비연소 반응의 연료와 산화제 액적의 분무유동을 해석하였다. 분무연소 모델로 DSF 모델과 Euler-Lagrange 방법을 사용하였다. 액적과 가스상 그리고 증발한 기체상의 커플효과는 PSIC 모델을 사용하여 계산하였으며 수치 해석 방법으로 SIMPLER 알고리즘과 QUICK Scheme을 사용하였다. 해석결과 연소실내 속도와 온도 분포를 알 수 있었고, 연료와 산화제의 몰분율 또한 알 수 있었다.

• 한국전산유체공학회지
• /
• 제2권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.

• 한국전산유체공학회지
• /
• 제2권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. "