In this paper, we discuss the a posteriori error estimates of the mixed finite element method for quadratic optimal control problems governed by linear parabolic equations. The state and the co-state are discretized by the high order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions. We derive a posteriori error estimates for both the state and the control approximation. Such estimates, which are apparently not available in the literature, are an important step towards developing reliable adaptive mixed finite element approximation schemes for the control problem.

In this paper, we give a general proof on convergence estimates for some regularization methods to solve a Cauchy problem for the Laplace equation in a rectangular domain. The regularization methods we considered are: a non-local boundary value problem method, a boundary Tikhonov regularization method and a generalized method. Based on the conditional stability estimates, the convergence estimates for various regularization methods are easily obtained under the simple verifications of some conditions. Numerical results for one example show that the proposed numerical methods are effective and stable.

In this paper we consider polynomials orthogonal with respect to the linear functional defined on the space of all algebraic polynomials by

where α,β> -1/2 are real numbers such that ℓ = |β - α| is a positive integer, and ζ∈ R{0}. We prove the existence of such orthogonal polynomials for some pairs of α and ζ and for all nonnegative integers ℓ. For such orthogonal polynomials we derive three-term recurrence relations and also some differential-difference relations. For such orthogonal polynomials the corresponding quadrature rules of Gaussian type are considered. Also, some numerical examples are included.

In this paper, the finite element approximation of a class of semilinear parabolic optimal control problems with pointwise control constraint is studied. We discretize the state and co-state variables by piecewise linear continuous functions, and the control variable is approximated by piecewise constant functions or piecewise linear discontinuous functions. Some a priori error estimates are derived for both the control and state approximations. The convergence orders are also obtained.

In this paper we study a class of multilevel high order time discretization procedures for the finite difference weighted essential non-oscillatory (WENO) schemes to solve the one-dimensional and two-dimensional shallow water equations with source terms. Multilevel time discretization methods can make full use of computed information by WENO spatial discretization and save CPU cost by holding the former computational values. Extensive simulations are performed, which indicate that, the finite difference WENO schemes with multilevel time discretization can achieve higher accuracy, and are more cost effective than WENO scheme with Runge-Kutta time discretization, while still maintaining nonoscillatory properties.

This paper explores the possibilities of very simple analysis on derivation of spiral regions for a single segment of cubic function matching positional, tangential, and curvature end conditions. Spirals are curves of monotone curvature with constant sign and have the potential advantage that the minimum and maximum curvature exists at their end points. Therefore, spirals are free from singularities, inflection points, and local curvature extrema. These properties make the study of spiral segments an interesting problem both in practical and aesthetic applications, like highway or railway designing or the path planning of non-holonomic mobile robots. Our main contribution is to simplify the procedure of existence methods while keeping it stable and providing flexile constraints for easy applications of spiral segments.

This paper is concerned with the numerical stability of implicit Runge-Kutta methods for nonlinear neutral Volterra delay-integro-differential equations with constant delay. Using a Halanay inequality generalized by Liz and Trofimchuk, we give two sufficient conditions for the stability of the true solution to this class of equations. Runge-Kutta methods with compound quadrature rule are considered. Nonlinear stability conditions for the proposed methods are derived. As an illustration of the application of these investigations, the asymptotic stability of the presented methods for Volterra delay-integro-differential equations are proved under some weaker conditions than those in the literature. An extension of the stability results to such equations with weakly singular kernel is also discussed.

We consider an optimal control problem with an 1D singularly perturbed differential state equation. For solving such problems one uses the enhanced system of the state equation and its adjoint form. Thus, we obtain a system of two convection-diffusion equations. Using linear finite elements on adapted grids we treat the effects of two layers arising at different boundaries of the domain. We proof uniform error estimates for this method on meshes of Shishkin type. We present numerical results supporting our analysis.

The immersed interface method is modified to compute Schrödinger equation with discontinuous potential. By building the jump conditions of the solution into the finite difference approximation near the interface, this method can give at least second order convergence rate for the numerical solution on uniform cartesian grids. The accuracy of this algorithm is tested via several numerical examples.