• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.771-779
• /
• 2004

Unsteady behavior of the early wake in the viscous flow field past an impulsively started semicircular cylinder is studied numerically. In this paper, we propose the hybrid diffusion scheme to simulate dynamic characteristics of wake such as a fishtail-like flapping and an alternate vortex-shedding more accurately. This diffusion scheme based on particle strength exchange is mixed with the stochastic nature of random walk method. Also, the viscous splitting algorithm which calculates convective and diffusion terms successively is applied in order to handle random walk method effectively. Consequently, the early behavior of wake due to the breakdown of symmetrical vortici balance is more practically simulated with the vortex particle method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1640-1647
• /
• 1993

The numerical diffusion arising from streamline-to-grid skewness produces a deteriorating effect on the numerical accuracy. The QUICKER scheme to reduce the numerical diffusion requires more computational effort than the HYBRID scheme. This paper deals with the relative computational efficiencies of adopting QUICKER scheme with a coarser grid system and of adopting HYBRID scheme with a denser grid system. Laminar driven cavity flow with Re=400, 1000 is used as a test problem. It is found that QUICKER scheme with a coarser grid system is more efficient than the HYBRID scheme with a denser grid system.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.2
• /
• pp.467-474
• /
• 1993

Viscous flow and heat transfer phenomena in an axisymmetric cylinder which models a diesel engine have been numerically studied. In order to search for a way to minimize numerical diffusion, the effectiveness and the appropriateness of two selected numerical schemes for convective terms in the governing equations have been tested. They are Linear Upwind Difference Scheme and Hybrid Scheme. Using a standard k-.epsilon. turbulence model, the calculation has been carried out basically up to 180.deg. of crank angle. As a result, it was shown from comparison with previous experimental data that Linear Upwind Difference Scheme is less influenced than Hybrid Scheme by the numerical diffusion and it was suggested that these effects of numerical diffusion can be more significant than those due to turbulence modeling.

• Journal of applied mathematics & informatics
• /
• v.36 no.1_2
• /
• pp.135-154
• /
• 2018

In this paper, an almost second order overlapping Schwarz method for singularly perturbed third order convection-diffusion type problem is constructed. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region we use the combination of classical finite difference scheme and central finite difference scheme on a uniform mesh while on the non-layer region we use the midpoint difference scheme on a uniform mesh. It is shown that the numerical approximations which converge in the maximum norm to the exact solution. We proved that, when appropriate subdomains are used, the method produces convergence of second order. Furthermore, it is shown that, two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantages of this method used with the proposed scheme are it reduce iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.

• Journal of applied mathematics & informatics
• /
• v.36 no.3_4
• /
• pp.181-203
• /
• 2018

In this paper, we have constructed an overlapping Schwarz method for singularly perturbed second order convection-diffusion equations. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region we use the central finite difference scheme on a uniform mesh while on the non-layer region we use the mid-point difference scheme on a uniform mesh. It is shown that the numerical approximations which converge in the maximum norm to the exact solution. When appropriate subdomains are used, the numerical approximations generated from the method are shown to be first order convergent. Furthermore it is shown that, two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantages of this method used with the proposed scheme is it reduces iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.

• Communications of the Korean Mathematical Society
• /
• v.20 no.1
• /
• pp.179-193
• /
• 2005

We present results of fully nonlinear time-dependent simulations of a thin liquid film flowing up an inclined plane. Equations of the type $h_t+f_y(h) = -{\in}^3{\nabla}{\cdot}(M(h){\nabla}{\triangle}h)$ arise in the context of thin liquid films driven by a thermal gradient with a counteracting gravitational force, where h = h(x, t) is the fluid film height. A hybrid scheme is constructed for the solution of two-dimensional higher-order nonlinear diffusion equations. Problems in the fluid dynamics of thin films are solved to demonstrate the accuracy and effectiveness of the hybrid scheme.

• Journal of applied mathematics & informatics
• /
• v.27 no.5_6
• /
• pp.1001-1015
• /
• 2009

In this paper, two hybrid difference schemes on the Shishkin mesh are constructed for solving a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a small parameter multiplying the highest derivative. We prove that the schemes are almost second order convergence in the supremum norm independent of the diffusion parameter. Error bounds for the numerical solution and its derivative are established. Numerical results are provided to illustrate the theoretical results.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.10a
• /
• pp.29-34
• /
• 1995

The good features of the analytic function expansion nodal (AFEN) method are utilized to develop a practical scheme jot the multigroup diffusion problems, in combination with the polynomial expansion nodal (PEN) method. The thermal group fluxes exhibiting strong gradients are solved by the AFEN method[1-6], while the fast group fluxes that are smoother than the thermal group fuzes are solved by the PEN method[7-9]. The scheme is applied to a MOX-fuel loaded core with good results.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1035-1054
• /
• 2010

We consider a mixed type singularly perturbed one dimensional elliptic problem with discontinuous source term. The domain under consideration is partitioned into two subdomains. A convection-diffusion and a reaction-diffusion type equations are posed on the first and second subdomains respectively. Two hybrid difference schemes on Shishkin mesh are constructed and we prove that the schemes are almost second order convergence in the maximum norm independent of the diffusion parameter. Error bounds for the numerical solution and its numerical derivative are established. Numerical results are presented which support the theoretical results.

• Kyungpook Mathematical Journal
• /
• v.62 no.4
• /
• pp.699-713
• /
• 2022

This paper proposes a hybrid successive over-relaxation iterative method for the numerical solution of a nonlinear diffusion in a one-dimensional porous medium. The considered mathematical model is discretized using a computational complexity reduction scheme called half-sweep finite differences. The local truncation error and the analysis of the stability of the scheme are discussed. The proposed iterative method, which uses explicit group technique and modified successive over-relaxation, is formulated systematically. This method improves the efficiency of obtaining the solution in terms of total iterations and program elapsed time. The accuracy of the proposed method, which is measured using the magnitude of absolute errors, is promising. Numerical convergence tests of the proposed method are also provided. Some numerical experiments are delivered using initial-boundary value problems to show the superiority of the proposed method against some existing numerical methods.