A main research field is the development, analysis, improvement and application of numerical methods for equations coming from CFD. These methods are based on finite element or finite volume methods as spatial discretization and implicit temporal discretizations.
Convectiondiffusionreaction equations model, for instance, transport processes which are driven by diffusion or convection. A challenging situation from the point of view of simulations is the c ase of dominant convection, because the solution has scales which cannot be resolved by the used grids (layers). Standard stabilized methods lead to smeared discrete solutions or to solutions possessing spurious oscillations. Spurious oscillations are not acceptable in many applications. On research topic is the development of oscillationfree numerical methods for convectiondominated equations. In addition, numerical methods for nonlinear transport equations with exponential nonlinearities are developed and analyzed.
The NavierStokes equations are the fundamental equations of fluid mechanics. A main topic of research are methods for the simulation of turbulent flows, in particular large eddy simulation (LES) and variational multiscale (VMS) methods. The basic difficulty in the simulation of turbulent flows consists also in the impossibility of resolving small flow structures (scales) on the used grids. LES and VMS methods are two approaches for the simulation of large scales, which model the effect of the unresolved small scales onto the large scales via a turbulence model. Another topic of research is the development of discretely massconserving discretizations. This property is essential for some of the considered applications.
In these applications, coupled systems of equations have to be solved that contain (among others) the types of equations described above. An example are population balance systems. These are coupled systems consisting of the NavierStokes equations (flow), convectiondiffusionreaction equations (chemical reaction, transport of energy) and transport equations for the particle size distribution. These transport equations are defined in a higher dimension than the other equations because the particle size distribution depends also on properties of the particles (internal coordinates). Another example are van Roosbroeck systems which describe the transport of carriers in a semiconductor crystal lattice. The equations of such systems contain exponential nonlinearities. Numerical methods for StokesDarcy systems, which model the coupling of free flows and flows in porous media, are investigated, too. In addition, it is studied to which extend reduced order models are able to increase the efficiency of simulations with only an acceptable loss of accuracy.
Flow through an aortic arch
Turbulent flow around a cylinder, velocity (left) and vorticity (right)
Publications
Monographs

R. Ahrens, Z. Lakdawala, A. Voigt, V. Wiedmeyer, V. John, S. Le Borne, K. Sundmacher, Numerical Methods for Coupled Population Balance Systems Applied to the Dynamical Simulation of Crystallization Processes, S. Heinrich, ed., Dynamic Flowsheet Simulation of Solids Processes, Springer, Cham, 2020, pp. 475518, (Chapter Published), DOI 10.1007/9783030451684_14 .

V. John, P. Knobloch, U. Wilbrandt, Finite Element Pressure Stabilizations for Incompressible Flow Problems, T. Bodnár, G. Galdi, Š. Nečasová, eds., Fluids Under Pressure, Birkhäuser, Cham, 2020, pp. pp 483573, (Chapter Published), DOI 10.1007/9783030396398_6 .
Abstract
Discretizations of incompressible flow problems with pairs of finite element spaces that do not satisfy a discrete infsup condition require a socalled pressure stabilization. This paper gives an overview and systematic assessment of stabilized methods, including the respective error analysis. 
U. Wilbrandt, StokesDarcy Equations  Analytic and Numerical Analysis, D. Bresch, V. John, M. Hieber, I. Kukavica, J. Robinson, Y. Shibata, eds., Lecture Notes in Mathematical Fluid Mechanics, Birkhäuser, Basel, 2019, 212 pages, (Monograph Published), DOI 10.1007/9783030029043 .

M. Hintermüller, M. Hinze, Ch. Kahle, T. Keil, Chapter 13: Fully Adaptive and Integrated Numerical Methods for the Simulation and Control of Variable Density Multiphase Flows Governed by Diffuse Interface Models, in: Transport Processes at Fluidic Interfaces, D. Bothe, A. Reusken, eds., Advances in Mathematical Fluid Mechanics, Birkhäuser, Springer International Publishing AG, Cham, 2017, pp. 305353, (Chapter Published), DOI 10.1007/9783319566023 .

V. John, Finite Element Methods for Incompressible Flow Problems, 51 of Springer Series in Computational Mathematics, Springer International Publishing AG, Cham, 2016, xiii+812 pages, (Monograph Published).

K. Gärtner, H. Si, A. Rand, N. Walkington, Chapter 11: 3D Delaunay Mesh Generation, in: Combinatorial Scientific Computing, U. Naumann, O. Schenk, eds., Computational Science Series, CRC Computational Science/Chapman & Hall, Boca Raton, 2012, pp. 299319, (Chapter Published).
Articles in Refereed Journals

M. Akbas, Th. Gallouët, A. Gassmann, A. Linke, Ch. Merdon, A gradientrobust wellbalanced scheme for the compressible isothermal Stokes problem, Computer Methods in Applied Mechanics and Engineering, 367 (2020), 113069, DOI 10.1016/j.cma.2020.113069 .
Abstract
A novel notion for constructing a wellbalanced scheme  a gradientrobust scheme  is introduced and a showcase application for a steady compressible, isothermal Stokes equations is presented. Gradientrobustness means that arbitrary gradient fields in the momentum balance are wellbalanced by the discrete pressure gradient  if there is enough mass in the system to compensate the force. The scheme is asymptoticpreserving in the sense that it degenerates for low Mach numbers to a recent infsup stable and pressurerobust discretization for the incompressible Stokes equations. The convergence of the coupled FEMFVM scheme for the nonlinear, isothermal Stokes equations is proved by compactness arguments. Numerical examples illustrate the numerical analysis, and show that the novel approach can lead to a dramatically increased accuracy in nearlyhydrostatic low Mach number flows. Numerical examples also suggest that a straightforward extension to barotropic situations with nonlinear equations of state is feasible. 
N. Ahmed, V. John, An assessment of two classes of variational multiscale methods for the simulation of incompressible turbulent flows, Computer Methods in Applied Mechanics and Engineering, 365 (2020), pp. 112997/1112997/20, DOI https://doi.org/10.1016/j.cma.2020.112997 .
Abstract
A numerical assessment of two classes of variational multiscale (VMS) methods for the simulation of incompressible flows is presented. Two types of residualbased VMS methods and two types of projectionbased VMS methods are included in this assessment. The numerical simulations are performed at turbulent channel flow problems with various friction Reynolds numbers. It turns out the the residualbased VMS methods, in particular when used with a pair of infsup stable finite elements, give usually the most accurate results for second order statistics. For this pair of finite element spaces, a flexible GMRES method with a Least Squares Commutator (LSC) preconditioner proved to be an efficient solver. 
D. Janke, A. Caiazzo, N. Ahmed, N. Alia, O. Knoth, B. Moreau, U. Wilbrandt, D. Willink, Th. Amon, V. John, On the feasibility of using open source solvers for the simulation of a turbulent air flow in a dairy barn, Computers and Electronics in Agriculture, 175 (2020), pp. 105546/1105546/16, DOI 10.1016/j.compag.2020.105546 .
Abstract
Two transient open source solvers, OpenFOAM and ParMooN, are assessed with respect to the simulation of the turbulent air flow inside and around a dairy barn. For this purpose, data were obtained in an experimental campaign at a 1:100 scaled wind tunnel model. Both solvers used different meshes, discretization schemes, and turbulence models. The experimental data and numerical results agree well for timeaveraged streamwise and verticalwise velocities. In particular, the air exchange was predicted with high accuracy by both solvers with relative errors less than 5 % compared to the experimental results. With respect to the turbulent quantities, good agreements at the second (downwind) half of the barn inside and especially outside the barn could be achieved, where both codes accurately predicted the flow separation and the rootmeansquare velocities. Deviations between simulations and experimental results regarding turbulent quantities could be observed in the first part of the barn, due to different inlet conditions between the experimental setup and the numerical simulations. Both solvers proved to be promising tools for the accurate prediction of timedependent phenomena in an agricultural context, e.g., like the transport of particulate matter or pathogenladen aerosols in and around agricultural buildings. 
C.K. Macnamara, A. Caiazzo, I. RamisConde, M.A.J. Chaplain, Computational modelling and simulation of cancer growth and migration within a 3D heterogeneous tissue: The effects of fibre and vascular structure, Journal of Computational Science, 40 (2020), 101067, DOI 10.1016/j.jocs.2019.101067 .
Abstract
The term cancer covers a multitude of bodily diseases, broadly categorised by having cells which do not behave normally. Since cancer cells can arise from any type of cell in the body, cancers can grow in or around any tissue or organ making the disease highly complex. Our research is focused on understanding the specific mechanisms that occur in the tumour microenvironment via mathematical and computational modeling. We present a 3D individualbased model which allows one to simulate the behaviour of, and spatiotemporal interactions between, cells, extracellular matrix fibres and blood vessels. Each agent (a single cell, for example) is fully realised within the model and interactions are primarily governed by mechanical forces between elements. However, as well as the mechanical interactions we also consider chemical interactions, for example, by coupling the code to a finite element solver to model the diffusion of oxygen from blood vessels to cells. The current state of the art of the model allows us to simulate tumour growth around an arbitrary bloodvessel network or along the striations of fibrous tissue. 
A. Linke, Ch. Merdon, M. Neilan, Pressurerobustness in quasioptimal a priori estimates for the Stokes problem, Electronic Transactions on Numerical Analysis, 52 (2020), pp. 281294, DOI 10.1553/etna_vol52s281 .

D. Peschka, S. Haefner, L. Marquant, K. Jacobs, A. Münch, B. Wagner, Signatures of slip in dewetting polymer films, Proceedings of the National Academy of Sciences of the United States of America, 116 (2019), pp. 92759284, DOI 10.1073/pnas.1820487116 .

N.R. Gauger, A. Linke, P. Schroeder, On highorder pressurerobust space discretisations, their advantages for incompressible high Reynolds number generalised Beltrami flows and beyond, SMAI Journal of Computational Mathematics, 5 (2019), pp. 89129.
Abstract
Recently, highorder space discretisations were proposed for the numerical simulation of the incompressible NavierStokes equations at high Reynolds numbers, even for complicated domains from simulation practice. Although the overall spatial approximation order of the algorithms depends on the approximation quality of the boundary (often not better than third order), competitively accurate and efficient results were reported. In this contribution, first, a possible explanation for this somewhat surprising result is proposed: the velocity error of highorder space discretisations is more robust against quantitatively large and complicated pressure fields than loworder methods. Second, it is demonstrated that novel pressurerobust methods are significantly more accurate than comparable classical, nonpressurerobust space discretisations, whenever the quadratic, nonlinear convection term is a nontrivial gradient field like in certain generalised Beltrami flows at high Reynolds number. Then, pressurerobust methods even allow to halve the (formal) approximation order without compromising the accuracy. Third, classical highorder space discretisations are outperformed by pressurerobust methods whenever the boundary is not approximated with highorder accuracy. This improved accuracy of (loworder) pressurerobust mixed methods is explained in terms of a HelmholtzHodge projector, which cancels out the nonlinear convection term in any generalised Beltrami flow, since it is a gradient field. The numerical results are illustrated by a novel numerical analysis for pressurerobust and classical space discretisations. Further, the relevance of these results is discussed for flows that are not of Beltrami type. 
L. Heltai, A. Caiazzo, Multiscale modeling of vascularized tissues via nonmatching immersed methods, International Journal of Numerical Methods in Biomedical Engineering, 35 (2019), pp. 3264/13264/32, DOI 10.1002/cnm.3264 .
Abstract
We consider a multiscale approach based on immersed methods for the efficient computational modeling of tissues composed of an elastic matrix (in two or threedimensions) and a thin vascular structure (treated as a codimension two manifold) at a given pressure. We derive different variational formulations of the coupled problem, in which the effect of the vasculature can be surrogated in the elasticity equations via singular or hypersingular forcing terms. These terms only depends on information defined on codimension two manifolds (such as vessel center line, cross sectional area, and mean pressure over cross section), thus drastically reducing the complexity of the computational model. We perform several numerical tests, ranging from simple cases with known exact solutions to the modeling of materials with random distributions of vessels. In the latter case, we use our immersed method to perform an in silico characterization of the mechanical properties of the effective biphasic material tissue via statistical simulations. 
A. Jha, V. John, A study of solvers for nonlinear AFC discretizations of convectiondiffusion equations, Computational & Applied Mathematics, 78 (2019), pp. 31173138, DOI 10.1016/j.camwa.2019.04.020 .

P.L. Lederer, Ch. Merdon, J. Schöberl, Refined a posteriori error estimation for classical and pressurerobust Stokes finite element methods, Journal of Numerical Mathematics, 142 (2019), pp. 713748.
Abstract
Recent works showed that pressurerobust modifications of mixed finite element methods for the Stokes equations outperform their standard versions in many cases. This is achieved by divergencefree reconstruction operators and results in pressureindependent velocity error estimates which are robust with respect to small viscosities. In this paper we develop a posteriori error control which reflects this robustness. 
L.O. Müller, A. Caiazzo, P.J. Blanco, Reducedorder unscented Kalman filter with observations in the frequency domain: Application to computational hemodynamics, IEEE Transactions on Biomedical Engineering, 66 (2019), pp. 12691276, DOI 10.1109/TBME.2018.2872323 .
Abstract
Objective: The aim of this work is to assess the potential of the reduced order unscented Kalman filter (ROUKF) in the context of computational hemodynamics, in order to estimate cardiovascular model parameters when employing real patientspecific data. Methods: The approach combines an efficient blood flow solver for onedimensional networks (for the forward problem) with the parameter estimation problem cast in the frequency space. Namely, the ROUKF is used to correct model parameter after each cardiac cycle, depending on the discrepancies of model outputs with respect to available observations properly mapped into the frequency space. Results: First we validate the filter in frequency domain applying it in the context of a set of experimental measurements for an in vitro model. Second, we perform different numerical experiments aiming at parameter estimation using patientspecific data. Conclusion: Our results demonstrate that the filter in frequency domain allows a faster and more robust parameter estimation, when compared to its time domain counterpart. Moreover, the proposed approach allows to estimate parameters that are not directly related to the network but are crucial for targeting interindividual parameter variability (e.g., parameters that characterize the cardiac output). Significance: The ROUKF in frequency domain provides a robust and flexible tool for estimating parameters related to cardiovascular mathematical models using in vivo data. 
P.W. Schroeder, V. John, P.L. Lederer, Ch. Lehrenfeld, G. Lube, J. Schöberl, On reference solutions and the sensitivity of the 2D KelvinHelmholtz instability problem, Computers & Mathematics with Applications. An International Journal, 77 (2019), pp. 10101028, DOI 10.1016/j.camwa.2018.10.030 .

J. DE Frutos, B. Garc'iaArchilla, V. John, J. Novo, Error analysis of non infsup stable discretizations of the timedependent NavierStokes equations with local projection stabilization, IMA Journal of Numerical Analysis, 39 (2019), pp. 17471786, DOI 10.1093/imanum/dry044 .

C. Bartsch, V. John, R.I.A. Patterson, Simulations of an ASA flow crystallizer with a coupled stochasticdeterministic approach, Comput. Chem. Engng., 124 (2019), pp. 350363, DOI 10.1016/j.compchemeng.2019.01.012 .
Abstract
A coupled solver for population balance systems is presented, where the flow, temperature, and concentration equations are solved with finite element methods, and the particle size distribution is simulated with a stochastic simulation algorithm, a socalled kinetic MonteCarlo method. This novel approach is applied for the simulation of an axisymmetric model of a tubular flow crystallizer. The numerical results are compared with experimental data. 
J. Fuhrmann, C. Guhlke, Ch. Merdon, A. Linke, R. Müller, Induced charge electroosmotic flow with finite ion size and solvation effects, Electrochimica Acta, 317 (2019), pp. 778785, DOI 10.1016/j.electacta.2019.05.051 .

A. Linke, L.G. Rebholz, Pressureinduced locking in mixed methods for timedependent (Navier)Stokes equations, Journal of Computational Physics, 388 (2019), pp. 350356, DOI 10.1016/j.jcp.2019.03.010 .
Abstract
We consider infsup stable mixed methods for the timedependent incompressible Stokes and NavierStokes equations, extending earlier work on the steady (Navier)Stokes Problem. A locking phenomenon is identified for classical infsup stable methods like the TaylorHood or the CrouzeixRaviart elements by a novel, elegant and simple numerical analysis and corresponding numerical experiments, whenever the momentum balance is dominated by forces of a gradient type. More precisely, a reduction of the L^{2} convergence order for high order methods, and even a complete stall of the L^{2} convergence order for lowestorder methods on preasymptotic meshes is predicted by the analysis and practically observed. On the other hand, it is also shown that (structurepreserving) pressurerobust mixed methods do not suffer from this locking phenomenon, even if they are of lowestorder. A connection to wellbalanced schemes for (vectorial) hyperbolic conservation laws like the shallow water or the compressible Euler equations is made. 
N. Alia, V. John, S. Ollila, Revisiting the singlephase flow model for liquid steel ladle stirred by gas, Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems. Elsevier Science Inc., New York, NY. English, English abstracts., 67 (2019), pp. 549556 (published online on 21.11.2018), DOI 10.1016/j.apm.2018.11.005 .

M. Akbas, A. Linke, L.G. Rebholz, P.W. Schroeder, The analogue of graddiv stabilization in DG methods for incompressible flows: Limiting behavior and extension to tensorproduct meshes, Computer Methods in Applied Mechanics and Engineering, 341 (2018), pp. 917938, DOI 10.1016/j.cma.2018.07.019 .
Abstract
Graddiv stabilization is a classical remedy in conforming mixed finite element methods for incompressible flow problems, for mitigating velocity errors that are sometimes called poor mass conservation. Such errors arise due to the relaxation of the divergence constraint in classical mixed methods, and are excited whenever the spacial discretization has to deal with comparably large and complicated pressures. In this contribution, an analogue of graddiv stabilization is presented for nonconforming flow discretizations of Discontinuous Galerkin or nonconforming finite element type. Here the key is the penalization of the jumps of the normal velocities over facets of the triangulation, which controls the measurevalued part of the distributional divergence of the discrete velocity solution. Furthermore, we characterize the limit for arbitrarily large penalization parameters, which shows that the proposed nonconforming Discontinuous Galerkin methods remain robust and accurate in this limit. Several numerical examples illustrate the theory and show their relevance for the simulation of practical, nontrivial flows. 
G.R. Barrenechea, V. John, P. Knobloch, R. Rankin, A unified analysis of algebraic flux correction schemes for convectiondiffusion equations, SeMA Journal. Boletin de la Sociedad Espannola de Matematica Aplicada, 75 (2018), pp. 655685, DOI 10.1007/s4032401801606 .

