3,261 research outputs found
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 stable for any real number
, and commute with the differential operators ,
, and . 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
-, - and -conforming
finite element spaces, with and without homogeneous boundary conditions. These
operators commute with the differential operators , ,
and , are -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
A nonintrusive Reduced Basis Method applied to aeroacoustic simulations
The Reduced Basis Method can be exploited in an efficient way only if the
so-called affine dependence assumption on the operator and right-hand side of
the considered problem with respect to the parameters is satisfied. When it is
not, the Empirical Interpolation Method is usually used to recover this
assumption approximately. In both cases, the Reduced Basis Method requires to
access and modify the assembly routines of the corresponding computational
code, leading to an intrusive procedure. In this work, we derive variants of
the EIM algorithm and explain how they can be used to turn the Reduced Basis
Method into a nonintrusive procedure. We present examples of aeroacoustic
problems solved by integral equations and show how our algorithms can benefit
from the linear algebra tools available in the considered code.Comment: 28 pages, 7 figure
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
- …