New a posteriori L∞(L2) and L2(L2)-error estimates of mixed finite element methods for general nonlinear parabolic optimal control problems

We study new a posteriori error estimates of the mixed finite element methods for general optimal control problems governed by nonlinear 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 in L^∞(J; L²Ω)-norm and L²(J; L²Ω)-norm for both the state, the co-state and the control approximation. Such estimates, which seem to be new, are an important step towards developing a reliable adaptive mixed finite element approximation for optimal control problems. Finally, the performance of the posteriori error estimators is assessed by two numerical examples.     

The best <Emphasis Type="Italic">L</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript> approximation of the Laplace operator by linear bounded operators in the classes of functions of two and three variables

Close two-sided estimates are obtained for the best approximation in the space L _p (? ^m ), m = 2 and 3, 1 ≤ p ≤ ∞, of the Laplace operator by linear bounded operators in the class of functions for which the second power of the Laplace operator belongs to the space L _p (? ^m ). We estimate the best constant in the corresponding Kolmogorov inequality and the error of the optimal recovery of values of the Laplace operator on functions from this class given with an error. We present an operator whose deviation from the Laplace operator is close to the best.     

Stabilized finite element method for the viscoelastic Oldroyd fluid flows

In this paper, a new stabilized finite element method based on two local Gauss integrations is considered for the two-dimensional viscoelastic fluid motion equations, arising from the Oldroyd model for the non-Newtonian fluid flows. This new stabilized method presents attractive features such as being parameter-free, or being defined for non-edge-based data structures. It confirms that the lowest equal-order P ₁???P ₁ triangle element and Q ₁???Q ₁ quadrilateral element are compatible. Moreover, the long time stabilities and error estimates for the velocity in H ¹-norm and for the pressure in L ²-norm are obtained. Finally, some numerical experiments are performed, which show that the new method is applied to this model successfully and can save lots of computational cost compared with the standard ones.     

An adaptive <Emphasis Type="Italic">C</Emphasis><Superscript>0</Superscript>IPG method for the Helmholtz transmission eigenvalue problem

The interior penalty methods using C⁰ Lagrange elements (C⁰IPG) developed in the recent decade for the fourth order problems are an interesting topic in academia at present. In this paper, we discuss the adaptive fashion of C⁰IPG method for the Helmholtz transmission eigenvalue problem. We give the a posteriori error indicators for primal and dual eigenfunctions, and prove their reliability and efficiency. We also give the a posteriori error indicator for eigenvalues and design a C⁰IPG adaptive algorithm. Numerical experiments show that this algorithm is efficient and can get the optimal convergence rate.     

Analysis of a new finite difference/local discontinuous Galerkin method for the fractional Cattaneo equation

In this paper, we first present a new finite difference scheme to approximate the time fractional derivatives, which is defined in the sense of Caputo, and give a semidiscrete scheme in time with the truncation error O((Δt)^3?α ), where Δt is the time step size. Then a fully discrete scheme based on the semidiscrete scheme for the fractional Cattaneo equation in which the space direction is approximated by a local discontinuous Galerkin method is presented and analyzed. We prove that the method is unconditionally stable and convergent with order O(h ^k+1 + (Δt)^3?α ), where k is the degree of piecewise polynomial. Numerical examples are also given to confirm the theoretical analysis.     

A priori error estimates of finite volume method for nonlinear optimal control problem

In this paper, we study a priori error estimates for the finite volume element approximation of nonlinear optimal control problem. The schemes use discretizations based on a finite volume method. For the variational inequality, we use the method of the variational discretization concept to obtain the control. Under some reasonable assumptions, we obtain some optimal order error estimates. The approximate order for the state, costate and control variables is O(h ²) or \(O\left( {{h^2}\sqrt {\left| {\ln h} \right|} } \right)\) in the sense of L ²-norm or L ^∞-norm. A numerical experiment is presented to test the theoretical results. Finally, we give some conclusions and future works.     

A new mixed finite element method based on the Crank-Nicolson scheme for Burgers’ equation

In this paper, a new mixed finite element method is used to approximate the solution as well as the flux of the 2D Burgers’ equation. Based on this new formulation, we give the corresponding stable conforming finite element approximation for the P₀² ? P₁ pair by using the Crank-Nicolson time-discretization scheme. Optimal error estimates are obtained. Finally, numerical experiments show the efficiency of the new mixed method and justify the theoretical results.     