C. Bertoglio, A. Caiazzo, Y. Bazilevs, M. Braack, M. EsmailyMoghadam, V. Gravemeier, A.L. Marsden, O. Pironneau, I.E. VignonClementel, W.A. Wall, Benchmark problems for numerical treatment of backflow at open boundaries, International Journal of Numerical Methods in Biomedical Engineering, 34 (2018), pp. e2918/1e2918/34, DOI 10.1002/cnm.2918 .
Abstract
In computational fluid dynamics, incoming velocity at open boundaries, or backflow, often yields to unphysical instabilities already for moderate Reynolds numbers. Several treatments to overcome these backflow instabilities have been proposed in the literature. However, these approaches have not yet been compared in detail in terms of accuracy in different physiological regimes, in particular due to the difficulty to generate stable reference solutions apart from analytical forms. In this work, we present a set of benchmark problems in order to compare different methods in different backflow regimes (with a full reversal flow and with propagating vortices after a stenosis). The examples are implemented in FreeFem++ and the source code is openly available, making them a solid basis for future method developments. 
P.W. Schroeder, Ch. Lehrenfeld, A. Linke, G. Lube, Towards computable flows and robust estimates for infsup stable FEM applied to the timedependent incompressible NavierStokes equations, SeMA Journal. Boletin de la Sociedad Espannola de Matematica Aplicada, 75 (2018), pp. 629653, DOI 10.1007/s4032401801571 .
Abstract
Infsup stable FEM applied to timedependent incompressible NavierStokes flows are considered. The focus lies on robust estimates for the kinetic and dissipation energies in a twofold sense. Firstly, pressurerobustness ensures the fulfilment of a fundamental invariance principle and velocity error estimates are not corrupted by the pressure approximability. Secondly, Resemirobustness means that constants appearing on the righthand side of kinetic and dissipation energy error estimates (including Gronwall constants) do not explicitly depend on the Reynolds number. Such estimates rely on an essential regularity assumption for the gradient of the velocity, which is discussed in detail. In the sense of best practice, we review and establish pressure and Resemirobust estimates for pointwise divergencefree H1conforming FEM (like ScottVogelius pairs or certain isogeometric based FEM) and pointwise divergencefree H(div)conforming discontinuous Galerkin FEM. For convectiondominated problems, the latter naturally includes an upwind stabilisation for the velocity which is not gradientbased. 
J. DE Frutos, B. Garc'iaArchilla, V. John, J. Novo, Analysis of the graddiv stabilization for the timedependent NavierStokes equations with infsup stable finite elements, Advances in Computational Mathematics, 44 (2018), pp. 195225.

L. Blank, A. Caiazzo, F. Chouly, A. Lozinski, J. Mura, Analysis of a stabilized penaltyfree Nitsche method for the Brinkman, Stokes, and Darcy problems, ESAIM: Mathematical Modelling and Numerical Analysis, 52 (2018), pp. 21492185, DOI 10.1051/m2an/2018063 .

N. Ahmed, C. Bartsch, V. John, U. Wilbrandt, An assessment of solvers for some saddle point problems emerging from the incompressible NavierStokes equations, Computer Methods in Applied Mechanics and Engineering, 331 (2018), pp. 492513, DOI 10.1016/j.cma.2017.12.004 .
Abstract
Efficient incompressible flow simulations, using infsup stable pairs of finite element spaces, require the application of efficient solvers for the arising linear saddle point problems. This paper presents an assessment of different solvers: the sparse direct solver UMFPACK, the flexible GMRES (FGMRES) method with different coupled multigrid preconditioners, and FGMRES with Least Squares Commutator (LSC) preconditioners. The assessment is performed for steadystate and timedependent flows around cylinders in 2d and 3d. Several pairs of infsup stable finite element spaces with second order velocity and first order pressure are used. It turns out that for the steadystate problems often FGMRES with an appropriate multigrid preconditioner was the most efficient method on finer grids. For the timedependent problems, FGMRES with LSC preconditioners that use an inexact iterative solution of the velocity subproblem worked best for smaller time steps. 
N. Ahmed, V. John, G. Matthies, J. Novo, A local projection stabilization/continuous GalerkinPetrov method for incompressible flow problems, Applied Mathematics and Computation, 333 (2018), pp. 304324, DOI 10.1016/j.amc.2018.03.088 .
Abstract
The local projection stabilization (LPS) method in space is considered to approximate the evolutionary Oseen equations. Optimal error bounds independent of the viscosity parameter are obtained in the continuousintime case for the approximations of both velocity and pressure. In addition, the fully discrete case in combination with higher order continuous GalerkinPetrov (cGP) methods is studied. Error estimates of order k + 1 are proved, where k denotes the polynomial degree in time, assuming that the convective term is timeindependent. Numerical results show that the predicted order is also achieved in the general case of timedependent convective terms. 
N. Ahmed, A. Linke, Ch. Merdon, On really lockingfree mixed finite element methods for the transient incompressible Stokes equations, SIAM Journal on Numerical Analysis, 56 (2018), pp. 185209.
Abstract
Infsup stable mixed methods for the steady incompressible Stokes equations that relax the divergence constraint are often claimed to deliver lockingfree discretizations. However, this relaxation leads to a pressuredependent contribution in the velocity error, which is proportional to the inverse of the viscosity, thus giving rise to a (different) locking phenomenon. However, a recently proposed modification of the right hand side alone leads to a discretization that is really lockingfree, i.e., its velocity error converges with optimal order and is independent of the pressure and the smallness of the viscosity. In this contribution, we extend this approach to the transient incompressible Stokes equations, where besides the right hand side also the velocity time derivative requires an improved space discretization. Semidiscrete and fullydiscrete apriori velocity and pressure error estimates are derived, which show beautiful robustness properties. Two numerical examples illustrate the superior accuracy of pressurerobust space discretizations in the case of small viscosities. 
M. Hintermüller, M. Hinze, Ch. Kahle, T. Keil, A goaloriented dualweighted adaptive finite element approach for the optimal control of a nonsmooth CahnHilliardNavierStokes system, Optimization and Engineering. International Multidisciplinary Journal to Promote Optimization Theory & Applications in Engineering Sciences, 19 (2018), pp. 629662, DOI 10.1007/s1108101893936 .
Abstract
This paper is concerned with the development and implementation of an adaptive solution algorithm for the optimal control of a timediscrete CahnHilliardNavierStokes system with variable densities. The free energy density associated to the CahnHilliard system incorporates the doubleobstacle potential which yields an optimal control problem for a family of coupled systems in each time instant of a variational inequality of fourth order and the NavierStokes equation. A dualweighed residual approach for goaloriented adaptive finite elements is presented which is based on the concept of Cstationarity. The overall error representation depends on primal residual weighted by approximate dual quantities and vice versa as well as various complementary mismatch errors. Details on the numerical realization of the adaptive concept and a report on numerical tests are given. 
V. John, S. Kaya, J. Novo, Finite element error analysis of a mantle convection model, International Journal of Numerical Analysis and Modeling. Wuhan University, Wuhan and Institute for Scientific Computing and Information(ISCI), Edmonton, Alberta. English., 15 (2018), pp. 677698, DOI 10.20347/WIAS.PREPRINT.2403 .
Abstract
A mantle convection model consisting of the stationary Stokes equations and a timedependent convectiondiffusion equation for the temperature is studied. The Stokes problem is discretized with a conforming infsup stable pair of finite element spaces and the temperature equation is stabilized with the SUPG method. Finite element error estimates are derived which show the dependency of the error of the solution of one problem on the error of the solution of the other equation. The dependency of the error bounds on the coefficients of the problem is monitored. 
V. John, P. Knobloch, J. Novo, Finite elements for scalar convectiondominated equations and incompressible flow problems  A never ending story?, Computing and Visualization in Science, 19 (2018), pp. 4763.
Abstract
The contents of this paper is twofold. First, important recent results concerning finite element methods for convectiondominated problems and incompressible flow problems are described that illustrate the activities in these topics. Second, a number of, in our opinion, important problems in these fields are discussed. 
A. Linke, Ch. Merdon, M. Neilan, F. Neumann, Quasioptimality of a pressurerobust nonconforming finite element method for the Stokes problem, Mathematics of Computation, 87 (2018), pp. 15431566, DOI 10.1090/mcom/3344 .
Abstract
Nearly all classical infsup stable mixed finite element methods for the incompressible Stokes equations are not pressurerobust, i.e., the velocity error is dependent on the pressure. However, recent results show that pressurerobustness can be recovered by a nonstandard discretization of the right hand side alone. This variational crime introduces a consistency error in the method which can be estimated in a straightforward manner provided that the exact velocity solution is sufficiently smooth. The purpose of this paper is to analyze the pressurerobust scheme with low regularity. The numerical analysis applies divergencefree H¹conforming Stokes finite element methods as a theoretical tool. As an example, pressurerobust velocity and pressure apriori error estimates will be presented for the (first order) nonconforming CrouzeixRaviart element. A key feature in the analysis is the dependence of the errors on the Helmholtz projector of the right hand side data, and not on the entire data term. Numerical examples illustrate the theoretical results. 
P. Farrell, A. Linke, Uniform second order convergence of a complete flux scheme on unstructured 1D grids for a singularly perturbed advectiondiffusion equation and some multidimensional extensions, Journal of Scientific Computing, 72 (2017), pp. 373395, DOI 10.1007/s1091501703617 .
Abstract
The accurate and efficient discretization of singularly perturbed advectiondiffusion equations on arbitrary 2D and 3D domains remains an open problem. An interesting approach to tackle this problem is the complete flux scheme (CFS) proposed by G. D. Thiart and further investigated by J. ten Thije Boonkkamp. For the CFS, uniform second order convergence has been proven on structured grids. We extend a version of the CFS to unstructured grids for a steady singularly perturbed advectiondiffusion equation. By construction, the novel finite volume scheme is nodally exact in 1D for piecewise constant source terms. This property allows to use elegant continuous arguments in order to prove uniform second order convergence on unstructured onedimensional grids. Numerical results verify the predicted bounds and suggest that by aligning the finite volume grid along the velocity field uniform second order convergence can be obtained in higher space dimensions as well. 
A. Linke, M. Neilan, L.G. Rebholz, N. Wilson, A connection between coupled and penalty projection timestepping schemes with FE spatial discretization for the NavierStokes equations, Journal of Numerical Mathematics, 25 (2017), pp. 229248, DOI 10.1515/jnma20161024 .
Abstract
We prove that in finite element settings where the divergencefree subspace of the velocity space has optimal approximation properties, the solution of Chorin/Temam projection methods for NavierStokes equations equipped with graddiv stabilization with parameter gamma, converge to the associated coupled method solution with rate 1/gamma as gamma goes to infinity. We prove this first for backward Euler schemes, and then extend the results to BDF2 schemes, and finally to schemes with outflow boundary conditions. Several numerical experiments are given which verify the convergence rate, and show how using projection methods in this setting with large graddiv stabilization parameters can dramatically improve accuracy. 
U. Wilbrandt, C. Bartsch, N. Ahmed, N. Alia, F. Anker, L. Blank, A. Caiazzo, S. Ganesan, S. Giere, G. Matthies, R. Meesala, A. Shamim, J. Venkatensan, V. John, ParMooN  A modernized program package based on mapped finite elements, Computers & Mathematics with Applications. An International Journal, 74 (2017), pp. 7488, DOI 10.1016/j.camwa.2016.12.020 .

S. Giere, V. John, Towards physically admissible reducedorder solutions for convectiondiffusion problems, Applied Mathematics Letters, 73 (2017), pp. 7883, DOI 10.1016/j.aml.2017.03.022 .

H. Antil, M. Hintermüller, R.H. Nochetto, Th.M. Surowiec, D. Wegner, Finite horizon model predictive control of electrowetting on dielectric with pinning, Interfaces and Free Boundaries. Mathematical Modelling, Analysis and Computation, 19 (2017), pp. 130, DOI 10.4171/IFB/375 .

J. Bulling, V. John, P. Knobloch, Isogeometric analysis for flows around a cylinder, Applied Mathematics Letters, 63 (2017), pp. 6570.

