An equilibrated a posteriori error estimator for an Interior Penalty Discontinuous Galerkin approximation of the p-Laplace problem

We are concerned with an Interior Penalty Discontinuous Galerkin (IPDG) approximation of the p-Laplace equation and an equilibrated a posteriori error estimator. The IPDG method can be derived from a discretization of the associated minimization problem involving appropriately defined reconstruction operators. The equilibrated a posteriori error estimator provides an upper bound for the discretization error in the broken W 1,p norm and relies on the construction of an equilibrated flux in terms of a numerical flux function associated with the mixed formulation of the IPDG approximation. The relationship with a residual-type a posteriori error estimator is established as well. Numerical results illustrate the performance of both estimators.

Error Constants for the Semi-Discrete Galerkin Approximation of the Linear Heat Equation

In this paper, we propose $$L^2(J;H^1_0(\Omega ))$$ L 2 ( J ; H 0 1 ( Ω ) ) and $$L^2(J;L^2(\Omega ))$$ L 2 ( J ; L 2 ( Ω ) ) norm error estimates that provide the explicit values of the error constants for the semi-discrete Galerkin approximation of the linear heat equation. The derivation of these error estimates shows the convergence of the approximation to the weak solution of the linear heat equation. Furthermore, explicit values of the error constants for these estimates play an important role in the computer-assisted existential proofs of solutions to semi-linear parabolic partial differential equations. In particular, the constants provided in this paper are better than the existing constants and, in a sense, the best possible.

Weak solutions to the time-fractional g-Navier–Stokes equations and optimal control

In this paper, we introduce the g-Navier–Stokes equations with time-fractional derivative of order α ∈ ( 0 , 1 ) {\alpha\in(0,1)} in domains of ℝ 2 {\mathbb{R}^{2}} . We then study the existence and uniqueness of weak solutions by means of the Galerkin approximation. Finally, an optimal control problem is considered and solved.

Instability of Vertical Throughflows in Bidisperse Porous Media

In this paper, the instability of a vertical fluid motion, or throughflow, is investigated in a horizontal bidisperse porous layer that is uniformly heated from below. By means of the order-1 Galerkin approximation method, the critical Darcy–Rayleigh number for the onset of steady instability is determined in closed form. The coincidence between the linear instability threshold and the global nonlinear stability threshold, in the energy norm, is shown.

Spectral Galerkin Approximation of Space Fractional Optimal Control Problem with Integral State Constraint

In this paper spectral Galerkin approximation of optimal control problem governed by fractional advection diffusion reaction equation with integral state constraint is investigated. First order optimal condition of the control problem is discussed. Weighted Jacobi polynomials are used to approximate the state and adjoint state. A priori error estimates for control, state, adjoint state and Lagrangian multiplier are derived. Numerical experiment is carried out to illustrate the theoretical findings.

Connecting a direct and a Galerkin approach to slow manifolds in infinite dimensions

In this paper, we study slow manifolds for infinite-dimensional evolution equations. We compare two approaches: an abstract evolution equation framework and a finite-dimensional spectral Galerkin approximation. We prove that the slow manifolds constructed within each approach are asymptotically close under suitable conditions. The proof is based upon Lyapunov-Perron methods and a comparison of the local graphs for the slow manifolds in scales of Banach spaces. In summary, our main result allows us to change between different characterizations of slow invariant manifolds, depending upon the technical challenges posed by particular fast-slow systems.

Non-symmetric isogeometric FEM-BEM couplings

We present a coupling of the Finite Element and the Boundary Element Method in an isogeometric framework to approximate either two-dimensional Laplace interface problems or boundary value problems consisting of two disjoint domains. We consider the Finite Element Method in the bounded domains to simulate possibly non-linear materials. The Boundary Element Method is applied in unbounded or thin domains where the material behavior is linear. The isogeometric framework allows to combine different design and analysis tools: first, we consider the same type of NURBS parameterizations for an exact geometry representation and second, we use the numerical analysis for the Galerkin approximation. Moreover, it facilitates to perform h- and p-refinements. For the sake of analysis, we consider the framework of strongly monotone and Lipschitz continuous operators to ensure well-posedness of the coupled system. Furthermore, we provide a priori error estimates. We additionally show an improved convergence behavior for the errors in functionals of the solution that may double the rate under certain assumptions. Numerical examples conclude the work which illustrate the theoretical results.

Non-Isothermal Creeping Flows in a Pipeline Network: Existence Results

This paper deals with a 3D mathematical model for the non-isothermal steady-state flow of an incompressible fluid with temperature-dependent viscosity in a pipeline network. Using the pressure and heat flux boundary conditions, as well as the conjugation conditions to satisfy the mass balance in interior junctions of the network, we propose the weak formulation of the nonlinear boundary value problem that arises in the framework of this model. The main result of our work is an existence theorem (in the class of weak solutions) for large data. The proof of this theorem is based on a combination of the Galerkin approximation scheme with one result from the field of topological degrees for odd mappings defined on symmetric domains.

A FEM-Green Approach for Magnetic Field Problems with Open Boundaries

A new finite element method/boundary element method (FEM/BEM) scheme is proposed for the solution of the 2D magnetic static and quasi-static problems with unbounded domains. The novelty is an original approach in the treatment of the outer region. The related domain integral is eliminated at the discrete level by using the finite element approximation of the fundamental solutions (Green’s functions) at every node of the related mesh. This “FEM-Green” approach replaces the standard boundary element method. It is simpler to implement because no integration on the boundary of the domain is required. Then, the method leads to a substantially reduced computational burden. Moreover, the coupling with finite elements is more natural since it is based on the same Galerkin approximation. Some examples with open boundary and nonlinear materials are presented and compared with the standard finite element method.