18,119 research outputs found
Unique continuation for the Helmholtz equation using stabilized finite element methods
In this work we consider the computational approximation of a unique
continuation problem for the Helmholtz equation using a stabilized finite
element method. First conditional stability estimates are derived for which,
under a convexity assumption on the geometry, the constants grow at most
linearly in the wave number. Then these estimates are used to obtain error
bounds for the finite element method that are explicit with respect to the wave
number. Some numerical illustrations are given.Comment: corrected typos; included suggestions from reviewer
Stabilization of Unstable Procedures: The Recursive Projection Method
Fixed-point iterative procedures for solving nonlinear parameter dependent problems can converge for some interval of parameter values and diverge as the parameter changes. The Recursive Projection Method (RPM), which stabilizes such procedures by computing a projection onto the unstable subspace is presented. On this subspace a Newton or special Newton iteration is performed, and the fixed-point iteration is used on the complement. As continuation in the parameter proceeds, the projection is efficiently updated, possibly increasing or decreasing the dimension of the unstable subspace. The method is extremely effective when the dimension of the unstable subspace is small compared to the dimension of the system. Convergence proofs are given and pseudo-arclength continuation on the unstable subspace is introduced to allow continuation past folds. Examples are presented for an important application of the RPM in which a “black-box” time integration scheme is stabilized, enabling it to compute unstable steady states. The RPM can also be used to accelerate iterative procedures when slow convergence is due to a few slowly decaying modes
A stabilized finite element method for inverse problems subject to the convection-diffusion equation. I: diffusion-dominated regime
The numerical approximation of an inverse problem subject to the
convection--diffusion equation when diffusion dominates is studied. We derive
Carleman estimates that are on a form suitable for use in numerical analysis
and with explicit dependence on the P\'eclet number. A stabilized finite
element method is then proposed and analysed. An upper bound on the condition
number is first derived. Combining the stability estimates on the continuous
problem with the numerical stability of the method, we then obtain error
estimates in local - or -norms that are optimal with respect to the
approximation order, the problem's stability and perturbations in data. The
convergence order is the same for both norms, but the -estimate requires
an additional divergence assumption for the convective field. The theory is
illustrated in some computational examples.Comment: 21 pages, 6 figures; in v2 we added two remarks and an appendix on
psiDOs, and made some minor correction
Stabilized nonconforming finite element methods for data assimilation in incompressible flows
We consider a stabilized nonconforming finite element method for data
assimilation in incompressible flow subject to the Stokes' equations. The
method uses a primal dual structure that allows for the inclusion of
nonstandard data. Error estimates are obtained that are optimal compared to the
conditional stability of the ill-posed data assimilation problem
Fully discrete finite element data assimilation method for the heat equation
We consider a finite element discretization for the reconstruction of the
final state of the heat equation, when the initial data is unknown, but
additional data is given in a sub domain in the space time. For the
discretization in space we consider standard continuous affine finite element
approximation, and the time derivative is discretized using a backward
differentiation. We regularize the discrete system by adding a penalty of the
-semi-norm of the initial data, scaled with the mesh-parameter. The
analysis of the method uses techniques developed in E. Burman and L. Oksanen,
Data assimilation for the heat equation using stabilized finite element
methods, arXiv, 2016, combining discrete stability of the numerical method with
sharp Carleman estimates for the physical problem, to derive optimal error
estimates for the approximate solution. For the natural space time energy norm,
away from , the convergence is the same as for the classical problem with
known initial data, but contrary to the classical case, we do not obtain faster
convergence for the -norm at the final time
A "poor man's" approach for high-resolution three-dimensional topology optimization of natural convection problems
This paper treats topology optimization of natural convection problems. A
simplified model is suggested to describe the flow of an incompressible fluid
in steady state conditions, similar to Darcy's law for fluid flow in porous
media. The equations for the fluid flow are coupled to the thermal
convection-diffusion equation through the Boussinesq approximation. The coupled
non-linear system of equations is discretized with stabilized finite elements
and solved in a parallel framework that allows for the optimization of high
resolution three-dimensional problems. A density-based topology optimization
approach is used, where a two-material interpolation scheme is applied to both
the permeability and conductivity of the distributed material. Due to the
simplified model, the proposed methodology allows for a significant reduction
of the computational effort required in the optimization. At the same time, it
is significantly more accurate than even simpler models that rely on convection
boundary conditions based on Newton's law of cooling. The methodology discussed
herein is applied to the optimization-based design of three-dimensional heat
sinks. The final designs are formally compared with results of previous work
obtained from solving the full set of Navier-Stokes equations. The results are
compared in terms of performance of the optimized designs and computational
cost. The computational time is shown to be decreased to around 5-20% in terms
of core-hours, allowing for the possibility of generating an optimized design
during the workday on a small computational cluster and overnight on a high-end
desktop
Stabilised finite element methods for ill-posed problems with conditional stability
In this paper we discuss the adjoint stabilised finite element method
introduced in, E. Burman, Stabilized finite element methods for nonsymmetric,
noncoercive and ill-posed problems. Part I: elliptic equations, SIAM Journal on
Scientific Computing, and how it may be used for the computation of solutions
to problems for which the standard stability theory given by the Lax-Milgram
Lemma or the Babuska-Brezzi Theorem fails. We pay particular attention to
ill-posed problems that have some conditional stability property and prove
(conditional) error estimates in an abstract framework. As a model problem we
consider the elliptic Cauchy problem and provide a complete numerical analysis
for this case. Some numerical examples are given to illustrate the theory.Comment: Accepted in the proceedings from the EPSRC Durham Symposium Building
Bridges: Connections and Challenges in Modern Approaches to Numerical Partial
Differential Equation
- …