P.L. Lederer, A. Linke, Ch. Merdon, J. Schöberl, Divergencefree reconstruction operators for pressurerobust Stokes discretizations with continuous pressure finite elements, SIAM Journal on Numerical Analysis, 55 (2017), pp. 12911314.
Abstract
Classical infsup stable mixed finite elements for the incompressible (Navier)Stokes equations are not pressurerobust, i.e., their velocity errors depend on the continuous pressure. However, a modification only in the right hand side of a Stokes discretization is able to reestablish pressurerobustness, as shown recently for several infsup stable Stokes elements with discontinuous discrete pressures. In this contribution, this idea is extended to low and high order TaylorHood and mini elements, which have continuous discrete pressures. For the modification of the right hand side a velocity reconstruction operator is constructed that maps discretely divergencefree test functions to exactly divergencefree ones. The reconstruction is based on local H(div)conforming flux equilibration on vertex patches, and fulfills certain orthogonality properties to provide consistency and optimal apriori error estimates. Numerical examples for the incompressible Stokes and NavierStokes equations confirm that the new pressurerobust TaylorHood and mini elements converge with optimal order and outperform signicantly the classical versions of those elements when the continuous pressure is comparably large. 
N. Ahmed, S. Becher, G. Matthies, Higherorder discontinuous Galerkin time stepping and local projection stabilization techniques for the transient Stokes problem, Computer Methods in Applied Mechanics and Engineering, 313 (2017), pp. 2852.
Abstract
We introduce and analyze discontinuous Galerkin time discretizations coupled with continuous finite element methods based on equalorder interpolation in space for velocity and pressure in transient Stokes problems. Spatial stability of the pressure is ensured by adding a stabilization term based on local projection. We present error estimates for the semidiscrete problem after discretization in space only and for the fully discrete problem. The fully discrete pressure shows an instability in the limit of small time step length. Numerical tests are presented which confirm our theoretical results including the pressure instability. 
N. Ahmed, T.Ch. Rebollo, V. John, S. Rubino, A review of variational multiscale methods for the simulation of turbulent incompressible flows, Archives of Computational Methods in Engineering. State of the Art Reviews, 24 (2017), pp. 115164.
Abstract
Various realizations of variational multiscale (VMS) methods for simulating turbulent incompressible flows have been proposed in the past fifteen years. All of these realizations obey the basic principles of VMS methods: They are based on the variational formulation of the incompressible NavierStokes equations and the scale separation is defined by projections. However, apart from these common basic features, the various VMS methods look quite different. In this review, the derivation of the different VMS methods is presented in some detail and their relation among each other and also to other discretizations is discussed. Another emphasis consists in giving an overview about known results from the numerical analysis of the VMS methods. A few results are presented in detail to highlight the used mathematical tools. Furthermore, the literature presenting numerical studies with the VMS methods is surveyed and the obtained results are summarized. 
N. Ahmed, T.Ch. Rebollo, V. John, S. Rubino, Analysis of a full spacetime discretization of the NavierStokes equations by a local projection stabilization method, IMA Journal of Numerical Analysis, 37 (2017), pp. 14371467, DOI 10.1093/imanum/drw048 .
Abstract
A finite element error analysis of a local projection stabilization (LPS) method for the timedependent NavierStokes equations is presented. The focus is on the highorder termbyterm stabilization method that has one level, in the sense that it is defined on a single mesh, and in which the projectionstabilized structure of standard LPS methods is replaced by an interpolationstabilized structure. The main contribution is on proving, theoretically and numerically, the optimal convergence order of the arising fully discrete scheme. In addition, the asymptotic energy balance is obtained for slightly smooth flows. Numerical studies support the analytical results and illustrate the potential of the method for the simulation of turbulent flows. Smooth unsteady flows are simulated with optimal order of accuracy. 
N. Ahmed, A. Linke, Ch. Merdon, Towards pressurerobust mixed methods for the incompressible NavierStokes equations, Computational Methods in Applied Mathematics, 18 (2018), pp. 353372 (published online on 18.11.2017), DOI 10.1515/cmam20170047 .
Abstract
In this contribution, classical mixed methods for the incompressible NavierStokes equations that relax the divergence constraint and are discretely infsup stable, are reviewed. Though the relaxation of the divergence constraint was claimed to be harmless since the beginning of the 1970ies, Poisson locking is just replaced by another more subtle kind of locking phenomenon, which is sometimes called poor mass conservation. Indeed, divergencefree mixed methods and classical mixed methods behave qualitatively in a different way: divergencefree mixed methods are pressurerobust, which means that, e.g., their velocity error is independent of the continuous pressure. The lack of pressurerobustness in classical mixed methods can be traced back to a consistency error of an appropriately defined discrete Helmholtz projector. Numerical analysis and numerical examples reveal that really lockingfree mixed methods must be discretely infsup stable and pressurerobust, simultaneously. Further, a recent discovery shows that lockingfree, pressurerobust mixed methods do not have to be divergencefree. Indeed, relaxing the divergence constraint in the velocity trial functions is harmless, if the relaxation of the divergence constraint in some velocity test functions is repaired, accordingly. 
N. Ahmed, On the graddiv stabilization for the steady Oseen and NavierStokes equations, Calcolo. A Quarterly on Numerical Analysis and Theory of Computation, 54 (2017), pp. 471501, DOI 10.1007/s100920160194z .
Abstract
This paper studies the parameter choice in the graddiv stabilization applied to the generalized problems of Oseen type. Stabilization parameters based on minimizing the H^{1}(Ω) error of the velocity are derived which do not depend on the viscosity parameter. For the proposed parameter choices, the H^{1}(Ω) error of the velocity is derived that shows a direct dependence on the viscosity parameter. Differences and common features to the situation for the Stokes equations are discussed. Numerical studies are presented which confirm the theoretical results. Moreover, for the Navier Stokes equations, numerical simulations were performed on a twodimensional ow past a circular cylinder. It turns out, for the MINI element, that the best results can be obtained without graddiv stabilization. 
V. John, A. Linke, Ch. Merdon, M. Neilan, L.G. Rebholz, On the divergence constraint in mixed finite element methods for incompressible flows, SIAM Review, 59 (2017), pp. 492544, DOI 10.1137/15M1047696 .
Abstract
The divergence constraint of the incompressible NavierStokes equations is revisited in the mixed finite element framework. While many stable and convergent mixed elements have been developed throughout the past four decades, most classical methods relax the divergence constraint and only enforce the condition discretely. As a result, these methods introduce a pressuredependent consistency error which can potentially pollute the computed velocity. These methods are not robust in the sense that a contribution from the righthand side, which influences only the pressure in the continuous equations, impacts both velocity and pressure in the discrete equations. This paper reviews the theory and practical implications of relaxing the divergence constraint. Several approaches for improving the discrete mass balance or even for computing divergencefree solutions will be discussed: graddiv stabilization, higher order mixed methods derived on the basis of an exact de Rham complex, $bH(mathrmdiv)$conforming finite elements, and mixed methods with an appropriate reconstruction of the test functions. Numerical examples illustrate both the potential effects of using nonrobust discretizations and the improvements obtained by utilizing pressurerobust discretizations. 
A. Linke, Ch. Merdon, W. Wollner, Optimal L2 velocity error estimate for a modified pressurerobust CrouzeixRaviart Stokes element, IMA Journal of Numerical Analysis, 37 (2017), pp. 354374, DOI 10.1093/imanum/drw019 .
Abstract
Recently, a novel approach for the robust discretization of the incompressible Stokes equations was proposed that slightly modifies the nonconforming CrouzeixRaviart element such that its velocity error becomes pressureindependent. The modification results in an O(h) consistency error that allows straightforward proofs for the optimal convergence of the discrete energy norm of the velocity and of the L2 norm of the pressure. However, though the optimal convergence of the velocity in the L2 norm was observed numerically, it appeared to be nontrivial to prove. In this contribution, this gap is closed. Moreover, the dependence of the error estimates on the discrete infsup constant is traced in detail, which shows that classical error estimates are extremely pessimistic on domains with large aspect ratios. Numerical experiments in 2D and 3D illustrate the theoretical findings. 
N. Rotundo, T.Y. Kim, W. Jiang, L. Heltai, E. Fried, Error estimates of Bspline based finiteelement method for the winddriven ocean circulation, Journal of Scientific Computing, 69 (2016), pp. 430459.
Abstract
We present the error analysis of a Bspline based finiteelement approximation of the streamfunction formulation of the large scale winddriven ocean circulation. In particular, we derive optimal error estimates for hrefinement using a Nitschetype variational formulations of the two simplied linear models of the stationary quasigeostrophic equations, namely the Stommel and StommelMunk models. Numerical results on rectangular and embedded geometries confirm the error analysis. 
A. Linke, G. Matthies, L. Tobiska, Robust arbitrary order mixed finite element methods for the incompressible Stokes equations with pressure independent velocity errors, ESAIM: Mathematical Modelling and Numerical Analysis, 50 (2016), pp. 289309.
Abstract
Standard mixed finite element methods for the incompressible NavierStokes equations that relax the divergence constraint are not robust against large irrotational forces in the momentum balance and the velocity error depends on the continuous pressure. This robustness issue can be completely cured by using divergencefree mixed finite elements which deliver pressureindependent velocity error estimates. However, the construction of H1conforming, divergencefree mixed finite element methods is rather difficult. Instead, we present a novel approach for the construction of arbitrary order mixed finite element methods which deliver pressureindependent velocity errors. The approach does not change the trial functions but replaces discretely divergencefree test functions in some operators of the weak formulation by divergencefree ones. This modification is applied to infsup stable conforming and nonconforming mixed finite element methods of arbitrary order in two and three dimensions. Optimal estimates for the incompressible Stokes equations are proved for the H1 and L2 errors of the velocity and the L2 error of the pressure. Moreover, both velocity errors are pressureindependent, demonstrating the improved robustness. Several numerical examples illustrate the results. 
C. Bertoglio, A. Caiazzo, A Stokesresidual backflow stabilization method applied to physiological flows, Journal of Computational Physics, 313 (2016), pp. 260278.
Abstract
In computational fluid dynamics incoming flow at open boundaries, or emphbackflow, often yields to unphysical instabilities for high Reynolds numbers. It is widely accepted that this is due to the incoming energy arising from the convection term, which cannot be empha priori controlled when the velocity field is unknown at the boundary. In order to improve the robustness of the numerical simulations, we propose a stabilized formulation based on a penalization of the residual of a weak Stokes problem on the open boundary, whose viscous part controls the incoming convective energy, while the inertial term contributes to the kinetic energy. We also present different strategies for the approximation of the boundary pressure gradient, which is needed for defining the stabilization term. The method has the advantage that it does not require neither artificial modifications or extensions of the computational domain. Moreover, it is consistent with the Womersley solution. We illustrate our approach on numerical examples  both academic and reallife  relevant to blood and respiratory flows. The results also show that the stabilization parameter can be reduced with the mesh size. 
P. Bringmann, C. Carstensen, Ch. Merdon, Guaranteed error control for the pseudostress approximation of the Stokes equations, Numerical Methods for Partial Differential Equations. An International Journal, 32 (2016), pp. 14111432.
Abstract
The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in $L^2$. Any standard mixed finite element function space can be utilized for this mixed formulation, e.g. the RaviartThomas discretization which is related to the CrouzeixRaviart nonconforming finite element scheme in the lowestorder case. The effective and guaranteed a posteriori error control for this nonconforming velocityoriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local infsup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy. 
A. Ern, D. Di Pietro, A. Linke, F. Schieweck, A discontinuous skeletal method for the viscositydependent Stokes problem, Computer Methods in Applied Mechanics and Engineering, 306 (2016), pp. 175195.
Abstract
We devise and analyze arbitraryorder nonconforming methods for the discretization of the viscositydependent Stokes equations on simplicial meshes. We keep track explicitly of the viscosity and aim at pressurerobust schemes that can deal with the practically relevant case of body forces with large curlfree part in a way that the discrete velocity error is not spoiled by large pressures. The method is inspired from the recent Hybrid HighOrder (HHO) methods for linear elasticity. After elimination of the auxiliary variables by static condensation, the linear system to be solved involves only discrete facebased velocities, which are polynomials of degree k >=0, and cellwise constant pressures. Our main result is a pressureindependent energyerror estimate on the velocity of order (k+1). The main ingredient to achieve pressureindependence is the use of a divergencepreserving velocity reconstruction operator in the discretization of the body forces. We also prove an L2pressure estimate of order (k+1) and an L2velocity estimate of order (k+2), the latter under elliptic regularity. The local mass and momentum conservation properties of the discretization are also established. Finally, two and threedimensional numerical results are presented to support the analysis. 
M. Khodayari, P. Reinsberg, A.A. AbdElLatif, Ch. Merdon, J. Fuhrmann, H. Baltruschat, Determining solubility and diffusivity by using a flow cell coupled to a mass spectrometer, ChemPhysChem, 17 (2016), pp. 16471655.

J. DE Frutos, B. Garc'iaArchilla, V. John, J. Novo, Graddiv stabilization for the evolutionary Oseen problem with infsup stable finite elements, Journal of Scientific Computing, 66 (2016), pp. 9911024.
Abstract
The approximation of the timedependent Oseen problem using infsup stable mixed finite elements in a Galerkin method with graddiv stabilization is studied. The main goal is to prove that adding a graddiv stabilization term to the Galerkin approximation has a stabilizing effect for small viscosity. Both the continuousintime and the fully discrete case (backward Euler method, the twostep BDF, and CrankNicolson schemes) are analyzed. In fact, error bounds are obtained that do not depend on the inverse of the viscosity in the case where the solution is sufficiently smooth. The bounds for the divergence of the velocity as well as for the pressure are optimal. The analysis is based on the use of a specific Stokes projection. Numerical studies support the analytical results. 
J. DE Frutos, V. John, J. Novo, Projection methods for incompressible flow problems with WENO finite difference schemes, Journal of Computational Physics, 309 (2016), pp. 368386.
Abstract
Weighted essentially nonoscillatory (WENO) finite difference schemes have been recommended in a competitive study of discretizations for scalar evolutionary convectiondiffusion equations [20]. This paper explores the applicability of these schemes for the simulation of incompressible flows. To this end, WENO schemes are used in several nonincremental and incremental projection methods for the incompressible NavierStokes equations. Velocity and pressure are discretized on the same grid. A pressure stabilization PetrovGalerkin (PSPG) type of stabilization is introduced in the incremental schemes to account for the violation of the discrete infsup condition. Algorithmic aspects of the proposed schemes are discussed. The schemes are studied on several examples with different features. It is shown that the WENO finite difference idea can be transferred to the simulation of incompressible flows. Some shortcomings of the methods, which are due to the splitting in projection schemes, become also obvious. 
J. Fuhrmann, A. Linke, Ch. Merdon, F. Neumann, T. Streckenbach, H. Baltruschat, M. Khodayari, Inverse modeling of thin layer flow cells for detection of solubility, transport and reaction coefficients from experimental data, Electrochimica Acta, 211 (2016), pp. 110.
Abstract
Thin layer flow cells are used in electrochemical research as experimental devices which allow to perform investigations of electrocatalytic surface reactions under controlled conditions using reasonably small electrolyte volumes. The paper introduces a general approach to simulate the complete cell using accurate numerical simulation of the coupled flow, transport and reaction processes in a flow cell. The approach is based on a mass conservative coupling of a divergencefree finite element method for fluid flow and a stable finite volume method for mass transport. It allows to perform stable and efficient forward simulations that comply with the physical bounds namely mass conservation and maximum principles for the involved species. In this context, several recent approaches to obtain divergencefree velocities from finite element simulations are discussed. In order to perform parameter identification, the forward simulation method is coupled to standard optimization tools. After an assessment of the inverse modeling approach using known realistic data, first results of the identification of solubility and transport data for O2 dissolved in organic electrolytes are presented. A plausibility study for a more complex situation with surface reactions concludes the paper and shows possible extensions of the scope of the presented numerical tools. 
V. John, K. Kaiser, J. Novo, Finite element methods for the incompressible Stokes equations with variable viscosity, ZAMM. Zeitschrift für Angewandte Mathematik und Mechanik, 96 (2016), pp. 205216.
Abstract
Finite element error estimates are derived for the incompressible Stokes equations with variable viscosity. The ratio of the supremum and the infimum of the viscosity appears in the error bounds. Numerical studies show that this ratio can be observed sometimes. However, often the numerical results show a weaker dependency on the viscosity. 
A. Linke, Ch. Merdon, On velocity errors due to irrotational forces in the NavierStokes momentum balance, Journal of Computational Physics, 313 (2016), pp. 654661.
Abstract
This contribution studies the influence of the pressure on the velocity error in finite element discretisations of the NavierStokes equations. Three simple benchmark problems that are all close to realworld applications convey that the pressure can be comparably large and is not to be underestimated. For widely used finite element methods like the TaylorHood finite element method, such relatively large pressures can lead to spurious oscillations and arbitrarily large errors in the velocity, even if the exact velocity is in the ansatz space. Only mixed finite element methods, whose velocity error is pressureindependent, like the ScottVogelius finite element method can avoid this influence. 
A. Linke, Ch. Merdon, Pressurerobustness and discrete Helmholtz projectors in mixed finite element methods for the incompressible NavierStokes equations, Computer Methods in Applied Mechanics and Engineering, 311 (2016), pp. 304326.
Abstract
Recently, it was understood how to repair a certain L2orthogonality of discretelydivergencefree vector fields and gradient fields such that the velocity error of infsup stable discretizations for the incompressible Stokes equations becomes pressureindependent. These new 'pressurerobust' Stokes discretizations deliver a small velocity error, whenever the continuous velocity field can be well approximated on a given grid. On the contrary, classical infsup stable Stokes discretizations can guarantee a small velocity error only, when both the velocity and the pressure field can be approximated well, simultaneously.
In this contribution, 'pressurerobustness' is extended to the timedependent NavierStokes equations. In particular, steady and timedependent potential flows are shown to build an entire class of benchmarks, where pressurerobust discretizations can outperform classical approaches significantly. Speedups will be explained by a new theoretical concept, the 'discrete Helmholtz projector' of an infsup stable discretization. Moreover, different discrete nonlinear convection terms are discussed, and skewsymmetric pressurerobust discretizations are proposed. 
G.R. Barrenechea, V. John, P. Knobloch, Some analytical results for an algebraic flux correction scheme for a steady convectiondiffusion equation in one dimension, IMA Journal of Numerical Analysis, 34 (2015), pp. 17291756.
Abstract
Algebraic flux correction schemes are nonlinear discretizations of convection dominated problems. In this work, a scheme from this class is studied for a steadystate convectiondiffusion equation in one dimension. It is proved that this scheme satisfies the discrete maximum principle. Also, as it is a nonlinear scheme, the solvability of the linear subproblems arising in a Picard iteration is studied, where positive and negative results are proved. Furthermore, the nonexistence of solutions for the nonlinear scheme is proved by means of counterexamples. Therefore, a modification of the method, which ensures the existence of a solution, is proposed. A weak version of the discrete maximum principle is proved for this modified method. 
CH. Brennecke, A. Linke, Ch. Merdon, J. Schöberl, Optimal and pressureindependent $L^2$ velocity error estimates for a modified CrouzeixRaviart Stokes element with BDM reconstructions, Journal of Computational Mathematics, 33 (2015), pp. 191208.
Abstract
Nearly all infsup stable mixed finite elements for the incompressible Stokes equations relax the divergence constraint. The price to pay is that a priori estimates for the velocity error become pressuredependent, while divergencefree mixed finite elements deliver pressureindependent estimates. A recently introduced new variational crime using lowestorder RaviartThomas velocity reconstructions delivers a much more robust modified CrouzeixRaviart element, obeying an optimal pressureindependent discrete H^{1} velocity estimate. Refining this approach, a more sophisticated variational crime employing the lowestorder BDM element is proposed, which also allows proving an optimal pressure independent L^{2 } velocity error. Numerical examples confirm the analysis and demonstrate the improved robustness in the NavierStokes case. 
V. John, J. Novo, Analysis of the pressure stabilized PetrovGalerkin (PSPG) method stabilization for the evolutionary Stokes equations avoiding timestep restrictions, SIAM Journal on Numerical Analysis, 53 (2015), pp. 10051031.
Abstract
Optimal error estimates for the pressure stabilized PetrovGalerkin (PSPG) method for the continuousintime discretization of the evolutionary Stokes equations are proved in the case of regular solutions. The main result is applicable to higher order finite elements. The error bounds for the pressure depend on the error of the pressure at the initial time. An approach is suggested for choosing the discrete initial velocity in such a way that this error is bounded. The “instability of the discrete pressure for small time steps”, which is reported in the literature, is discussed on the basis of the analytical results. Numerical studies confirm the theoretical results, showing in particular that this instability does not occur for the proposed initial condition. 
S. Giere, T. Iliescu, V. John, D. Wells, SUPG reduced order models for convectiondominated convectiondiffusionreaction equations, Computer Methods in Applied Mechanics and Engineering, 289 (2015), pp. 454474.
Abstract
This paper presents a StreamlineUpwind PetrovGalerkin (SUPG) reduced order model (ROM) based on Proper Orthogonal Decomposition (POD). This ROM is investigated theoretically and numerically for convectiondominated convectiondiffusionreaction equations. The SUPG finite element method was used on realistic meshes for computing the snapshots, leading to some noise in the POD data. Numerical analysis is used to propose the scaling of the stabilization parameter for the SUPGROM. Two approaches are used: One based on the underlying finite element discretization and the other one based on the POD truncation. The resulting SUPGROMs and the standard Galerkin ROM (GROM) are studied numerically. For many settings, the results obtained with the SUPGROMs are more accurate. Finally, one of the choices for the stabilization parameter is recommended. 
A. Caiazzo, R. Guibert, Y. Boudjemline, I.E. VignonClementel, Efficient blood flow simulations for the design of stented valve reducer in enlarged ventricular outflow tracts, Cardiovascular Engineering and Technology, 6 (2015), pp. 485500.
Abstract
Tetralogy of Fallot is a congenital heart disease characterized over time, after the initial repair, by the absence of a functioning pulmonary valve, which causes regurgitation, and by progressive enlargement of the right ventricle and pulmonary arteries. Due to this pathological anatomy, available transcatheter valves are usually too small to be deployed in the enlarged right ventricular outflow tracts (RVOT). To avoid surgical valve replacement, an alternative consists in implanting a reducer prior to or in combination with a transcatheter valve. We describe a computational model to study the effect of a stented valve RVOT reducer on the hemodynamics in enlarged ventricular outflow tracts. To this aim, blood flow in the right ventricular outflow tract is modeled via the incompressible NavierStokes equations coupled to a simplified valve model, numerically solved with a standard finite element method and with a reduced order model based on Proper Orthogonal Decomposition (POD). Numerical simulations are based on a patient geometry obtained from medical imaging and boundary conditions tuned according to measurements of inlet flow rates and pressures. Different geometrical models of the reducer are built, varying its length and/or diameter, and compared with the initial devicefree state. Simulations thus investigate multiple device configurations and describe the effect of geometry on hemodynamics. Forces exerted on the valve and on the reducer are monitored, varying with geometrical parameters. Results support the thesis that the reducer does not introduce significant pressure gradients, as was found in animal experiments. Finally, we demonstrate how computational complexity can be reduced with POD. 
A. Caiazzo, G. Montecinos, L.O. Müller, E.M. Haacke, E.F. Toro, Computational haemodynamics in stenotic internal jugular veins, Journal of Mathematical Biology, 70 (2015), pp. 745772.
Abstract
Stenosis in internal jugular veins (IJVs) are frequently associated to pathological venous circulation and insufficient cerebral blood drainage. In this work, we set up a computational framework to assess the relevance of IJV stenoses through numerical simulation, combining medical imaging, patientspecific data and a mathematical model for venous occlusions. Coupling a threedimensional (3D) description of blood flow in IJVs with a reduced onedimesional model (1D) for major intracranial veins, we are able to model different anatomical configurations, an aspect of importance to understand the impact of IJV stenosis in intracranial venous haemodynamics. We investigate several stenotic configurations in a physiologic patientspecific regime, quantifying the effect of the stenosis in terms of venous pressure increase and wall shear stress patterns. Simulation results are in qualitative agreement with reported pressure anomalies in pathological cases. Moreover, they demonstrate the potential of the proposed multiscale framework for individualbased studies and computeraided diagnosis. 
A. Caiazzo, I. RamisConde, Multiscale modeling of palisade formation in glioblastoma multiforme, Journal of Theoretical Biology, 383 (2015), pp. 145156.
Abstract
Palisades are characteristic tissue aberrations that arise in glioblastomas. Observation of palisades is considered as a clinical indicator of the transition from a noninvasive to an invasive tumour. In this article we propose a computational model to study the influence of genotypic and phenotypic heterogeneity in palisade formation. For this we produced three dimensional realistic simulations, based on a multiscale hybrid model, coupling the evolution of tumour cells and the oxygen diffusion in tissue, that depict the shape of palisades during its formation. Our results can be summarized as the following: (1) we show that cell heterogeneity is a crucial factor in palisade formation and tumour growth; (2) we present results that can explain the observed fact that recursive tumours are more malignant than primary tumours; and (3) the presented simulations can provide to clinicians and biologists for a better understanding of palisades 3D structure as well as glioblastomas growth dynamics 
A. Linke, Ch. Merdon, Guaranteed energy error estimators for a modified robust CrouzeixRaviart Stokes element, Journal of Scientific Computing, 64 (2015), pp. 541558.
Abstract
This paper provides guaranteed upper energy error bounds for a modified lowestorder nonconforming CrouzeixRaviart finite element method for the Stokes equations. The modification from [A. Linke 2014, On the role of the Helmholtzdecomposition in mixed methods for incompressible flows and a new variational crime] is based on the observation that only the divergencefree part of the righthand side should balance the vector Laplacian. The new method has optimal energy error estimates and can lead to errors that are smaller by several magnitudes, since the estimates are pressureindependent. An efficient a posteriori velocity error estimator for the modified method also should involve only the divergencefree part of the righthand side. Some designs to approximate the Helmholtz projector are compared and verified by numerical benchmark examples. They show that guaranteed error control for the modified method is possible and almost as sharp as for the unmodified method. 
A. Linke, On the role of the Helmholtz decomposition in mixed methods for incompressible flows and a new variational crime, Computer Methods in Applied Mechanics and Engineering, 268 (2014), pp. 782800.
Abstract
According to the Helmholtz decomposition, the irrotational parts of the momentum balance equations of the incompressible NavierStokes equations are balanced by the pressure gradient. Unfortunately, nearly all mixed methods for incompressible flows violate this fundamental property, resulting in the wellknown numerical instability of poor mass conservation. The origin of this problem is the lack of L2orthogonality between discretely divergencefree velocities and irrotational vector fields. In order to cure this, a new variational crime using divergencefree velocity reconstructions is proposed. Applying lowest order RaviartThomas velocity reconstructions to the nonconforming CrouzeixRaviart element allows to construct a cheap flow discretization for general 2d and 3d simplex meshes that possesses the same advantageous robustness properties like divergencefree flow solvers. In the Stokes case, optimal apriori error estimates for the velocity gradients and the pressure are derived. Moreover, the discrete velocity is independent of the continuous pressure. Several detailed linear and nonlinear numerical examples illustrate the theoretical findings. 
E. Jenkins, V. John, A. Linke, L.G. Rebholz, On the parameter choice in graddiv stabilization for the Stokes equations, Advances in Computational Mathematics, 40 (2014), pp. 491516.
Abstract
Graddiv stabilization has been proved to be a very useful tool in discretizations of incompressible flow problems. Standard error analysis for infsup stable conforming pairs of finite element spaces predicts that the stabilization parameter should be optimally chosen to be O(1). This paper revisits this choice for the Stokes equations on the basis of minimizing the $H^1$ error of the velocity and the $L^2$ error of the pressure. It turns out, by applying a refined error analysis, that the optimal parameter choice is more subtle than known so far in the literature. It depends on the used norm, the solution, the family of finite element spaces, and the type of mesh. Depending on the situation, the optimal stabilization parameter might range from being very small to very large. The analytic results are supported by numerical examples. 
C. Bertoglio, A. Caiazzo, A tangential regularization method for backflow stabilization in hemodynamics, Journal of Computational Physics, 261 (2014), pp. 162171.
Abstract
In computational simulations of fluid flows, instabilities at the Neumann boundaries may appear during backflow regime. It is widely accepted that this is due to the incoming energy at the boundary, coming from the convection term, which cannot be controlled when the velocity field is unknown. We propose a stabilized formulation based on a local regularization of the fluid velocity along the tangential directions on the Neumann boundaries. The stabilization term is proportional to the amount of backflow, and does not require any further assumption on the velocity profile. The perfomance of the method is assessed on a two and threedimensional Womersley flows, as well as considering a hemodynamic physiological regime in a patientspecific aortic geometry. 
R. Eymard, J. Fuhrmann, A. Linke, On MAC schemes on triangular Delaunay meshes, their convergence and application to coupled flow problems, Numerical Methods for Partial Differential Equations. An International Journal, 30 (2014), pp. 13971424.
Abstract
We study two classical generalized MAC schemes on unstructured triangular Delaunay meshes for the incompressible Stokes and NavierStokes equations and prove their convergence for the first time. These generalizations use the duality between Voronoi and triangles of Delaunay meshes, in order to construct two staggered discretization schemes. Both schemes are especially interesting, since compatible finite volume discretizations for coupled convectiondiffusion equations can be constructed which preserve discrete maximum principles. In the first scheme, called tangential velocity scheme, the pressures are defined at the vertices of the mesh, and the discrete velocities are tangential to the edges of the triangles. In the second scheme, called normal velocity scheme, the pressures are defined in the triangles, and the discrete velocities are normal to the edges of the triangles. For both schemes, we prove the convergence in $L^2$ for the velocities and the discrete rotations of the velocities for the Stokes and the NavierStokes problem. Further, for the normal velocity scheme, we also prove the strong convergence of the pressure in $L^2$. Linear and nonlinear numerical examples illustrate the theoretical predictions. 
J. DE Frutos, B. Garc'iaArchilla, V. John, J. Novo, An adaptive SUPG method for evolutionary convectiondiffusion equations, Computer Methods in Applied Mechanics and Engineering, 273 (2014), pp. 219237.
Abstract
An adaptive algorithm for the numerical simulation of timedependent convectiondiffusionreaction equations will be proposed and studied. The algorithm allows the use of the natural extension of any error estimator for the steadystate problem for controlling local refinement and coarsening. The main idea consists in considering the SUPG solution of the evolutionary problem as the SUPG solution of a particular steadystate convectiondiffusion problem with data depending on the computed solution. The application of the error estimator is based on a heuristic argument by considering a certain term to be of higher order. This argument is supported in the onedimensional case by numerical analysis. In the numerical studies, particularly the residualbased error estimator from [18] will be applied, which has proved to be robust in the SUPG norm. The effectivity of this error estimator will be studied and the numerical results (accuracy of the solution, fineness of the meshes) will be compared with results obtained by utilizing the adaptive algorithm proposed in [6]. 
A. Caiazzo, T. Iliescu, V. John, S. Schyschlowa, A numerical investigation of velocitypressure reduced order models for incompressible flows, Journal of Computational Physics, 259 (2014), pp. 598616.
Abstract
This report has two main goals. First, it numerically investigates three velocitypressure reduced order models (ROMs) for incompressible flows. The proper orthogonal decomposition (POD) is used to generate the modes. One method computes the ROM pressure solely based on the velocity POD modes, whereas the other two ROMs use pressure modes as well. To the best of the authors' knowledge, one of the latter methods is novel. The second goal is to numerically investigate the impact of the snapshot accuracy on the ROMs accuracy. Numerical studies are performed on a twodimensional laminar flow past a circular obstacle. It turns out that, both in terms of accuracy and efficiency, the two ROMs that utilize pressure modes are clearly superior to the ROM that uses only velocity modes. The numerical results also show a strong correlation of the accuracy of the snap shots with the accuracy of the ROMs. 
A. Caiazzo, V. John, U. Wilbrandt, On classical iterative subdomain methods for the StokesDarcy problem, Computer & Geosciences, 18 (2014), pp. 711728.
Abstract
Iterative subdomain methods for the StokesDarcy problem that use Robin boundary conditions on the interface are reviewed. Their common underlying structure and their main differences are identified. In particular, it is clarified that there are different updating strategies for the interface conditions. For small values of fluid viscosity and hydraulic permeability, which are relevant in applications from geosciences, it is shown in numerical studies that only one of these updating strategies leads to an efficient numerical method, if this strategy is used in combination with appropriate parameters in the Robin boundary conditions. In particular, it is observed that the values of appropriate parameters are larger than those proposed so far. Not only the size but also the ratio of appropriate Robin parameters depends on the coefficients of the problem. 
V. John, L. Schumacher, A study of isogeometric analysis for scalar convectiondiffusion equations, Applied Mathematics Letters, 27 (2014), pp. 4348.

