58 research outputs found
Each H
We consider the solution of second order elliptic PDEs in Rd
with inhomogeneous Dirichlet data by means of an hâadaptive FEM with
fixed polynomial order p â N. As model example serves the Poisson
equation with mixed DirichletâNeumann boundary conditions, where the inhomogeneous
Dirichlet data are discretized by use of an H1Â /Â 2âstable
projection, for instance, the L2âprojection for
p = 1 or the ScottâZhang projection for general p â„ 1.
For error estimation, we use a residual error estimator which includes the Dirichlet data
oscillations. We prove that each H1Â /Â 2âstable projection
yields convergence of the adaptive algorithm even with quasiâoptimal convergence rate.
Numerical experiments with the ScottâZhang projection conclude the work
A study on spline quasi-interpolation based quadrature rules for the isogeometric Galerkin BEM
Two recently introduced quadrature schemes for weakly singular integrals
[Calabr\`o et al. J. Comput. Appl. Math. 2018] are investigated in the context
of boundary integral equations arising in the isogeometric formulation of
Galerkin Boundary Element Method (BEM). In the first scheme, the regular part
of the integrand is approximated by a suitable quasi--interpolation spline. In
the second scheme the regular part is approximated by a product of two spline
functions. The two schemes are tested and compared against other standard and
novel methods available in literature to evaluate different types of integrals
arising in the Galerkin formulation. Numerical tests reveal that under
reasonable assumptions the second scheme convergences with the optimal order in
the Galerkin method, when performing -refinement, even with a small amount
of quadrature nodes. The quadrature schemes are validated also in numerical
examples to solve 2D Laplace problems with Dirichlet boundary conditions
Multiscale modeling in micromagnetics : existence of solutions and numerical integration
Various applications ranging from spintronic devices, giant magnetoresistance sensors, and magnetic storage devices, include magnetic parts on very different length scales. Since the consideration of the Landau-Lifshitz-Gilbert equation (LLG) constrains the maximum element size to the exchange length within the media, it is numerically not attractive to simulate macroscopic parts with this approach. On the other hand, the magnetostatic Maxwell equations do not constrain the element size, but cannot describe the short-range exchange interaction accurately. A combination of both methods allows one to describe magnetic domains within the micromagnetic regime by use of LLG and also considers the macroscopic parts by a nonlinear material law using the Maxwell equations. In our work, we prove that under certain assumptions on the nonlinear material law, this multiscale version of LLG admits weak solutions. Our proof is constructive in the sense that we provide a linear-implicit numerical integrator for the multiscale model such that the numerically computable finite element solutions admit weak H1-convergence (at least for a subsequence) towards a weak solution
Wavelet boundary element methods â Adaptivity and goal-oriented error estimation
This article is dedicated to the adaptive wavelet boundary element method. It computes an approximation to the unknown solution of the boundary integral equation under consideration with a rate , whenever the solution can be approximated with this rate in the setting determined by the underlying wavelet basis. The computational cost scale linearly in the number of degrees of freedom. Goal-oriented error estimation for evaluating linear output functionals of the solution is also considered. An algorithm is proposed that approximately evaluates a linear output functional with a rate , whenever the primal solution can be approximated with a rate and the dual solution can be approximated with a rate , while the cost still scale linearly in . Numerical results for an acoustic scattering problem and for the point evaluation of the potential in case of the Laplace equation are reported to validate and quantify the approach
Linear elastic fracture simulation directly from CAD: 2D NURBS-based implementation and role of tip enrichment
Adaptive FEM with Optimal Convergence Rates for a Certain Class of Nonsymmetric and Possibly Nonlinear Problems
Optimal convergence for adaptive IGA boundary element methods for weakly-singular integral equations
Stability of the derivative of a canonical product
With each sequence α=(αn)nâN of pairwise distinct and non-zero points which are such that the canonical product  Pα(z):=limrââââŁÎ±nâŁâ€r(1âz/αn) converges, the sequence  αâČ:=(Pα'(αn))nâN is associated. We give conditions on the difference ÎČâα of two sequences which ensure that ÎČ' and α' are comparable in the sense that  âc,C>0: c|α'n|â€|ÎČ'n|â€C|α'n|, nâN. The values α'n play an important role in various contexts. As a selection of applications we present: an inverse spectral problem, a class of entire functions and a continuation problem
- âŠ