We study the Rayleigh-Stokes problem for a generalized second-grade fluid which involves a Riemann-Liouville fractional derivative in time, and present an analysis of the problem in the continuous, space semidiscrete and fully discrete formulations. We establish the Sobolev regularity of the homogeneous problem for both smooth and nonsmooth initial data [Formula: see text], including [Formula: ...

This work develops adaptive time stepping algorithms for the approximation of a functional of a diffusion with jumps based on a jump augmented Monte Carlo Euler–Maruyama method, which achieve a prescribed precision. The main result is the derivation of new expansions for the time discretization error, with computable leading order term in a posteriori form, which are based on stochastic flows a...

Abstract. An adaptive algorithm, based on residual type a posteriori indicators of errors measured in L∞(L2) and L2(L2) norms, for a numerical scheme consisting of implicit Euler method in time and discontinuous Galerkin method in space for linear parabolic fourth order problems is presented. The a posteriori analysis is performed for convex domains in two and three space dimensions for local s...

In this paper we present the efficient parallel implementation of elastoplastic problems based on the TFETI (Total Finite Element Tearing and Interconnecting) domain decomposition method. This approach allow us to use parallel solution and compute this nonlinear problem on the supercomputers and decrease the solution time and compute problems with millions of DOFs. In our approach we consider a...

A second-order singularly perturbed parabolic equation in one space dimension is considered. For this equation, we give computable a posteriori error estimates in the maximum norm for two semidiscretisations in time and a full discretisation using P1 FEM in space. Both the Backward-Euler method and the Crank-Nicolson method are considered, and certain critical details of the implementation are ...

By applying the Euler–Rayleigh method to a specific representation of Jacobi polynomials as hypergeometric functions, we obtain new bounds for their largest zeros. In particular, derive upper and lower bound $$1-x_{nn}^2(\lambda )$$ , with $$x_{nn}(\lambda being zero n-th ultraspherical polynomial $$P_n^{(\lambda )}$$ . For every fixed $$\lambda >-1/2$$ limit ratio our does not exceed 1.6. This...

Abstract. We derive optimal order, residual-based a posteriori error estimates for time discretizations by the two–step BDF method for linear parabolic equations. Appropriate reconstructions of the approximate solution play a key role in the analysis. To utilize the BDF method we employ one step by both the trapezoidal method or the backward Euler scheme. Our a posteriori error estimates are of...

In standard deformable object simulation in computer animation, all the mesh elements or vertices are timestepped synchronously, i.e., under the same timestep. Previous asynchronous methods have been largely limited to explicit integration. We demonstrate how to perform spatially-varying timesteps for the widely popular implicit backward Euler integrator. Spatiallyvarying timesteps are useful w...

H(div) conforming and discontinuous Galerkin (DG) methods are designed for incompressible Euler’s equation in two and three dimension. Error estimates are proved for both the semi-discrete method and fully-discrete method using backward Euler time stepping. Numerical examples exhibiting the performance of the methods are given.

This paper studies the chaotic behavior of a series ferroresonant circuit. At first, we derive a state variable equation containing a nonlinear hysteretic circuit element represented by a Chua-type magnetization model. Second, the backward Euler method with automatic modification gives transient response of the state variable equation. The characteristic values of the state transition matrix ar...