R. Bordás, V. John, E. Schmeyer, D. Thévenin, Numerical methods for the simulation of a coalescencedriven droplet size distribution, Theoretical and Computational Fluid Dynamics. SpringerVerlag, Berlin., 27 (2013), pp. 253271.
Abstract
A droplet size distribution in a turbulent flow field is considered and modeled by means of a population balance system. This paper studies different numerical methods for the 4D population balance equation and their impact on an output of interest, the timespaceaveraged droplet size distribution at the outlet which is known from experiments. These methods include different interpolations of the experimental data at the inlet, various discretizations in time and space, and different schemes for computing the aggregation integrals. It will be shown that notable changes in the output of interest might occur. In addition, the efficiency of the studied methods is discussed. 
G.R. Barrenechea, V. John, P. Knobloch, A local projection stabilization finite element method with nonlinear crosswind diffusion for convectiondiffusionreaction equations, ESAIM: Mathematical Modelling and Numerical Analysis, 47 (2013), pp. 13351366.
Abstract
An extension of the local projection stabilization (LPS) finite element method for convectiondiffusionreaction equations is presented and analyzed, both in the steadystate and the transient setting. In addition to the standard LPS method, a nonlinear crosswind diffusion term is introduced that accounts for the reduction of spurious oscillations. The existence of a solution can be proved and, depending on the choice of the stabilization parameter, also its uniqueness. Error estimates are derived which are supported by numerical studies. These studies demonstrate also the reduction of the spurious oscillations. 
A. Bradji, J. Fuhrmann, Some abstract error estimates of a finite volume scheme for a nonstationary heat equation on general nonconforming multidimensional spatial meshes, Applications of Mathematics, 58 (2013), pp. 138.
Abstract
A general class of nonconforming meshes has been recently studied for stationary anisotropic heterogeneous diffusion problems by R. Eymard and coworkers. Thanks to these basic ideas developed for stationary problems, we derive a new discretization scheme in order to approximate the nonstationary heat problem. The unknowns of this scheme are the values at the centre of the control volumes, at some internal interfaces, and at the mesh points of the time discretization. Although the numerical scheme stems from the finite volume method, its formulation is based on the discrete version for the weak formulation defined for the heat problem. We derive error estimates for the solution in discrete norm, and an error estimate for an approximation of the gradient, in a general framework in which the discrete bilinear form is satisfying ellipticity. We prove in particular, that, when the discrete flux is calculated using a stabilized discrete gradient, the convergence order is h+k , where h (resp. k) is the mesh size of the spatial (resp. time) discretization. This estimate is valid under the regularity assumption that the exact solution is twice continuously differentiable in time and space. These error estimates are useful because they allow us to get error estimates for the approximations of the exact solution and its first derivatives. 
B. Cousins, S. Le Borne, A. Linke, Z. Wang, Efficient linear solvers for incompressible flow simulations using ScottVogelius finite elements, Numerical Methods for Partial Differential Equations. An International Journal, 29 (2013), pp. 12171237.
Abstract
Recent research has shown that in some practically relevant situations like multiphysics flows (Galvin et al., Comput Methods Appl Mech Eng, 2012) divergencefree mixed finite elements may have a significantly smaller discretization error than standard nondivergencefree mixed finite elements. To judge the overall performance of divergencefree mixed finite elements, we investigate linear solvers for the saddle point linear systems arising in ScottVogelius finite element implementations of the incompressible NavierStokes equations. We investigate both direct and iterative solver methods. Due to discontinuous pressure elements in the case of ScottVogelius (SV) elements, considerably more solver strategies seem to deliver promising results than in the case of standard mixed finite elements such as TaylorHood elements. For direct methods, we extend recent preliminary work using sparse banded solvers on the penalty method formulation to finer meshes and discuss extensions. For iterative methods, we test augmented Lagrangian and H LU preconditioners with GMRES, on both full and statically condensed systems. Several numerical experiments are provided that show these classes of solvers are well suited for use with SV elements and could deliver an interesting overall performance in several applications. 
V. John, J. Novo, A robust SUPG norm a posteriori error estimator for stationary convectiondiffusion equations, Computer Methods in Applied Mechanics and Engineering, 255 (2013), pp. 289305.
Abstract
A robust residualbased a posteriori estimator is proposed for the SUPG finite element method applied to stationary convectiondiffusionreaction equations. The error in the natural SUPG norm is estimated. The main concern of this paper is the consideration of the convectiondominated regime. A global upper bound and a local lower bound for the error are derived, where the global upper estimate relies on some hypotheses. Numerical studies demonstrate the robustness of the estimator and the fulfillment of the hypotheses. A comparison to other residualbased estimators with respect to the adaptive grid refinement is also provided. 
A. Linke, L. Rebholz, On a reduced sparsity stabilization of graddiv type for incompressible flow problems, Computer Methods in Applied Mechanics and Engineering, 261262 (2013), pp. 142153.
Abstract
We introduce a new operator for stabilizing error that arises from the weak enforcement of mass conservation in finite element simulations of incompressible flow problems. We show this new operator has a similar positive effect on velocity error as the wellknown and very successful graddiv stabilization operator, but the new operator is more attractive from an implementation standpoint because it yields a sparser block structure matrix. That is, while graddiv produces fully coupled block matrices (i.e. blockfull), the matrices arising from the new operator are blockupper triangular in two dimensions, and in three dimensions the 2,1 and 3,1 blocks are empty. Moreover, the diagonal blocks of the new operator's matrices are identical to those of graddiv. We provide error estimates and numerical examples for finite element simulations with the new operator, which reveals the significant improvement in accuracy it can provide. Solutions found using the new operator are also compared to those using usual graddiv stabilization, and in all cases, solutions are found to be very similar. 
R. Bordás, V. John, E. Schmeyer, D. Thévenin, Measurement and simulation of a droplet population in a turbulent flow field, Computers and Fluids. Pergamon Press, Oxford., 66 (2012), pp. 5262.
Abstract
The interaction of a disperse droplet population (spray) in a turbulent flow field is studied by combining wind tunnel experiments with simulations based on the model of a population balance system. The behavior of the droplets is modeled numerically by a population balance equation. Velocities of the air and of the droplets are determined by nonintrusive measurements. A direct discretization of the 4D equation for the droplet size distribution is used in the simulations. Important components of the numerical algorithm are a variational multiscale method for turbulence modeling, an upwind scheme for the 4D equation and a preprocessing approach to evaluate the aggretation integrals. The simulations of this system accurately predict the modifications of the droplet size distribution from the inlet to the outlet of the measurement section. Since the employed configuration is simple and considering that all measurement data are freely available thanks to an Internetbased repository, the considered experiment is proposed as a benchmark problem for the simulation of disperse twophase turbulent flows. 
R. Eymard, Th. Gallouët, R. Herbin, A. Linke, Finite volume schemes for the biharmonic problem on general meshes, Mathematics of Computation, 81 (2012), pp. 20192048.
Abstract
We propose a finite volume scheme for the approximation of a biharmonic problem with Dirichlet boundary conditions. We prove that the piecewise constant approximate solution converges to the exact solution, as well as the discrete approximate of the gradient and the discrete approximate of the Laplacian of the exact solution. These results are confirmed by numerical results. 
K. Galvin, A. Linke, L. Rebholz, N. Wilson, Stabilizing poor mass conservation in incompressible flow problems with large irrotational forcing and application to thermal convection, Computer Methods in Applied Mechanics and Engineering, 237240 (2012), pp. 166176.
Abstract
We consider the problem of poor mass conservation in mixed finite element algorithms for flow problems with large rotationfree forcing in the momentum equation. We provide analysis that suggests for such problems, obtaining accurate solutions necessitates either the use of pointwise divergencefree finite elements (such as ScottVogelius), or heavy graddiv stabilization of weakly divergencefree elements. The theory is demonstrated in numerical experiments for a benchmark natural convection problem, where large irrotational forcing occurs with high Rayleigh numbers. 
W. Hackbusch, V. John, A. Khachatryan, C. Suciu, A numerical method for the simulation of an aggregationdriven population balance system, International Journal for Numerical Methods in Fluids, 69 (2012), pp. 16461660.
Abstract
A population balance system which models the synthesis of urea is studied in this paper. The equations for the flow field, the mass and the energy balances are given in a threedimensional domain and the equation for the particle size distribution (PSD) in a fourdimensional domain. This problem is convectiondominated and aggregationdriven. Both features require the application of appropriate numerical methods. This paper presents a numerical approach for simulating the population balance system which is based on finite element schemes, a finite difference method and a modern method to evaluate convolution integrals that appear in the aggregation term. Two experiments are considered and the numerical results are compared with experimental data. Unknown parameters in the aggregation kernel have to be calibrated. For appropriately chosen parameters, good agreements are achieved of the experimental data and the numerical results computed with the proposed method. A detailed study of the computational results reveals the influence of different parts of the aggregation kernel. 
M. Augustin, A. Caiazzo, A. Fiebach, J. Fuhrmann, V. John, A. Linke, R. Umla, An assessment of discretizations for convectiondominated convectiondiffusion equations, Computer Methods in Applied Mechanics and Engineering, 200 (2011), pp. 33953409.
Abstract
The performance of several numerical schemes for discretizing convectiondominated convectiondiffusion equations will be investigated with respect to accuracy and efficiency. Accuracy is considered in measures which are of interest in applications. The study includes an exponentially fitted finite volume scheme, the StreamlineUpwind PetrovGalerkin (SUPG) finite element method, a spurious oscillations at layers diminishing (SOLD) finite element method, a finite element method with continuous interior penalty (CIP) stabilization, a discontinuous Galerkin (DG) finite element method, and a total variation diminishing finite element method (FEMTVD). A detailed assessment of the schemes based on the Hemker example will be presented. 
M. Case, V. Ervin, A. Linke, L. Rebholz, N. Wilson, Stable computing with an enhanced physics based scheme for the 3D NavierStokes equations, International Journal of Numerical Analysis and Modeling. Wuhan University, Wuhan and Institute for Scientific Computing and Information(ISCI), Edmonton, Alberta. English., 8 (2011), pp. 118136.
Abstract
We study extensions of an earlier developed energy and helicity preserving scheme for the 3D NavierStokes equations and apply them to a more general class of problems. The scheme is studied together with stabilizations of graddiv type in order to mitigate the effect of the Bernoulli pressure error on the velocity error. We prove stability, convergence, discuss conservation properties, and present numerical experiments that demonstrate the advantages of the scheme. 
M. Case, V. Ervin, A. Linke, L. Rebholz, A connection between ScottVogelius and graddiv stabilized TaylorHood FE approximations of the NavierStokes equations, SIAM Journal on Numerical Analysis, 49 (2011), pp. 14611481.
Abstract
This article studies two methods for obtaining excellent mass conservation in finite element computations of the NavierStokes equations using continuous velocity fields. Under mild restrictions, the ScottVogelius element pair has recently been shown to be infsup stable and have optimal approximation properties, while also providing pointwise mass conservation. We present herein the first numerical tests of this element pair for the time dependent NavierStokes equations. We also prove that, again under these mild restrictions, the limit of the graddiv stabilized TaylorHood solutions to the NavierStokes problem converges to the ScottVogelius solution as the stabilization parameter tends to infinity. That is, in this setting, we provide theoretical justification that choosing the parameter large does not destroy the solution. A limiting result is also proven for the general case. Numerical tests are provided which verify the theory, and show how both ScottVogelius and graddiv stabilized TaylorHood (with large stabilization parameter) elements can provide accurate results with excellent mass conservation for NavierStokes approximations. 
A. Caiazzo, D. Evans, J.L. Falcone, J. Hegewald, E. Lorenz, ET AL., A complex automata approach for instent restenosis: Twodimensional multiscale modeling and simulations, Journal of Computational Science, 2 (2011), pp. 917.

