A remark on least-squares mixed element methods for reaction–diffusion problems

In this paper, we propose a least-squares mixed element procedure for a reaction–diffusion problem based on the first-order system. By selecting the least-squares functional properly, the resulting procedure can be split into two independent symmetric positive definite schemes, one of which is for the unknown variable and the other of which is for the unknown flux variable, which lead to the optimal order H¹(Ω)$H^1(\Omega )$ and L²(Ω)$L^2(\Omega )$ norm error estimates for the primal unknown and optimal H(div;Ω)$H(\mathrm{div};\Omega )$ norm error estimate for the unknown flux. Finally, we give some numerical examples.     

A remark on least-squares Galerkin procedures for pseudohyperbolic equations

In this paper, we introduce two split least-squares Galerkin finite element procedures for pseudohyperbolic equations arising in the modelling of nerve conduction process. By selecting the least-squares functional properly, the procedures can be split into two sub-procedures, one of which is for the primitive unknown variable and the other is for the flux. The convergence analysis shows that both the two methods yield the approximate solutions with optimal accuracy in L²(Ω)$L^2\left(\Omega \right)$ norm for u$u$ and _ut$u_t$ and (L²(Ω))²${\left(L^2\left(\Omega \right)\right)}^2$ norm for the flux σ$\sigma$. Moreover, the two methods get approximate solutions with first-order and second-order accuracy in time increment, respectively. A numerical example is given to show the efficiency of the introduced schemes.     

A priori error analysis of a discontinuous Galerkin approximation for a kind of compressible miscible displacement problems

A kind of compressible miscible displacement problems which include molecular diffusion and dispersion in porous media are investigated.A symmetric interior penalty discontinuous Galerkin (SIPG) method is applied to the coupled system of flow and transport.Using the induction hypotheses instead of the cut-off operator and the interpolation projection properties,a priori hp error estimates are presented.The error bounds in L2(H1) norm for concentration and in L∞(L2) norm for velocity are optimal in h and suboptimal in p with a loss of power 1/2.     

Superconvergence of Legendre projection methods for the eigenvalue problem of a compact integral operator

In this paper, we consider the Galerkin and collocation methods for the eigenvalue problem of a compact integral operator with a smooth kernel using the Legendre polynomials of degree ≤n. We prove that the error bounds for eigenvalues are of the order O(n^−2r) and the gap between the spectral subspaces are of the orders O(n^−r) in L₂-norm and O(n^1/2−r) in the infinity norm, where r denotes the smoothness of the kernel. By iterating the eigenvectors we show that the iterated eigenvectors converge with the orders of convergence O(n^−2r) in both L₂-norm and infinity norm. We illustrate our results with numerical examples.     

Long-time behavior of reaction-diffusion equations with dynamical boundary condition

In this paper, we study the long-time behavior of the reaction-diffusion equation with dynamical boundary condition, where the nonlinear terms f and g satisfy the polynomial growth condition of arbitrary order. Some asymptotic regularity of the solution has been proved. As an application of the asymptotic regularity results, we can not only obtain the existence of a global attractor A in (H¹(Ω)∩L^p(Ω))×L^q(Γ) immediately, but also can show further that A attracts every L²(Ω)×L²(Γ)-bounded subset with (H¹(Ω)∩L^p+δ(Ω))×L^q+κ(Γ)-norm for any δ,κ∈[0,∞).     

Asymptotic behaviour of nonlocal reaction-diffusion equations

The existence of a global attractor in L²(Ω) is established for a reaction-diffusion equation on a bounded domain Ω in R^d with Dirichlet boundary conditions, where the reaction term contains an operator F:L²(Ω)→L²(Ω) which is nonlocal and possibly nonlinear. Existence of weak solutions is established, but uniqueness is not required. Compactness of the multivalued flow is obtained via estimates obtained from limits of Galerkin approximations. In contrast with the usual situation, these limits apply for all and not just for almost all time instants.     

A coupling method based on new MFE and FE for fourth-order parabolic equation

