scholarly journals Adaptive Optimal -Stage Runge-Kutta Methods for Solving Reaction-Diffusion-Chemotaxis Systems

2011 ◽  
Vol 2011 ◽  
pp. 1-25
Author(s):  
Jui-Ling Yu
Keyword(s):  
Method Of Lines ◽  
Reaction Diffusion ◽  
Stability And Convergence ◽  
Optimal Time ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Time Stepping ◽  
Application Problem ◽  
Runge Kutta

We present a class of numerical methods for the reaction-diffusion-chemotaxis system which is significant for biological and chemistry pattern formation problems. To solve reaction-diffusion-chemotaxis systems, efficient and reliable numerical algorithms are essential for pattern generations. Along with the implementation of the method of lines, implicit or semi-implicit schemes are typical time stepping solvers to reduce the effect on time step constrains due to the stability condition. However, these two schemes are usually difficult to employ. In this paper, we propose an adaptive optimal time stepping strategy for the explicit -stage Runge-Kutta method to solve reaction-diffusion-chemotaxis systems. Instead of relying on empirical approaches to control the time step size, variable time step sizes are given explicitly. Yet, theorems about stability and convergence of the algorithm are provided in analyzing robustness and efficiency. Numerical experiment results on a testing problem and a real application problem are shown.

scholarly journals A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations

Complexity ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-15
Author(s):  
Philku Lee ◽  
George V. Popescu ◽  
Seongjai Kim
Keyword(s):  
Numerical Solution ◽  
Reaction Diffusion ◽  
Second Order ◽  
Reaction Diffusion Equations ◽  
Step Size ◽  
Time Step ◽  
Potential Oscillations ◽  
Time Step Size ◽  
Time Stepping ◽  
Spurious Oscillations

After a theory of morphogenesis in chemical cells was introduced in the 1950s, much attention had been devoted to the numerical solution of reaction-diffusion (RD) partial differential equations (PDEs). The Crank–Nicolson (CN) method has been a common second-order time-stepping procedure. However, the CN method may introduce spurious oscillations for nonsmooth data unless the time step size is sufficiently small. This article studies a nonoscillatory second-order time-stepping procedure for RD equations, called a variable-θmethod, as a perturbation of the CN method. In each time level, the new method detects points of potential oscillations to implicitly resolve the solution there. The proposed time-stepping procedure is nonoscillatory and of a second-order temporal accuracy. Various examples are given to show effectiveness of the method. The article also performs a sensitivity analysis for the numerical solution of biological pattern forming models to conclude that the numerical solution is much more sensitive to the spatial mesh resolution than the temporal one. As the spatial resolution becomes higher for an improved accuracy, the CN method may produce spurious oscillations, while the proposed method results in stable solutions.

scholarly journals Analysis of a Decoupled Time-Stepping Scheme for Evolutionary Micropolar Fluid Flows

2016 ◽  
Vol 2016 ◽  
pp. 1-13 ◽  
Author(s):  
S. S. Ravindran
Keyword(s):  
Micropolar Fluid ◽  
Stokes Equations ◽  
Fluid Model ◽  
Navier Stokes ◽  
Optimal Order ◽  
Step Size ◽  
Time Step ◽  
Order Error ◽  
Time Step Size ◽  
Time Stepping

Micropolar fluid model consists of Navier-Stokes equations and microrotational velocity equations describing the dynamics of flows in which microstructure of fluid is important. In this paper, we propose and analyze a decoupled time-stepping algorithm for the evolutionary micropolar flow. The proposed method requires solving only one uncoupled Navier-Stokes and one microrotation subphysics problem per time step. We derive optimal order error estimates in suitable norms without assuming any stability condition or time step size restriction.

scholarly journals A Minimizing Movement approach to a class of scalar reaction-diffusion equations

2020 ◽  
Author(s):  
Florentine Catharina Fleißner
Keyword(s):  
General Setting ◽  
Gradient Flow ◽  
Reaction Diffusion ◽  
Reaction Diffusion Equations ◽  
Diffusion Equations ◽  
Step Size ◽  
Time Step ◽  
Flux Boundary Condition ◽  
Time Step Size ◽  
Kantorovich Distance