J. Fuhrmann, H. Langmach, A. Linke, A numerical method for mass conservative coupling between fluid flow and solute transport, Applied Numerical Mathematics. An IMACS Journal, 61 (2011), pp. 530553.
Abstract
We present a new coupled discretization approach for species transport in an incompressible fluid. The NavierStokes equations for the flow are discretized by the divergencefree ScottVogelius element on barycentrically refined meshes guaranteeing LBB stability. The convectiondiffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape. 
A.L. Bowers, B.R. Cousins, A. Linke, L.G. Rebholz, New connections between finite element formulations of the NavierStokes equations, Journal of Computational Physics, 229 (2010), pp. 90209025.
Abstract
We show the velocity solutions to the convective, skewsymmetric, and rotational Galerkin finite element formulations of the NavierStokes equations are identical if ScottVogelius elements are used, and thus all three formulations will the same pointwise divergence free solution velocity. A connection is then established between the formulations for graddiv stabilized TaylorHood elements: under mild restrictions, the formulations' velocity solutions converge to each other (and to the ScottVogelius solution) as the stabilization parameter tends to infinity. Thus the benefits of using ScottVogelius elements can be obtained with the less expensive TaylorHood elements, and moreover the benefits of all the formulations can be retained if the rotational formulation is used. Numerical examples are provided that confirm the theory. 
A. Glitzky, J.A. Griepentrog, Discrete SobolevPoincaré inequalities for Voronoi finite volume approximations, SIAM Journal on Numerical Analysis, 48 (2010), pp. 372391.
Abstract
We prove a discrete SobolevPoincare inequality for functions with arbitrary boundary values on Voronoi finite volume meshes. We use Sobolev's integral representation and estimate weakly singular integrals in the context of finite volumes. We establish the result for star shaped polyhedral domains and generalize it to the finite union of overlapping star shaped domains. In the appendix we prove a discrete Poincare inequality for space dimensions greater or equal to two. 
V. John, A. Kindl, O.C. Suciu, Finite element LES and VMS methods on tetrahedral meshes, Journal of Computational and Applied Mathematics, (2010), pp. 30953102.

V. John, A. Kindl, A variational multiscale method for turbulent flow simulation with adaptive large scale space, Journal of Computational Physics, 229 (2010), pp. 301312.

V. John, A. Kindl, Numerical studies of finite element variational multiscale methods for turbulent flow simulations, Computer Methods in Applied Mechanics and Engineering, 199 (2010), pp. 841852.

V. John, J. Rang, Adaptive time step control for the incompressible NavierStokes equations, Computer Methods in Applied Mechanics and Engineering, 199 (2010), pp. 514524.

V. John, M. Roland, On the impact of the scheme for solving the higherdimensional equation in coupled population balance systems, International Journal for Numerical Methods in Engineering, 82 (2010), pp. 14501474.

S. Bartels, M. Jensen, R. Müller, Discontinuous Galerkin finite element convergence for incompressible miscible displacement problems of low regularity, SIAM Journal on Numerical Analysis, 47 (2009), pp. 37203743.

P. Rudolph, M. Czupalla, N. Dropka, Ch. FrankRotsch, F.M. Kiessling, O. Klein, B. Lux, W. Miller, U. Rehse, O. Root, Crystal growth from melt in combined heatermagnet modules, Journal of the Korean Crystal Growth and Crystal Technology, 19 (2009), pp. 215222.

A. Glitzky, K. Gärtner, Energy estimates for continuous and discretized electroreactiondiffusion systems, Nonlinear Analysis. Theory, Methods & Applications. An International Multidisciplinary Journal. Series A: Theory and Methods, 70 (2009), pp. 788805.
Abstract
We consider electroreactiondiffusion systems consisting of continuity equations for a finite number of species coupled with a Poisson equation. We take into account heterostructures, anisotropic materials and rather general statistic relations.
We investigate thermodynamic equilibria and prove for solutions to the evolution system the monotone and exponential decay of the free energy to its equilibrium value. Here the essential idea is an estimate of the free energy by the dissipation rate which is proved indirectly.
The same properties are shown for an implicit time discretized version of the problem. Moreover, we provide a space discretized scheme for the electroreactiondiffusion system which is dissipative (the free energy decays monotonously). On a fixed grid we use for each species different Voronoi boxes which are defined with respect to the anisotropy matrix occurring in the flux term of this species. 
A. Glitzky, Energy estimates for electroreactiondiffusion systems with partly fast kinetics, Discrete and Continuous Dynamical Systems, 25 (2009), pp. 159174.
Abstract
We start from a basic model for the transport of charged species in heterostructures containing the mechanisms diffusion, drift and reactions in the domain and at its boundary. Considering limit cases of partly fast kinetics we derive reduced models. This reduction can be interpreted as some kind of projection scheme for the weak formulation of the basic electroreactiondiffusion system. We verify assertions concerning invariants and steady states and prove the monotone and exponential decay of the free energy along solutions to the reduced problem and to its fully implicit discretetime version by means of the results of the basic problem. Moreover we make a comparison of prolongated quantities with the solutions to the basic model. 
V. John, T. Mitkova, M. Roland, K. Sundmacher, L. Tobiska, A. Voigt, Simulations of population balance systems with one internal coordinate using finite element methods, Chemical Engineering Sciences, 64 (2009), pp. 733741.

O. Klein, Ch. Lechner, P.É. Druet, P. Philip, J. Sprekels, Ch. FrankRotsch, F.M. Kiessling, W. Miller, U. Rehse, P. Rudolph, Numerical simulations of the influence of a traveling magnetic field, generated by an internal heatermagnet module, on liquid encapsulated Czochralski crystal growth, Magnetohydrodynamics. Consultants Bureau, New York (US). Consultants Bureau, New York. Translation from: Magnitnaya Gidrodinamika., 45 (2009), pp. 557567.

A. Glitzky, Exponential decay of the free energy for discretized electroreactiondiffusion systems, Nonlinearity, 21 (2008), pp. 19892009.
Abstract
Our focus are electroreactiondiffusion systems consisting of continuity equations for a finite number of species coupled with a Poisson equation. We take into account heterostructures, anisotropic materials and rather general statistical relations. We introduce a discretization scheme (in space and fully implicit in time) using a fixed grid but for each species different Voronoi boxes which are defined with respect to the anisotropy matrix occurring in the flux term of this species. This scheme has the special property that it preserves the main features of the continuous systems, namely positivity, dissipativity and flux conservation. For the discretized electroreactiondiffusion system we investigate thermodynamic equilibria and prove for solutions to the evolution system the monotone and exponential decay of the free energy to its equilibrium value. The essential idea is an estimate of the free energy by the dissipation rate which is proved indirectly. 
O. Klein, Ch. Lechner, P.É. Druet, P. Philip, J. Sprekels, Ch. FrankRotsch, F.M. Kiessling, W. Miller, U. Rehse, P. Rudolph, Numerical simulation of Czochralski crystal growth under the influence of a traveling magnetic field generated by an internal heatermagnet module (HMM), Journal of Crystal Growth, 310 (2008), pp. 15231532.

H. Stephan, A. Khrabustovskyi, Positivity and time behavior of a linear reactiondiffusion system, nonlocal in space and time, Mathematical Methods in the Applied Sciences, 31 (2008), pp. 18091834.
Abstract
We consider a general linear reactiondiffusion system in three dimensions and time, containing diffusion (local interaction), jumps (nonlocal interaction) and memory effects. We prove a maximum principle, and positivity of the solution, and investigate its asymptotic behavior. Moreover, we give an explicite expression of the limit of the solution for large times. In order to obtain these results we use the following method: We construct a Riemannian manifold with complicated microstructure depending on a small parameter. We study the asymptotic behavior of the solution of a simple diffusion equation on this manifold as the small parameter tends to zero. It turns out that the homogenized system coincides with the original reactiondiffusion system what allows us to investigate its properties. 
CH. Lechner, O. Klein, P.É. Druet, Development of a software for the numerical simulation of VCz growth under the influence of a traveling magnetic field, Journal of Crystal Growth, 303 (2007), pp. 161164.

J. Geiser, O. Klein, P. Philip, Influence of anisotropic thermal conductivity in the apparatus insulation for sublimation growth of SiC: Numerical investigation of heat transfer, Crystal Growth & Design, 6 (2006), pp. 20212028.

J. Geiser, O. Klein, P. Philip, Transient numerical study of temperature gradients during sublimation growth of SiC: Dependence on apparatus design, Journal of Crystal Growth, 297 (2006), pp. 2032.

O. Klein, F. Luterotti, R. Rossi, Existence and asymptotic analysis of a phase field model for supercooling, Quarterly of Applied Mathematics, 64 (2006), pp. 291319.

P. Philip, O. Klein, Transient conductiveradiative heat transfer: Discrete existence and uniqueness for a finite volume scheme, Mathematical Models & Methods in Applied Sciences, 15 (2005), pp. 227258.

A. Glitzky, R. Hünlich, Global existence result for pair diffusion models, SIAM Journal on Mathematical Analysis, 36 (2005), pp. 12001225.

A. Glitzky, R. Hünlich, Stationary energy models for semiconductor devices with incompletely ionized impurities, ZAMM. Zeitschrift für Angewandte Mathematik und Mechanik, 85 (2005), pp. 778792.

J. Griepentrog, On the unique solvability of a nonlocal phase separation problem for multicomponent systems, Banach Center Publications, 66 (2004), pp. 153164.

H. Gajewski, I.V. Skrypnik, On unique solvability of nonlocal driftdiffusiontype problems, Nonlinear Analysis. Theory, Methods & Applications. An International Multidisciplinary Journal. Series A: Theory and Methods, 56 (2004), pp. 803830.

H. Gajewski, I.V. Skrypnik, To the uniqueness problem for nonlinear parabolic equations, Discrete and Continuous Dynamical Systems, 10 (2004), pp. 315336.

A. Glitzky, R. Hünlich, Stationary solutions of twodimensional heterogeneous energy models with multiple species, Banach Center Publications, 66 (2004), pp. 135151.

A. Glitzky, Electroreactiondiffusion systems with nonlocal constraints, Mathematische Nachrichten, 277 (2004), pp. 1446.

O. Klein, Asymptotic behaviour for a phasefield model with hysteresis in onedimensional thermoviscoplasticity, Applications of Mathematics, 49 (2004), pp. 309341.

H. Gajewski, K. Zacharias, On a nonlocal phase separation model, Journal of Mathematical Analysis and Applications, 286 (2003), pp. 1131.

G. Albinus, H. Gajewski, R. Hünlich, Thermodynamic design of energy models of semiconductor devices, Nonlinearity, 15 (2002), pp. 367383.

H. Gajewski, On a nonlocal model of nonisothermal phase separation, Advances in Mathematical Sciences and Applications, 12 (2002), pp. 569586.

A. Glitzky, R. Hünlich, Global properties of pair diffusion models, Advances in Mathematical Sciences and Applications, 11 (2001), pp. 293321.

A. Glitzky, R. Hünlich, Electroreactiondiffusion systems including cluster reactions of higher order, Mathematische Nachrichten, 216 (2000), pp. 95118.

J.A. Griepentrog, An application of the Implicit Function Theorem to an energy model of the semiconductor theory, ZAMM. Zeitschrift für Angewandte Mathematik und Mechanik, 79 (1999), pp. 4351.
Abstract
In this paper we deal with a mathematical model for the description of heat conduction and carrier transport in semiconductor heterostructures. We solve a coupled system of nonlinear elliptic differential equations consisting of the heat equation with Joule heating as a source, the Poisson equation for the electric field an driftdiffusion equations with temperature dependent coefficients describing the charge and current conservation, subject to general thermal and electrical boundary conditions. We prove the existence and uniqueness of Holder continuous weak solutions near thermodynamic equilibria points using the Implicit Function Theorem. To show the differentiability of maps corresponding to the weak formulation of the problem we use regularity results from the theory of nonsmooth linear elliptic boundary value problems in SobolevCampanato spaces.
Contributions to Collected Editions

A. Linke, Ch. Merdon, On the significance of pressurerobustness for the space discretization of incompressible high reynolds number flows, in: Proceedings ofFinite Volumes for Complex Applications IX  Methods, Theoretical Aspects, Examples  FVCA 9, Bergen, Norway, June 2020, R. Klöfkorn, E. Keilegavlen, A.F. Radu, J. Fuhrmann, eds., Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham et al., 2020, pp. 103112.

A. Linke, Ch. Merdon, Wellbalanced discretisation for the compressible Stokes problem by gradientrobustness, in: Proceedings ofFinite Volumes for Complex Applications IX  Methods, Theoretical Aspects, Examples  FVCA 9, Bergen, Norway, June 2020, R. Klöfkorn, E. Keilegavlen, A.F. Radu, J. Fuhrmann, eds., Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham et al., 2020, pp. 113121.

A. Jha, V. John, On basic iteration schemes for nonlinear AFC discretizations, in: Boundary and Interior Layers, Computational and Asymptotic Methods, BAIL 2018, G.N. Barrenechea, J. Mackenzie, eds., 135 of Lecture Notes in Computational Science and Engineering, Springer, Cham, 2020, pp. 113128, DOI https://doi.org/10.1007/9783030418007_7 .
Abstract
Algebraic flux correction (AFC) finite element discretizations of steadystate convectiondiffusionreaction equations lead to a nonlinear problem. This paper presents first steps of a systematic study of solvers for these problems. Two basic fixed point iterations and a formal Newton method are considered. It turns out that the fixed point iterations behave often quite differently. Using a sparse direct solver for the linear problems, one of them exploits the fact that only one matrix factorization is needed to become very efficient in the case of convergence. For the behavior of the formal Newton method, a clear picture is not yet obtained. 
J. Fuhrmann, C. Guhlke, A. Linke, Ch. Merdon, R. Müller, Models and numerical methods for electrolyte flows, in: Topics in Applied Analysis and Optimisation, M. Hintermüller, J.F. Rodrigues, eds., CIM Series in Mathematical Sciences, Springer Nature Switzerland AG, Cham, 2019, pp. 183209.

J. Fuhrmann, C. Guhlke, A. Linke, Ch. Merdon, R. Müller, Voronoi finite volumes and pressure robust finite elements for electrolyte models with finite ion sizes, in: Numerical Geometry, Grid Generation and Scientific Computing. Proceedings of the 9th International Conference, NUMGRID 2018 / Voronoi 150, V.A. Garanzha, L. Kamenski, H. Si, eds., 131 of Lecture Notes in Computational Science and Engineering, Springer Nature Switzerland AG, Cham, 2019, pp. 7383, DOI 10.1007/9783030234362 .

N. Ahmed, A. Linke, Ch. Merdon, Towards pressurerobust mixed methods for the incompressible NavierStokes equations, in: Finite Volumes for Complex Applications VIII  Methods and Theoretical Aspects, FVCA 8, Lille, France, June 2017, C. Cancès, P. Omnes, eds., 199 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing AG, Cham, 2017, pp. 351359.

P. Farrell, A. Linke, Uniform second order convergence of a complete flux scheme on nonuniform 1D grids, in: Finite Volumes for Complex Applications VIII  Methods and Theoretical Aspects, FVCA 8, Lille, France, June 2017, C. Cancès, P. Omnes, eds., 199 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing AG, Cham, 2017, pp. 303310.

S. Ganesan, V. John, G. Matthies, R. Meesala, A. Shamim, U. Wilbrandt, An object oriented parallel finite element scheme for computations of PDEs: Design and implementation, in: 2016 IEEE 23rd International Conference on High Performance Computing Workshops (PDF only), pp. 106115, DOI 10.1109/HiPCW.2016.19 .

