• Title/Summary/Keyword: Finite Differential Method

Search Result 511, Processing Time 0.032 seconds

STABILITY OF POSITIVE PERIODIC NUMERICAL SOLUTION OF AN EPIDEMIC MODEL

  • Kim, Mi-Young
    • Korean Journal of Mathematics
    • /
    • v.13 no.2
    • /
    • pp.149-159
    • /
    • 2005
  • We study an age-dependent s-i-s epidemic model with spatial diffusion. The model equations are described by a nonlinear and nonlocal system of integro-differential equations. Finite difference methods along the characteristics in age-time domain combined with finite elements in the spatial variable are applied to approximate the solution of the model. Stability of the discrete periodic solution is investigated.

  • PDF

A PETROV-GALERKIN METHOD FOR A SINGULARLY PERTURBED ORDINARY DIFFERENTIAL EQUATION WITH NON-SMOOTH DATA

  • Zheng T.;Liu F.
    • Journal of applied mathematics & informatics
    • /
    • v.22 no.1_2
    • /
    • pp.317-329
    • /
    • 2006
  • In this paper, a singularly perturbed ordinary differential equation with non-smooth data is considered. The numerical method is generated by means of a Petrov-Galerkin finite element method with the piecewise-exponential test function and the piecewise-linear trial function. At the discontinuous point of the coefficient, a special technique is used. The method is shown to be first-order accurate and singular perturbation parameter uniform convergence. Finally, numerical results are presented, which are in agreement with theoretical results.

UNIFORMLY CONVERGENT NUMERICAL SCHEME FOR A SINGULARLY PERTURBED DIFFERENTIAL-DIFFERENCE EQUATIONS ARISING IN COMPUTATIONAL NEUROSCIENCE

  • DABA, IMIRU TAKELE;DURESSA, GEMECHIS FILE
    • Journal of applied mathematics & informatics
    • /
    • v.39 no.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.

Nonlocal integral elasticity analysis of beam bending by using finite element method

  • Taghizadeh, M.;Ovesy, H.R.;Ghannadpour, S.A.M.
    • Structural Engineering and Mechanics
    • /
    • v.54 no.4
    • /
    • pp.755-769
    • /
    • 2015
  • In this study, a 2-D finite element formulation in the frame of nonlocal integral elasticity is presented. Subsequently, the bending problem of a nanobeam under different types of loadings and boundary conditions is solved based on classical beam theory and also 3-D elasticity theory using nonlocal finite elements (NL-FEM). The obtained results are compared with the analytical and numerical results of nonlocal differential elasticity. It is concluded that the classical beam theory and the nonlocal differential elasticity can separately lead to significant errors for the problem under consideration as distinct from 3-D elasticity and nonlocal integral elasticity respectively.

Semi-analytical elastostatic analysis of two-dimensional domains with similar boundaries

  • Deeks, Andrew J.
    • Structural Engineering and Mechanics
    • /
    • v.14 no.1
    • /
    • pp.99-118
    • /
    • 2002
  • The scaled-boundary finite element method is a novel semi-analytical technique, combining the advantages of the finite element and the boundary element methods with unique properties of its own. The method works by weakening the governing differential equations in one coordinate direction through the introduction of shape functions, then solving the weakened equations analytically in the other (radial) coordinate direction. These coordinate directions are defined by the geometry of the domain and a scaling centre. This paper presents a general development of the scaled boundary finite-element method for two-dimensional problems where two boundaries of the solution domain are similar. Unlike three-dimensional and axisymmetric problems of the same type, the use of logarithmic solutions of the weakened differential equations is found to be necessary. The accuracy and efficiency of the procedure is demonstrated through two examples. The first of these examples uses the standard finite element method to provide a comparable solution, while the second combines both solution techniques in a single analysis. One significant application of the new technique is the generation of transition super-elements requiring few degrees of freedom that can connect two regions of vastly different levels of discretisation.

