• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.43-52
• /
• 2007

A comparative study on flux functions for the 2-dimensional Euler equations has been conducted. Explicit 4-stage Runge-Kutta method is used to integrate the equations. Flux functions used in the study are Steger-Warming's, van Leer's, Godunov's, Osher's(physical order and natural order), Roe's, HLLE, AUSM, AUSM+, AUSMPW+ and M-AUSMPW+. The performance of MUSCL limiters and MLP limiters in conjunction with flux functions are compared extensively for steady and unsteady problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.36-40
• /
• 2006

A comparative study on flux functions for the 2-dimensional Euler equations has been conducted. Explicit 4-stage Runge-Kutta method is used to integrate the equations. Flux functions used in the study are Steger-Warming's, van Leer's. Godunov's, Osher's(physical order and natural order), Roe's, HILE, AUSM, AUSM+ and AUSMPW+. The performance of MUSCL limiters and MLP limiters in conjunction with flux functions are compared extensively for steady and unsteady problems.

• Bulletin of the Korean Mathematical Society
• /
• v.47 no.1
• /
• pp.29-37
• /
• 2010

We establish the existence of weak solutions in an infinite subsonic channel in the self-similar plane to the two-dimensional Burgers system. We consider a boundary value problem in a fixed domain such that a part of the domain is degenerate, and the system becomes a second order elliptic equation in the channel. The problem is motivated by the study of the weak shock reflection problem and 2-D Riemann problems. The two-dimensional Burgers system is obtained through an asymptotic reduction of the 2-D full Euler equations to study weak shock reflection by a ramp.

• Convergence Security Journal
• /
• v.5 no.3
• /
• pp.93-97
• /
• 2005

MacCormack's method is an explicit, second order finite difference scheme that is widely used in the solution of hyperbolic partial differential equations. Apparently, however, it has shown entropy violations under small discontinuity. This non-physical shock grows fast and eventually all the meaningful information of the solution disappears. Some modifications of MacCormack's methods follow ideas of central schemes with an advantage of second order accuracy for space and conserve the high order accuracy for time step also. Numerical results are shown to perform well for the one-dimensional Burgers' equation and Euler equations gas dynamic.

• Journal of applied mathematics & informatics
• /
• v.26 no.1_2
• /
• pp.15-27
• /
• 2008

We present and discuss an algorithm for the numerical solution of initial value problems of the form $D_*^\alpha$y(t) = f(t, y(t)), y(0) = y0, where $D_*^\alpha$y is the derivative of y of order $\alpha$ in the sense of Caputo and 0<${\alpha}{\leq}1$. The algorithm is based on the fractional Euler's method which can be seen as a generalization of the classical Euler's method. Numerical examples are given and the results show that the present algorithm is very effective and convenient.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.1-11
• /
• 2016

Considerable efforts are required to develop a monotone, robust and stable high-order numerical scheme for solving the hyperbolic system. The discontinuous Galerkin(DG) method is a natural choice, but elimination of the spurious oscillations from the high-order solutions demands a new development of proper limiters for the DG method. There are several available limiters for controlling or removing unphysical oscillations from the high-order approximate solution; however, very few studies were directed to analyze the exact role of the limiters in the hyperbolic systems. In this study, the performance of the several well-known limiters is examined by comparing the high-order($p^1$, $p^2$, and $p^3$) approximate solutions with the exact solutions. It is shown that the accuracy of the limiter is in general problem-dependent, although the Hermite WENO limiter and maximum principle limiter perform better than the TVD and generalized moment limiters for most of the test cases. It is also shown that application of the troubled cell indicators may improve the accuracy of the limiters under some specific conditions.

• Journal of Integrative Natural Science
• /
• v.6 no.1
• /
• pp.53-56
• /
• 2013

We are interested in the rate of convergence of solutions of 2D Navier-Stokes equations in a smooth bounded domain as the viscosity tends to zero under Navier friction condition. If the initial velocity is smooth enough($u{\in}W^{2,p}$, p>2), it is known that the rate of convergence is linearly propotional to the viscosity. Here, we consider the rate of convergence for nonsmooth velocity fields when the gradient of the corresponding solution of the Euler equations belongs to certain Orlicz spaces. As a corollary, if the initial vorticity is bounded and small enough, we obtain a sublinear rate of convergence.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

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

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is studied for the compressible flow calculations. Fourier analysis on the local time stepping and block-Jacobi preconditioned residual operators is performed using the linearized 2-D Navier-Stokes equations. It fumed out that block-Jacobi preconditioner has better performance in eigenvalue clustering. They are implemented in the 2-D compressible Euler and Wavier-Stokes calculations with multigrid methods to verify that the block-Jacobi preconditioned multistage time stepping shows better performance in convergence acceleration.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.9 no.1
• /
• pp.79-90
• /
• 2005

In this paper we study the numerical solutions of the stochastic differential equations of the form $$du(x,\;t)=f(x,\;t,\;u)dt\;+\;g(x,\;t,\;u)dW(t)\;+\;\sum\limits_{|q|\leq2m}\;A_q(x,\;t)D^qu(x,\;t)dt$$ where $0\;{\leq}\;t\;{\leq}\;T,\;x\;{\in}\;R^{\nu}$, ($R^{nu}$ is the $\nu$-dimensional Euclidean space). Here $u\;{\in}\;R^n$, W(t) is an n-dimensional Brownian motion, $$f\;:\;R^{n+\nu+1}\;{\rightarrow}\;R^n,\;g\;:\;R^{n+\nu+1}\;{\rightarrow}\;R^{n{\times}n},$$, and $$A_q\;:\;R^{\nu}\;{\times}\;[0,\;T]\;{\rightarrow}\;R^{n{\times}n}$$ where ($A_q,\;|\;q\;|{\leq}\;2m$) is a family of square matrices whose elements are sufficiently smooth functions on $R^{\nu}\;{\times}\;[0,\;T]\;and\;D^q\;=\;D^{q_1}_1_{\ldots}_{\ldots}D^{q_{\nu}}_{\nu},\;D_i\;=\;{\frac{\partial}{\partial_{x_i}}}$.