CH. Brennecke, A. Linke, Ch. Merdon, J. Schöberl, Optimal and pressureindependent $L^2$ velocity error estimates for a modified CrouzeixRaviart element with BDM reconstructions, in: Finite Volumes for Complex Applications VII  Methods and Theoretical Aspects  FVCA 7, Berlin, June 2014, J. Fuhrmann, M. Ohlberger, Ch. Rohde, eds., 77 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham et al., 2014, pp. 159167.
Abstract
Nearly all infsup stable mixed finite elements for the incompressible Stokes equations relax the divergence constraint. The price to pay is that a priori estimates for the velocity error become pressuredependent, while divergencefree mixed finite elements deliver pressureindependent estimates. A recently introduced new variational crime using lowestorder RaviartThomas velocity reconstructions delivers a much more robust modified CrouzeixRaviart element, obeying an optimal pressureindependent discrete H^{1} velocity estimate. Refining this approach, a more sophisticated variational crime employing the lowestorder BDM element is proposed, which also allows proving an optimal pressureindependent L^{2 } velocity error. Numerical examples confirm the analysis and demonstrate the improved robustness in the NavierStokes case. 
J. Fuhrmann, A. Linke, Ch. Merdon, Coupling of fluid flow and solute transport using a divergencefree reconstruction of the CrouzeixRaviart element, in: Finite Volumes for Complex Applications VII  Elliptic, Parabolic and Hyperbolic Problems  FVCA 7, Berlin, June 2014, J. Fuhrmann, M. Ohlberger, Ch. Rohde, eds., 78 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham et al., 2014, pp. 587595.

A. Glitzky, J.A. Griepentrog, On discrete SobolevPoincaré inequalities for Voronoi finite volume approximations, in: Finite volumes for complex applications VI: Problems and perspectives, J. Fořt, J. Fürst, J. Halama, R. Herbin, F. Hubert, eds., Springer Proceedings in Mathematics 4, Springer, Heidelberg, 2011, pp. 533541.

S. Bartels, R. Müller, Die kalte Zunge, in: Besser als Mathe  Moderne angewandte Mathematik aus dem MATHEON zum Mitmachen, K. Biermann, M. Grötschel, B. LutzWestphal, eds., Reihe: Populär, Vieweg+Teubner, Wiesbaden, 2010, pp. 227235.

M. Jensen, R. Müller, Stable CrankNicolson discretisation for incompressible miscible displacement problems of low regularity, in: Numerical Mathematics and Advanced Applications 2009, Part 2, G. Kreiss, P. Lötstedt, A. Målqvist, M. Neytcheva, eds., Springer, Heidelberg et al., pp. 469477.
Abstract
In this article we study the numerical approximation of incompressible miscible displacement problems with a linearised CrankNicolson time discretisation, combined with a mixed finite element and discontinuous Galerkin method. At the heart of the analysis is the proof of convergence under low regularity requirements. Numerical experiments demonstrate that the proposed method exhibits secondorder convergence for smooth and robustness for rough problems. 
V. John, E. Schmeyer, On finite element methods for 3D timedependent convectiondiffusionreaction equations with small diffusion, in: BAIL 2008  Boundary and Interior Layers, A. Hegarty, N. Kopteva, E. O'Riordan, M. Stynes, eds., 69 of Lecture Notes in Computational Science and Engineering, Springer, Berlin/Heidelberg, 2009, pp. 173181.

H. Gajewski, H.Chr. Kaiser, H. Langmach, R. Nürnberg, R.H. Richter, Mathematical modelling and numerical simulation of semiconductor detectors, in: Mathematics  Key Technology for the Future. Joint Projects Between Universities and Industry, W. Jäger, H.J. Krebs, eds., Springer, Berlin [u.a.], 2003, pp. 355364.

R. Hünlich, G. Albinus, H. Gajewski, A. Glitzky, W. Röpke, J. Knopke, Modelling and simulation of power devices for highvoltage integrated circuits, in: Mathematics  Key Technology for the Future. Joint Projects Between Universities and Industry, W. Jäger, H.J. Krebs, eds., Springer, Berlin [u.a.], 2003, pp. 401412.

U. Bandelow, H. Gajewski, H.Chr. Kaiser, Modeling combined effects of carrier injection, photon dynamics and heating in Strained MultiQuantumWell Laser, in: Physics and Simulation of Optoelectronic Devices VIII, R.H. Binder, P. Blood, M. Osinski, eds., 3944 of Proceedings of SPIE, SPIE, Bellingham, WA, 2000, pp. 301310.

G. Schwarz, E. Schöll, R. Nürnberg, H. Gajewski, Simulation of current filamentation in an extended driftdiffusion model, in: EQUADIFF 99: International Conference on Differential Equations, Berlin 1999, B. Fiedler, K. Gröger, J. Sprekels, eds., 2, World Scientific, Singapore [u. a.], 2000, pp. 13341336.

H. Gajewski, K. Zacharias, On a reactiondiffusion system modelling chemotaxis, in: EQUADIFF 99: International Conference on Differential Equations, Berlin 1999, B. Fiedler, K. Gröger, J. Sprekels, eds., 2, World Scientific, Singapore [u. a.], 2000, pp. 10981103.

R. Hünlich, A. Glitzky, On energy estimates for electrodiffusion equations arising in semiconductor technology, in: Partial differential equations. Theory and numerical solution, W. Jäger, J. Nečas, O. John, K. Najzar, eds., 406 of Chapman & Hall Research Notes in Mathematics, Chapman & Hall, Boca Raton, FL, 2000, pp. 158174.

H.Chr. Kaiser, J. Rehberg, About some mathematical questions concerning the embedding of SchrödingerPoisson systems into the driftdiffusion model of semiconductor devices, in: EQUADIFF 99: International Conference on Differential Equations, Berlin 1999, B. Fiedler, K. Gröger, J. Sprekels, eds., 2, World Scientific, Singapore [u. a.], 2000, pp. 13281333.
Preprints, Reports, Technical Reports

D. Abdel, P. Vágner, J. Fuhrmann, P. Farrell, Modelling charge transport in perovskite solar cells: Potentialbased and limiting ion vacancy depletion, Preprint no. 2780, WIAS, Berlin, 2020.
Abstract, PDF (1067 kByte)
From MaxwellStefan diffusion and general electrostatics, we derive a driftdiffusion model for charge transport in perovskite solar cells (PSCs) where any ion in the perovskite layer may flexibly be chosen to be mobile or immobile. Unlike other models in the literature, our model is based on quasi Fermi potentials instead of densities. This allows to easily include nonlinear diffusion (based on for example FermiDirac, GaussFermi or Blakemore statistics) as well as limit the ion vacancy depletion (via the FermiDirac integral of order 1). The latter will be motivated by a grandcanonical formalism of ideal lattice gas. Furthermore, our model allows to use different statistics for different species. We discuss the thermodynamic equilibrium, electroneutrality as well as generation/recombination. Finally, we present numerical finite volume simulations to underline the importance of limiting ion vacancy depletion. 
U. Wilbrandt, N. Alia, V. John, Optimal control of a buoyancydriven liquid steel stirring modeled with singlephase NavierStokes equations, Preprint no. 2776, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2776 .
Abstract, PDF (1469 kByte)
Gas stirring is an important process used in secondary metallurgy. It allows to homogenize the temperature and the chemical composition of the liquid steel and to remove inclusions which can be detrimental for the endproduct quality. In this process, argon gas is injected from two nozzles at the bottom of the vessel and rises by buoyancy through the liquid steel thereby causing stirring, i.e., a mixing of the bath. The gas flow rates and the positions of the nozzles are two important control parameters in practice. A continuous optimization approach is pursued to find optimal values for these control variables. The effect of the gas appears as a volume force in the singlephase incompressible NavierStokes equations. Turbulence is modeled with the Smagorinsky Large Eddy Simulation (LES) model. An objective functional based on the vorticity is used to describe the mixing in the liquid bath. Optimized configurations are compared with a default one whose design is based on a setup from industrial practice. 
D. Frerichs, V. John, On reducing spurious oscillations in discontinuous Galerkin (DG) methods for steadystate convectiondiffusionreaction equations, Preprint no. 2769, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2769 .
PDF (1663 kByte) 
P.L. Lederer, Ch. Merdon, Guaranteed upper bounds for the velocity error of pressurerobust Stokes discretisations, Preprint no. 2750, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2750 .
Abstract, PDF (422 kByte)
This paper improves guaranteed error control for the Stokes problem with a focus on pressurerobustness, i.e. for discretisations that compute a discrete velocity that is independent of the exact pressure. A PragerSynge type result relates the errors of divergencefree primal and H(div)conforming dual mixed methods (for the velocity gradient) with an equilibration constraint that needs special care when discretised. To relax the constraints on the primal and dual method, a more general result is derived that enables the use of a recently developed mass conserving mixed stress discretisation to design equilibrated fluxes that yield pressureindependent guaranteed upper bounds for any pressurerobust (but not necessarily divergencefree) primal discretisation. Moreover, a provably efficient local design of the equilibrated fluxes is presented that reduces the numerical costs of the error estimator. All theoretical findings are verified by numerical examples which also show that the efficiency indices of our novel guaranteed upper bounds for the velocity error are close to 1. 
N. Ahmed, G.R. Barrenechea, E. Burman, J. Guzmán, A. Linke, Ch. Merdon, A pressurerobust discretization of Oseen's equation using stabilization in the vorticity equation, Preprint no. 2740, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2740 .
Abstract, PDF (843 kByte)
Discretization of NavierStokes' equations using pressurerobust finite element methods is considered for the high Reynolds number regime. To counter oscillations due to dominating convection we add a stabilization based on a bulk term in the form of a residualbased least squares stabilization of the vorticity equation supplemented by a penalty term on (certain components of) the gradient jump over the elements faces. Since the stabilization is based on the vorticity equation, it is independent of the pressure gradients, which makes it pressurerobust. Thus, we prove pressureindependent error estimates in the linearized case, known as Oseen's problem. In fact, we prove an O(h^{k}+^{1}/2) error estimate in the L^{2}norm that is known to be the best that can be expected for this type of problem. Numerical examples are provided that, in addition to confirming the theoretical results, show that the present method compares favorably to the classical residualbased SUPG stabilization. 
L. Baňas, R. Lasarzik, A. Prohl, Numerical analysis for nematic electrolytes, Preprint no. 2717, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2717 .
Abstract, PDF (8587 kByte)
We consider a system of nonlinear PDEs modeling nematic electrolytes, and construct a dissipative solution with the help of its implementable, structureinheriting spacetime discretization. Computational studies are performed to study the mutual effects of electric, elastic, and viscous effects onto the molecules in a nematic electrolyte. 
TH. Apel, V. Kempf, A. Linke, Ch. Merdon, A nonconforming pressurerobust finite element method for the Stokes equations on anisotropic meshes, Preprint no. 2702, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2702 .
Abstract, PDF (509 kByte)
Most classical finite element schemes for the (Navier)Stokes equations are neither pressurerobust, nor are they infsup stable on general anisotropic triangulations. A lack of pressurerobustness may lead to large velocity errors, whenever the Stokes momentum balance is dominated by a strong and complicated pressure gradient. It is a consequence of a method, which does not exactly satisfy the divergence constraint. However, infsup stable schemes can often be made pressurerobust just by a recent, modified discretization of the exterior forcing term, using H(div)conforming velocity reconstruction operators. This approach has so far only been analyzed on shaperegular triangulations. The novelty of the present contribution is that the reconstruction approach for the CrouzeixRaviart method, which has a stable Fortin operator on arbitrary meshes, is combined with results on the interpolation error on anisotropic elements for reconstruction operators of RaviartThomas and BrezziDouglasMarini type, generalizing the method to a large class of anisotropic triangulations. Numerical examples confirm the theoretical results in a 2D and a 3D test case. 
D. Frerichs, Ch. Merdon, Divergencepreserving reconstructions on polygons and a really pressurerobust virtual element method for the Stokes problem, Preprint no. 2683, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2683 .
Abstract, PDF (314 kByte)
Non divergencefree discretisations for the incompressible Stokes problem may suffer from a lack of pressurerobustness characterised by large discretisations errors due to irrotational forces in the momentum balance. This paper argues that also divergencefree virtual element methods (VEM) on polygonal meshes are not really pressurerobust as long as the righthand side is not discretised in a careful manner. To be able to evaluate the righthand side for the testfunctions, some explicit interpolation of the virtual testfunctions is needed that can be evaluated pointwise everywhere. The standard discretisation via an L^{2} bestapproximation does not preserve the divergence and so destroys the orthogonality between divergencefree testfunctions and possibly eminent gradient forces in the righthand side. To repair this orthogonality and restore pressurerobustness another divergencepreserving reconstruction is suggested based on RaviartThomas approximations on local subtriangulations of the polygons. All findings are proven theoretically and are demonstrated numerically in two dimensions. The construction is also interesting for hybrid highorder methods on polygonal or polyhedral meshes. 
L. Blank, E. Meneses Rioseco, U. Wilbrandt, A. Caiazzo, Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers, Preprint no. 2656, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2656 .
Abstract, PDF (20 MByte)
Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermohydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradientfree optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and investigate numerically the optimal positioning and spacing of multiwell systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multiwell arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geologocal configurations, and can be thus flexibly used in the context of multivariable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multiwell systems. 
C. Grässle, M. Hintermüller, M. Hinze, T. Keil, Simulation and control of a nonsmooth CahnHilliard NavierStokes system with variable fluid densities, Preprint no. 2617, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2617 .
Abstract, PDF (11 MByte)
We are concerned with the simulation and control of a two phase flow model governed by a coupled CahnHilliard NavierStokes system involving a nonsmooth energy potential.We establish the existence of optimal solutions and present two distinct approaches to derive suitable stationarity conditions for the bilevel problem, namely C and strong stationarity. Moreover, we demonstrate the numerical realization of these concepts at the hands of two adaptive solution algorithms relying on a specifically developed goaloriented error estimator.In addition, we present a model order reduction approach using proper orthogonal decomposition (PODMOR) in order to replace highfidelity models by low order surrogates. In particular, we combine POD with spaceadapted snapshots and address the challenges which are the consideration of snapshots with different spatial resolutions and the conservation of a solenoidal property. 
M. Hintermüller, T. Keil, Optimal control of geometric partial differential equations, Preprint no. 2612, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2612 .
Abstract, PDF (15 MByte)
Optimal control problems for geometric (evolutionary) partial differential inclusions are considered. The focus is on problems which, in addition to the nonlinearity due to geometric evolution, contain optimization theoretic challenges because of nonsmoothness. The latter might stem from energies containing nonsmooth constituents such as obstacletype potentials or terms modeling, e.g., pinning phenomena in microfluidics. Several techniques to remedy the resulting constraint degeneracy when deriving stationarity conditions are presented. A particular focus is on Yosidatype mollifications approximating the original degenerate problem by a sequence of nondegenerate nonconvex optimal control problems. This technique is also the starting point for the development of numerical solution schemes. In this context, also dualweighted residual based error estimates are addressed to facilitate an adaptive mesh refinement. Concerning the underlying state model, sharp and diffuse interface formulations are discussed. While the former always allows for accurately tracing interfacial motion, the latter model may be dictated by the underlying physical phenomenon, where near the interface mixed phases may exist, but it may also be used as an approximate model for (sharp) interface motion. In view of the latter, (sharp interface) limits of diffuse interface models are addressed. For the sake of presentation, this exposition confines itself to phase field type diffuse interface models and, moreover, develops the optimal control of either of the two interface models along model applications. More precisely, electrowetting on dielectric is used in the sharp interface context, and the control of multiphase fluids involving spinodal decomposition highlights the phase field technique. Mathematically, the former leads to a HeleShaw flow with geometric boundary conditions involving a complementarity system due to contact line pinning, and the latter gives rise to a CahnHilliard NavierStokes model including a nonsmooth obstacle type potential leading to a variational inequality constraint.
Talks, Poster

P. Farrell, D. Peschka, Challenges in driftdiffusion semiconductor simulations, Finite Volumes for Complex Applications IX (Online Event), Belgien, Norway, June 15  19, 2020.

J. Fuhrmann, D.H. Doan, A. Glitzky, M. Liero, G. Nika, Unipolar driftdiffusion simulation of Sshaped currentvoltage relations for organic semiconductor devices, Finite Volumes for Complex Applications IX (Online Event), Bergen, Norway, June 15  19, 2020.

V. John, On the provable convergence order for the kinetic energy of FEMs for the incompressible NavierStokes equations, 7th European Seminar on Computing (ESCO 2020) (Online Event), Pilsen, Czech Republic, June 8  12, 2020.

R. Lasarzik, Dissipative solutions in the context of the numerical approximation of nematic electrolytes (online talk), Oberseminar Numerik, Universität Bielefeld, Fakultät für Mathematik, June 23, 2020.

A. Linke, Ch. Merdon, On the significance of pressurerobustness for the space discretization of incompressible high reynolds number flows, Finite Volumes for Complex Applications IX (Online Event), Belgien, Norway, June 15  19, 2020.

H. Si, Adaptive exponential time integration of the NavierStokes equations, 3rd AIAA Sonic Boom Prediction Workshop, January 5  10, 2020, American Institute of Aeronautics and Astronautics SciTech Forum, Orlando, Florida, USA, January 10, 2020.

D. Peschka, Mathematical modeling and simulation of substrateflow interaction using generalized gradient flow, Begutachtungskolloquium für die Anträge des SPP 2171 ``Dynamische Benetzung flexibler, adaptiver und schaltbarer Oberflächen'', Mainz, February 7  8, 2019.

A. Caiazzo, Data assimilation in onedimensional hemodynamics, European Conference on Numerical Mathematics and Advanced Applications (ENUMATH 2019), Minisymposium 36 ``DataDriven Computational Fluid Dynamics (Part 2)'', September 30  October 4, 2019, Eindhoven University of Technology, Netherlands, October 1, 2019.

A. Caiazzo, Geothermal reservoir: Modeling, simulation and optimization for district heating in hot sedimentary acquires, Leibniz MMS Days 2019, March 20  22, 2019, Universität Rostock , LeibnizInstitut für Atmosphärenphysik, Kühlungsborn, March 22, 2019.

V. John, Algebraic finite element stabilizations for convectiondiffusion equations, Workshop on Computational Modeling and Numerical Analysis (WCMNA 2019), February 25  28, 2019, Laboratório Nacional de Computação Científica, Petrópolis, Brazil, February 26, 2019.

