In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(...In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order tim...Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.展开更多
A group of asymmetric difference schemes to approach the Korteweg-de Vries (KdV) equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segme...A group of asymmetric difference schemes to approach the Korteweg-de Vries (KdV) equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.展开更多
In this paper, two different numerical schemes, namely the Runge-Kutta fourth order method and the implicit Euler method with perturbation method of the second degree, are applied to solve the nonlinear thermal wave i...In this paper, two different numerical schemes, namely the Runge-Kutta fourth order method and the implicit Euler method with perturbation method of the second degree, are applied to solve the nonlinear thermal wave in one and two dimensions using the differential quadrature method. The aim of this paper is to make comparison between previous numerical schemes and detect which is more efficient and more accurate by comparing the obtained results with the available analytical ones and computing the computational time.展开更多
In this paper, a new kind of alternating direction implicit (ADI) Crank-Nicolson-type orthogonal spline collocation (OSC) method is formulated for the two-dimensional frac-tional evolution equation with a weakly s...In this paper, a new kind of alternating direction implicit (ADI) Crank-Nicolson-type orthogonal spline collocation (OSC) method is formulated for the two-dimensional frac-tional evolution equation with a weakly singular kernel arising in the theory of linear viscoelas-ticity. The novel OSC method is used for the spatial discretization, and ADI Crank-Nicolson-type method combined with the second order fractional quadrature rule are considered for thetemporal component. The stability of proposed scheme is rigourously established, and nearlyoptimal order error estimate is also derived. Numerical experiments are conducted to supportthe predicted convergence rates and also exhibit expected super-convergence phenomena.展开更多
In this research article, two finite difference implicit numerical schemes are described to approximate the numerical solution of the two-dimension modified reaction diffusion Fisher’s system which exists in coupled ...In this research article, two finite difference implicit numerical schemes are described to approximate the numerical solution of the two-dimension modified reaction diffusion Fisher’s system which exists in coupled form. Finite difference implicit schemes show unconditionally stable and second-order accurate nature of computational algorithm also the validation and comparison of analytical solution, are done through the examples having known analytical solution. It is found that the numerical schemes are in excellent agreement with the analytical solution. We found, second-implicit scheme is much faster than the first with good rate of convergence also we used NVIDA devices to accelerate the computations and efficiency of the algorithm. Numerical results show our proposed schemes with use of HPC (High performance computing) are very efficient and reliable.展开更多
In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two ...In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.展开更多
Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly im...Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.展开更多
We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equati...We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation (image selective smoothing model) given by Alvarez et al. (Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer. Anal., 1992, 29: 845-866). Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes, unconditionally stable in L∞ and W1'2 (W1'1) sense in isotropic (anisotropic) diffu- sion domain.展开更多
In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and di...In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and discuss the iteratively B-convergence of the Newton iterative process for solving the algebraic equations of the scheme, secondly we present a strategy providing initial values parallelly for the iterative process. Finally, some numerical results show that our parallel scheme is higher efficient as N is not so large.展开更多
Referring to the construction way of Lax-Wendroff scheme, new IRS (Implicit Residual Smoothing) schemes have been developed for hyperbolic, parabolic and hyperparabolic equations. These IRS schemes have 2nd-order or 3...Referring to the construction way of Lax-Wendroff scheme, new IRS (Implicit Residual Smoothing) schemes have been developed for hyperbolic, parabolic and hyperparabolic equations. These IRS schemes have 2nd-order or 3 rd-order time accuracy, and can extend the stability region of basic explicit time-stepping scheme greatly and thus can permit higher CFL number in the calculation of flow field. The, central smoothing and upwind-bias smoothing techniques have been developed tao. Based on one-dimensional linear model equation, it has been found that the scheme is unconditionally stable according to the von-Neumann analysis. The limitation of Dawes' method, which has been applied in turbomachinery widespreadly, has been discussed on solving steady flow and viscous flow. It is shown that stable solution of this method is not completely independent with the value of time step. In the end, numerical results by using IRS schemes and Dawes' method as well as TVD ( total variation diminishing) scheme and four-stage Runge-Kutta technique are presented to verify, the analytical conclusions.展开更多
An EBCOT-based method is proposed to reduce the priority of background coefficients in the ROI code block without compromising algorithm complexity. The region of interest is encoded to a higher quality level than bac...An EBCOT-based method is proposed to reduce the priority of background coefficients in the ROI code block without compromising algorithm complexity. The region of interest is encoded to a higher quality level than background, and the target searching time in video-guided penetrating missile can be shortened. Three kinds of coding schemes based on EBCOT are discussed. Experimental results demonstrate that the proposed method shows higher compression efficiency, lower complexity, and good reconstructed ROI image quality in the lower channel capacity.展开更多
In this paper, a modified implicit Kirk-multistep iteration scheme and a strong convergence result for a general class of maps in a normed linear space was established. It was also shown that the convergence of this i...In this paper, a modified implicit Kirk-multistep iteration scheme and a strong convergence result for a general class of maps in a normed linear space was established. It was also shown that the convergence of this iteration scheme is equivalent to the convergency of some other implicit Kirk-type iteration (implicit Kirk-Noor, implicit Kirk-Ishikawa and implicit Kirk-Mann iterations) for the same class of maps. Some numerical examples were considered to show that the equivalence of convergence results to the fixed point is true. The results unify most equivalence results in literature.展开更多
In order to investigate parameters of FAE (fuel air explosive) explosion, the two-phase gas-droplet conservation equations with two-dimensional axial symmetry in the Euler coordinate were used. High-resolution implici...In order to investigate parameters of FAE (fuel air explosive) explosion, the two-phase gas-droplet conservation equations with two-dimensional axial symmetry in the Euler coordinate were used. High-resolution implicit TVD ( total variation diminishing) schemes were applied to gas phase equations and MacCormack schemes to liquid equations. The formation and propagation of gas-droplet detonation wave were simulated numerically. The simulation results and the others are compared with a good agreement.展开更多
An implicit numerical scheme is developed based on the simplified marker and cell (SMAC) method to solve Reynolds-averaged equations in general curvilinear coordinates for three-dimensional (3-D) unsteady incompre...An implicit numerical scheme is developed based on the simplified marker and cell (SMAC) method to solve Reynolds-averaged equations in general curvilinear coordinates for three-dimensional (3-D) unsteady incompressible turbulent flow. The governing equations include the Reynolds-averaged momentum equations, in which contravariant velocities are unknown variables, pressure-correction Poisson equation and k- s turbulent equations. The governing equations are discretized in a 3-D MAC staggered grid system. To improve the numerical stability of the implicit SMAC scheme, the higherorder high-resolution Chakravarthy-Osher total variation diminishing (TVD) scheme is used to discretize the convective terms in momentum equations and k- e equations. The discretized algebraic momentum equations and k- s equations are solved by the time-diversion multiple access (CTDMA) method. The algebraic Poisson equations are solved by the Tschebyscheff SLOR (successive linear over relaxation) method with alternating computational directions. At the end of the paper, the unsteady flow at high Reynolds numbers through a simplified cascade made up of NACA65-410 blade are simulated with the program written according to the implicit numerical scheme. The reliability and accuracy of the implicit numerical scheme are verified through the satisfactory agreement between the numerical results of the surface pressure coefficient and experimental data. The numerical results indicate that Reynolds number and angle of attack are two primary factors affecting the characteristics of unsteady flow.展开更多
The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part o...The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part of the operator and an implicit Euler discretization of the linear part. Finite difference schemes are used for the spatial part. This finally leads to the numerical solution of a sparse linear system that can be solved efficiently.展开更多
A mimetic finite difference scheme for the transient heat equation under Robin’s conditions is presented. The scheme uses second order gradient and divergence mimetic operators, on a staggered grid, to approximate th...A mimetic finite difference scheme for the transient heat equation under Robin’s conditions is presented. The scheme uses second order gradient and divergence mimetic operators, on a staggered grid, to approximate the space derivatives. The temporal derivative is replaced by a first order backward difference approximation to obtain an implicit formulation. The resulting scheme contains nonstandard finite difference stencils. An original convergence analysis by the matrix’s method shows that the proposed scheme is unconditionally stable. A comparative study against standard finite difference schemes, based on central difference or first order one side approximations, reveals the advantages of our scheme without being its implementation more expensive or difficult to achieve.展开更多
This work aimed to construct an epidemic model with fuzzy parameters.Since the classical epidemic model doesnot elaborate on the successful interaction of susceptible and infective people,the constructed fuzzy epidemi...This work aimed to construct an epidemic model with fuzzy parameters.Since the classical epidemic model doesnot elaborate on the successful interaction of susceptible and infective people,the constructed fuzzy epidemicmodel discusses the more detailed versions of the interactions between infective and susceptible people.Thenext-generation matrix approach is employed to find the reproduction number of a deterministic model.Thesensitivity analysis and local stability analysis of the systemare also provided.For solving the fuzzy epidemic model,a numerical scheme is constructed which consists of three time levels.The numerical scheme has an advantage overthe existing forward Euler scheme for determining the conditions of getting the positive solution.The establishedscheme also has an advantage over existing non-standard finite difference methods in terms of order of accuracy.The stability of the scheme for the considered fuzzy model is also provided.From the plotted results,it can beobserved that susceptible people decay by rising interaction parameters.展开更多
A mathematical model comprising of nonlinear reaction, diffusion, and convection mechanisms seen in natural and anthropogenic processes is numerically investigated here. It is proposed that a higher order numerical sc...A mathematical model comprising of nonlinear reaction, diffusion, and convection mechanisms seen in natural and anthropogenic processes is numerically investigated here. It is proposed that a higher order numerical scheme of finite difference method be used in conjunction with an iterative approach in order to solve the nonlinear one dimensional convection-diffusion-reaction equation. To account for the wide variety of physical characteristics and boundary conditions, an iterative approach is presented that yields a reliable and precise solution every time. We examined the accuracy and operational efficiency of two distinct finite difference approaches. The efficiency of the system is determined by comparing the estimated results to the appropriate analytical solution by adhering to established norms. Coherence and convergence were analyzed for each approach. The simulation results demonstrate the efficacy and accuracy of these methods in solving nonlinear convection- diffusion-reaction equations. Convection-diffusion-reaction equation modeling is critical for employing the offered results in heat and mass transport processes.展开更多
文摘In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金Supported by the Discipline Construction and Teaching Research Fund of LUTcte(20140089)
文摘Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.
基金Project supported by the National Natural Science Foundation of China(No.10671113)the Natural Science Foundation of Shandong Province of China(No.Y2003A04)
文摘A group of asymmetric difference schemes to approach the Korteweg-de Vries (KdV) equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.
文摘In this paper, two different numerical schemes, namely the Runge-Kutta fourth order method and the implicit Euler method with perturbation method of the second degree, are applied to solve the nonlinear thermal wave in one and two dimensions using the differential quadrature method. The aim of this paper is to make comparison between previous numerical schemes and detect which is more efficient and more accurate by comparing the obtained results with the available analytical ones and computing the computational time.
基金supported by National Nature Science Foundation of China(11701168,11601144 and 11626096)Hunan Provincial Natural Science Foundation of China(2018JJ3108,2018JJ3109 and 2018JJ4062)+1 种基金Scientific Research Fund of Hunan Provincial Education Department(16K026 and YB2016B033)China Postdoctoral Science Foundation(2018M631403)
文摘In this paper, a new kind of alternating direction implicit (ADI) Crank-Nicolson-type orthogonal spline collocation (OSC) method is formulated for the two-dimensional frac-tional evolution equation with a weakly singular kernel arising in the theory of linear viscoelas-ticity. The novel OSC method is used for the spatial discretization, and ADI Crank-Nicolson-type method combined with the second order fractional quadrature rule are considered for thetemporal component. The stability of proposed scheme is rigourously established, and nearlyoptimal order error estimate is also derived. Numerical experiments are conducted to supportthe predicted convergence rates and also exhibit expected super-convergence phenomena.
文摘In this research article, two finite difference implicit numerical schemes are described to approximate the numerical solution of the two-dimension modified reaction diffusion Fisher’s system which exists in coupled form. Finite difference implicit schemes show unconditionally stable and second-order accurate nature of computational algorithm also the validation and comparison of analytical solution, are done through the examples having known analytical solution. It is found that the numerical schemes are in excellent agreement with the analytical solution. We found, second-implicit scheme is much faster than the first with good rate of convergence also we used NVIDA devices to accelerate the computations and efficiency of the algorithm. Numerical results show our proposed schemes with use of HPC (High performance computing) are very efficient and reliable.
文摘In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.
文摘Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.
文摘We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation (image selective smoothing model) given by Alvarez et al. (Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer. Anal., 1992, 29: 845-866). Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes, unconditionally stable in L∞ and W1'2 (W1'1) sense in isotropic (anisotropic) diffu- sion domain.
基金national natural science foundation natural science foundation of Gansu province.
文摘In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and discuss the iteratively B-convergence of the Newton iterative process for solving the algebraic equations of the scheme, secondly we present a strategy providing initial values parallelly for the iterative process. Finally, some numerical results show that our parallel scheme is higher efficient as N is not so large.
文摘Referring to the construction way of Lax-Wendroff scheme, new IRS (Implicit Residual Smoothing) schemes have been developed for hyperbolic, parabolic and hyperparabolic equations. These IRS schemes have 2nd-order or 3 rd-order time accuracy, and can extend the stability region of basic explicit time-stepping scheme greatly and thus can permit higher CFL number in the calculation of flow field. The, central smoothing and upwind-bias smoothing techniques have been developed tao. Based on one-dimensional linear model equation, it has been found that the scheme is unconditionally stable according to the von-Neumann analysis. The limitation of Dawes' method, which has been applied in turbomachinery widespreadly, has been discussed on solving steady flow and viscous flow. It is shown that stable solution of this method is not completely independent with the value of time step. In the end, numerical results by using IRS schemes and Dawes' method as well as TVD ( total variation diminishing) scheme and four-stage Runge-Kutta technique are presented to verify, the analytical conclusions.
基金the National Defense Foundation Research (A2320060093)the National Natural Science Foundation of China (60572039)
文摘An EBCOT-based method is proposed to reduce the priority of background coefficients in the ROI code block without compromising algorithm complexity. The region of interest is encoded to a higher quality level than background, and the target searching time in video-guided penetrating missile can be shortened. Three kinds of coding schemes based on EBCOT are discussed. Experimental results demonstrate that the proposed method shows higher compression efficiency, lower complexity, and good reconstructed ROI image quality in the lower channel capacity.
文摘In this paper, a modified implicit Kirk-multistep iteration scheme and a strong convergence result for a general class of maps in a normed linear space was established. It was also shown that the convergence of this iteration scheme is equivalent to the convergency of some other implicit Kirk-type iteration (implicit Kirk-Noor, implicit Kirk-Ishikawa and implicit Kirk-Mann iterations) for the same class of maps. Some numerical examples were considered to show that the equivalence of convergence results to the fixed point is true. The results unify most equivalence results in literature.
文摘In order to investigate parameters of FAE (fuel air explosive) explosion, the two-phase gas-droplet conservation equations with two-dimensional axial symmetry in the Euler coordinate were used. High-resolution implicit TVD ( total variation diminishing) schemes were applied to gas phase equations and MacCormack schemes to liquid equations. The formation and propagation of gas-droplet detonation wave were simulated numerically. The simulation results and the others are compared with a good agreement.
文摘An implicit numerical scheme is developed based on the simplified marker and cell (SMAC) method to solve Reynolds-averaged equations in general curvilinear coordinates for three-dimensional (3-D) unsteady incompressible turbulent flow. The governing equations include the Reynolds-averaged momentum equations, in which contravariant velocities are unknown variables, pressure-correction Poisson equation and k- s turbulent equations. The governing equations are discretized in a 3-D MAC staggered grid system. To improve the numerical stability of the implicit SMAC scheme, the higherorder high-resolution Chakravarthy-Osher total variation diminishing (TVD) scheme is used to discretize the convective terms in momentum equations and k- e equations. The discretized algebraic momentum equations and k- s equations are solved by the time-diversion multiple access (CTDMA) method. The algebraic Poisson equations are solved by the Tschebyscheff SLOR (successive linear over relaxation) method with alternating computational directions. At the end of the paper, the unsteady flow at high Reynolds numbers through a simplified cascade made up of NACA65-410 blade are simulated with the program written according to the implicit numerical scheme. The reliability and accuracy of the implicit numerical scheme are verified through the satisfactory agreement between the numerical results of the surface pressure coefficient and experimental data. The numerical results indicate that Reynolds number and angle of attack are two primary factors affecting the characteristics of unsteady flow.
文摘The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part of the operator and an implicit Euler discretization of the linear part. Finite difference schemes are used for the spatial part. This finally leads to the numerical solution of a sparse linear system that can be solved efficiently.
文摘A mimetic finite difference scheme for the transient heat equation under Robin’s conditions is presented. The scheme uses second order gradient and divergence mimetic operators, on a staggered grid, to approximate the space derivatives. The temporal derivative is replaced by a first order backward difference approximation to obtain an implicit formulation. The resulting scheme contains nonstandard finite difference stencils. An original convergence analysis by the matrix’s method shows that the proposed scheme is unconditionally stable. A comparative study against standard finite difference schemes, based on central difference or first order one side approximations, reveals the advantages of our scheme without being its implementation more expensive or difficult to achieve.
基金the support of Prince Sultan University for paying the article processing charges(APC)of this publication.
文摘This work aimed to construct an epidemic model with fuzzy parameters.Since the classical epidemic model doesnot elaborate on the successful interaction of susceptible and infective people,the constructed fuzzy epidemicmodel discusses the more detailed versions of the interactions between infective and susceptible people.Thenext-generation matrix approach is employed to find the reproduction number of a deterministic model.Thesensitivity analysis and local stability analysis of the systemare also provided.For solving the fuzzy epidemic model,a numerical scheme is constructed which consists of three time levels.The numerical scheme has an advantage overthe existing forward Euler scheme for determining the conditions of getting the positive solution.The establishedscheme also has an advantage over existing non-standard finite difference methods in terms of order of accuracy.The stability of the scheme for the considered fuzzy model is also provided.From the plotted results,it can beobserved that susceptible people decay by rising interaction parameters.
文摘A mathematical model comprising of nonlinear reaction, diffusion, and convection mechanisms seen in natural and anthropogenic processes is numerically investigated here. It is proposed that a higher order numerical scheme of finite difference method be used in conjunction with an iterative approach in order to solve the nonlinear one dimensional convection-diffusion-reaction equation. To account for the wide variety of physical characteristics and boundary conditions, an iterative approach is presented that yields a reliable and precise solution every time. We examined the accuracy and operational efficiency of two distinct finite difference approaches. The efficiency of the system is determined by comparing the estimated results to the appropriate analytical solution by adhering to established norms. Coherence and convergence were analyzed for each approach. The simulation results demonstrate the efficacy and accuracy of these methods in solving nonlinear convection- diffusion-reaction equations. Convection-diffusion-reaction equation modeling is critical for employing the offered results in heat and mass transport processes.