Search Results

Now showing 1 - 10 of 16
  • Item
    A review of variational multiscale methods for the simulation of turbulent incompressible flows
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 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.
  • Item
    A pressure-robust discretization of Oseen's equation using stabilization in the vorticity equation
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Ahmed, Naveed; Barrenechea, Gabriel R.; Burman, Erik; Guzmán, Johnny; Linke, Alexander; Merdon, Christian
    Discretization of Navier--Stokes' equations using pressure-robust 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 residual-based 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 pressure-robust. Thus, we prove pressureindependent error estimates in the linearized case, known as Oseen's problem. In fact, we prove an O(hk+1/2) error estimate in the L2-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 residual-based SUPG stabilization.
  • Item
    ParMooN - a modernized program package based on mapped finite elements
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 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.
  • Item
    Adaptive time step control for higher order variational time discretizations applied to convection-diffusion equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 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.
  • Item
    An assessment of solvers for algebraically stabilized discretizations of convection-diffusion-reaction equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2021) Jha, Abhinav; Pártl, Ondřej; Ahmed, Naveed; Kuzmin, Dmitri
    We consider flux-corrected finite element discretizations of 3D convection-dominated transport problems and assess the computational efficiency of algorithms based on such approximations. The methods under investigation include flux-corrected transport schemes and monolithic limiters. We discretize in space using a continuous Galerkin method and P1 or Q1 finite elements. Time integration is performed using the Crank-Nicolson method or an explicit strong stability preserving Runge-Kutta method. Nonlinear systems are solved using a fixed-point iteration method, which requires solution of large linear systems at each iteration or time step. The great variety of options in the choice of discretization methods and solver components calls for a dedicated comparative study of existing approaches. To perform such a study, we define new 3D test problems for time dependent and stationary convection-diffusion-reaction equations. The results of our numerical experiments illustrate how the limiting technique, time discretization and solver impact on the overall performance.
  • Item
    Detection of electromagnetic inclusions using topological sensitivity
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Wahab, Abdul; Abbas, Tasawar; Ahmed, Naveed; Zia, Qazi Muhammad Zaigham
    In this article a topological sensitivity framework for far field detection of a diametrically small electromagnetic inclusion is established. The cases of single and multiple measurements of the electric far field scattering amplitude at a fixed frequency are taken into account. The performance of the algorithm is analyzed theoretically in terms of its resolution and sensitivity for locating an inclusion. The stability of the framework with respect to measurement and medium noises is discussed. Moreover, the quantitative results for signal-to-noise ratio are presented. A few numerical results are presented to illustrate the detection capabilities of the proposed framework with single and multiple measurements.
  • Item
    A local projection stabilization/continuous Galerkin-Petrov method for incompressible flow problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 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.
  • Item
    Analysis of a full space-time discretization of the Navier-Stokes equations by a local projection stabilization
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Ahmed, Naveed; Rebollo, Tomás Chacón; John, Volker; Rubino, Samuele
    A finite element error analysis of a local projection stabilization (LPS) method for the time-dependent Navier-Stokes equations is presented. The focus is on the highorder term-by-term stabilization method that has one level, in the sense that it is defined on a single mesh, and in which the projection-stabilized structure of standard LPS methods is replaced by an interpolation-stabilized 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 ows. Smooth unsteady flows are simulated with optimal order of accuracy.
  • Item
    An assessment of two classes of variational multiscale methods for the simulation of incompressible turbulent flows
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Ahmed, Naveed; John, Volker
    A numerical assessment of two classes of variational multiscale (VMS) methods for the simulation of incompressible flows is presented. Two types of residual-based VMS methods and two types of projection-based 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 residual-based VMS methods, in particular when used with a pair of inf-sup 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.
  • Item
    On the grad-div stabilization for the steady Oseen and Navier-Stokes equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Ahmed, Naveed
    This paper studies the parameter choice in the grad-div stabilization applied to the generalized problems of Oseen type. Stabilization parameters based on minimizing the H1 (Omega) error of the velocity are derived which do not depend on the viscosity parameter. For the proposed parameter choices, the H1 (Omega) 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 two-dimensional flow past a circular cylinder. It turns out, for the MINI element, that the best results can be obtained without grad-div stabilization.