NUMERICAL DISCRETIZATION OF A POPULATION DIFFUSION EQUATION

  • Cho, Sung-Min;Kim, Dong-Ho;Kim, Mi-Young;Park, Eun-Jae
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.14 no.3
    • /
    • pp.189-200
    • /
    • 2010
  • A numerical method is proposed and analyzed to approximate a mathematical model of age-dependent population dynamics with spatial diffusion. The model takes a form of nonlinear and nonlocal system of integro-differential equations. A finite difference method along the characteristic age-time direction is considered and primal mixed finite elements are used in the spatial variable. A priori error estimates are derived for the relevant variables.

Comparison of elastic buckling loads for liquid storage tanks

  • Mirfakhraei, P.;Redekop, D.
    • Steel and Composite Structures
    • /
    • v.2 no.3
    • /
    • pp.161-170
    • /
    • 2002
  • The problem of the elastic buckling of a cylindrical liquid-storage tank subject to horizontal earthquake loading is considered. An equivalent static loading is used to represent the dynamic effect. A theoretical solution based on the nonlinear Fl$\ddot{u}$gge shell equations is developed, and numerical results are found using the new differential quadrature method. A second solution is obtained using the finite element package ADINA. A major motivation of the study was to show that the new method can serve to verify finite element solutions for cylindrical shell buckling problems. For this purpose the paper concludes with a comparison of buckling results for a number of cases covering a wide range in tank geometry.

Numerical analysis of the differential pressure venturi-cone flowmeter (차압식 벤튜리콘 유량계에 대한 유동해석)

  • 윤준용;맹주성;이정원
    • Korean Journal of Air-Conditioning and Refrigeration Engineering
    • /
    • v.10 no.6
    • /
    • pp.714-720
    • /
    • 1998
  • The differential pressure venturi-cone flowmeter is an advanced flowmeter which has many advantages such as wide range of measurement, high accuracy, excellent flow turn-down ratio, low headless, short installation pipe length requirement, and etc. Like other differential pressure flowmeters, the venturi-cone flowmeter uses the law of energy conservation, but its shape and position make it perform better than others. The cone acts as its own flow conditioner and mixer, fully conditioning and mixing the flow prior to measurement. For the analysis, we used Reynolds-averaged Wavier-Stokes equations and k-$\omega$ turbulence model. The equations were fully transformed into the computational domain, the pressure-velocity coupling was made through SIMPLER algorithm, and the equations were discretized using finite analytic solutions of the liberalized equations(Finite Analytic Method). To control the separation phenomenon on the cone surface, we proposed a new shape of cone, and analyzed the flowfield in the new flowmeter system, and found the improvement on the performance of the new cone flowmeter.

  • PDF

Finite Element Solution of Ordinary Differential Equation by the Discontinuous Galerkin Method (불연속 갤러킨 방법에 의한 상미분방정식의 유한요소해석)

  • 김지경
    • Computational Structural Engineering
    • /
    • v.6 no.4
    • /
    • pp.83-88
    • /
    • 1993
  • A time-discontinuous Galerkin method based upon using a finite element formulation in time has evolved. This method, working from the differential equation viewpoint, is different from those which have been generally used. They admit discontinuities with respect to the time variable at each time step. In particular, the elements can be chosen arbitrarily at each time step with no connection with the elements corresponding to the previous step. Interpolation functions and weighting functions are taken to be discontinuous across inter-element boundaries. These methods lead to a unconditional stable higher-order accurate ordinary differential equation solver.

  • PDF

THE h × p FINITE ELEMENT METHOD FOR OPTIMAL CONTROL PROBLEMS CONSTRAINED BY STOCHASTIC ELLIPTIC PDES

  • LEE, HYUNG-CHUN;LEE, GWOON
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.19 no.4
    • /
    • pp.387-407
    • /
    • 2015
  • This paper analyzes the $h{\times}p$ version of the finite element method for optimal control problems constrained by elliptic partial differential equations with random inputs. The main result is that the $h{\times}p$ error bound for the control problems subject to stochastic partial differential equations leads to an exponential rate of convergence with respect to p as for the corresponding direct problems. Numerical examples are used to confirm the theoretical results.