|
| |
|
|
The effect of irregular interfaces on the BDDC method for the Navier-Stokes equations
Hanek, M. ; Šístek, Jakub ; Burda, P.
We investigate the effect of interface irregularity on the convergence of the BDDC method for Navier-Stokes equations. A benchmark problem of a sequence of contracting channels is proposed to evaluate the robustness of the iterative solver with respect to element aspect ratios at the interface. Partitioners based on graph of the mesh and the geometry of the domain are compared. It is shown, that the convergence is significantly improved by avoiding irregular interfaces for the benchmark problem as well as for an industrial problem of oil flow in hydrostatic bearing.
|
|
|
A particular smooth interpolation that generates splines
Segeth, Karel
There are two grounds the spline theory stems from -- the algebraic one (where splines are understood as piecewise smooth functions satisfying some continuity conditions) and the variational one (where splines are obtained via minimization of some quadratic functionals with constraints). We use the general variational approach called $it smooth interpolation$ introduced by Talmi and Gilat and show that it covers not only the cubic spline and its 2D and 3D analogues but also the well known tension spline (called also spline with tension). We present the results of a 1D numerical example that characterize some properties of the tension spline.
|
|
|
On the quality of local flux reconstructions for guaranteed error bounds
Vejchodský, Tomáš
In this contribution we consider elliptic problems of a reaction-diffucion type discretized by the finite element method and study the quality of guaranteed upper bounds of the error. In particular, we concentrate on complementary error bounds whose values are determined by suitable flux reconstructions. We present numerical experiments comparing the performance of the local flux reconstruction of Ainsworth and Vejchodský [2] and the reconstruction of Braess and Schröberl [5]. We evaluate the efficiency of these flux reconstructions by their comparison with the optimal flux reconstruction computed as a global minimization problem.
|
|
|
A note on tension spline
Segeth, Karel
Spline theory is mainly grounded on two approaches: the algebraic one (where splines are understood as piecewise smooth functions) and the variational one (where splines are obtained via minimization of quadratic functionals with constraints). We show that the general variational approach called smooth interpolation introduced by Talmi and Gilat covers not only the cubic spline but also the well known tension spline (called also spline in tension or spline with tension). We present the results of a 1D numerical example that show the advantages and drawbacks of the tension spline.
|
|
|
Smooth approximation spaces based on a periodic system
Segeth, Karel
A way of data approximation called smooth was introduced by Talmi and Gilat in 1977. Such an approach employs a (possibly infinite) linear combination of smooth basis functions with coefficients obtained as the unique solution of a minimization problem. While the minimization guarantees the smoothness of the approximant and its derivatives, the constraints represent the interpolating or smoothing conditions at nodes. In the contribution, a special attention is paid to the periodic basis system $exp(-ii kx)$. A 1D numerical example is presented.
|
|
|
An application of the BDDC method to the Navier-Stokes equations in 3-D cavity
Hanek, M. ; Šístek, Jakub ; Burda, P.
We deal with numerical simulation of incompressible flow governed by the Navier-Stokes equations. The problem is discretised using the finite element method, and the arising system of nonlinear equations is solved by Picard iteration. We explore the applicability of the Balancing Domain Decomposition by Constraints (BDDC) method to nonsymmetric problems arising from such linearisation. One step of BDDC is applied as the preconditioner for the stabilized variant of the biconjugate gradient (BiCGstab) method. We present results for a 3-D cavity problem computed on 32 cores of a parallel supercomputer.
|
|
|
Different approaches to interface weights in the BDDC method in 3D
Čertíková, M. ; Šístek, Jakub ; Burda, P.
In this paper, we discuss the choice of weights in averaging of local (subdomain) solutions on the interface for the BDDC method (Balancing Domain Decomposition by Constraints). We try to find relations among different choices of the interface weights and compare them numerically on model problems of the Poisson equation and linear elasticity in 3D. Problems with jumps in coefficients of material properties are considered and both regular and irregular interfaces between subdomains are tested.
|
|
|
A parallel finite element solver for unsteady incompressible Navier-Stokes equations
Šístek, Jakub
A parallel solver for unsteady incompressible Navier-Stokes equations is presented. It is based on the finite element method combined with the pressure-correction approach. Semi-implicit treatment of the convective term is considered, leading to five systems of linear algebraic equations to be solved in each time step. Krylov subspace iterative methods are employed for the solution of these systems with a particular emphasis on efficient parallel preconditioners. A simulation of a benchmark problem of incompressible viscous flow around a sphere at Reynolds number 300 is presented and compared with literature.
|
host ::
RSS.