Geometric partial differential equations: Theory, numerics and applications

This workshop concentrated on partial differential equations involving stationary and evolving surfaces in which geometric quantities play a major role. Mutual interest in this emerging field stimulated the interaction between analysis, numerical solution, and applications

Space-time Methods for Time-dependent Partial Differential Equations

Modern discretizations of time-dependent PDEs consider the full problem in the space-time cylinder and aim to overcome limitations of classical approaches such as the method of lines (first discretize in space and then solve the resulting ODE) and the Rothe method (first discretize in time and then solve the PDE). A main advantage of a holistic space-time method is the direct access to space-time adaptivity and to the backward problem (required for the dual problem in optimization or error control). Moreover, this allows for parallel solution strategies simultaneously in time and space. Several space-time concepts where proposed (different conforming and nonconforming space-time finite elements, the parareal method, wavefront relaxation etc.) but this topic has become a rapidly growing field in numerical analysis and scientific computing. In this workshop the focus is the development of adaptive and flexible space-time discretization methods for solving parabolic and hyperbolic space-time partial differential equations

Spectral methods for partial differential equations

Origins of spectral methods, especially their relation to the Method of Weighted Residuals, are surveyed. Basic Fourier, Chebyshev, and Legendre spectral concepts are reviewed, and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic, and mixed type. Fluid dynamical applications are emphasized

Higher-order implicit-explicit multi-domain compressible Navier-Stokes solvers

This paper presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional multi-domains. Building up on the recent single-domain ADI-based high-order Navier-Stokes solvers (Bruno and Cubillos, Journal of Computational Physics 307 (2016) 476-495) this article presents multi-domain implicit-explicit methods of high-order of temporal accuracy. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of high-order backward differentiation formulae (BDF) and an alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Nearly dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. As demonstrated via a variety of numerical experiments in two and three dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, robust stability properties, limited dispersion, and high parallel efficiency

Solution of 3-dimensional time-dependent viscous flows. Part 2: Development of the computer code

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow equations to aid in the design of helicopter rotors. The development of a computer code to solve a three dimensional unsteady approximate form of the Navier-Stokes equations employing a linearized block emplicit technique in conjunction with a QR operator scheme is described. Results of calculations of several Cartesian test cases are presented. The computer code can be applied to more complex flow fields such as these encountered on rotating airfoils

Constraint-consistent Runge-Kutta methods for one-dimensional incompressible multiphase flow

New time integration methods are proposed for simulating incompressible multiphase flow in pipelines described by the one-dimensional two-fluid model. The methodology is based on 'half-explicit' Runge-Kutta methods, being explicit for the mass and momentum equations and implicit for the volume constraint. These half-explicit methods are constraint-consistent, i.e., they satisfy the hidden constraints of the two-fluid model, namely the volumetric flow (incompressibility) constraint and the Poisson equation for the pressure. A novel analysis shows that these hidden constraints are present in the continuous, semi-discrete, and fully discrete equations. Next to constraint-consistency, the new methods are conservative: the original mass and momentum equations are solved, and the proper shock conditions are satisfied; efficient: the implicit constraint is rewritten into a pressure Poisson equation, and the time step for the explicit part is restricted by a CFL condition based on the convective wave speeds; and accurate: achieving high order temporal accuracy for all solution components (masses, velocities, and pressure). High-order accuracy is obtained by constructing a new third order Runge-Kutta method that satisfies the additional order conditions arising from the presence of the constraint in combination with time-dependent boundary conditions. Two test cases (Kelvin-Helmholtz instabilities in a pipeline and liquid sloshing in a cylindrical tank) show that for time-independent boundary conditions the half-explicit formulation with a classic fourth-order Runge-Kutta method accurately integrates the two-fluid model equations in time while preserving all constraints. A third test case (ramp-up of gas production in a multiphase pipeline) shows that our new third order method is preferred for cases featuring time-dependent boundary conditions

A High-Order Radial Basis Function (RBF) Leray Projection Method for the Solution of the Incompressible Unsteady Stokes Equations

A new projection method based on radial basis functions (RBFs) is presented for discretizing the incompressible unsteady Stokes equations in irregular geometries. The novelty of the method comes from the application of a new technique for computing the Leray-Helmholtz projection of a vector field using generalized interpolation with divergence-free and curl-free RBFs. Unlike traditional projection methods, this new method enables matching both tangential and normal components of divergence-free vector fields on the domain boundary. This allows incompressibility of the velocity field to be enforced without any time-splitting or pressure boundary conditions. Spatial derivatives are approximated using collocation with global RBFs so that the method only requires samples of the field at (possibly scattered) nodes over the domain. Numerical results are presented demonstrating high-order convergence in both space (between 5th and 6th order) and time (up to 4th order) for some model problems in two dimensional irregular geometries.Comment: 34 pages, 8 figure