No abstract

In this article, we propose a novel discontinuous Galerkin method for convectiondiffusion-reaction problems, characterized by three main properties. The first is that the method is hybridizable; this renders it efficiently implementable and competitive with the main existing methods for these problems. The second is that, when the method uses polynomial approximations of the same degree for both the total flux and the scalar variable, optimal convergence properties are obtained for both variables; this is in sharp contrast with all other discontinuous methods for this problem. The third is that the method exhibits superconvergence properties of the approximation to the scalar variable; this allows us to postprocess the approximation in an element-by-element fashion to obtain another approximation to the scalar variable which converges faster than the original one. In this paper, we focus on the efficient implementation of the method and on the validation of its computational performance. With this aim, extensive numerical tests are devoted to explore the convergence properties of the novel scheme, to compare it with other methods in the diffusiondominated regime, and to display its stability and accuracy in the convection-dominated case.

No abstract

We consider the initial and boundary value problem for a system of partial differential equations describing the motion of a fluid-solid mixture under the assumption of full saturation. The ability of the fluid phase to flow within the solid skeleton is described by the permeability tensor, which is assumed here to be a multiple of the identity and to depend nonlinearly on the volumetric solid strain. In particular, we study the problem of existence of weak solutions in bounded domains, accounting for non-zero volumetric and boundary forcing terms. We investigate the influence of viscoelasticity on the solution functional setting and on the regularity requirements for the forcing terms. The theoretical analysis shows that different time regularity requirements are needed for the volumetric source of linear momentum and the boundary source of traction depending on whether or not viscoelasticity is present. The theoretical results are further investigated via numerical simulations based on a novel dual mixed hybridized finite element discretization. When the data are sufficiently regular, the simulations show that the solutions satisfy the energy estimates predicted by the theoretical analysis. Interestingly, the simulations also show that, in the purely elastic case, the Darcy velocity and the related fluid energy might become unbounded if indeed the data do not enjoy the time regularity required by the theory.

No abstract