Bivariate polynomial interpolation on the square at new nodal sets

As known, the problem of choosing ``good'' nodes is a central one in polynomial interpolation. While the problem is essentially solved in one dimension (all good nodal sequences are asymptotically equidistributed with respect to the arc-cosine metric), in several variables it still represents a substantially open question. In this work we consider new nodal sets for bivariate polynomial interpolation on the square. First, we consider fast Leja points for tensor-product interpolation. On the other hand, for classical polynomial interpolation on the square we experiment four families of points which are (asymptotically) equidistributed with respect to the Dubiner metric, which extends to higher dimension the arc-cosine metric. One of them, nicknamed Padua points, gives numerically a Lebesgue constant growing like log square of the degree

Tensor Product Gauss-Lobatto Points are Fekete Points for the Cube

We show that the tensor product of Gauss-Lobatto points are Fekete points (for the tensor product basis) for the cube

hp-FEM for Two-component Flows with Applications in Optofluidics

This thesis is concerned with the application of hp-adaptive finite element methods to a mathematical model of immiscible two-component flows. With the aim of simulating the flow processes in microfluidic optical devices based on liquid-liquid interfaces, we couple the time-dependent incompressible Navier-Stokes equations with a level set method to describe the flow of the fluids and the evolution of the interface between them

Nodale Spektralelemente und unstrukturierte Gitter - Methodische Aspekte und effiziente Algorithmen

Die Dissertation behandelt methodische und algorithmische Aspekte der Spektralelementemethode zur räumlichen Diskretisierung partieller Differentialgleichungen. Die Weiterentwicklung einer symmetriebasierten Faktorisierung ermöglicht effiziente Operatoren für Tetraederelemente. Auf Grundlage einer umfassenden Leistungsanalyse werden Engpässe in der Implementierung der Operatoren identifiziert und durch algorithmische Modifikationen der Methode eliminiert

Spectral element approximation of the incompressible Navier-Stokes equations in a moving domain and applications

In this thesis we address the numerical approximation of the incompressible Navier-Stokes equations evolving in a moving domain with the spectral element method and high order time integrators. First, we present the spectral element method and the basic tools to perform spectral discretizations of the Galerkin or Galerkin with Numerical Integration (G-NI) type. We cover a large range of possibilities regarding the reference elements, basis functions, interpolation points and quadrature points. In this approach, the integration and differentiation of the polynomial functions is done numerically through the help of suitable point sets. Regarding the differentiation, we present a detailed numerical study of which points should be used to attain better stability (among the choices we present). Second, we introduce the incompressible steady/unsteady Stokes and Navier-Stokes equations and their spectral approximation. In the unsteady case, we introduce a combination of Backward Differentiation Formulas and an extrapolation formula of the same order for the time integration. Once the equations are discretized, a linear system must be solved to obtain the approximate solution. In this context, we consider the solution of the whole system of equations combined with a block type preconditioner. The preconditioner is shown to be optimal in terms of number of iterations used by the GMRES method in the steady case, but not in the unsteady one. Another alternative presented is to use algebraic factorization methods of the Yosida type and decouple the calculation of velocity and pressure. A benchmark is also presented to access the numerical convergence properties of this type of methods in our context. Third, we extend the algorithms developed in the fixed domain case to the Arbitrary Lagrangian Eulerian framework. The issue of defining a high order ALE map is addressed. This allows to construct a computational domain that is described with curved elements. A benchmark using a direct method to solve the linear system or the Yosida-q methods is presented to show the convergence orders of the method proposed. Finally, we apply the developed method with an implicit fully coupled and semi-implicit approach, to solve a fluid-structure interaction problem for a simple 2D hemodynamics example