• Journal of applied mathematics & informatics
• /
• 제37권5_6호
• /
• pp.357-382
• /
• 2019

In this paper, we consider boundary value problem for a weakly coupled system of two singularly perturbed differential equations of convection diffusion type with discontinuous source term. In general, solution of this type of problems exhibits interior and boundary layers. A numerical method based on streamline diffusiom finite element and Shishkin meshes is presented. We derive an error estimate of order $O(N^{-2}\;{\ln}^2\;N$) in the maximum norm with respect to the perturbation parameters. Numerical experiments are also presented to support our theoritical results.

• Journal of applied mathematics & informatics
• /
• 제8권3호
• /
• pp.811.1-818
• /
• 2001

• 한국전산유체공학회지
• /
• 제7권1호
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• 한국전산유체공학회지
• /
• 제4권1호
• /
• pp.8-18
• /
• 1999

In devising a numerical approximation for the convective spatial transport of a fluid mechanical quantity, it is noted that the convective motion of a scalar quantity occurs in one-way, or from upstream to downstream. This consideration leads to a new scheme termed a pure upwind difference scheme (PUDS) in which an estimated value for a fluid mechanical quantity at a control surface is not influenced from downstream values. The formal accuracy of the proposed scheme is third order accurate. Two typical benchmark problems of a wall-driven fluid flow in a square cavity and a buoyancy-driven natural convection in a tall cavity are computed to evaluate performance of the proposed method. for comparison, the widely used simple upwind scheme, power-law scheme, and QUICK methods are also considered. Computation results are encouraging: the proposed PUDS sensitized to the convection direction produces the least numerical diffusion among tested convection schemes, and, notable improvements in representing recirculation of fluid stream and spatial change of a scalar. Although the formal accuracy of PUDS and QUICK are the same, the accuracy difference of approximately a single order is observed from the revealed results.

• 설비공학논문집
• /
• 제3권2호
• /
• pp.97-102
• /
• 1991

Finite element method has been developed recently for the solution of the convection-diffusion problems. Finite element method has several advantages over finite difference method, but its requirement of the larger memory size of the computer has prevented from wide application. In the present study, line-by-line technique has been implemented to finite element method to overcome this disadvantage. Two dimensional laminar natural convection in square cavity was chosen as an example in this study. The numerical result shows good agreement with bench mark solution and the size of the coefficient marix has been reduced drastically.

• Journal of applied mathematics & informatics
• /
• 제33권5_6호
• /
• pp.635-648
• /
• 2015

We consider a weakly coupled system of singularly perturbed convection-diffusion equations with discontinuous source term. The diffusion term of each equation is associated with a small positive parameter of different magnitude. Presence of discontinuity and different parameters creates boundary and weak interior layers that overlap and interact. A numerical method is constructed for this problem which involves an appropriate piecewise uniform Shishkin mesh. The numerical approximations are proved to converge to the continuous solutions uniformly with respect to the singular perturbation parameters. Numerical results are presented which illustrates the theoretical results.

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

• Journal of applied mathematics & informatics
• /
• 제26권5_6호
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• 대한기계학회논문집B
• /
• 제25권7호
• /
• pp.989-996
• /
• 2001

High-speed and low-speed flows are simulated numerically by flowfield-dependent mixed explicit-implicit (FDMEI) method. This algorithm depends on implicitness parameters of convection, diffusion, diffusion gradients, and source terms which are calculated from the changes of local Mach, Reynolds, Peclet, and Damkohler numbers between adjacent nodes. Convection phenomena or shock waves are resolved from Mach number-dependent implicitness parameters whereas diffusion or viscous actions are simulated by Reynolds number or Peclet number-dependent implicitness parameters. Fluctuation components of all variables are properly accommodated spatially and temporally in the FDMEI procedure. To illustrate, some benchmark example problems are presented for comparisons of the FDMEI results with other available data. These results appear to be encouraging and point toward the need for further investigations of the FDMEI theory.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 1994년도 Proceedings of the Korea Automatic Control Conference, 9th (KACC) ; Taejeon, Korea; 17-20 Oct. 1994
• /
• pp.747-752
• /
• 1994

Ussing's flux ratio theorem (1978) reflects a reciprocal relationship behavior between the unidirectional fluxes in asymmetric steady diffusion-convection in a membrane slab. This surprising result has led to many subsequent studies in a wide range of applications, in particular involving linear models of time dependent problems in biology and physiology. Ussing's theorem and its extensions are inherently linear in character. It is of considerable interest to ask to what extent these results apply, if at all, in situations involving, for example, nonlinear reaction. A physiologically interesting situation has been considered by Weisiger et at. (1989, 1991, 1992) and by McNabb et al. (1990, 1991) who studied the role of albumin in the transport of ligands across aqueous diffusion barriers in a liver membrane slab. The results are that there exist reciprocal relationships between unidirectional fluxes in the steady state, although albumin is chemically interacting in a nonlinear way of the diffusion processes. However, the results do not hold in general at early times. Since this type of study first started, it has been speculated about when and how the Ussing's flux ratio theorem fails in a general diffusion-convection-reaction system. In this paper we discuss the validity of Ussing-type theorems in time-dependent situations, and consider the limiting time behavior of a general nonlinear diffusion system with interaction.