• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.156-157
• /
• 2003

A simple method is proposed to split the flux vector of the Euler equations by introducing two artificial wave speeds. The direction of wave propagation can be adjusted by these two wave speeds. This idea greatly simplifies the upwinding, and leads to a new family of upwind schemes. Numerical flux function for multi-dimensional Euler equations is formulated for any grid system, structured or unstructured. A remarkable simplicity of the scheme is that it successfully achieves one-sided approximation for all waves without recourse to any matrix operation. Moreover, its accuracy is comparable with the exact Riemann solver. For 1-D Euler equations, the scheme actually surpasses the exact solver in avoiding expansion shocks without any additional entropy fix. The scheme can exactly resolve stationary contact discontinuities, and it is also freed of the carbuncle problem in multi­dimensional computations.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.49-55
• /
• 2003

A comprehensive study has been made for the investigation of the convergence characteristics of the LU scheme for the Euler equations using von Neumann stability analysis. The stability results indicate that the convergence rate is governed by a specific parameter combination. Based on this insight it is shown that the LU scheme will not suffer convergence deterioration at any grid aspect ration if the local time step is defined using appropriate parameter combination. The numerical results demonstrate that this time step definition gives uniform convergence for grid aspect ratios from one to $1\times10^4$.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.9
• /
• pp.1-11
• /
• 2004

A comprehensive study has been made for the investigation of the convergence and stability characteristics of the LU scheme for solving the Euler equations on unstructured meshes. The von Neumann stability analysis technique was initially applied to a scalar model equation, and then the analysis was extended to the Euler equations. The results indicated that the convergence rate is governed by a specific combination of flow parameters. Based on this insight, it was shown that the LU scheme does not suffer any convergence deterioration at all grid aspect ratios, as long as the local time step is defined using an appropriate parameter combination.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.175-177
• /
• 2003

The convergence characteristics of the LV scheme for the Euler equations have been investigated by using the Von Neumann stability analysis. The results indicated that the convergence rate is governed by a specific combination of CFD parameters. Based on this insight, it is shown that the convergence characteristics of the LV scheme is not deteriorated at any grid aspect-ratio as long as the local time step is defined based on the parameter combination. The numerical results demonstrated that this time step definition provide a uniform convergence for grid aspect-ratios between one to$1{\times}10^{4}$.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.2
• /
• pp.115-122
• /
• 2008

The effects of characteristic condition number on the convergence of preconditioned Euler equations were investigated. The two-dimensional preconditioned Euler equations adopting Choi and Merkle's preconditioning and the temperature preconditioning are considered. Preconditioned Roe's FDS scheme was adopted for spatial discretization and preconditioned LU-SGS scheme was used for time integration. It is shown that the convergence characteristics of the Euler equations are strongly affected by the characteristic condition number, and there is an optimal characteristic condition number for a problem. The optimal characteristic condition numbers for the Choi and Merkle's preconditioning and temperature preconditioning are different.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• v.11B no.4
• /
• pp.121-128
• /
• 2001

Numerical schemes for the simulation of the cold gas flow in the SF6 puffer type circuit breaker is presented. The governing equation is axisymmetric compressible Euler Equation and FVM is used to analyze the behavior of flow. The upwind scheme is used to avoid numerical instability and MUSCL is used to obtain high order accuracy. For the efficient calculation, AF-ADI scheme is used. The simulation result shows good agreement with the experimental data.

• Journal of the Korean Statistical Society
• /
• v.33 no.4
• /
• pp.459-470
• /
• 2004

We study the approximation of the solution of linear stochastic evolution equations driven by infinite-dimensional fractional Brownian motion with Hurst parameter H > 1/2 through discretization of space and time. The rate of convergence of an approximation for Euler scheme is established.

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