Fourth order time-stepping for low dispersion Korteweg-de Vries and nonlinear Schrödinger equations.
We introduce fractional-order Bessel functions (FBFs) to obtain an approximate solution for various kinds of differential equations. Our main aim is to consider the new functions based on Bessel polynomials to the fractional calculus. To calculate derivatives and integrals, we use Caputo fractional derivatives and Riemann-Liouville fractional integral definitions. Then, operational matrices of fractional-order derivatives and integration for FBFs are derived. Also, we discuss an error estimate between...
In this article, we analyze the stability of various numerical schemes for differential models of viscoelastic fluids. More precisely, we consider the prototypical Oldroyd-B model, for which a free energy dissipation holds, and we show under which assumptions such a dissipation is also satisfied for the numerical scheme. Among the numerical schemes we analyze, we consider some discretizations based on the log-formulation of the Oldroyd-B system proposed by Fattal and Kupferman in [J. Non-Newtonian...
The error analysis of preconditioned waveform relaxation iterations for differential systems is presented. This analysis extends and refines previous results by Burrage, Jackiewicz, Nørsett and Renaut by incorporating all terms in the expansion of the error of waveform relaxation iterations in the Laplace transform domain. Lower bounds for the size of the window of rapid convergence are also obtained. The theory is illustrated for waveform relaxation methods applied to differential systems resulting...
This paper aims at the development of numerical schemes for nonlinear reaction diffusion problems with a convection that blows up in a finite time. A full discretization of this problem that preserves the blow - up property is presented as well as a numerical simulation. Efficiency of the method is derived via a numerical comparison with a classical scheme based on the Runge Kutta scheme.
We present a fully adaptive multiresolution scheme for spatially one-dimensional quasilinear strongly degenerate parabolic equations with zero-flux and periodic boundary conditions. The numerical scheme is based on a finite volume discretization using the Engquist-Osher numerical flux and explicit time stepping. An adaptive multiresolution scheme based on cell averages is then used to speed up the CPU time and the memory requirements of the underlying finite volume scheme, whose first-order...
A posteriori error estimates for a nonlinear parabolic problem are introduced. A fully discrete scheme is studied. The space discretization is based on a concept of hierarchical finite element basis functions. The time discretization is done using singly implicit Runge-Kutta method (SIRK). The convergence of the effectivity index is proven.
The Poisson-Boltzmann (PB) model is an effective approach for the electrostatics analysis of solvated biomolecules. The nonlinearity associated with the PB equation is critical when the underlying electrostatic potential is strong, but is extremely difficult to solve numerically. In this paper, we construct two operator splitting alternating direction implicit (ADI) schemes to efficiently and stably solve the nonlinear PB equation in a pseudo-transient continuation approach. The operator splitting...
We consider an initial and Dirichlet boundary value problem for a fourth-order linear stochastic parabolic equation, in one space dimension, forced by an additive space-time white noise. Discretizing the space-time white noise a modelling error is introduced and a regularized fourth-order linear stochastic parabolic problem is obtained. Fully-discrete approximations to the solution of the regularized problem are constructed by using, for discretization in space, a Galerkin finite element method...