2,603 research outputs found
Space-time adaptive finite elements for nonlocal parabolic variational inequalities
This article considers the error analysis of finite element discretizations
and adaptive mesh refinement procedures for nonlocal dynamic contact and
friction, both in the domain and on the boundary. For a large class of
parabolic variational inequalities associated to the fractional Laplacian we
obtain a priori and a posteriori error estimates and study the resulting
space-time adaptive mesh-refinement procedures. Particular emphasis is placed
on mixed formulations, which include the contact forces as a Lagrange
multiplier. Corresponding results are presented for elliptic problems. Our
numerical experiments for -dimensional model problems confirm the
theoretical results: They indicate the efficiency of the a posteriori error
estimates and illustrate the convergence properties of space-time adaptive, as
well as uniform and graded discretizations.Comment: 47 pages, 20 figure
A multiscale method for heterogeneous bulk-surface coupling
In this paper, we construct and analyze a multiscale (finite element) method
for parabolic problems with heterogeneous dynamic boundary conditions. As
origin, we consider a reformulation of the system in order to decouple the
discretization of bulk and surface dynamics. This allows us to combine
multiscale methods on the boundary with standard Lagrangian schemes in the
interior. We prove convergence and quantify explicit rates for low-regularity
solutions, independent of the oscillatory behavior of the heterogeneities. As a
result, coarse discretization parameters, which do not resolve the fine scales,
can be considered. The theoretical findings are justified by a number of
numerical experiments including dynamic boundary conditions with random
diffusion coefficients
Discontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints
We consider time discretization methods for abstract parabolic problems with
inhomogeneous linear constraints. Prototype examples that fit into the general
framework are the heat equation with inhomogeneous (time dependent) Dirichlet
boundary conditions and the time dependent Stokes equation with an
inhomogeneous divergence constraint. Two common ways of treating such linear
constraints, namely explicit or implicit (via Lagrange multipliers) are
studied. These different treatments lead to different variational formulations
of the parabolic problem. For these formulations we introduce a modification of
the standard discontinuous Galerkin (DG) time discretization method in which an
appropriate projection is used in the discretization of the constraint. For
these discretizations (optimal) error bounds, including superconvergence
results, are derived. Discretization error bounds for the Lagrange multiplier
are presented. Results of experiments confirm the theoretically predicted
optimal convergence rates and show that without the modification the (standard)
DG method has sub-optimal convergence behavior.Comment: 35 page
Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods
In this paper, we discuss a general multiscale model reduction framework
based on multiscale finite element methods. We give a brief overview of related
multiscale methods. Due to page limitations, the overview focuses on a few
related methods and is not intended to be comprehensive. We present a general
adaptive multiscale model reduction framework, the Generalized Multiscale
Finite Element Method. Besides the method's basic outline, we discuss some
important ingredients needed for the method's success. We also discuss several
applications. The proposed method allows performing local model reduction in
the presence of high contrast and no scale separation
Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity
Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV).Peer ReviewedPostprint (author's final draft
Improving Newton's method performance by parametrization: the case of Richards equation
The nonlinear systems obtained by discretizing degenerate parabolic equations
may be hard to solve, especially with Newton's method. In this paper, we apply
to Richards equation a strategy that consists in defining a new primary unknown
for the continuous equation in order to stabilize Newton's method by
parametrizing the graph linking the pressure and the saturation. The resulting
form of Richards equation is then discretized thanks to a monotone Finite
Volume scheme. We prove the well-posedness of the numerical scheme. Then we
show under appropriate non-degeneracy conditions on the parametrization that
Newton\^as method converges locally and quadratically. Finally, we provide
numerical evidences of the efficiency of our approach
- …