I³MS - Lindgren Seminar

Location: AICES Seminar Room 115, 1st floor, Schinkelstr. 2, 52062 Aachen

Dr. Eric Lindgren - Modelling of Electrostatic Self-Assembly in Many-Body Dielectric Systems

AICES Graduate School, RWTH Aachen University


A numerical method based on a Galerkin approximation of an integral equation formulation to compute electrostatic interactions between many dielectric particles will be introduced. The method is sufficiently general, as it is able to treat systems embedded in a homogeneous dielectric medium, containing an arbitrary number of spherical particles of arbitrary size, charge, dielectric constant and position in the three-dimensional space. Simple numerical examples will be presented to illustrate the capabilities of the model, and special focus will be given to the influence of non-additive mutual polarization between particles in an electrostatic interaction. Calculations that successfully reproduce many of the observed patterns of behaviour of two experimental studies relating to electrostatic self-assembly will also be presented. The first study relates to experiments on the assembly of polymer particles that have been subjected to tribocharging, and the second study explores events observed following collisions between single particles and small clusters composed of charged particles derived from a metal oxide composite. Finally, current developments relating to the continuum treatment of ionic species in the medium will be briefly addressed.

I³MS - Rozza Seminar

Location: AICES Seminar Room 115, 1st floor, Schinkelstr. 2, 52062 Aachen

Prof. Gianluigi Rozza, Ph.D - Reduced Order Methods: State of the Art and Perspectives with a Special Focus on Computational Fluid Dynamics

SISSA, International School for Advanced Studies, Mathematics Area, mathLab, Trieste, Italy


In this talk, we provide the state of the art of Reduced Order Methods (ROM) for parametric Partial Differential Equations (PDEs), and we focus on some perspectives in their current trends and developments, with a special interest in parametric problems arising in offline-online Computational Fluid Dynamics (CFD). Systems modelled by PDEs are depending by several complex parameters in need of being reduced, even before the computational phase in a pre-processing step, in order to reduce parameter space.
Efficient parametrizations (random inputs, geometry, physics) are very important to be able to properly address an offline-online decoupling of the computational procedures and to allow competitive computational performances. Current ROM developments in CFD include: a better use of stable high fidelity methods, considering also spectral element method, to enhance the quality of the reduced model too; more efficient sampling techniques to reduce the number of the basis functions, retained as snapshots, as well as the dimension of online systems; the improvements of the certification of accuracy based on residual based error bounds and of the stability factors, as well as the the guarantee of the stability of the approximation with proper space enrichments. For nonlinear systems, also the investigation on bifurcations of parametric solutions are crucial and they may be obtained thanks to a reduced eigenvalue analysis of the linearised operator. All the previous aspects are very important in CFD problems to be able to focus in real time on complex parametric industrial and biomedical flow problems, or even in a control flow setting, and to couple viscous flows -velocity, pressure, as well as thermal field - with a structural field or a porous medium, thus requiring also an efficient reduced parametric treatment of interfaces between different physics.

I³MS - Kirchhart Seminar

Location: AICES Seminar Room 115, 1st floor, Schinkelstr. 2, 52062 Aachen

Dr. Matthias Kirchhart - Vortex Methods for Incompressible Flows

AICES Graduate School, RWTH Aachen University


Vortex methods are numerical schemes for solving the incompressible
Navier–Stokes equations. These equations accurately describe the motion
of both gases and liquids as we encounter them in everyday life, i.e.,
at velocities far below the speed of sound and not subject to extreme
temperatures or pressures. It is hard to overestimate their importance
in engineering applications, where they can for example be used to
minimise air-resistance and thereby fuel consumption of cars. However,
current numerical schemes for these equations face severe problems when
applied to turbulent flows: stringent time-step constraints,
instabilities, or the introduction of significant amounts of artificial,
spurious viscosity make their application infeasible or render the
results unusable.

Vortex methods, on the other hand, are particle methods that are based
on the vorticity formulation of the Navier–Stokes equations. This
formulation comes with two main benefits: the pressure variable is
eliminated and the system consists of separate dynamic and kinematic
parts, which can be treated independently with semi-analytical schemes.
The dynamic part is discretised using particles, which are then
convected with the flow. This natural treatment of convection renders
the method virtually free of artificial viscosity. The kinematic part of
the equations is solved using a solver based on the Biot–Savart law,
which guarantees incompressibility in the strong, point-wise sense. In
addition, in the two-dimensional case, the resulting schemes can be
shown to also conserve circulation, linear momentum, angular momentum,
and energy. These properties make vortex methods an interesting
alternative to current, mesh-based alternatives.

In this talk we will first describe a simple vortex method in the
unbounded, two-dimensional setting to illustrate the intuition of
vortex methods. We then present basic results from their analysis,
before moving on to discuss some of the problems vortex methods are
facing in three-dimensional bounded domains. We present recent research
results on one of these problems and conclude with an outlook to further
research opportunities.

