The time fractional order differential equations are fundamental tools that are used for modeling neuronal dynamics. These equations are obtained by substituting the time derivative of order  where , in the standard equation with the Caputo fractional formula. In this paper, two implicit difference schemes: the linearly Euler implicit and the Crank-Nicolson (CN) finite difference schemes, are employed in solving a one-dimensional time-fractional semilinear equation with Dirichlet boundary conditions. Moreover, the consistency, stability and convergence of the proposed schemes are investigated. We prove that the IEM is unconditionally stable, while CNM is conditionally stable. Furthermore, a comparative study between these two s

This paper is concerned with the numerical blow-up solutions of semi-linear heat equations, where the nonlinear terms are of power type functions, with zero Dirichlet boundary conditions. We use explicit linear and implicit Euler finite difference schemes with a special time-steps formula to compute the blow-up solutions, and to estimate the blow-up times for three numerical experiments. Moreover, we calculate the error bounds and the numerical order of convergence arise from using these methods. Finally, we carry out the numerical simulations to the discrete graphs obtained from using these methods to support the numerical results and to confirm some known blow-up properties for the studied problems.

In this paper, the blow-up solutions for a parabolic problem, defined in a bounded domain, are studied. Namely, we consider the upper blow-up rate estimate for heat equation with a nonlinear Neumann boundary condition defined on a ball in Rn.

This paper is concerned with the blow-up solutions of a system of two reaction-diffusion equations coupled in both equations and boundary conditions. In order to understand how the reaction terms and the boundary terms affect the blow-up properties, the lower and upper blow-up rate estimates are derived. Moreover, the blow-up set under some restricted assumptions is studied.

This paper deals with the blow-up properties of positive solutions to a parabolic system of two heat equations, defined on a ball in  associated with coupled Neumann boundary conditions of exponential type. The upper bounds of blow-up rate estimates are derived. Moreover, it is proved that the blow-up in this problem can only occur on the boundary.

    This paper is concerned with a Coupled Reaction-diffusion system defined in a ball with homogeneous Dirichlet boundary conditions. Firstly, we studied the blow-up set showing that, under some conditions, the blow-up in this problem occurs only at a single point. Secondly, under some restricted assumptions on the reaction terms, we established the upper (lower) blow-up rate estimates. Finally, we considered the Ignition system in general dimensional space as an application to our results.

This paper is concerned with the numerical solutions of the vorticity transport equation (VTE) in two-dimensional space with homogenous Dirichlet boundary conditions. Namely, for this problem, the Crank-Nicolson finite difference equation is derived.  In addition, the consistency and stability of the Crank-Nicolson method are studied. Moreover, a numerical experiment is considered to study the convergence of the Crank-Nicolson scheme and to visualize the discrete graphs for the vorticity and stream functions. The analytical result shows that the proposed scheme is consistent, whereas the numerical results show that the solutions are stable with small space-steps and at any time levels.

In this paper, a numerical approximation for a time fractional one-dimensional bioheat equation (transfer paradigm) of temperature distribution in tissues is introduced. It deals with the Caputo fractional derivative with order for time fractional derivative and new mixed nonpolynomial spline for second order of space derivative. We also analyzed the convergence and stability by employing Von Neumann method for the present scheme.