V. John, Algebraic finite element stabilizations for convectiondiffusion equations, Workshop ''Towards Computable Flows'', April 26  27, 2019, GeorgAugustUniversität Göttingen, Institut für Numerische und Angewandte Mathematik, April 26, 2019.

V. John, Finite element methods for incompressible flows, Workshop on Computational Modeling and Numerical Analysis (WCMNA 2019), February 25  28, 2019, Laboratório Nacional de Computação Científica, Petrópolis, Brazil.

V. John, Finite elements for scalar convectiondominated equations and incompressible flow problems  A never ending story?, Conference on Applied Mathematics, August 19  21, 2019, Lahore University of Management Sciences, Pakistan, August 19, 2019.

V. John, Finite elements for scalar convectiondominated equations and incompressible flow problems  A never ending story?, IndoGerman Conference on Computational Mathematics (IGCM), December 2  4, 2019, Indian Institute of Science, Department of Computer and Data Sciences, Bangalore, India, December 3, 2019.

V. John, On $L^2(Omega)$ estimates for finite element methods for evolutionary convectiondominated problems, PIMSGermany Workshop on Discretization of Variational Eigenvalue and Flow Problems, June 24  26, 2019, Universität Heidelberg, June 25, 2019.

V. John, Variational Multiscale (VMS) methods for the simulation of turbulent incompressible flows, University of Groningen, Bernoulli Institute, Computational Mechanics & Numerical Mathematics, Netherlands, September 23, 2019.

V. John, Variational Multiscale (VMS) methods for the simulation of turbulent incompressible flows, Indian Institute of Science, Department of Computational and Data Science, Bangalore, India, November 28, 2019.

A. Linke, On pressurerobustness, coherent structures and vortexdominated flows, Leibniz MMS Days 2019, March 20  22, 2019, Universität Rostock , LeibnizInstitut für Atmosphärenphysik, Kühlungsborn, March 22, 2019.

N. Alia, Modeling and optimization of a gasstirred liquid flow for steelmaking processes, The 20th European Conference on Mathematics for Industry (ECMI), MS27: MSO for steel production and manufacturing, June 18  22, 2018, University Budapest, Institute of Mathematics at Eötvös Loránd, Hungary, June 19, 2018.

R. Ahrens, F. Anker, C. Bartsch, A. Voigt, V. Wiedmeyer, K. Sundmacher, V. John, S. Le Borne, Advanced numerical methods for the simulation of population balance systems, 6th International Conference on Population Balance Modelling (PBM2018), Ghent University, Belgium, May 7  9, 2018.

L. Blank, A robust finite element method for the Brinkman problem, 13th International Workshop on Variational Multiscale and Stabilized Finite Elemements, WeierstraßInstitut, Berlin, December 5, 2018.

L. Blank, An unconditionally stable, low order, and robust finite element method for the numerical simulation of porous media flow, 39th Northern German Colloquium on Applied Analysis and Numerical Mathematics (NoKo 2018), June 1  2, 2018, Technische Universität Braunschweig, June 2, 2018.

C. Bartsch, V. John, R.I.A. Patterson, A new mixed stochasticdeterministic simulation approach to particle populations in fluid flows, 6th International Conference on Population Balance Modelling (PBM2018), Ghent University, Belgium, May 7  9, 2018.

A. Caiazzo, A penaltyfree Nitsche method for the Stokes, Darcy and Brinkman problems, Universität Augsburg, Lehrstuhl für Numerische Mathematik, May 15, 2018.

A. Caiazzo, Mathematical modeling and simulations of geothermal reservoirs, 16th European Finite Element Fair (EFEF 2018), June 8  9, 2018, Universität Heidelberg,, June 8, 2018.

A. Caiazzo, Mathematical modeling and simulations of geothermal reservoirs, LeibnizInstitut für Angewandte Geophysik, Hannover, November 7, 2018.

A. Caiazzo, Multiscale and reducedorder modeling of biphasic materials with application to tissue elastography, Rheinische FriedrichWilhelmsUniversität Bonn, Institut für Numerische Simulation, November 23, 2018.

A. Caiazzo, Multiscale modeling of individualbased cancer models with arbitrary vasculature and fiber structure, 6th European Conference on Computational Mechanics, 7th European Conference on Computational Fluid Dynamics (ECCMECFD 2018), June 11  15, 2018, University of Glasgow, UK, June 12, 2018.

A. Caiazzo, Robust open boundary conditions and efficient data assimilation in multiscale hemodynamics, International Symposium ``Modeling, Simulation and Optimization of the Cardiovascular System'', October 22  24, 2018, Universität Magdeburg, October 22, 2018.

A. Caiazzo, Towards the personalization of (1D) bloodflow simulations, University of Amsterdam, Computational Science Lab, Netherlands, September 21, 2018.

V. John, A new mixed stochasticdeterministic simulation approach for particle populations in fluid flows, 6th European Seminar on Computing (ESCO), June 3  8, 2018, University of West Bohemia, Pilsen, Czech Republic, June 6, 2018.

V. John, Angewandte Mathematik für Strömungssimulationen, MartinLutherUniversität Halle, Fachbereich Mathematik, October 27, 2018.

V. John, Finite elements for scalar convectiondominated equations and incompressible flow problems  A never ending story?, International Conference of Boundary and Interior Layers (BAIL 2018), June 18  22, 2018, University of Strathclyde, Glasgow, Scotland, UK, June 18, 2018.

V. John, Variational Multiscale (VMS) methods for the simulation of turbulent incompressible flows, University of Strathclyde, Department of Mathematics and Statistics, Glasgow, UK, November 14, 2018.

A. Linke, On highorder pressurerobust space discretisations, their advantages for incompressible high Reynolds number generalised Beltrami flows and beyond, 13th International Workshop on Variational Multiscale and Stabilized Finite Elemements, WeierstraßInstitut, Berlin, December 5, 2018.

A. Linke, On the role of the HelmholtzHodge projector for a novel pressurerobust discretization theory for the incompressible NavierStokes equations, The University of Texas at Austin, Institute for Computational Engineering and Scienses, USA, September 18, 2018.

A. Linke, On the role of the HelmholtzLeray projector for a novel pressurerobust discretization theory for the incompressible NavierStokes equations, Workshop ``Finite Element Exterior Calculus (FEEC) and High Order Methods'', June 4  6, 2018, University of Oslo, Faculty of Mathematics and Natural Sciences, Norway, June 6, 2018.

A. Linke, On the role of the HelmholtzLeray projector for a novel pressurerobust discretization theory for the incompressible NavierStokes equations, GeorgAugustUniversität Göttingen, Institut für Numerische und Angewandte Mathematik, January 23, 2018.

A. Linke, On the role of the HelmholtzLeray projector for a novel pressurerobust discretization theory for the incompressible NavierStokes equations, Aix Marseille Université, Institut de Mathématiques, France, April 3, 2018.

A. Linke, On the role of the HelmholtzLeray projector for a novel pressurerobust discretization theory for the incompressible NavierStokes equations, Clemson University, Department of Mathematical Sciences, South Carolina, USA, August 17, 2018.

A. Linke, On the role of the HelmholtzLeray projector in the space discretization of the NavierStokes equations, 3rd Leibniz MMS Days 2018, February 28  March 2, 2018, Wissenschaftszentrum Leipzig, March 1, 2018.

A. Linke, Towards pressurerobust mixed methods for the incompressible NavierStokes equations, 89th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2018), March 19  23, 2018, Technische Universität München, March 22, 2018.

A. Linke, User talk: Pressurerobust space discretization for the incompressible NavierStokes equations, 2nd NGSolve User Meeting, July 4  6, 2018, GeorgAugustUniversität Göttingen, Institut für Numerische und Angewandte Mathematik, July 5, 2018.

U. Wilbrandt, Iterative subdomain methods for the StokesDarcy coupling, 6th European Conference on Computational Mechanics, 7th European Conference on Computational Fluid Dynamics (ECCMECFD 2018), June 11  15, 2018, University of Glasgow, UK, June 11, 2018.

T. Keil, Simulation and control of a nonsmooth CahnHilliard NavierStokes system with variable fluid densities (with Carmen Graessle), Annual Meeting of the DFG Priority Programme 1962, October 9  11, 2017, Kremmen (Sommerfeld), October 11, 2017.

T. Keil, Strong stationarity conditions for the optimal control of a CahnHilliardNavierStokes system, 14th International Conference on Free Boundary Problems: Theory and Applications, Theme Session 8 ``Optimization and Control of Interfaces'', July 9  14, 2017, Shanghai Jiao Tong University, China, July 10, 2017.

CH. Merdon, A novel concept for the discretisation of the coupled NernstPlanckPoissonNavierStokes system, 14th Symposium on Fuel Cell Modelling and Experimental Validation (MODVAL 14), March 2  3, 2017, Karlsruher Institut für Technologie, Institut für Angewandte Materialien, Karlsruhe, Germany, March 3, 2017.

CH. Merdon, Druckrobuste FiniteElementeMethoden für die NavierStokesGleichungen, Universität Paderborn, Institut für Mathematik, April 25, 2017.

CH. Merdon, Pressurerobustness in mixed finite element discretisations for the NavierStokes equations, Universität des Saarlandes, Fakultät für Mathematik und Informatik, July 12, 2017.

D. Peschka, Modelling and simulation of suspension flow, Graduate Seminar PDE in the Sciences, Universität Bonn, Institut für Angewandte Mathematik, January 20, 2017.

N. Kumar, J.H.M. Ten Thije Boonkkamp, B. Koren, A. Linke, A nonlinear flux approximation scheme for the viscous Burgers equation, 8th International Symposium on Finite Volumes for Complex Applications (FVCA 8), Université Lille 1, Villeneuve d'Ascq, France, June 12  16, 2017.

N. Ahmed, Higherorder discontinuous Galerkin time discretizations for the evolutionary NavierStokes equations, Technische Universität Dresden, Institut für Numerische Mathematik, March 9, 2017.

N. Ahmed, On really lockingfree mixed finite element methods for the transient incompressible Stokes equations, CASM International Conference on Applied Mathematics, May 22  24, 2017, Lahore University of Management Sciences, Centre for Advanced Studies in Mathematics, Pakistan, May 22, 2017.

C. Bartsch, ParMooN  A parallel finite element solver, Part I, Indian Institute of Science, Supercomputer Education and Research Centre, Bangalore, India, March 16, 2017.

M. Hintermüller, Optimal control of multiphase fluids based on non smooth models, 14th International Conference on Free Boundary Problems: Theory and Applications, Theme Session 8 ``Optimization and Control of Interfaces'', July 9  14, 2017, Shanghai Jiao Tong University, China, July 10, 2017.

M. Hintermüller, Total variation diminishing RungeKutta methods for the optimal control of conservation laws: Stability and orderconditions, SIAM Conference on Optimization, Minisymposium MS111 ``Optimization with Balance Laws on Graphs'', May 22  25, 2017, Vancouver, British Columbia, Canada, May 25, 2017.

V. John, Analytical and numerical results for algebraic flux correction schemes, 12th International Workshop on Variational Multiscale and Stabilization Methods (VMS2017), April 26  28, 2017, Edificio Celestino Mutis, Campus Reina Mercedes, Sevilla, Spain, April 26, 2017.

V. John, Finite element methods for incompressible flow problems, May 14  18, 2017, Beijing Computational Science Research Center, Applied and Computational Mathematics, China.

V. John, Finite elements for scalar convectiondominated equations and incompressible flow problems  A never ending story?, 30th Chemnitz FEM Symposium, September 25  27, 2017, Bundesinstitut für Erwachsenenbildung, St. Wolfgang / Strobl, Austria, September 27, 2017.

V. John, Variational multiscale (VMS) methods for the simulation of turbulent incompressible flows, CDS: Computational Science Symposium, March 16  18, 2017, Indian Institute of Science, Department of Computer and Data Sciences, Bangalore, India, March 16, 2017.

V. John, Variational multiscale (VMS) methods for the simulation of turbulent incompressible flows, Mahindra École Centrale, School of Natural Sciences, Hyderabad, India, March 9, 2017.

V. John, Variational multiscale (VMS) methods for the simulation of turbulent incompressible flows, Chinese Academy of Sciences, Academy of Mathematics and Systems Science, Beijing, May 10, 2017.

V. John, Variational multiscale (VMS) methods for the simulation of turbulent incompressible flows, Peking University, School of Mathematical Sciences, Beijing, China, May 11, 2017.

A. Linke, Towards pressurerobust mixed methods for the incompressible NavierStokes equations, 12th International Workshop on Variational Multiscale and Stabilization Methods (VMS2017), April 26  28, 2017, Edificio Celestino Mutis, Campus Reina Mercedes, Sevilla, Spain, April 27, 2017.

CH. Merdon, Pressurerobust finite element methods for the NavierStokes equations, GAMM Workshop on Numerical Analysis, November 1  2, 2017, RheinischWestfälische Technische Hochschule Aachen, November 2, 2017.

CH. Merdon, Pressurerobust mixed finite element methods for the NavierStokes equations, scMatheon Workshop RMMM 8  Berlin 2017, Reliable Methods of Mathematical Modeling, July 31  August 3, 2017, HumboldtUniversität zu Berlin, August 2, 2017.

U. Wilbrandt, ParMooN  A parallel finite element solver, Part II, Indian Institute of Science, Supercomputer Education and Research Centre, Bangalore, India, March 16, 2017.

N. Alia, V. John, Optimal control of ladle stirring, 1st Leibniz MMS Mini Workshop on CFD & GFD, WIAS Berlin, September 8  9, 2016.

CH. Merdon, J. Fuhrmann, A. Linke, A.A. AbdElLatif, M. Khodayari, P. Reinsberg, H. Baltruschat, Inverse modelling of thin layer flow cells and RRDEs, The 67th Annual Meeting of the International Society of Electrochemistry, Den Haag, Netherlands, August 21  26, 2016.

D. Peschka, A free boundary problem for the flow of viscous liquid bilayers, ERC Workshop on Modeling Materials and Fluids using Variational Methods, February 22  26, 2016, WIAS Berlin, Berlin, February 26, 2016.

D. Peschka, Multiphase flows with contact lines: Solid vs liquid substrates, Industrial and Applied Mathematics Seminar, University of Oxford, Mathematical Institute, UK, October 27, 2016.

N. Ahmed, A review of VMS methods for the simulation of turbulent incompressible flows, International Conference on Differential Equations and Applications, May 26  28, 2016, Lahore University of Management Sciences, Pakistan, May 27, 2016.

N. Ahmed, On the graddiv stabilization for the steady Oseen and NavierStokes evaluations, International Conference of Boundary and Interior Layers (BAIL 2016), August 15  19, 2016, Beijing Computational Science Research Center, Beijing, China, August 15, 2016.

V. John, A review of variational multiscale methods for the simulation of turbulent incompressible flows, 15th Mathematics of Finite Elements and Applications, June 14  17, 2016, Brunel University London, London, UK, June 17, 2016.

V. John, A review of variational multiscale methods for the simulation of turbulent incompressible flows, ``Variational Multiscale and Stabilization Methods'' (VMS2016), March 16  18, 2016, OttovonGuericke Universität Magdeburg, Magdeburg, March 17, 2016.

V. John, A review of variational multiscale methods for the simulation of turbulent incompressible flows, International Conference of Boundary and Interior Layers (BAIL 2016), August 15  19, 2016, Beijing Computational Science Research Center, Beijing, China, August 15, 2016.

V. John, A review of variational multiscale methods for the simulation of turbulent incompressible flows, 7th American Mathematical Society Meeting, Special Session on Above and Beyond Fluid Flow Studies, October 8  9, 2016, Denver, Colorado, USA, October 9, 2016.

V. John, A survey on the analysis and numerical analysis of some turbulence models, Technische Universität Darmstadt, Fachbereich Mathematik, January 20, 2016.

V. John, Ein weites Feld  Wissenschaftliche Beiträge von Prof. Dr. Lutz Tobiska, Festkolloquium aus Anlass des 65. Geburtstags von Prof. Dr. Lutz Tobiska, Universität Magdeburg, Institut für Analysis und Numerik, March 31, 2016.

V. John, On the divergence constraint in mixed finite element methods for incompressible flows, 5th European Seminar on Computing (ESCO 2016), June 5  10, 2016, Pilsen, Czech Republic, June 7, 2016.

V. John, On the divergence constraint in mixed finite element methods for incompressible flows, Beijing Computational Science Research Center, China, August 23, 2016.

V. John, The role of the pressure in finite element methods for incompressible flow problems, Summer School 2016 ``Fluids under Pressure'' and Workshop, August 29  September 2, 2016, Nečas Center for Mathematical Modeling, Prague, Czech Republic.

CH. Merdon, Inverse modeling of thin layer flow cells for detection of solubility transport and reaction coefficients from experimental data, 17th Topical Meeting of the International Society of Electrochemistry Multiscale Analysis of Electrochemical Systems, May 31  June 3, 2015, Saint Malo Congress Center, France, June 1, 2015.

N. Ahmed, Adaptive time step control with variational time stepping schemes for convectiondiffusionreaction equations, 26th Biennial Numerical Analysis Conference, June 23  26, 2015, University of Strathclyde, Glasgow, UK, June 23, 2015.

N. Ahmed, Higher order time stepping schemes, 13th European Finite Element Fair, June 5  6, 2015, Charles University in Prague, Praha, Czech Republic, June 5, 2015.

A. Caiazzo, Assessment of Kalman filtering for parameter identification in onedimensional blood flow model, 4th International Conference on Computational & Mathematical Biomedical Engineering, June 29  July 1, 2015, Ecole Normale Supérieure de Cachan, Cachan, France, June 29, 2015.

V. John, S. Rockel, S. Sobolev, Viscoelastic mantle flow  Numerical modeling of geological phenomena, 4th Annual GEOSIM Workshop, Potsdam, November 17  18, 2015.

U. Wilbrandt, Classical iterative subdomain methods for the StokesDarcy problem, RheinischWestfälische Technische Hochschule Aachen, Institut für Geometrie und Praktische Mathematik, Aachen, June 11, 2015.

A. Linke, On the role of the Helmholtz decomposition in mixed methods for incompressible flows, International Conference of Boundary and Interior Layers (BAIL 2014), September 15  19, 2014, Carles University Prague, Czech Republic, September 17, 2014.

