Asymptotic analysis and error estimates for an optimal control problem with oscillating boundaries

In this article, we consider a distributed optimal control problem associated with the Laplacian in a domain with rapidly oscillating boundary. For simplicity, we consider a rectangular region in 2d with oscillations on one part of the boundary. We consider two types of functionals, namely a functional involving the L ²-norm of the state variable and another one involving its H ¹-norm. The homogenization of the optimality system is obtained and then we derive appropriate error estimates in both cases.     

Sampling of Discontinuous Dirac Systems

Homogenization and error analysis of an optimal interior control problem in the framework of Stokes’ system, on a domain with rapidly oscillating boundary, are the subject matters of this article. We consider a three dimensional domain constituted of a parallelepiped with a large number of rectangular cylinders at the top of it. An interior control is applied in a proper subdomain of the parallelepiped, away from the oscillating volume. We consider two types of functionals, namely a functional involving the L ²-norm of the state variable and another one involving its H ¹-norm. The asymptotic analysis of optimality systems for both cases, when the cross sectional area of the rectangular cylinders tends to zero, is done here. Our major contribution is to derive error estimates for the state, the co-state and the associated pressures, in appropriate functional spaces.     

On the convergence and superconvergence for a class of two-dimensional time fractional reaction-subdiffusion equations

This paper analyzes a class of two-dimensional (2-D) time fractional reaction-subdiffusion equations with variable coefficients. The high-order L2-1_σ time-stepping scheme on graded meshes is presented to deal with the weak singularity at the initial time t = 0, and the bilinear finite element method (FEM) on anisotropic meshes is used for spatial discretization. Using the modified discrete fractional Grönwall inequality, and combining the interpolation operator and the projection operator, the L²-norm error estimation and H¹-norm superclose results are rigorously proved. The superconvergence result in the H¹-norm is derived by applying the interpolation postprocessing technique. Finally, numerical examples are presented to verify the validation of our theoretical analysis.     

Correctors for some asymptotic problems

In the theory of anisotropic singular perturbation boundary value problems, the solution u _ɛ does not converge, in the H ¹-norm on the whole domain, towards some u ₀. In this paper we construct correctors to have good approximations of u _ɛ in the H ¹-norm on the whole domain. Since the anisotropic singular perturbation problems can be connected to the study of the asymptotic behaviour of problems defined in cylindrical domains becoming unbounded in some directions, we transpose our results for such problems.     

Stationary solutions of the non-linear Boltzmann equation in a bounded spatial domain

A contraction mapping (or, alternatively, an implicit function theory) argument is applied in combination with the Fredholm alternative to prove the existence of a unique stationary solution of the non-linear Boltzmann equation on a bounded spatial domain under a rather general reflection law at the piecewise C¹ boundary. The boundary data are to be small in a weighted L_∞-norm.     

Weak Galerkin finite element methods for electric interface model with nonhomogeneous jump conditions

In this paper, the weak Galerkin finite element method (WG-FEM) is applied to a pulsed electric model arising in biological tissue when a biological cell is exposed to an electric field. A fitted WG-FEM is proposed to approximate the voltage of the pulsed electric model across the physical media involving an electric interface (surface membrane), and heterogeneous permittivity and a heterogeneous conductivity. This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes. Optimal pointwise-in-time error estimates in L²-norm and H¹-norm are shown to hold for the semidiscrete scheme even if the regularity of the solution is low on the whole domain. Furthermore, a fully discrete approximation based on backward Euler scheme is analyzed and related optimal error estimates are derived.     

Mapped Finite Elements on Hexahedra. Necessary and Sufficient Conditions for Optimal Interpolation Errors

This paper considers finite elements which are defined on hexahedral cells via a reference transformation which is in general trilinear. For affine reference mappings, the necessary and sufficient condition for an interpolation order O(h ^k+1) in the L ²-norm and O(h ^k ) in the H ¹-norm is that the finite dimensional function space on the reference cell contains all polynomials of degree less than or equal to k. The situation changes in the case of a general trilinear reference transformation. We will show that on general meshes the necessary and sufficient condition for an optimal order for the interpolation error is that the space of polynomials of degree less than or equal to k in each variable separately is contained in the function space on the reference cell. Furthermore, we will show that this condition can be weakened on special families of meshes. These families which are obtained by applying usual refinement techniques can be characterized by the asymptotic behaviour of the semi-norms of the reference mapping.     

