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 \(v\), including \(v\in L^2(\Omega )\). A space semidiscrete Galerkin scheme using continuous piecewise linear finite elements is developed, and optimal with respect to initial data regularity error estimates for the finite element approximations are derived. Further, two fully discrete schemes based on the backward Euler method and second-order backward difference method and the related convolution quadrature are developed, and optimal error estimates are derived for the fully discrete approximations for both smooth and nonsmooth initial data. Numerical results for one- and two-dimensional examples with smooth and nonsmooth initial data are presented to illustrate the efficiency of the method, and to verify the convergence theory.

The abstract Cauchy problem for the distributed order fractional evolution equation in the Caputo and in the Riemann–Liouville sense is studied for operators generating a strongly continuous one-parameter semigroup on a Banach space. Continuous as well as discrete distribution of fractional time-derivatives of order less than one is considered. The problem is reformulated as an abstract Volterra integral equation. It is proven that its kernel satisfies certain complete monotonicity properties. Based on these properties, the well-posedness of the problem is established and a series expansion of the solution is obtained. In case of ordered Banach space this representation implies positivity of the solution operator. In addition, a subordination formula is obtained.

The fractional order differential equation \(u'(t)=Au(t)+\gamma D_t^{\alpha} Au(t)+f(t), \ t>0\), \(u(0)=a\in X\) is studied, where \(A\) is an operator generating a strongly continuous one-parameter semigroup on a Banach space \(X\), \(D_t^{\alpha}\) is the Riemann–Liouville fractional derivative of order \(\alpha \in (0,1)\), \(\gamma>0\) and \(f\) is an \(X\)-valued function. Equations of this type appear in the modeling of unidirectional viscoelastic flows. Well-posedness is proven, and a subordination identity is obtained relating the solution operator of the considered problem and the \(C_{0}\)-semigroup, generated by the operator \(A\). As an example, the Rayleigh–Stokes problem for a generalized second-grade fluid is considered.