In this article, a coupling method of new mixed finite element (MFE) and finite element (FE) is proposed and analyzed for fourth-order parabolic partial differential equation. First, the fourth-order parabolic equation is split into the coupled system of second-order equations. Then, an equation is solved by finite element method, the other equation is approximated by the new mixed finite element method, whose flux belongs to the square integrable space replacing the classical H(div;Ω) space. The stability for fully discrete scheme is derived, and both semi-discrete and fully discrete error estimates are obtained. Moreover, the optimal a priori error estimates in L ² and H ¹-norm for both the scalar unknown u and the diffusion term γ and a priori error estimate in (L ²)²-norm for its flux σ are derived. Finally, some numerical results are provided to validate our theoretical analysis.     

ОБ ЁкВИВАлЕНтНых НОР МИРОВкАх пОДпРОстРА НстВL 2

Denote byH(Ω) the subspace of functionsp fromL ₂(Ω) for which \(\mathop \smallint \limits_\Omega p dx = 0\) , and define the norm inH(Ω) as follows It is proved that this norm is equivalent to the norm ofL ₂(Ω) assuming that the boundary ofΩ is piecewise Lipschitzian. As a consequence it is derived that for linear elasticity problems the stresses are uniformly bounded inL ₂ with respect to the Poisson constant.     

Systems of p-Laplacian differential inclusions with large diffusion

In this paper we consider coupled systems of p-Laplacian differential inclusions and we prove, under suitable conditions, that a homogenization process occurs when diffusion parameters become arbitrarily large. In fact we obtain that the attractors are continuous at infinity on L²(Ω)×L²(Ω) topology, with respect to the diffusion coefficients, and the limit set is the attractor of an ordinary differential problem.     

New error estimates of nonconforming mixed finite element methods for the Stokes problem

In this paper, we revisit the classical error estimates of nonconforming Crouzeix–Raviart type finite elements for the Stokes equations. By introducing some quasi‐interpolation operators and using the special properties of these nonconforming elements, it is proved that their consistency errors can be bounded by their approximation errors together with a high‐order term, especially which can be of arbitrary order provided that f in the right‐hand side is piecewise smooth enough. Furthermore, we show an interesting result that both in the energy norm and L² norm the consistency errors are dominated by the approximation errors of their finite element spaces. As byproducts, we derive the error estimates in both energy and L² norms under the regularity assumption ( u ,p) ∈ H ^{1 + s}(Ω) × H^s(Ω) with any s ∈ (0,1], which fills the gap in the a priori error estimate of these nonconforming elements with low regularity . Furthermore, a robust convergence is proved with minimal regularity assumption s = 0. These results seem to be missing in the literature. Numerical tests are provided, confirming the analysis, especially the new results on the L² convergence. Copyright © 2013 John Wiley & Sons, Ltd.     

Elliptic boundary value problems with λ-dependent boundary conditions

In this paper second order elliptic boundary value problems on bounded domains Ω⊂Rn with boundary conditions on ∂Ω depending nonlinearly on the spectral parameter are investigated in an operator theoretic framework. For a general class of locally meromorphic functions in the boundary condition a solution operator of the boundary value problem is constructed with the help of a linearization procedure. In the special case of rational Nevanlinna or Riesz-Herglotz functions on the boundary the solution operator is obtained in an explicit form in the product Hilbert space L²(Ω)⊕(L²m(∂Ω)), which is a natural generalization of known results on λ-linear elliptic boundary value problems and λ-rational boundary value problems for ordinary second order differential equations.     

The convergence of non-Newtonian fluids to Navier-Stokes equations

In this paper, the convergence of solutions for incompressible dipolar viscous non-Newtonian fluids is investigated. We obtain the conclusion that the solutions of non-Newtonian fluids converge to the solutions of Navier-Stokes equations in the sense of L²-norm (resp. H¹-norm), as the viscosities tend to zero and the initial data belong to H¹(Ω) (resp. H²(Ω)). Moreover, we obtain L^∞-norm convergence of solutions if the initial data belong to H²(Ω).     

Local mesh refinement for the discretization of Neumann boundary control problems on polyhedra

This paper deals with optimal control problems constrained by linear elliptic partial differential equations. The case where the right‐hand side of the Neumann boundary is controlled, is studied. The variational discretization concept for these problems is applied, and discretization error estimates are derived. On polyhedral domains, one has to deal with edge and corner singularities, which reduce the convergence rate of the discrete solutions, that is, one cannot expect convergence order two for linear finite elements on quasi‐uniform meshes in general. As a remedy, a local mesh refinement strategy is presented, and a priori bounds for the refinement parameters are derived such that convergence with optimal rate is guaranteed. As a by‐product, finite element error estimates in the H¹(Ω)‐norm, L²(Ω)‐norm and L²(Γ)‐norm for the boundary value problem are obtained, where the latter one turned out to be the main challenge. Copyright © 2015 John Wiley & Sons, Ltd.     

