• 한국가시화정보학회지
• /
• 제20권3호
• /
• pp.74-85
• /
• 2022

We propose a Cubic-Interpolated Pseudo-Particle Lattice Boltzmann method (CIP-LBM) for the convection-diffusion equation (CDE) based on the Bhatnagar-Gross-Krook (BGK) scheme equation. The CIP-LBM relies on an accurate numerical lattice equilibrium particle distribution function on the advection term and the use of a splitting technique to solve the Lattice Boltzmann equation. Different schemes of lattice spaces such as D1Q3, D2Q5, and D2Q9 have been used for simulating a variety of problems described by the CDE. All simulations were carried out using the BGK model, although another LB scheme based on a collision term like two-relation time or multi-relaxation time can be easily applied. To show quantitative agreement, the results of the proposed model are compared with an analytical solution.

• Journal of applied mathematics & informatics
• /
• 제28권1_2호
• /
• pp.109-130
• /
• 2010

In this paper, a numerical method for a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with the mixed type boundary conditions is presented. Parameter-uniform error bounds for the numerical solution and also to numerical derivative are established. Numerical results are provided to illustrate the theoretical results.

• Journal of applied mathematics & informatics
• /
• 제39권5_6호
• /
• pp.655-676
• /
• 2021

A parameter uniform numerical scheme is proposed for solving singularly perturbed parabolic partial differential-difference convection-diffusion equations with a small delay and advance parameters in reaction terms and spatial variable. Taylor's series expansion is applied to approximate problems with the delay and advance terms. The resulting singularly perturbed parabolic convection-diffusion equation is solved by utilizing the implicit Euler method for the temporal discretization and finite difference method for the spatial discretization on a uniform mesh. The proposed numerical scheme is shown to be an ε-uniformly convergent accurate of the first order in time and second-order in space directions. The efficiency of the scheme is proved by some numerical experiments and by comparing the results with other results. It has been found that the proposed numerical scheme gives a more accurate approximate solution than some available numerical methods in the literature.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제7권2호
• /
• pp.75-81
• /
• 2003

Compact schemes are shown to be effective for a class of problems including convection-diffusion equations when combined with multigrid algorithms [7, 8] and V-cycle convergence is proved[5]. We apply the multigrid algorithm for an semianalytic finite difference scheme, which is desinged to preserve high order accuracy despite of singularities.

• 한국농공학회지
• /
• 제36권1호
• /
• pp.128-140
• /
• 1994

Coastal areas, especially embayments are apt to be polluted easily and many embayments in Korea are already suffering from pollution problems. To manage such pollution, it is strongly needed to develop technique to trace movements of pollution. Such technique cove- ring the embayment affected by the tidal influence, should take account both of the convection and the diffusion motions which cause lots of problems in numerical calculation. In this study, a Eulerian-Lagrangian Analysis(ELA) model was applied to Young Il bay and tested for its applicablity, which was developed by using the Eulerian-Lagrangian Method that reduce the numerical disperison and oscillation by way of solving convection and diffusion terrns separately. Concentration of Chemical Oxygen Demand(COD) and Suspend Solid(SS) of the embay- ment area were estimated by the model and compared with the observed values and the sound results were obtained. At the same time the diffsion coefficient and decay coefficient for Chemical Oxygen Demand in the Young II Bay were found as Dx = Dy = 20m$^2$/sec, kd=2.5 ${\times}$ 10-5/sec respectively, and for Suspend Solid, Dx =Dy = 30m$^2$/sec, kd=5.0${\times}$ 10-5/sec

• 한국수자원학회논문집
• /
• 제31권5호
• /
• pp.599-612
• /
• 1998

수송방정식의 양면적은 특성으로 인하여 이송항이 지배적인 흐름에 있어서 수송방정식의 수채해석은 매우 난해하다. 특히 유한요소법을 사용하여 수치해석할 때, 상류방향으로 더 많은 가중치를 두기 위하여 변화된 시행함수를 사용한다. 이러한 방법을 Petrov-Galerkin 방법이라고 한다. 본 논문에서는 N+1 과 N+2 Petrov-Galerkin 방법을 격자 Peclet 수가 큰 수송문제에 적용하였다. 주파수맞춤 기법을 사용하여 N+2 Petrov- Galerkin 방법을 격자 Peclet 수가 큰 소송문제에 적용하였다. 주파수맞춤 기법을 사용하여 N+2 Petrov-Galerkin 방법의 적정 풍상정도를 찾아내었고, 그 결과를 토의하였다. 이 기법의 시행함수는 이송항과 확산항의 상대적 크기에 따라 그 모양이 변화된다. 수치실험을 통하여 제시된 새로운 수치해석기법의 우수성을 설명하였다.

• Journal of applied mathematics & informatics
• /
• 제23권1_2호
• /
• pp.141-152
• /
• 2007

In this paper a numerical method is presented to solve singularly perturbed two points boundary value problems for second order ordinary differential equations consisting a discontinuous source term. First, in this method, an asymptotic expansion approximation of the solution of the boundary value problem is constructed using the basic ideas of a well known perturbation method WKB. Then some initial value problems and terminal value problems are constructed such that their solutions are the terms of this asymptotic expansion. These initial value problems are happened to be singularly perturbed problems and therefore fitted mesh method (Shishkin mesh) are used to solve these problems. Necessary error estimates are derived and examples provided to illustrate the method.

• Journal of applied mathematics & informatics
• /
• 제28권1_2호
• /
• pp.31-48
• /
• 2010

In this article, we consider singularly perturbed Boundary Value Problems(BVPs) for second order Ordinary Differential Equations (ODEs) with Neumann boundary condition and discontinuous source term. A parameter-uniform error bound for the solution is established using the Streamline-Diffusion Finite Element Method (SDFEM) on a piecewise uniform meshes. We prove that the method is almost second order of convergence in the maximum norm, independently of the perturbation parameter. Further we derive superconvergence results for scaled derivatives of solution of the same problem. Numerical results are provided to substantiate the theoretical results.

• Journal of applied mathematics & informatics
• /
• 제28권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.

• Journal of applied mathematics & informatics
• /
• 제22권1_2호
• /
• pp.49-65
• /
• 2006

A robust numerical method for a singularly perturbed second-order ordinary differential equation having two parameters with a discontinuous source term is presented in this article. Theoretical bounds are derived for the derivatives of the solution and its smooth and singular components. An appropriate piecewise uniform mesh is constructed, and classical upwind finite difference schemes are used on this mesh to obtain the discrete system of equations. Parameter-uniform error bounds for the numerical approximations are established. Numerical results are provided to illustrate the convergence of the numerical approximations.