A. Linke, Ch. Merdon, Optimal and pressureindependent $L^2$ velocity error estimates for a modified CrouzeixRaviart element with BDM reconstructions, The International Symposium of Finite Volumes for Complex Applications VII (FVCA 7), BerlinBrandenburgische Akademie der Wissenschaften, June 15  20, 2014.

A. Linke, Ein neues Konstruktionsprinzip zur divergenzfreien Diskretisierung der inkompressiblen NavierStokesGleichungen, RuhrUniversität Bochum, Fakultät für Mathematik, July 17, 2014.

A. Linke, On the role of the Helmholtz decomposition in incompressible flows and a new variational crime, NonLinear PDE and Applications: Theoretical and Numerical Study, May 5  7, 2014, Abdelmalek Essadi University, Tanger, Morocco, May 6, 2014.

A. Linke, On the role of the Helmholtz decomposition in mixed methods for incompressible flows and a new variational crime, Technische Universität Wien, Institut für Analysis und Scientific Computing, Austria, April 2, 2014.

A. Linke, On the role of the Helmholtz decomposition in mixed methods for incompressible flows and a new variational crime, Technische Universität HamburgHarburg, Institut für Mathematik, January 7, 2014.

A. Linke, On the role of the Helmholtz decomposition in mixed methods for incompressible flows and a new variational crime, FriedrichAlexanderUniversität ErlangenNürnberg, Fachbereich Mathematik, November 20, 2014.

A. Linke, On the role of the Helmholtz decomposition in mixed methods for incompressible flows and a new variational crime, GeorgAugustUniversität Göttingen, Institut für Numerische und Angewandte Mathematik, December 9, 2014.

V. John, On isogeometric analysis for convectiondiffusion and NavierStokes equations, International Conference of Boundary and Interior Layers (BAIL 2014), September 15  19, 2014, Carles University Prague, Czech Republic, September 15, 2014.

J. Fuhrmann, A. Linke, Ch. Merdon, M. Khodayari, H. Baltruschat, Detection of solubility, transport and reaction coefficients from experimental data by inverse modeling of thin layer flow cells, 65th Annual Meeting of the International Society of Electrochemistry, Lausanne, Switzerland, August 31  September 5, 2014.

J. Fuhrmann, A. Linke, Ch. Merdon, Coupling of fluid flow and solute transport using a divergencefree reconstruction of the CrouzeixRaviart element, The International Symposium of Finite Volumes for Complex Applications VII (FVCA 7), BerlinBrandenburgische Akademie der Wissenschaften, June 15  20, 2014.

J. Fuhrmann, Activity based finite volume methods for generalised NernstPlanckPoisson systems, The International Symposium of Finite Volumes for Complex Applications VII (FVCA 7), BerlinBrandenburgische Akademie der Wissenschaften, June 15  20, 2014.

J. Fuhrmann, Voronoi finite volume methods, XVI. MathematicaTag, Berlin, December 9, 2014.

V. John, S. Rockel, S. Sobolev, Viscoelastic mantle flow, 3rd Annual GEOSIM Workshop, Potsdam, November 13  14, 2014.

V. John, Finite element methods for incompressible flow problems, Indian National Program on Differential Equations  Theory, Computation, and Applications (NPDETCA) Advanced Workshop on Finite Element Methods for NavierStokes Equations, September 8  12, 2014, Indian Institute of Science, Numerical Mathematics & Scientific Computing, Bangalore.

V. John, Numerical methods for convectiondominated equations, 2nd Workshop ``Populationsbilanzen'' of the DFG Priority Program SPP 1679 ``DynSimFP  Dynamic Simulation of Interconnected Solids Processes'', May 6  7, 2014, Technische Universität HamburgHarburg, May 6, 2014.

V. John, On the simulation of mantle convection, Symposium on Simulation and Optimization of Extreme Fluids, November 10  12, 2014, Internationales Wissenschaftsforum Heidelberg, November 12, 2014.

V. John, Turbulent flows and their numerical simulation, Humboldt Kolleg on Interdisciplinary Science: Catalyst for Sustainable Progress, September 4  6, 2014, Indian Institute of Science, Numerical Mathematics & Scientific Computing, Bangalore, September 5, 2014.

U. Wilbrandt, Iterative subdomain methods for StokesDarcy problems, Norddeutsches Kolloquium über Angewandte Analysis und Numerische Mathematik (NoKo), May 9  10, 2014, ChristianAlbrechtsUniversität zu Kiel, May 9, 2014.

V. John, On recent analytical and numerical investigations of the SUPG method, Advances in Computational Mechanics (ACM 2013)  A Conference Celebrating the 70th Birthday of Thomas J.R. Hughes, February 24  27, 2013, San Diego, USA, February 27, 2013.

A. Linke, Coupled flows and poor mass conservation, Workshop ``Complex grids and fluid flows, conclusion of VFSitCom, National Research Project'', April 2  4, 2012, RhôneAlpes, Lyon, France, April 3, 2012.

V. John, On reduced order modeling methods for incompressible flows based proper orthogonal decomposition, 6th Variational Multiscale Methods Workshop (VMS 2012), June 27  29, 2012, ChristianAlbrechtsUniversität zu Kiel, June 28, 2012.

V. John, On the analysis and numerical analysis of some turbulence models, Workshop ``Connections Between Regularized and LargeEddy Simulation Methods for Turbulence'', May 14  17, 2012, Banff International Research Station for Mathematical Innovation and Discovery, Canada, May 15, 2012.

A. Glitzky, J.A. Griepentrog, Discrete SobolevPoincaré inequalities for Voronoi finite volume approximations, Finite Volumes for Complex Applications VI (FVCA 6), Prague, Czech Republic, June 6  10, 2011.

O.C. Suciu, A numerical method for the simulation of population balance systems with one internal coordinate, 8th International Conference on LargeScale Scientific Computations, June 6  10, 2011, Bulgarian Academy of Sciences, Sozopol, June 7, 2011.

A. Caiazzo, J. Fuhrmann, V. John, Finite elementfinite volume coupling for simulations of thin porous layer fuel cells, Workshop on Simulation of Flow in Porous Media and Applications in Waste Management and CO2 Sequestration, John Radon Institute for Computational and Applied Mathematics, Linz, Austria, October 3  7, 2011.

A. Caiazzo, DarcyStokes coupling for simulation of solute transport in thinlayer channel fuel cells, Workshop on Fluid Dynamics and Porous Media, September 12  14, 2011, University of Coimbra, Department of Mathematics, Portugal, September 13, 2011.

A. Caiazzo, Model reduction approaches for simulation of cardiovascular stents and pulmonary valve, Laboratory of Modeling and Scientific Computing, Department of Mathematics, Milan, Italy, October 27, 2011.

A. Caiazzo, Physical and mathematicalbased reduced order modeling in computational hemodynamics, Wroclaw University of Technology, Institute of Mathematics and Computer Science, Poland, November 16, 2011.

E. Schmeyer, V. John, R. Bordás, D. Thévenin, Referenzexperimente im Windkanal, numerische Simulation und Validierung, Annual Meeting 2010 of DFG Priority Program 1276 MetStröm, Zuse Institute Berlin, October 28  29, 2010.

O.C. Suciu, On numerical methods for simulation of 3D/4D of population balance system, SimParTurS Meeting, MaxPlanckInstitut Leipzig, February 17, 2010.

V. John, A posteriori optimization of parameters in stabilized methods for convectiondiffusion problems, BIRS Workshop ``Nonstandard Discretizations for Fluid Flows'', November 21  26, 2010, Banff International Research Station, Canada.

V. John, Error analysis of the SUPG finite element discretization of evolutionary convectiondiffusionreaction equations, 6th Variational Multiscale Methods Workshop (VMS 2010), May 27  28, 2010, Université de Pau et des Pays de l'Adour, France, May 27, 2010.

V. John, Error analysis of the SUPG finite element discretization of evolutionary convectiondiffusionreaction equations, International Conference of Boundary and Interior Layers (BAIL 2010), July 5  9, 2010, University of Zaragoza, Spain, July 8, 2010.

V. John, On the numerical simulation of population balance systems, Technische Universität Dresden, Fakultät Mathematik und Naturwissenschaften, April 20, 2010.

V. John, Variational multiscale methods for the simulation of turbulent flows, Volkswagen Wolfsburg, August 27, 2010.

M. Ehrhardt, J. Fuhrmann, A. Linke, Finite volume methods for the simulation of flow cell experiments, Workshop ``New Trends in Model Coupling  Theory, Numerics & Applications'' (NTMC'09), Paris, France, September 2  4, 2009.

M. Ehrhardt, A high order finite element method for waves in periodic structures, 9th International Conference on Spectral and High Order Methods (ICOSAHOM09), Minisymposium ``Highorder Methods for Linear and Nonlinear Wave Equations'', June 22  26, 2009, Norwegian University of Science and Technology, Trondheim, June 24, 2009.

V. John, A variational multiscale method for turbulent flow simulation with adaptive large scale space, Workshop on Computational Multiscale Modeling, November 26  27, 2009, University of Twente, Enschede, Netherlands, November 27, 2009.

V. John, On the numerical simulation of population balance systems, Karlsruher Institut für Technologie, Fakultät für Mathematik, December 9, 2009.

H.Chr. Kaiser, Transient KohnSham theory, Jubiläumssymposium ``Licht  Materialien  Modelle'' (100 Jahre Innovation aus Adlershof), BerlinAdlershof, September 7  8, 2009.

K. Hoke, Numerical treatment of the KohnSham system for semiconductor devices, Workshop on Mathematical Aspects of Transport in Mesoscopic Systems, Dublin, Ireland, December 4  7, 2008.

E. Bänsch, H. Berninger, U. Böhm, A. Bronstert, M. Ehrhardt, R. Forster, J. Fuhrmann, R. Klein, R. Kornhuber, A. Linke, A. Owinoh, J. Volkholz, Pakt für Forschung und Innovation: Das Forschungsnetzwerk ``Gekoppelte Strömungsprozesse in Energie und Umweltforschung'', Show of the Leibniz Association ``Exzellenz durch Vernetzung. Kooperationsprojekte der deutschen Wissenschaftsorganisationen mit Hochschulen im Pakt für Forschung und Innovation'', Berlin, November 12, 2008.

M. Ehrhardt, O. Gloger, Th. Dietrich, O. Hellwich, K. Graf, E. Nagel, Level Set Methoden zur Segmentierung von kardiologischen MRBildern, 22. Treffpunkt Medizintechnik: Fortschritte in der medizinischen Bildgebung, Charité, Campus Virchow Klinikum Berlin, May 22, 2008.

A. Glitzky, Analysis of spinpolarized driftdiffusion models, 79th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2008), session ``Applied Analysis'', March 31  April 4, 2008, University of Bremen, April 1, 2008.

A. Glitzky, Energy estimates for continuous and discretized reactiondiffusion systems in heterostructures, Annual Meeting of the Deutsche MathematikerVereinigung 2008, minisymposium ``Analysis of ReactionDiffusion Systems with Internal Interfaces'', September 15  19, 2008, FriedrichAlexanderUniversität ErlangenNürnberg, September 15, 2008.

A. Glitzky, Energy estimates for space and time discretized electroreactiondiffusion systems, Conference on Differential Equations and Applications to Mathematical Biology, June 23  27, 2008, Université Le Havre, France, June 26, 2008.

H.Chr. Kaiser, A driftdiffusion model for semiconductors with internal interfaces, Annual Meeting of the Deutsche MathematikerVereinigung 2008, Minisymposium ``Analysis of ReactionDiffusion Systems with Internal Interfaces'', September 15  19, 2008, FriedrichAlexanderUniversität ErlangenNürnberg, September 15, 2008.

H.Chr. Kaiser, A thermodynamic approach to transient KohnSham theory, 100th Statistical Mechanics Conference, December 13  18, 2008, Rutgers, The State University of New Jersey, New Brunswick, USA, December 16, 2008.

H.Chr. Kaiser, On driftdiffusion KohnSham theory, 79th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2008), session ``Applied Analysis'', March 31  April 4, 2008, University of Bremen, April 1, 2008.

A. Glitzky, Energy estimates for reactiondiffusion processes of charged species, 6th International Congress on Industrial and Applied Mathematics (ICIAM), July 16  20, 2007, ETH Zürich, Switzerland, July 16, 2007.

H.Chr. Kaiser, A driftdiffusion model of transient KohnSham theory, First Joint International Meeting between the American Mathematical Society and the Polish Mathematical Society, Special Session ``Mathematics of Large Quantum Systems'', July 31  August 3, 2007, University of Warsaw, Poland, August 3, 2007.

K. Gärtner, A. Glitzky, Th. Koprucki, Analysis and simulation of spinpolarized driftdiffusion models, Evaluation Colloquium of the DFG Priority Program SPP 1285 ``Semiconductor Spintronics'', Bad Honnef, December 14  15, 2006.

A. Glitzky, Energy estimates for electroreactiondiffusion systems with partly fast kinetics, 6th AIMS International Conference on Dynamical Systems, Differential Equations & Applications, June 25  28, 2006, Université de Poitiers, France, June 27, 2006.

A. Glitzky, An application of the Implicit Function Theorem to stationary energy models for semiconductor devices, International Workshop ``Regularity for nonlinear and linear PDEs in nonsmooth domains'', September 4  7, 2005, Universität Stuttgart, Hirschegg, Austria, September 5, 2005.

A. Glitzky, Stationary energy models for semicoductor devices with incompletely ionized impurities, 2nd Joint Meeting of AMS, DMV, ÖMG, June 16  19, 2005, Johannes Gutenberg Universität, Mainz, June 19, 2005.

H. Gajewski, R. Hünlich, H.Chr. Kaiser, M. Baro, Quantum mechanical and macroscopic models for optoelectronic devices, DFG Research Center sc Matheon, Technische Universität Berlin, July 19, 2004.

H. Gajewski, Zur Numerik des Ladungsträgertransports in Halbleiterbauelementen, Technische Universität München, Institut fär Technische Elektrophysik, February 5, 2004.

M. Baro, H. Gajewski, R. Hünlich, H.Chr. Kaiser, Optoelektronische Bauelemente: mikroskopische & makroskopische Modelle, MathInside  Überall ist Mathematik, event of the DFG Research Center ``Mathematics for Key Technologies'' on the occasion of the Open Day of Urania, Berlin, September 13, 2003  December 3, 2004.

A. Glitzky, R. Hünlich, Stationary solutions of twodimensional heterogeneous energy models with multiple species, Nonlocal Elliptic and Parabolic Problems, September 9  11, 2003, Bedlewo, Poland, September 10, 2003.

H.Chr. Kaiser, Classical solutions of van Roosbroeck's equations with discontinuous coefficients and mixed boundary conditions on twodimensional space domains, 19th GAMM Seminar Leipzig on Highdimensional problems  Numerical treatment and applications, January 23  25, 2003, MaxPlanckInstitut für Mathematik in den Naturwissenschaften, Leipzig, January 25, 2003.
External Preprints

J. Fuhrmann, C. Guhlke, Ch. Merdon, A. Linke, R. Müller, Induced charge electroosmotic flow with finite ion size and solvation effects, Preprint no. arXiv:1901.06941, Cornell University Library, 2019, DOI 10.1016/j.electacta.2019.05.051 .

A. Linke, Ch. Merdon, M. Neilan, Pressurerobustness in quasioptimal a priori estimates for the Stokes problem, Preprint no. arXiv:1906.03009, Cornell University Library, arXiv.org, 2019.

N. Alia, V. John, S. Ollila, Revisiting the singlephase flow model for liquid steel ladle stirred by gas, Preprint no. arXiv.1811.11535, Cornell University Library, arXiv.org, 2018, DOI 10.1016/j.apm.2018.11.005 .
Abstract
Ladle stirring is an important step of the steelmaking process to homogenize the temperature and the chemical composition of the liquid steel and to remove inclusions before casting. Gas is injected from the bottom of the bath to induce a turbulent flow of the liquid steel. Multiphase modeling of ladle stirring can become computationally expensive, especially when used within optimal flow control problems. This paper focuses therefore on singlephase flow models. It aims at improving the existing models from the literature. Simulations in a 2d axialsymmetrical configuration, as well as, in a real 3d laboratoryscale ladle, are performed. The results obtained with the present model are in a relative good agreement with experimental data and suggest that it can be used as an efficient model in optimal flow control problems. 
N.R. Gauger, A. Linke, P. Schroeder, On highorder pressurerobust space discretisations, their advantages for incompressible high Reynolds number generalised Beltrami flows and beyond, Preprint no. arXiv:1808.10711, Cornell University Library, arXiv.org, 2018.
Abstract
Recently, highorder space discretisations were proposed for the numerical simulation of the incompressible NavierStokes equations at high Reynolds numbers, even for complicated domains from simulation practice. Although the overall spatial approximation order of the algorithms depends on the approximation quality of the boundary (often not better than third order), competitively accurate and efficient results were reported. In this contribution, first, a possible explanation for this somewhat surprising result is proposed: the velocity error of highorder space discretisations is more robust against quantitatively large and complicated pressure fields than loworder methods. Second, it is demonstrated that novel pressurerobust methods are significantly more accurate than comparable classical, nonpressurerobust space discretisations, whenever the quadratic, nonlinear convection term is a nontrivial gradient field like in certain generalised Beltrami flows at high Reynolds number. Then, pressurerobust methods even allow to halve the (formal) approximation order without compromising the accuracy. Third, classical highorder space discretisations are outperformed by pressurerobust methods whenever the boundary is not approximated with highorder accuracy. This improved accuracy of (loworder) pressurerobust mixed methods is explained in terms of a HelmholtzHodge projector, which cancels out the nonlinear convection term in any generalised Beltrami flow, since it is a gradient field. The numerical results are illustrated by a novel numerical analysis for pressurerobust and classical space discretisations. Further, the relevance of these results is discussed for flows that are not of Beltrami type.
Contact
Contributing Groups of WIAS
Projects/Grants
 Dynamic wetting and dewetting of viscous liquid droplets films on viscoelastic substrates
 Mathematical modeling and simulation of substrateflow interaction using generalized gradient flows
 Modeling and Analysis of Suspension Flows
 Multicriteria optimization subject to equilibrium constraints at the example of gas markets