Optimal Control of the Obstacle for a Parabolic Variational Inequality

An optimal control problem for a parabolic obstacle variational inequality is considered. The obstacle in L²(0, T; H²(Ω) ∩ H¹₀(Ω)) with ψ_t ∈ L²(Q) is taken as the control, and the solution to the obstacle problem is taken as the state. The goal is to find the optimal control so that the state is close to the desired profile while the norm of the obstacle is not too large. Existence and necessary conditions for the optimal control are established.     

On the function S(T) in the theory of the Riemann zeta-function

The function S(T) is the error term in the formula for the number of zeros of the Riemann zeta-function above the real axis and up to height T in the complex plane. We assume the Riemann hypothesis, and examine how well S(T) can be approximated by a Dirichlet polynomial in the L² norm.     

Boundary trace embedding theorems for variable exponent Sobolev spaces

We study boundary trace embedding theorems for variable exponent Sobolev space W1,p(⋅)(Ω). Let Ω be an open (bounded or unbounded) domain in RN satisfying strong local Lipschitz condition. Under the hypotheses that p∈L^∞(Ω), 1?infp(x)?supp(x)<N, |∇p|∈Lγ(⋅)(Ω), where γ∈L^∞(Ω) and infγ(x)>N, we prove that there is a continuous boundary trace embedding W1,p(⋅)(Ω)→Lq(⋅)(∂Ω) provided q(⋅), a measurable function on ∂Ω, satisfies condition for x∈∂Ω.     

A uniformly convergent continuous-discontinuous Galerkin method for singularly perturbed problems of convection-diffusion type

In this paper, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving singularly perturbed convection-diffusion problems. On Shishkin mesh with linear elements, a rate O(N^-1lnN) in an associated norm is established, where N is the number of elements. Numerical experiments complement the theoretical results. Moreover, a rate O(N^-2ln²N) in a discrete L^∞ norm, and O(N^-2) in L² norm, are observed numerically on the Shishkin mesh.     

The polynomial Fourier transform with minimum mean square error for noisy data

In 2006, Naoki Saito proposed a Polyharmonic Local Fourier Transform (PHLFT) to decompose a signal f∈L²(Ω) into the sum of a polyharmonic componentu and a residualv, where Ω is a bounded and open domain in R^d. The solution presented in PHLFT in general does not have an error with minimal energy. In resolving this issue, we propose the least squares approximant to a given signal in L²([−1,1]) using the combination of a set of algebraic polynomials and a set of trigonometric polynomials. The maximum degree of the algebraic polynomials is chosen to be small and fixed. We show in this paper that the least squares approximant converges uniformly for a Hölder continuous function. Therefore Gibbs phenomenon will not occur around the boundary for such a function. We also show that the PHLFT converges uniformly and is a near least squares approximation in the sense that it is arbitrarily close to the least squares approximant in L² norm as the dimension of the approximation space increases. Our experiments show that the proposed method is robust in approximating a highly oscillating signal. Even when the signal is corrupted by noise, the method is still robust. The experiments also reveal that an optimum degree of the trigonometric polynomial is needed in order to attain the minimal l² error of the approximation when there is noise present in the data set. This optimum degree is shown to be determined by the intrinsic frequency of the signal. We also discuss the energy compaction of the solution vector and give an explanation to it.     

On non-effective weights in Orlicz spaces

Given a weight w in Ω ⊂ ∝^N, |Ω| < ∞ and a Young function φ, we consider the weighted modular ∫_Ω ω(f(x))w(x)dx and the resulting weighted Orlicz space L_ω(w). For a Young function Ω ∉ Δ₂(∞) we present a necessary and sufficient conditions in order that L_ω(w) = L_ω(X_Ω) up to the equivalence of norms. We find a necessary and sufficient condition for ω in order that there exists an unbounded weight w such that the above equality of spaces holds. By way of applications we simplify criteria from [5] for continuity of the composition operator from L_ω into itself when ω Δ₂(∞) and obtain necessary and sufficient condition in order that the composition operator maps L_ω. continuously onto L_ω.     

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.