2015, Ahmed, Naveed, Rebollo, Tomás Chacón, John, Volker, Rubino, Samuele

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 Navier-Stokes 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.

2017, John, Volker, Knobloch, Petr, Novo, Julia

The contents of this paper is twofold. First, important recent results concerning finite element methods for convection-dominated 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.

2015, John, Volker, Linke, Alexander, Merdon, Christian, Neilan, Michael, Rebholz, Leo G.

The divergence constraint of the incompressible Navier-Stokes 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 pressure-dependent consistency error which can potentially pollute the computed velocity. These methods are not robust in the sense that a contribution from the right-hand side, which in fluences 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 divergence-free solutions will be discussed: grad-div stabilization, higher order mixed methods derived on the basis of an exact de Rham complex, H(div)-conforming finite elements, and mixed methods with an appropriate reconstruction of the test functions. Numerical examples illustrate both the potential effects of using non-robust discretizations and the improvements obtained by utilizing pressure-robust discretizations.

2016, Ahmed, Naveed, John, Volker, Matthies, Gunar, Novo, Julia

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 continuous-in-time case for the approximations of both velocity and pressure. In addition, the fully discrete case in combination with higher order continuous Galerkin-Petrov (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 time-independent. Numerical results show that the predicted order is also achieved in the general case of time-dependent convective terms.

2019, John, Volker, Knobloch, Petr, Wilbrandt, Ulrich

Discretizations of incompressible flow problems with pairs of finite element spaces that do not satisfy a discrete inf-sup condition require a so-called pressure stabilization. This paper gives an overview and systematic assessment of stabilized methods, including the respective error analysis.

2014, Ahmed, Naveed, John, Volker

Higher order variational time stepping schemes allow an efficient post-processing for computing a higher order solution. This paper presents an adaptive algorithm whose time step control utilizes the post-processed solution. The algorithm is applied to convection-dominated convection-diffusion equations. It is shown that the length of the time step properly reflects the dynamics of the solution. The numerical costs of the adaptive algorithm are discussed.

2012, Barrenechea, Gabriel R., John, Volker, Knobloch, Petr

An extension of the local projection stabilization (LPS) finite element method for convection-diffusion-reaction equations is presented and analyzed, both in the steady-state 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.

2016, Wilbrandt, Ulrich, Bartsch, Clemens, Ahmed, Naveed, Alia, Najib, Anker, Felix, Blank, Laura, Caiazzo, Alfonso, Ganesa, Sashikumaar, Giere, Swetlana, Matthies, Gunar, Meesala, Raviteja, Shamim, Abdus, Venkatesan, Jagannath, John, Volker

PARMOON is a program package for the numerical solution of elliptic and parabolic partial differential equations. It inherits the distinct features of its predecessor MOONMD [28]: strict decoupling of geometry and finite element spaces, implementation of mapped finite elements as their definition can be found in textbooks, and a geometric multigrid preconditioner with the option to use different finite element spaces on different levels of the multigrid hierarchy. After having presented some thoughts about in-house research codes, this paper focuses on aspects of the parallelization, which is the main novelty of PARMOON. Numerical studies, performed on compute servers, assess the efficiency of the parallelized geometric multigrid preconditioner in comparison with parallel solvers that are available in external libraries. The results of these studies give a first indication whether the cumbersome implementation of the parallelized geometric multigrid method was worthwhile or not.

2013, John, Volker, Suciu, Carina

The accurate and efficient simulation of bi-variate population balance systems is nowadays a great challenge since the domain spanned by the external and internal coordinates is five-dimensional. This report considers direct discretizations of this equation in tensorproduct domains. In this situation, finite difference methods can be applied. The studied model includes the transport of dissolved potassium dihydrogen phosphate (KDP) and of energy (temperature) in a laminar flow field as well as the nucleation and growth of KDP particles. Two discretizations of the coupled model will be considered which differ only in the discretization of the population balance equation: a first order monotone upwind scheme and a third order essentially on-oscillatory (ENO) scheme. The Dirac term on the right-hand side of this equation is discretized with a finite volume method. The numerical results show that much different results are obtained even in the class of direct discretizations.

2010, John, Volker, Novo, Julia

Conditions on the stabilization parameters are explored for different approaches in deriving error estimates for the SUPG finite element stabilization of time-dependent convection-diffusion-reaction equations that is combined with the backward Euler method. Standard energy arguments lead to estimates for stabilization parameters that depend on the length of the time step. The stabilization vanishes in the time-continuous limit. However, based on numerical experiences, this seems not to be the correct behavior. For this reason, the time-continuous case is analyzed under certain conditions on the coefficients of the equation and the finite element method. An error estimate with the standard order of convergence is derived for stabilization parameters of the same form that is optimal for the steady-state problem. Numerical studies support the analytical results.