BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alexandre Ern

Alexandre Ern appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and Science

Catalog footprint

What is connected

17works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

Hybrid high-order and weak Galerkin methods for the biharmonic problem

We devise and analyze two hybrid high-order (HHO) methods for the numerical approximation of the biharmonic problem. The methods support polyhedral meshes, rely on the primal formulation of the problem, and deliver $O(h^{k+1})$ $H^2$-error estimates when using polynomials of order $k\ge0$ to approximate the normal derivative on the mesh (inter)faces. Both HHO methods hinge on a stabilization in the spirit of Lehrenfeld--Schöberl for second-order PDEs. The cell unknowns are polynomials of order $(k+2)$ that can be eliminated locally by means of static condensation. The face unknowns approximating the trace of the solution on the mesh (inter)faces are polynomials of order $(k+1)$ in the first HHO method which is valid in dimension two and uses an original stabilization involving the canonical hybrid finite element, and they are of order $(k+2)$ for the second HHO method which is valid in arbitrary dimension and uses only $L^2$-orthogonal projections in the stabilization. A comparative discussion with the weak Galerkin methods from the literature is provided, highlighting the close connections and the improvements proposed herein. Additionally, we show how the two HHO methods can be combined with a Nitsche-like boundary-penalty technique to weakly enforce the boundary conditions. An originality in the devised Nitsche's technique is to avoid any penalty parameter that must be large enough. Finally, numerical results showcase the efficiency of the proposed methods and indicate that the HHO methods can generally outperform discontinuous Galerkin methods and even be competitive with $C^0$-interior penalty methods on triangular meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local decay rates of best-approximation errors using vector-valued finite elements for fields with low regularity and integrable curl or divergence

We estimate best-approximation errors using vector-valued finite elements for fields with low regularity in the scale of fractional-order Sobolev spaces. By assuming additionally that the target field has a curl or divergence property, we establish upper bounds on these errors that can be localized to the mesh cells. These bounds are derived using the quasi-interpolation errors with or without boundary prescription derived in [A. Ern and J.-L. Guermond, ESAIM Math. Model. Numer. Anal., 51 (2017), pp.~1367--1385]. By using the face-to-cell lifting operators analyzed in [A. Ern and J.-L. Guermond, Found. Comput. Math., (2021)], and exploiting the additional assumption made on the curl or the divergence of the target field, a localized upper bound on the quasi-interpolation error is derived. As an illustration, we show how to apply these results to the error analysis of the curl-curl problem associated with Maxwell's equations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quasi-static crack propagation with a Griffith criterion using a variational discrete element method

A variational discrete element method is applied to simulate quasi-static crack propagation. Cracks are considered to propagate between the mesh cells through the mesh facets. The elastic behaviour is parametrized by the continuous mechanical parameters (Young modulus and Poisson ratio). A discrete energetic cracking criterion coupled to a discrete kinking criterion guide the cracking process. Two-dimensional numerical examples are presented to illustrate the robustness and versatility of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A variational discrete element method for quasi-static and dynamic elasto-plasticity

We propose a new discrete element method supporting general polyhedral meshes. The method can be understood as a lowest-order discontinuous Galerkin method parametrized by the continuous mechanical parameters (Young's modulus and Poisson's ratio). We consider quasi-static and dynamic elasto-plasticity, and in the latter situation, a pseudo-energy conserving time-integration method is employed. The computational cost of the time-stepping method is moderate since it is explicit and used with a naturally diagonal mass matrix. Numerical examples are presented to illustrate the robustness and versatility of the method for quasi-static and dynamic elasto-plastic evolutions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An explicit pseudo-energy conserving time-integration scheme for Hamiltonian dynamics

We propose a new explicit pseudo-energy and momentum conserving scheme for the time integration of Hamiltonian systems. The scheme, which is formally second-order accurate, is based on two key ideas: the integration during the time-steps of forces between free-flight particles and the use of momentum jumps at the discrete time nodes leading to a two-step formulation for the acceleration. The pseudo-energy conservation is established under exact force integration, whereas it is valid to second-order accuracy in the presence of quadrature errors. Moreover, we devise an asynchronous version of the scheme that can be used in the framework of slow-fast time-stepping strategies. The scheme is validated against classical benchmarks and on nonlinear or inhomogeneous wave propagation problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Equivalence of local-and global-best approximations, a simple stable local commuting projector, and optimal $hp$ approximation estimates in $H(\mathrm{div})$

Given an arbitrary function in H(div), we show that the error attained by the global-best approximation by H(div)-conforming piecewise polynomial Raviart-Thomas-Nédélec elements under additional constraints on the divergence and normal flux on the boundary, is, up to a generic constant, equivalent to the sum of independent local-best approximation errors over individual mesh elements, without constraints on the divergence or normal fluxes. The generic constant only depends on the shape-regularity of the underlying simplicial mesh, the space dimension, and the polynomial degree of the approximations. The analysis also gives rise to a stable, local, commuting projector in H(div), delivering an approximation error that is equivalent to the local-best approximation. We next present a variant of the equivalence result, where robustness of the constant with respect to the polynomial degree is attained for unbalanced approximations. These two results together further enable us to derive rates of convergence of global-best approximations that are fully optimal in both the mesh size h and the polynomial degree p, for vector fields that only feature elementwise the minimal necessary Sobolev regularity. We finally show how to apply our findings to derive optimal a priori hp-error estimates for mixed and least-squares finite element methods applied to a model diffusion problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Polynomial-degree-robust H(curl)-stability of discrete minimization in a tetrahedron

