Implicit-explicit Runge–Kutta schemes and finite elements with symmetric stabilization for advection-diffusion equations

We analyze a two-stage implicit-explicit Runge–Kutta scheme for time discretization of advection-diffusion equations. Space discretization uses continuous, piecewise affine finite elements with interelement gradient jump penalty; discontinuous Galerkin methods can be considered as well. The advective and stabilization operators are treated explicitly, whereas the diffusion operator is treated implicitly. Our analysis hinges on L 2 -energy estimates on discrete functions in physical space. Our main results are stability and quasi-optimal error estimates for smooth solutions under a standard hyperbolic CFL restriction on the time step, both in the advection-dominated and in the diffusion-dominated regimes. The theory is illustrated by numerical examples

The motion of an axisymmetric body falling in a tube at moderate Reynolds numbers

This study concerns the rectilinear and periodic paths of an axisymmetric solid body (short-length cylinder and disk of diameter d and thickness h) falling in a vertical tube of diameter D. We investigated experimentally the influence of the confinement ratio (S=d/D<0.8) on the motion of the body, for different aspect ratios (χ=d/h=3, 6 and 10), Reynolds numbers (80<Re<320) and a density ratio between the fluid and the body close to unity. For a given body, the Reynolds number based on its mean vertical velocity is observed to decrease when S increases. The critical Reynolds number for the onset of the periodic motion decreases with S in the case of thin bodies (χ=10), whereas it appears unaffected by S for thicker bodies (χ=3 and 6). The characteristics of the periodic motion are also strongly modified by the confinement ratio. A thick body (χ=3) tends to go back to a rectilinear path when S increases, while a thin body (χ=10) displays oscillations of growing amplitude with S until it touches the tube (at about S=0.5). For a given aspect ratio, however, the amplitudes of the oscillations follow a unique curve for all S, which depends only on the relative distance of the Reynolds number to the threshold of path instability. In parallel, numerical simulations of the wake of a body held fixed in a uniform confined flow were carried out. The simulations allowed us to determine in this configuration the effect of the confinement ratio on the thresholds for wake instability (loss of axial symmetry at Rec₁ and loss of stationarity at Rec₂) and on the maximal velocity Vw in the recirculating region of the stationary axisymmetric wake. The evolution with χ and S of Vw at Rec₁ was used to define a Reynolds number Re*. Remarkably, for a freely moving body, Re* remains almost constant when S varies, regardless of the nature of the path

Analysis of Compatible Discrete Operator Schemes for the Stokes Equations on Polyhedral Meshes

Compatible Discrete Operator schemes preserve basic properties of the continuous model at the discrete level. They combine discrete differential operators that discretize exactly topological laws and discrete Hodge operators that approximate constitutive relations. We devise and analyze two families of such schemes for the Stokes equations in curl formulation, with the pressure degrees of freedom located at either mesh vertices or cells. The schemes ensure local mass and momentum conservation. We prove discrete stability by establishing novel discrete Poincar\'e inequalities. Using commutators related to the consistency error, we derive error estimates with first-order convergence rates for smooth solutions. We analyze two strategies for discretizing the external load, so as to deliver tight error estimates when the external load has a large irrotational or divergence-free part. Finally, numerical results are presented on three-dimensional polyhedral meshes

Mollification in strongly Lipschitz domains with application to continuous and discrete De Rham complex

We construct mollification operators in strongly Lipschitz domains that do not invoke non-trivial extensions, are $L^p$ stable for any real number $p\in[1,\infty]$, and commute with the differential operators $\nabla$, $\nabla{\times}$, and $\nabla{\cdot}$. We also construct mollification operators satisfying boundary conditions and use them to characterize the kernel of traces related to the tangential and normal trace of vector fields. We use the mollification operators to build projection operators onto general $H^1$-, $\mathbf{H}(\text{curl})$- and $\mathbf{H}(\text{div})$-conforming finite element spaces, with and without homogeneous boundary conditions. These operators commute with the differential operators $\nabla$, $\nabla{\times}$, and $\nabla{\cdot}$, are $L^p$-stable, and have optimal approximation properties on smooth functions

A converse to Fortin's Lemma in Banach spaces

The converse of Fortin's Lemma in Banach spaces is established in this Note

Equilibrated tractions for the Hybrid High-Order method

We show how to recover equilibrated face tractions for the hybrid high-order method for linear elasticity recently introduced in [D. A. Di Pietro and A. Ern, A hybrid high-order locking-free method for linear elasticity on general meshes, Comput. Meth. Appl. Mech. Engrg., 2015, 283:1-21], and prove that these tractions are optimally convergent

Variants of the Empirical Interpolation Method: symmetric formulation, choice of norms and rectangular extension

The Empirical Interpolation Method (EIM) is a greedy procedure that constructs approximate representations of two-variable functions in separated form. In its classical presentation, the two variables play a non-symmetric role. In this work, we give an equivalent definition of the EIM approximation, in which the two variables play symmetric roles. Then, we give a proof for the existence of this approximation, and extend it up to the convergence of the EIM, and for any norm chosen to compute the error in the greedy step. Finally, we introduce a way to compute a separated representation in the case where the number of selected values is different for each variable. In the case of a physical field measured by sensors, this is useful to discard a broken sensor while keeping the information provided by the associated selected field.Comment: 7 page

Accurate and efficient evaluation of the a posteriori error estimator in the reduced basis method

The reduced basis method is a model reduction technique yielding substantial savings of computational time when a solution to a parametrized equation has to be computed for many values of the parameter. Certification of the approximation is possible by means of an a posteriori error bound. Under appropriate assumptions, this error bound is computed with an algorithm of complexity independent of the size of the full problem. In practice, the evaluation of the error bound can become very sensitive to round-off errors. We propose herein an explanation of this fact. A first remedy has been proposed in [F. Casenave, Accurate \textit{a posteriori} error evaluation in the reduced basis method. \textit{C. R. Math. Acad. Sci. Paris} \textbf{350} (2012) 539--542.]. Herein, we improve this remedy by proposing a new approximation of the error bound using the Empirical Interpolation Method (EIM). This method achieves higher levels of accuracy and requires potentially less precomputations than the usual formula. A version of the EIM stabilized with respect to round-off errors is also derived. The method is illustrated on a simple one-dimensional diffusion problem and a three-dimensional acoustic scattering problem solved by a boundary element method.Comment: 26 pages, 10 figures. ESAIM: Mathematical Modelling and Numerical Analysis, 201