[1] G.-H. Cottet and P. D. Koumoutsakos. Vortex Methods. Theory and
Practice. Cambridge University Press, 2000. ISBN: 0521621860.

[2] A. J. Majda and A. L. Bertozzi. Vorticity and Incompressible Flow.
Cambridge University Press, Nov. 2001. ISBN: 0521630576.

[3] M. Kirchhart and S. Obi. ‘A Smooth Partition of Unity Finite Element
Method for Vortex Particle Regularization’. In: SIAM Journal on
Scientific Computing 39.5 (Oct. 2017), pp. A2345–A2364. ISSN: 1064–8275.
DOI: 10.1137/17M1116258.

I³MS - Lockerby Seminar

Location: AICES Seminar Room 115, 1st floor, Schinkelstr. 2, 52062 Aachen

Prof. Duncan Lockerby, Ph.D. - Simulating Low-Speed Micro-Scale Gas Flows: A Method of Fundamental Solutions

School of Engineering, University of Warwick, UK


Creeping flow (also known as Stokes flow) is flow at vanishingly small Reynolds number. This classic topic has many applications, particularly in micro- and nano-flow technology, owing to both the small scale and the low speeds encountered. In these conditions the inertia of the fluid can be neglected, and (when local-equilibirum assumptions are valid) the momentum and continuity equation reduce to the Stokes equations. A common mathematical tool for the analysis of creeping flows is a fundamental singularity solution to the Stokes equations, known as the Stokeslet (first derived by Lorentz in 1897). This fundamental solution – in essence a Green’s function for the Stokes equations – is the flow response to a Dirac delta forcing term applied to the momentum equation.

In the most straightforward use, a flow field is represented by a superposition of Stokeslets (positioned outside of the domain) that are given a combination of strengths chosen to satisfy the same number of conditions at nodes on the boundary. This approach is known as the method of fundamental solutions (MFS) or the superposition method.  The approach has the advantage of having a flow domain that is meshless, and a dimensionality that is reduced in order by one (the boundary is discretized rather than the volume). The Stokeslet is also in itself of interest, as a fundamental solution to the Stokes equations, and can be used to conveniently derive certain analytical results.

In dilute gas flows departed far from local thermodynamic equilibrium (i.e. at non negligible Knudsen numbers, such as in rarefied conditions or at the micro/nano scale) the Stokes constitutive law becomes invalid/inaccurate. As such, the Stokeslet has limited applicability in the analysis of creeping gas flows at the micro/nano scale. However, constitutive closures exist that extend the applicability of the continuum treatment to higher Knudsen numbers. Notably, Grad’s family of moment equations(and particularly their ‘regularised’ counterparts; e.g. the R13 equations due to Struchtrup and Torrilhon) have attracted significant attention in recent years.

In this talk we introduce fundamental solutions to the the linearised R13 equations for very low Reynolds number flows; equivalent to the Stokeslet, but applicable to higher Knudsen numbers. A simple numerical implementation of the method of fundamental solutions (MFS) is presented for some three-dimensional creeping flows. Incorporation of these new fundamental solutions into the method of fundamental solutions (MFS) allows for efficient computation of three-dimensional gas microflows at remarkably low computational cost. The R13-MFS approach accurately recovers analytic solutions for low-speed flow around a stationary sphere and heat transfer from a hot sphere, capturing non-equilibrium flow phenomena missing from lower-order solutions. To demonstrate the potential of the new approach, the influence of kinetic effects on the hydrodynamic interaction between approaching solid microparticles is calculated.

Willcox Seminar

Location: AICES Seminar Room 115, 1st floor, Schinkelstr. 2, 52062

Prof. Karen Willcox, Ph.D. - Data-Driven Operator Inference for Learning Physics-Based Low-Dimensional Models

Department of Aeronautics and Astronautics, Massachusetts Institute of Technology, USA


This talk presents a non-intrusive data-driven approach for learning low-dimensional models for systems governed by time-dependent partial differential equations. Projection-based model reduction constructs the operators of a reduced model by projecting the equations of a full high-fidelity model onto a reduced space. Traditionally, this projection is intrusive, which means that the full-model operators are required explicitly in assembled form or implicitly through a routine that returns the action of the operators on a given vector; however, in many situations these full-model operators may be inaccessible. Our non-intrusive operator inference approach solves an optimization problem to infer approximations of the reduced operators directly from input and state data, without requiring the full model itself. The inferred operators are the solution of a least-squares problem and converge, with sufficient state trajectory data, in the Frobenius norm to the reduced operators that would be obtained via an intrusive projection of the full-model operators. Thus, while the approach is data-driven, it also embeds the physical constraints associated with the underlying system governing equations. Joint work with Benjamin Peherstorfer (U. Wisconsin Madison).