Insensitizing controls for a nonlinear parabolic equation with a nonlinear term involving the state and the gradient in unbounded domains

In this paper, we consider the existence of insensitizing control for a semilinear heat equation involving gradient terms in unbounded domain Ω. In this case, we prove the existence of controls insensitizing the L²-norm of the observation of the solution in an open subset of the domain. The proofs of the main results in this paper involve such inequalities and rely on the study of these linear problems and appropriate fixed point arguments.     

Optimal convergence analysis of an immersed interface finite element method

We analyze an immersed interface finite element method based on linear polynomials on noninterface triangular elements and piecewise linear polynomials on interface triangular elements. The flux jump condition is weakly enforced on the smooth interface. Optimal error estimates are derived in the broken H ¹-norm and L ²-norm.     

The pointwise estimates of a conservative difference scheme for Burgers' equation

In this article, we are concerned with the numerical analysis of a nonlinear implicit difference scheme for Burgers' equation. A priori estimation of the analytical solution is provided in the sense of L^∞ -norm when the initial value is bounded in H¹-norm. Conservation, boundedness, and unique solvability are proved at length. Inspired by the method of the priori estimation for the analytical solution, we prove the convergence and stability of the difference scheme in L^∞ -norm. Finally, numerical examples are carried out to verify our theoretical results.     

Extrapolation cascadic multigrid method on piecewise uniform grid

The triangular linear fnite elements on piecewise uniform grid for an elliptic problem in convex polygonal domain are discussed.Global superconvergence in discrete H1-norm and global extrapolation in discrete L2-norm are proved.Based on these global estimates the conjugate gradient method(CG)is efective,which is applied to extrapolation cascadic multigrid method(EXCMG).The numerical experiments show that EXCMG is of the global higher accuracy for both function and gradient.     

Characterization of functions as radiation patterns in linear elasticity

In this paper a necessary and sufficient condition for a pair of vector functions to be radiation patterns is presented. More precisely, it is proved that two vector functions, the first in the radial direction and the second in the tangential one, are radiation patterns if and only if there are two entire harmonic vector functions whose radial and tangential projections, respectively, are identical with the previous functions on the unit sphere and whose L²-norm over a sphere of radius R is a function of exponential type in the variable R.     

Convergence and superconvergence analysis for nonlinear delay reaction–diffusion system with nonconforming finite element

In this article, we propose and analyze several numerical methods for the nonlinear delay reaction–diffusion system with smooth and nonsmooth solutions, by using Quasi-Wilson nonconforming finite element methods in space and finite difference methods (including uniform and nonuniform L1 and L2-1_σ schemes) in time. The optimal convergence results in the senses of L²-norm and broken H¹-norm, and H¹-norm superclose results are derived by virtue of two types of fractional Grönwall inequalities. Then, the interpolation postprocessing technique is used to establish the superconvergence results. Moreover, to improve computational efficiency, fast algorithms by using sum-of-exponential technique are built for above proposed numerical schemes. Finally, we present some numerical experiments to confirm the theoretical correctness and show the effectiveness of the fast algorithms.     

On the k-sample Behrens-Fisher problem for high-dimensional data

For several decades, much attention has been paid to the two-sample Behrens-Fisher (BF) problem which tests the equality of the means or mean vectors of two normal populations with unequal variance/covariance structures. Little work, however, has been done for the k-sample BF problem for high dimensional data which tests the equality of the mean vectors of several high-dimensional normal populations with unequal covariance structures. In this paper we study this challenging problem via extending the famous Scheffe’s transformation method, which reduces the k-sample BF problem to a one-sample problem. The induced one-sample problem can be easily tested by the classical Hotelling’s T ² test when the size of the resulting sample is very large relative to its dimensionality. For high dimensional data, however, the dimensionality of the resulting sample is often very large, and even much larger than its sample size, which makes the classical Hotelling’s T ² test not powerful or not even well defined. To overcome this difficulty, we propose and study an L ²-norm based test. The asymptotic powers of the proposed L ²-norm based test and Hotelling’s T ² test are derived and theoretically compared. Methods for implementing the L ²-norm based test are described. Simulation studies are conducted to compare the L ²-norm based test and Hotelling’s T ² test when the latter can be well defined, and to compare the proposed implementation methods for the L ²-norm based test otherwise. The methodologies are motivated and illustrated by a real data example. The work was supported by the National University of Singapore Academic Research Grant (Grant No. R-155-000-085-112)     