The purpose of this paper is to introduce a Minimizing Movement approach to scalar reaction-diffusion equations of the form \partial_t u \ = \ \Lambda\cdot \mathrm{div}[u(\nabla F'(u) + \nabla V)] \ - \ \Sigma\cdot (F'(u) + V) u, \quad \text{ in } (0, +\infty)\times\Omega, with parameters $\Lambda, \Sigma > 0$ and no-flux boundary condition u(\nabla F'(u) + \nabla V)\cdot {\sf n} \ = \ 0, \quad \text{ on } (0, +\infty)\times\partial\Omega, which is built on their gradient-flow-like structure in the space $\mathcal{M}(\bar{\Omega})$ of finite nonnegative Radon measures on $\bar{\Omega}\subset\xR^d$, endowed with the recently introduced Hellinger-Kantorovich distance $\HK_{\Lambda, \Sigma}$. It is proved that, under natural general assumptions on $F: [0, +\infty)\to\xR$ and $V:\bar{\Omega}\to\xR$, the Minimizing Movement scheme \mu_\tau^0:=u_0\mathscr{L}^d \in\mathcal{M}(\bar{\Omega}), \quad \mu_\tau^n \text{ is a minimizer for } \mathcal{E}(\cdot)+\frac{1}{2\tau}\HK_{\Lambda, \Sigma}(\cdot, \mu_\tau^{n-1})^2, \ n\in\xN, for \mathcal{E}: \mathcal{M}(\bar{\Omega}) \to (-\infty, +\infty], \ \mathcal{E}(\mu):= \begin{cases} \int_\Omega{[F(u(x))+V(x)u(x)]\xdif x} &\text{ if } \mu=u\mathscr{L}^d, \\ +\infty &\text{ else}, \end{cases} yields weak solutions to the above equation as the discrete time step size $\tau\downarrow 0$. Moreover, a superdifferentiability property of the Hellinger-Kantorovich distance $\HK_{\Lambda, \Sigma}$, which will play an important role in this context, is established in the general setting of a separable Hilbert space.

The Dispersion in One- and Two-Dimensional Finite Element Solutions to the Heat Equation

2002 ◽  
Author(s):  
W. Dauksher ◽  
A. F. Emery
Keyword(s):  
Finite Element ◽  
Heat Equations ◽  
Two Dimensional ◽  
Step Size ◽  
Matrix Formulation ◽  
Time Step ◽  
Time Step Size ◽  
Time Stepping ◽  
Dimensional Finite Element ◽  
The One

The dispersive errors in the finite element solution to the one- and two-dimensional heat equations are examined as a function of element type and size, capacitance matrix formulation, time stepping scheme and time step size.

scholarly journals Error Estimates of a Continuous Galerkin Time Stepping Method for Subdiffusion Problem

2021 ◽  
Vol 88 (3) ◽  
Author(s):  
Yuyuan Yan ◽  
Bernard A. Egwu ◽  
Zongqi Liang ◽  
Yubin Yan
Keyword(s):  
Piecewise Linear ◽  
Time Discretization ◽  
Convergence Order ◽  
Step Size ◽  
Time Step ◽  
Transform Method ◽  
Nonsmooth Data ◽  
Time Step Size ◽  
Time Stepping ◽  
Continuous Galerkin

AbstractA continuous Galerkin time stepping method is introduced and analyzed for subdiffusion problem in an abstract setting. The approximate solution will be sought as a continuous piecewise linear function in time t and the test space is based on the discontinuous piecewise constant functions. We prove that the proposed time stepping method has the convergence order $$O(\tau ^{1+ \alpha }), \, \alpha \in (0, 1)$$ O ( τ 1 + α ) , α ∈ ( 0 , 1 ) for general sectorial elliptic operators for nonsmooth data by using the Laplace transform method, where $$\tau $$ τ is the time step size. This convergence order is higher than the convergence orders of the popular convolution quadrature methods (e.g., Lubich’s convolution methods) and L-type methods (e.g., L1 method), which have only $$O(\tau )$$ O ( τ ) convergence for the nonsmooth data. Numerical examples are given to verify the robustness of the time discretization schemes with respect to data regularity.

scholarly journals Energy-corrected FEM and explicit time-stepping for parabolic problems

2019 ◽  
Vol 53 (6) ◽  
pp. 1893-1914
Author(s):  
Piotr Swierczynski ◽  
Barbara Wohlmuth
Keyword(s):  
Finite Element ◽  
Finite Element Approximation ◽  
Regularity Of Solutions ◽  
Parabolic Problems ◽  
Element Approximation ◽  
Computational Domain ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Time Stepping

The presence of corners in the computational domain, in general, reduces the regularity of solutions of parabolic problems and diminishes the convergence properties of the finite element approximation introducing a so-called “pollution effect”. Standard remedies based on mesh refinement around the singular corner result in very restrictive stability requirements on the time-step size when explicit time integration is applied. In this article, we introduce and analyse the energy-corrected finite element method for parabolic problems, which works on quasi-uniform meshes, and, based on it, create fast explicit time discretisation. We illustrate these results with extensive numerical investigations not only confirming the theoretical results but also showing the flexibility of the method, which can be applied in the presence of multiple singular corners and a three-dimensional setting. We also propose a fast explicit time-stepping scheme based on a piecewise cubic energy-corrected discretisation in space completed with mass-lumping techniques and numerically verify its efficiency.

scholarly journals Super-Time-Stepping Scheme for Option Pricing

2020 ◽  
Vol 13 (13) ◽  
pp. 51-54
Author(s):  
Kedar Nath Uprety ◽  
Harithar Khanal ◽  
Ananta Upreti
Keyword(s):  
Option Pricing ◽  
Euler Method ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Time Stepping ◽  
Black Scholes ◽  
Explicit Finite Difference ◽  
The Stability ◽  
Explicit Finite Difference Method

We solve the Black - Scholes equation for option pricing numerically using an Explicit finite difference method. To overcome the stability restriction of the explicit scheme for parabolic partial differential equations in the time step size Courant-Friedrichs-Lewy (CFL) condition, we employ a Super Time Stepping (STS) strategy based on modified Chebyshev polynomial. The numerical results show that the STS scheme boasts of large efficiency gains compared to the standard explicit Euler method.

scholarly journals A decoupled Crank-Nicolson time-stepping scheme for thermally coupled magneto-hydrodynamic system

2017 ◽  
Vol 8 (1) ◽  
pp. 43-62
Author(s):  
S. S. Ravindran
Keyword(s):  
Heat Equation ◽  
Mixed Finite Element ◽  
Mhd Equations ◽  
Optimal Order ◽  
Step Size ◽  
Time Step ◽  
Order Error ◽  
Time Step Size ◽  
Time Stepping ◽  
Crank Nicolson

Thermally coupled magneto-hydrodynamics (MHD) studies the dynamics of electro-magnetically and thermally driven flows,involving MHD equations coupled with heat equation. We introduce a partitioned method that allows one to decouplethe MHD equations from the heat equation at each time step and solve them separately. The extrapolated Crank-Nicolson time-stepping scheme is used for time discretizationwhile mixed finite element method is used for spatial discretization. We derive optimal order error estimates in suitable norms without assuming any stability condition or restrictions on the time step size. We prove the unconditional stability of the scheme. Numerical experiments are used to illustrate the theoretical results.

On the Behavior of a Numerical Approximation to the Rotatory Inertia and Transverse Shear Plate

1973 ◽  
Vol 40 (4) ◽  
pp. 977-982 ◽  
Author(s):  
R. D. Krieg
Keyword(s):  
Transverse Shear ◽  
Plate Thickness ◽  
Critical Time ◽  
Stability And Convergence ◽  
Numerical Verification ◽  
Rotatory Inertia ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Shear Plate

A set of finite-difference equations which approximate the plane-strain motion of a rotatory inertia and transverse shear plate are examined in detail. The equations are stated in such a form that the usual norm is an energy expression. Proper posedness, consistency, stability, and convergence are all examined. The von Neumann stability analysis is supplemented by a numerical verification. The critical time step size is shown to fall into three different regimes depending on the ratio of mesh size to plate thickness. The largest allowable time step size is found to be dictated not by mesh spacing, but rather by physical dimensions of the plate.

scholarly journals Implicit–explicit (IMEX) Runge–Kutta methods for non-hydrostatic atmospheric models

2018 ◽  
Vol 11 (4) ◽  
pp. 1497-1515 ◽  
Author(s):  
David J. Gardner ◽  
Jorge E. Guerra ◽  
François P. Hamon ◽  
Daniel R. Reynolds ◽  
Paul A. Ullrich ◽  
...  
Keyword(s):  
Acoustic Waves ◽  
Time Integration ◽  
Coupled System ◽  
Atmospheric Model ◽  
Step Size ◽  
Time Step ◽  
Solution Quality ◽  
Time Step Size ◽  
Runge Kutta ◽  
The Impact

Abstract. The efficient simulation of non-hydrostatic atmospheric dynamics requires time integration methods capable of overcoming the explicit stability constraints on time step size arising from acoustic waves. In this work, we investigate various implicit–explicit (IMEX) additive Runge–Kutta (ARK) methods for evolving acoustic waves implicitly to enable larger time step sizes in a global non-hydrostatic atmospheric model. The IMEX formulations considered include horizontally explicit – vertically implicit (HEVI) approaches as well as splittings that treat some horizontal dynamics implicitly. In each case, the impact of solving nonlinear systems in each implicit ARK stage in a linearly implicit fashion is also explored.The accuracy and efficiency of the IMEX splittings, ARK methods, and solver options are evaluated on a gravity wave and baroclinic wave test case. HEVI splittings that treat some vertical dynamics explicitly do not show a benefit in solution quality or run time over the most implicit HEVI formulation. While splittings that implicitly evolve some horizontal dynamics increase the maximum stable step size of a method, the gains are insufficient to overcome the additional cost of solving a globally coupled system. Solving implicit stage systems in a linearly implicit manner limits the solver cost but this is offset by a reduction in step size to achieve the desired accuracy for some methods. Overall, the third-order ARS343 and ARK324 methods performed the best, followed by the second-order ARS232 and ARK232 methods.

Sign in / Sign up

Export Citation Format

Share Document