Approximation of Gaussian Random Fields: General Results and Optimal Wavelet Representation of the Lévy Fractional Motion

We investigate the approximation rate for certain centered Gaussian fields by a general approach. Upper estimates are proved in the context of so–called Hölder operators and lower estimates follow from the eigenvalue behavior of some related self–adjoint integral operator in a suitable L ₂(μ)–space. In particular, we determine the approximation rate for the Lévy fractional Brownian motion X _H with Hurst parameter H∈(0,1), indexed by a self–similar set T?? ^N of Hausdorff dimension D. This rate turns out to be of order n ^?H/D (log?n)^1/2. In the case T=[0,1] ^N we present a concrete wavelet representation of X _H leading to an approximation of X _H with the optimal rate n ^?H/N (log?n)^1/2.     

Some estimates for the number of vertices of integer polyhedra

A polyhedron is called integer if its every vertex has integer coordinates. We consider integer polyhedra P _I = conv(P ∩ ? ^d ) defined implicitly; that is, no system of linear inequalities is known for P _I but some is known for P. Some estimates are given for the number of vertices of P _I .     

Error estimates of numerical methods for the nonlinear Dirac equation in the nonrelativistic limit regime

We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation(NLDE) in the nonrelativistic limit regime, involving a small dimensionless parameter 0 ε≤ 1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O(ε~2) and O(1) in time and space,respectively. We begin with the conservative Crank-Nicolson finite difference(CNFD) method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ as well as the small parameter 0 ε≤ 1. Based on the error bound, in order to obtain ‘correct' numerical solutions in the nonrelativistic limit regime, i.e., 0 ε■ 1, the CNFD method requests the ε-scalability: τ = O(ε~3) and h= O(ε~(1/2)). Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and timesplitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O(ε~2) and h = O(1) when 0 ε■1. Extensive numerical results are reported to confirm our error estimates.     

Hitting Times of Points and Intervals for Symmetric Lévy Processes

For one-dimensional symmetric Lévy processes, which hit every point with positive probability, we give sharp bounds for the tail function P ^x (T _B >t), where T _B is the first hitting time of B which is either a single point or an interval. The estimates are obtained under some weak type scaling assumptions on the characteristic exponent of the process. We apply these results to prove sharp two-sided estimates of the transition density of the process killed after hitting B.     

Modulation space estimates for damped fractional wave equation

Instead of the L~p estimates,we study the modulation space estimates for the solution to the damped wave equation.Decay properties for both the linear and semilinear equations are obtained.The estimates in modulation space differ in many aspects from those in L~p space.     

Superconvergence and recovery type a posteriori error estimation for hybrid stress finite element method

Superconvergence and recovery type a posteriori error estimators are analyzed for Pian and Sumihara's 4-node hybrid stress quadrilateral finite element method for linear elasticity problems. Superconvergence of order O(h~(1+min){α,1}) is established for both the displacement approximation in H~1-norm and the stress approximation in L~2-norm under a mesh assumption, where α 0 is a parameter characterizing the distortion of meshes from parallelograms to quadrilaterals. Recovery type approximations for the displacement gradients and the stress tensor are constructed, and a posteriori error estimators based on the recovered quantities are shown to be asymptotically exact. Numerical experiments confirm the theoretical results.     

Estimates of Riesz Constants for the Dirac System with an Integrable Potential

We consider the Dirac operator on the interval [0, π] with an integrable potential P = (p_ij (x)) _i,j=1 ² and strongly regular boundary conditions U. It is well known that for any integrable potential P the system {y_n}_n∈Z of root functions of the strongly regular operator L_P,U is a Riesz basis in the space H = L₂[0, π] × L₂[0, π]. We obtain estimates, uniform on every compact set of potentials, of the Riesz constants ||T||||T^?1||, where T is the operator of transition to an orthonormal basis.     

On stabilizability of evolution systems of partial differential equations on ℝ<Superscript><Emphasis Type="Italic">n</Emphasis></Superscript>×[0,+∞) by time-delayed feedback controlsby time-delayed feedback controls

In this work we obtain sufficient conditions for stabilizability by time-delayed feedback controls for the system

$\frac{{\partial w\left( {x,t} \right)}}{{\partial t}} = A(D_x )w(x,t) - A(D_x )u(x,t), x \in \mathbb{R}^n , t > h, $

