Preconditioning of weighted H(div)-norm and applications to numerical simulation of highly heterogeneous media

In this paper we propose and analyze a preconditioner for a system arising from a finite element approximation of second order elliptic problems describing processes in highly het- erogeneous media. Our approach uses the technique of multilevel methods and the recently proposed preconditioner based on additive Schur complement approximation by J. Kraus (see [8]). The main results are the design and a theoretical and numerical justification of an iterative method for such problems that is robust with respect to the contrast of the media, defined as the ratio between the maximum and minimum values of the coefficient (related to the permeability/conductivity).Comment: 28 page

Robust Methods for Multiscale Coarse Approximations of Diffusion Models in Perforated Domains

For the Poisson equation posed in a domain containing a large number of polygonal perforations, we propose a low-dimensional coarse approximation space based on a coarse polygonal partitioning of the domain. Similarly to other multiscale numerical methods, this coarse space is spanned by locally discrete harmonic basis functions. Along the subdomain boundaries, the basis functions are piecewise polynomial. The main contribution of this article is an error estimate regarding the H1-projection over the coarse space which depends only on the regularity of the solution over the edges of the coarse partitioning. For a specific edge refinement procedure, the error analysis establishes superconvergence of the method even if the true solution has a low general regularity. Combined with domain decomposition (DD) methods, the coarse space leads to an efficient two-level iterative linear solver which reaches the fine-scale finite element error in few iterations. It also bodes well as a preconditioner for Krylov methods and provides scalability with respect to the number of subdomains. Numerical experiments showcase the increased precision of the coarse approximation as well as the efficiency and scalability of the coarse space as a component of a DD algorithm.Comment: 32 pages, 14 figures, submitted to Journal of Computational Physic

On the linear convergence of additive Schwarz methods for the pp-Laplacian

We consider additive Schwarz methods for boundary value problems involving the $p$-Laplacian. While the existing theoretical estimates for the convergence rate of additive Schwarz methods for the $p$-Laplacian are sublinear, the actual convergence rate observed by numerical experiments is linear. In this paper, we bridge the gap between these theoretical and numerical results by analyzing the linear convergence rate of additive Schwarz methods for the $p$-Laplacian. In order to estimate the linear convergence rate of the methods, we present two essential components. Firstly, we present a new abstract convergence theory of additive Schwarz methods written in terms of a quasi-norm. This quasi-norm exhibits behavior similar to the Bregman distance of the convex energy functional associated to the problem. Secondly, we provide a quasi-norm version of the Poincar'{e}--Friedrichs inequality, which is essential for deriving a quasi-norm stable decomposition for a two-level domain decomposition setting.Comment: 23 pages, 2 figure

Multiscale, Multiphysics Modelling of Coastal Ocean Processes: Paradigms and Approaches

This Special Issue includes papers on physical phenomena, such as wind-driven flows, coastal flooding, and turbidity currents, and modeling techniques, such as model comparison, model coupling, parallel computation, and domain decomposition. These papers illustrate the need for modeling coastal ocean flows with multiple physical processes at different scales. Additionally, these papers reflect the current status of such modeling of coastal ocean flows, and they present a roadmap with numerical methods, data collection, and artificial intelligence as future endeavors

THE HP VERSION OF THE VIRTUAL ELEMENT METHOD

The interest in Galerkin methods for the approximation of partial differential equations on polytopal meshes has recently grown. The virtual element method (VEM) is one of the most successful approaches enabling computation on such meshes. So far, only the h version of the method has been investigated; here, the convergence of the error is obtained by keeping fixed the dimension of local spaces, while refining the mesh. Contrarily, the p version of a Galerkin method consists in achieving convergence by keeping fixed the decomposition of the domain and increasing the dimension of local spaces. The combination of the h and p versions goes under the name of hp version of the method under consideration. The present thesis aims to dovetail the technology of VEM with the p and the hp refinement strategies. Particular emphasis is offered to the approximation of Poisson and Laplace problems, a priori and a posteriori error analysis, multigrid solvers, stabilization of the method and ill-conditioning of the stiffness matrix

Proceedings of SIRM 2023 - The 15th European Conference on Rotordynamics

It was our great honor and pleasure to host the SIRM Conference after 2003 and 2011 for the third time in Darmstadt. Rotordynamics covers a huge variety of different applications and challenges which are all in the scope of this conference. The conference was opened with a keynote lecture given by Rainer Nordmann, one of the three founders of SIRM “Schwingungen in rotierenden Maschinen”. In total 53 papers passed our strict review process and were presented. This impressively shows that rotordynamics is relevant as ever. These contributions cover a very wide spectrum of session topics: fluid bearings and seals; air foil bearings; magnetic bearings; rotor blade interaction; rotor fluid interactions; unbalance and balancing; vibrations in turbomachines; vibration control; instability; electrical machines; monitoring, identification and diagnosis; advanced numerical tools and nonlinearities as well as general rotordynamics. The international character of the conference has been significantly enhanced by the Scientific Board since the 14th SIRM resulting on one hand in an expanded Scientific Committee which meanwhile consists of 31 members from 13 different European countries and on the other hand in the new name “European Conference on Rotordynamics”. This new international profile has also been emphasized by participants of the 15th SIRM coming from 17 different countries out of three continents. We experienced a vital discussion and dialogue between industry and academia at the conference where roughly one third of the papers were presented by industry and two thirds by academia being an excellent basis to follow a bidirectional transfer what we call xchange at Technical University of Darmstadt. At this point we also want to give our special thanks to the eleven industry sponsors for their great support of the conference. On behalf of the Darmstadt Local Committee I welcome you to read the papers of the 15th SIRM giving you further insight into the topics and presentations