372 research outputs found
Optimal-order isogeometric collocation at Galerkin superconvergent points
In this paper we investigate numerically the order of convergence of an
isogeometric collocation method that builds upon the least-squares collocation
method presented in [1] and the variational collocation method presented in
[2]. The focus is on smoothest B-splines/NURBS approximations, i.e, having
global continuity for polynomial degree . Within the framework of
[2], we select as collocation points a subset of those considered in [1], which
are related to the Galerkin superconvergence theory. With our choice, that
features local symmetry of the collocation stencil, we improve the convergence
behaviour with respect to [2], achieving optimal -convergence for odd
degree B-splines/NURBS approximations. The same optimal order of convergence is
seen in [1], where, however a least-squares formulation is adopted. Further
careful study is needed, since the robustness of the method and its
mathematical foundation are still unclear.Comment: 21 pages, 20 figures (35 pdf images
Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
We present a class of spline finite element methods for time-domain wave
propagation which are particularly amenable to explicit time-stepping. The
proposed methods utilize a discontinuous Galerkin discretization to enforce
continuity of the solution field across geometric patches in a multi-patch
setting, which yields a mass matrix with convenient block diagonal structure.
Over each patch, we show how to accurately and efficiently invert mass matrices
in the presence of curved geometries by using a weight-adjusted approximation
of the mass matrix inverse. This approximation restores a tensor product
structure while retaining provable high order accuracy and semi-discrete energy
stability. We also estimate the maximum stable timestep for spline-based finite
elements and show that the use of spline spaces result in less stringent CFL
restrictions than equivalent piecewise continuous or discontinuous finite
element spaces. Finally, we explore the use of optimal knot vectors based on L2
n-widths. We show how the use of optimal knot vectors can improve both
approximation properties and the maximum stable timestep, and present a simple
heuristic method for approximating optimal knot positions. Numerical
experiments confirm the accuracy and stability of the proposed methods
- …