Convergence and superconvergence analysis of finite element methods on graded meshes for singularly and semisingularly perturbed reaction–diffusion problems

The bilinear finite element methods on appropriately graded meshes are considered both for solving singular and semisingular perturbation problems. In each case, the quasi-optimal order error estimates are proved in the ε-weighted H¹-norm uniformly in singular perturbation parameter ε, up to a logarithmic factor. By using the interpolation postprocessing technique, the global superconvergent error estimates in ε-weighted H¹-norm are obtained. Numerical experiments are given to demonstrate validity of our theoretical analysis.     

An improvement and an extension of the Elzinga & Hearn's algorithm to the 1-center problem in ℝn withl 2b-normswithl 2b-norms

Summary In this paper we deal with the 1-center problem in ℝⁿ when the distance is measured by anyl _2b-norm. This type of norm generalizes the Euclidean norm (l ₂-norm) and can be used to estimate road distances or travel times in Locational Analysis, and to measure dissimilarities between data in Cluster Analysis. The problem is to find the smallestb-ellipsoid containing a given finite setA of points in ℝⁿ, which generalizes the one of finding the smallest sphere containingA (1-center problem with thel ₂-norm). We show that this problem has a unique optimal solution. For thel ₂-norm, we use the Elzinga-Hearn algorithm. New starting rules are proposed to initialize and to improve the algorithm. On the other hand, the Elzinga-Hearn algorithm is extended to solve the problem withl _2b-norms. Computational results are given for six differentl _2b-norms, when these new starting rules are used in order to show which is the best starting rule. Problems of up to 5.000 points in ℝⁿ,n=2,4,6,8 and 10, are solved in a few seconds.     

A new stabilized finite volume method for the stationary Stokes equations

In this paper we develop and study a new stabilized finite volume method for the two-dimensional Stokes equations. This method is based on a local Gauss integration technique and the conforming elements of the lowest-equal order pair (i.e., the P ₁–P ₁ pair). After a relationship between this method and a stabilized finite element method is established, an error estimate of optimal order in the H ¹-norm for velocity and an estimate in the L ²-norm for pressure are obtained. An optimal error estimate in the L ²-norm for the velocity is derived under an additional assumption on the body force. This work is supported in part by the NSF of China 10701001 and by the US National Science Foundation grant DMS-0609995 and CMG Chair Funds in Reservoir Simulation.     

Weak solutions for a diffuse interface model for two-phase flows of incompressible fluids with different densities and nonlocal free energies

We consider a diffuse interface model for the flow of two viscous incompressible Newtonian fluids with different densities in a bounded domain in two and three space dimensions and prove existence of weak solutions for it. In contrast to earlier contributions, we study a model with a singular nonlocal free energy, which controls the H^α/2-norm of the volume fraction. We show existence of weak solutions for large times with the aid of an implicit time discretization.     

On the existence of square-integrable solutions for systems with weak nonlinearity

The paper considers the problem of structural stability of systems under disturbance of coefficients having small L ²(ℝ)-norm. We derive conditions which guarantee that for every solution of the perturbed system there exists a solution of the original system which is close to the former in L ²(ℝ)-norm.     

On Morrey’s estimate of the Sobolev norms of solutions of elliptic equations

We give a complete proof of Morrey’s estimate for the W ^1,p -norm of a solution of a second-order elliptic equation on a domain in terms of the L ₁-norm of this solution. The dependence of the constant in this estimate on the coefficients of the equation is investigated.