Automatic time-step adaptation of the forward Euler method in simulation of models of ion channels and excitable cells and tissue

2008 ◽  
Vol 16 (6) ◽  
pp. 639-644 ◽  
Author(s):  
Topi Korhonen ◽  
Pasi Tavi
2015 ◽  
Vol 25 (11) ◽  
pp. 2015-2042
Author(s):  
Erik Burman

We propose an error analysis for a shock capturing finite element method for the Burgers' equation using the duality theory due to Tadmor. The estimates use a one-sided Lipschitz stability (Lip+-stability) estimate on the discrete solution and are obtained in a weak norm, but thanks to a total variation a priori bound on the discrete solution and an interpolation inequality, error estimates in Lp-norms (1 ≤ p < ∞) are deduced. Both first-order artificial viscosity and a nonlinear shock capturing term that formally is of second order are considered. For the discretization in time we use the forward Euler method. In the numerical section we verify the convergence order of the nonlinear scheme using the forward Euler method and a second-order strong stability preserving Runge–Kutta method. We also study the Lip+-stability property numerically and give some examples of when it holds strictly and when it is violated.


Mathematics ◽  
2019 ◽  
Vol 7 (10) ◽  
pp. 935 ◽  
Author(s):  
Simone Fiori

The present paper recalls a formulation of non-conservative system dynamics through the Lagrange–d’Alembert principle expressed through a generalized Euler–Poincaré form of the system equation on a Lie group. The paper illustrates applications of the generalized Euler–Poincaré equations on the rotation groups to a gyrostat satellite and a quadcopter drone. The numerical solution of the dynamical equations on the rotation groups is tackled via a generalized forward Euler method and an explicit Runge–Kutta integration method tailored to Lie groups.


Author(s):  
R. Vigneswaran ◽  
S. Thilaganathan

We consider a phase space stability error control for numerical simulation of dynamical systems. Standard adaptive algorithm used to solve the linear systems perform well during the finite time of integration with fixed initial condition and performs poorly in three areas. To overcome the difficulties faced the Phase Space Error control criterion was introduced. A new error control was introduced by R. Vigneswaran and Tony Humbries which is generalization of the error control first proposed by some other researchers. For linear systems with a stable hyperbolic fixed point, this error control gives a numerical solution which is forced to converge to the fixed point. In earlier, it was analyzed only for forward Euler method applied to the linear system whose coefficient matrix has real negative eigenvalues. In this paper we analyze forward Euler method applied to the linear system whose coefficient matrix has complex eigenvalues with negative large real parts. Some theoretical results are obtained and numerical results are given.


2016 ◽  
Vol 2016 ◽  
pp. 1-16
Author(s):  
Tomoyuki Suzuki ◽  
Keisuke Takasao ◽  
Noriaki Yamazaki

We consider a one-dimensional Allen-Cahn equation with constraint from the viewpoint of numerical analysis. The constraint is provided by the subdifferential of the indicator function on the closed interval, which is the multivalued function. Therefore, it is very difficult to perform a numerical experiment of our equation. In this paper we approximate the constraint by the Yosida approximation. Then, we study the approximating system of the original model numerically. In particular, we give the criteria for the standard forward Euler method to give the stable numerical experiments of the approximating equation. Moreover, we provide the numerical experiments of the approximating equation.


2018 ◽  
Vol 2018 ◽  
pp. 1-12
Author(s):  
Yingying Zhang ◽  
Yicang Zhou

A discrete population model integrated using the forward Euler method is investigated. The qualitative bifurcation analysis indicates that the model exhibits rich dynamical behaviors including the existence of the equilibrium state, the flip bifurcation, the Neimark-Sacker bifurcation, and two invariant closed curves. The conditions for existence of these bifurcations are derived by using the center manifold and bifurcation theory. Numerical simulations and bifurcation diagrams exhibit the complex dynamical behaviors, especially the occurrence of two invariant closed curves.


2020 ◽  
Vol 20 (4) ◽  
pp. 717-725 ◽  
Author(s):  
Vidar Thomée

AbstractFor a spatially periodic convection-diffusion problem, we analyze a time stepping method based on Lie splitting of a spatially semidiscrete finite element solution on time steps of length k, using the backward Euler method for the diffusion part and a stabilized explicit forward Euler approximation on {m\geq 1} intervals of length {k/m} for the convection part. This complements earlier work on time splitting of the problem in a finite difference context.


2021 ◽  
Vol 10 (6) ◽  
pp. 1239
Author(s):  
Alexandru Cojocaru ◽  
Emilia Burada ◽  
Adrian-Tudor Bălșeanu ◽  
Alexandru-Florian Deftu ◽  
Bogdan Cătălin ◽  
...  

As the average age and life expectancy increases, the incidence of both acute and chronic central nervous system (CNS) pathologies will increase. Understanding mechanisms underlying neuroinflammation as the common feature of any neurodegenerative pathology, we can exploit the pharmacology of cell specific ion channels to improve the outcome of many CNS diseases. As the main cellular player of neuroinflammation, microglia play a central role in this process. Although microglia are considered non-excitable cells, they express a variety of ion channels under both physiological and pathological conditions that seem to be involved in a plethora of cellular processes. Here, we discuss the impact of modulating microglia voltage-gated, potential transient receptor, chloride and proton channels on microglial proliferation, migration, and phagocytosis in neurodegenerative diseases.


2020 ◽  
Vol 20 (4) ◽  
pp. 769-782
Author(s):  
Amiya K. Pani ◽  
Vidar Thomée ◽  
A. S. Vasudeva Murthy

AbstractWe analyze a second-order in space, first-order in time accurate finite difference method for a spatially periodic convection-diffusion problem. This method is a time stepping method based on the first-order Lie splitting of the spatially semidiscrete solution. In each time step, on an interval of length k, of this solution, the method uses the backward Euler method for the diffusion part, and then applies a stabilized explicit forward Euler approximation on {m\geq 1} intervals of length {\frac{k}{m}} for the convection part. With h the mesh width in space, this results in an error bound of the form {C_{0}h^{2}+C_{m}k} for appropriately smooth solutions, where {C_{m}\leq C^{\prime}+\frac{C^{\prime\prime}}{m}}. This work complements the earlier study [V. Thomée and A. S. Vasudeva Murthy, An explicit-implicit splitting method for a convection-diffusion problem, Comput. Methods Appl. Math. 19 2019, 2, 283–293] based on the second-order Strang splitting.


Sign in / Sign up

Export Citation Format

Share Document