Goal of this paper is to suitably combine a model with an anisotropic mesh adaptation for the numerical simulation of nonlinear advection-diffusion-reaction systems and incompressible flows in ecological and environmental applications. Using the reduced-basis method terminology, the proposed approach leads to a noticeable computational saving of the online phase with respect to the resolution of the reference model on nonadapted grids. The search of a suitable adapted model/mesh pair is to be meant, instead, in an offline fashion.

In this paper, we investigate the superconvergence property and the L∞-error estimates of mixed finite element methods for a semilinear elliptic control problem. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. We derive some superconvergence results for the control variable. Moreover, we derive L∞-error estimates both for the control variable and the state variables. Finally, a numerical example is given to demonstrate the theoretical results.

An inverse problem of reconstructing the initial condition for a time fractional diffusion equation is investigated. On the basis of the optimal control framework, the uniqueness and first order necessary optimality condition of the minimizer for the objective functional are established, and a time-space spectral method is proposed to numerically solve the resulting minimization problem. The contribution of the paper is threefold: 1) a priori error estimate for the spectral approximation is derived; 2) a conjugate gradient optimization algorithm is designed to efficiently solve the inverse problem; 3) some numerical experiments are carried out to show that the proposed method is capable to find out the optimal initial condition, and that the convergence rate of the method is exponential if the optimal initial condition is smooth.

In this paper, we propose the Laguerre spectral method for high order problems with mixed inhomogeneous boundary conditions. It is also available for approximated solutions growing fast at infinity. The spectral accuracy is proved. Numerical results demonstrate its high effectiveness.

In this paper we extend the source transfer domain decomposition method (STDDM) introduced by the authors to solve the Helmholtz problems in two-layered media, the Helmholtz scattering problems with bounded scatterer, and Helmholtz problems in 3D unbounded domains. The STDDM is based on the decomposition of the domain into non-overlapping layers and the idea of source transfer which transfers the sources equivalently layer by layer so that the solution in the final layer can be solved using a PML method defined locally outside the last two layers. The details of STDDM is given for each extension. Numerical results are presented to demonstrate the efficiency of STDDM as a preconditioner for solving the discretization problem of the Helmholtz problems considered in the paper.

The z-transform is introduced to analyze a full discretization method for a partial integro-differential equation (PIDE) with a weakly singular kernel. In this method, spectral collocation is used for the spatial discretization, and, for the time stepping, the finite difference method combined with the convolution quadrature rule is considered. The global stability and convergence properties of complete discretization are derived and numerical experiments are reported.