where D _x =(-i?/?x ₁,...-i?/?x _n ), A(σ) and B(σ) are polynomial matrices (m×m), det B(σ)≡0 on ? ⁿ , w is an unknown function, u(·,t)=P(D _x )w(·,t?h) is a control, h>0. Here P is an infinite differentiable matrix (m×m), and the norm of each of its derivatives does not exceed Γ(1+|σ|²)^γ for some Γ, γ∈? depending on the order of this derivative. Necessary conditions for stabilizability of this system are also obtained. In particular, we study the stabilizability problem for the systems corresponding to the telegraph equation, the wave equation, the heat equation, the Schrödinger equation and another model equation. To obtain these results we use the Fourier transform method, the Lojasiewicz inequality and the Tarski—Seidenberg theorem and its corollaries. To choose an appropriate P and stabilize this system, we also prove some estimates of the real parts of the zeros of the quasipolynomial det {Iλ-A(σ)+B(σ)P(σ)e ^-hλ.

     

On the <Emphasis Type="Italic">σ</Emphasis>-countable compactness of spaces of continuous functions with the set-open topology

For a Tychonoff space X, we obtain a criterion of the σ-countable compactness of the space of continuous functions C(X) with the set-open topology. In particular, for the class of extremally disconnected spaces X, we prove that the space C _λ(X) is σ-countably compact if and only if X is a pseudocompact space, the set X(P) of all P-points of the space X is dense in X, and the family λ consists of finite subsets of the set X(P).     

Relative widths of Sobolev classes in the uniform and integral metrics

Let W_p^r be the Sobolev class consisting of 2π-periodic functions f such that ‖f^(r)‖_p ≤ 1. We consider the relative widths d_n(W_p^r, MW_p^r, L_p), which characterize the best approximation of the class W_p^r in the space L_p by linear subspaces for which (in contrast to Kolmogorov widths) it is additionally required that the approximating functions g should lie in MW_p^r, i.e., ‖g^(r)‖_p ≤ M. We establish estimates for the relative widths in the cases of p = 1 and p = ∞; it follows from these estimates that for almost optimal (with error at most Cn^?r, where C is an absolute constant) approximations of the class W_p^r by linear 2n-dimensional spaces, the norms of the rth derivatives of some approximating functions are not less than cln min(n, r) for large n and r.     

Paranormal Elements in Normed Algebra

For a normed algebra A and natural numbers k we introduce and investigate the ∥ · ∥ closed classes P _k (A). We show that P₁(A) is a subset of P _k (A) for all k. If T in P₁(A), then Tⁿ lies in P₁(A) for all natural n. If A is unital, U, V ∈ A are such that ∥U∥ = ∥V∥ = 1, VU = I and T lies in P _k (A), then UTV lies in P _k (A) for all natural k. Let A be unital, then 1) if an element T in P₁(A) is right invertible, then any right inverse element T^?1 lies in P₁(A); 2) for ßßIßß = 1 the class P₁(A) consists of normaloid elements; 3) if the spectrum of an element T, T ∈ P₁(A) lies on the unit circle, then ∥TX∥ = ∥X∥ for all X ∈ A. If A = B(H), then the class P₁(A) coincides with the set of all paranormal operators on a Hilbert space H.     

Quantum quaternion spheres

For the quantum symplectic group SP _q (2n), we describe the C ^?-algebra of continuous functions on the quotient space S P _q (2n)/S P _q (2n?2) as an universal C ^?-algebra given by a finite set of generators and relations. The proof involves a careful analysis of the relations, and use of the branching rules for representations of the symplectic group due to Zhelobenko. We then exhibit a set of generators of the K-groups of this C ^?-algebra in terms of generators of the C ^?-algebra.     

Boolean Lattices: Ramsey Properties and Embeddings

A subposet Q ^′ of a poset Q is a copy of a poset P if there is a bijection f between elements of P and Q ^′ such that x ≤ y in P iff f(x) ≤ f(y) in Q. For posets P, P ^′, let the poset Ramsey number R(P, P ^′) be the smallest N such that no matter how the elements of the Boolean lattice Q _N are colored red and blue, there is a copy of P with all red elements or a copy of P ^′ with all blue elements. We provide some general bounds on R(P, P ^′) and focus on the situation when P and P ^′ are both Boolean lattices. In addition, we give asymptotically tight bounds for the number of copies of Q _n in Q _N and for a multicolor version of a poset Ramsey number.