2,241 research outputs found
A numerical method to solve the Stokes problem with a punctual force in source term
The aim of this note is to present a numerical method to solve the Stokes
problem in a bounded domain with a Dirac source term, which preserves
optimality for any approximation order by the finite-element method. It is
based on the knowledge of a fundamental solution to the associated operator
over the whole space. This method is motivated by the modeling of the movement
of active thin structures in a viscous fluid.Comment: Comptes Rendus M{\'e}canique, Elsevier, 2015,
http://www.sciencedirect.com/science/article/pii/S163107211400181
Splitting method for elliptic equations with line sources
In this paper, we study the mathematical structure and numerical
approximation of elliptic problems posed in a (3D) domain when the
right-hand side is a (1D) line source . The analysis and approximation
of such problems is known to be non-standard as the line source causes the
solution to be singular. Our main result is a splitting theorem for the
solution; we show that the solution admits a split into an explicit, low
regularity term capturing the singularity, and a high-regularity correction
term being the solution of a suitable elliptic equation. The splitting
theorem states the mathematical structure of the solution; in particular, we
find that the solution has anisotropic regularity. More precisely, the solution
fails to belong to in the neighbourhood of , but exhibits
piecewise -regularity parallel to . The splitting theorem can
further be used to formulate a numerical method in which the solution is
approximated via its correction function . This approach has several
benefits. Firstly, it recasts the problem as a 3D elliptic problem with a 3D
right-hand side belonging to , a problem for which the discretizations and
solvers are readily available. Secondly, it makes the numerical approximation
independent of the discretization of ; thirdly, it improves the
approximation properties of the numerical method. We consider here the Galerkin
finite element method, and show that the singularity subtraction then recovers
optimal convergence rates on uniform meshes, i.e., without needing to refine
the mesh around each line segment. The numerical method presented in this paper
is therefore well-suited for applications involving a large number of line
segments. We illustrate this by treating a dataset (consisting of
line segments) describing the vascular system of the brain
- …