In this paper we construct a new H(div)-conforming projection-based p-interpolation operator that assumes only Hr(K) -1/2(div, K)-regularity (r > 0) on the reference element (either triangle or square) K. We show that this operator is stable with...
We consider a coupled PDE model of various fluid-structure interactions seen in nature. It has recently been shown by the authors [Contemp. Math. 440, 2007] that this model admits of an explicit semigroup generator representation 𝓐:D(𝓐)⊂ H → H, where H is the associated space of fluid-structure initial data. However, the argument for the maximality criterion was indirect, and did not provide for an explicit solution Φ ∈ D(𝓐) of the equation (λI-𝓐)Φ =F for given F ∈ H and λ > 0. The present...
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 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.
We introduce a new stable MINI-element pair for incompressible Stokes equations on quadrilateral meshes, which uses the smallest number of bubbles for the velocity. The pressure is discretized with the P1-midpoint-edge-continuous elements and each component of the velocity field is done with the standard Q1-conforming elements enriched by one bubble a quadrilateral. The superconvergence in the pressure of the proposed pair is analyzed on uniform rectangular meshes, and tested numerically on uniform...
This work is concerned with the introduction of a new numerical scheme based on the discontinuous Galerkin (DG) method. We propose to follow the methodology of higher order finite volume schemes and introduce a reconstruction operator into the DG scheme. This operator constructs higher order piecewise polynomial reconstructions from the lower order DG scheme. Such a procedure was proposed already in [2] based on heuristic arguments, however we provide a rigorous derivation, which justifies the increased...
This study is mainly dedicated to the development and analysis of non-overlapping domain decomposition methods for solving continuous-pressure finite element formulations of the Stokes problem. These methods have the following special features. By keeping the equations and unknowns unchanged at the cross points, that is, points shared by more than two subdomains, one can interpret them as iterative solvers of the actual discrete problem directly issued from the finite element scheme. In this way,...
This study is mainly dedicated to the development and analysis of non-overlapping domain decomposition methods for solving continuous-pressure finite element formulations of the Stokes problem. These methods have the following special features. By keeping the equations and unknowns unchanged at the cross points, that is, points shared by more than two subdomains, one can interpret them as iterative solvers of the actual discrete problem directly issued from the finite element scheme. In this way,...
The proof of the Friedrichs' inequality on the class of finite dimensional spaces used in the finite element method is given. In particular, the approximate spaces generated by simplicial isoparametric elements are considered.
When solving parabolic problems by the so-called Rothe method (see K. Rektorys, Czech. Math. J. 21 (96), 1971, 318-330 and other authors), some difficulties of theoretical nature are encountered in the case of nonhomogeneous initial and boundary conditions. As a rule, these difficulties lead to rather unnatural additional conditions imposed on the corresponding bilinear form and the initial and boundary functions. In the present paper, it is shown how to remove such additional assumptions in the...
For domains which are star-shaped w.r.t. at least one point, we give new bounds on the constants in Jackson-inequalities in Sobolev spaces. For convex domains, these bounds do not depend on the eccentricity. For non-convex domains with a re-entrant corner, the bounds are uniform w.r.t. the exterior angle. The main tool is a new projection operator onto the space of polynomials.