We prove that the minimizer in the Nédélec polynomial space of some degree p > 0 of a discrete minimization problem performs as well as the continuous minimizer in H(curl), up to a constant that is independent of the polynomial degree p. The minimization problems are posed for fields defined on a single non-degenerate tetrahedron in R^3 with polynomial constraints enforced on the curl of the field and its tangential trace on some faces of the tetrahedron. This result builds upon [L. Demkowicz, J. Gopalakrishnan, J. Schöberl SIAM J. Numer. Anal. 47 (2009), 3293--3324] and [M. Costabel, A. McIntosh, Math. Z. 265 (2010), 297--320] and is a fundamental ingredient to build polynomial-degree-robust a posteriori error estimators when approximating the Maxwell equations in several regimes leading to a curl-curl problem.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An adaptive stabilized conforming finite element method via residual minimization on dual discontinuous Galerkin norms

We design and analyze a new adaptive stabilized finite element method. We construct a discrete approximation of the solution in a continuous trial space by minimizing the residual measured in a dual norm of a discontinuous test space that has inf-sup stability. We formulate this residual minimization as a stable saddle-point problem which delivers a stabilized discrete solution and a residual representation that drives the adaptive mesh refinement. Numerical results on an advection-reaction model problem show competitive error reduction rates when compared to discontinuous Galerkin methods on uniformly refined meshes and smooth solutions. Moreover, the technique leads to optimal decay rates for adaptive mesh refinement and solutions having sharp layers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite element quasi-interpolation and best approximation

This paper introduces a quasi-interpolation operator for scalar- and vector-valued finite element spaces constructed on affine, shape-regular meshes with some continuity across mesh interfaces.This operator gives optimal estimates of the best approximation error in any $L^p$-norm assuming regularity in the fractional Sobolev spaces $W^{r,p}$, where $p\in [1,\infty]$ and the smoothness index $r$ can be arbitrarily close to zero. The operator is stable in $L^1$, leaves the corresponding finite element space point-wise invariant whether homogeneous boundary conditions are imposed or not. The theory is illustrated on $H^1$-, $\mathbf{H}(\text{curl})$- and $\mathbf{H}(\text{div})$-conforming spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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é 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A nonintrusive method to approximate linear systems with nonlinear parameter dependence

We consider a family of linear systems $A_μα=C$ with system matrix $A_μ$ depending on a parameter $μ$ and for simplicity parameter-independent right-hand side $C$. These linear systems typically result from the finite-dimensional approximation of a parameter-dependent boundary-value problem. We derive a procedure based on the Empirical Interpolation Method to obtain a separated representation of the system matrix in the form $A_μ\approx\sum_{m}β_m(μ)A_{μ_m}$ for some selected values of the parameter. Such a separated representation is in particular useful in the Reduced Basis Method. The procedure is called nonintrusive since it only requires to access the matrices $A_{μ_m}$. As such, it offers a crucial advantage over existing approaches that instead derive separated representations requiring to enter the code at the level of assembly. Numerical examples illustrate the performance of our new procedure on a simple one-dimensional boundary-value problem and on three-dimensional acoustic scattering problems solved by a boundary element method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coupled BEM-FEM for the convected Helmholtz equation with non-uniform flow in a bounded domain

We consider the convected Helmholtz equation modeling linear acoustic propagation at a fixed frequency in a subsonic flow around a scattering object. The flow is supposed to be uniform in the exterior domain far from the object, and potential in the interior domain close to the object. Our key idea is the reformulation of the original problem using the Prandtl--Glauert transformation on the whole flow domain, yielding (i) the classical Helmholtz equation in the exterior domain and (ii) an anisotropic diffusive PDE with skew-symmetric first-order perturbation in the interior domain such that its transmission condition at the coupling boundary naturally fits the Neumann condition from the classical Helmholtz equation. Then, efficient off-the-shelf tools can be used to perform the BEM-FEM coupling, leading to two novel variational formulations for the convected Helmholtz equation. The first formulation involves one surface unknown and can be affected by resonant frequencies, while the second formulation avoids resonant frequencies and involves two surface unknowns. Numerical simulations are presented to compare the two formulations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Mass conservative BDF-discontinuous Galerkin/explicit finite volume schemes for coupling subsurface and overland flows

Robust and accurate schemes are designed to simulate the coupling between subsurface and overland flows. The coupling conditions at the interface enforce the continuity of both the normal flux and the pressure. Richards' equation governing the subsurface flow is discretized using a Backward Differentiation Formula and a symmetric interior penalty Discontinuous Galerkin method. The kinematic wave equation governing the overland flow is discretized using a Godunov scheme. Both schemes individually are mass conservative and can be used within single-step or multi-step coupling algorithms that ensure overall mass conservation owing to a specific design of the interface fluxes in the multi-step case. Numerical results are presented to illustrate the performances of the proposed algorithms.

0 citations0 reviewsNo code linkedOpen access