WIAS Preprints

Permanent URI for this collection

https://oa.tib.eu/renate/handle/123456789/14629

Browse

Now showing 1 - 16 of 16

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.
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.
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.
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.
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.
Higher order continuous Galerkin-Petrov time stepping schemes for transient convection-diffusion-reaction equations
(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Ahmed, Naveed; Matties, Gunar
We present the analysis for the higher order continuous Galerkin-Petrov (cGP) time discretization schemes in combination with the one-level local projection stabilization in space applied to time-dependent convection-diffusion-reaction problems. Optimal a-priori error estimates will be proved. Numerical studies support the theoretical results. Furthermore, a numerical comparison between continuous Galerkin-Petrov and discontinuous Galerkin time discretization schemes will be given.
Higher-order discontinuous Galerkin time stepping and local projection stabilization techniques for the transient Stokes problem
(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Ahmed, Naveed; Becher, Simon; Matthies, Gunar
We introduce and analyze discontinuous Galerkin time discretizations coupled with continuous finite element methods based on equal-order 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 semi-discrete 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.
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.
Numerical studies of higher order variational time stepping schemes for evolutionary Navier-Stokes equations
(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Ahmed, Naveed; Matthies, Gunar
We present in this paper numerical studies of higher order variational time stepping schemes combined with finite element methods for simulations of the evolutionary Navier-Stokes equations. In particular, conforming inf-sup stable pairs of finite element spaces for approximating velocity and pressure are used as spatial discretization while continuous GalerkinPetrov methods (cGP) and discontinuous Galerkin (dG) methods are applied as higher order variational time discretizations. Numerical results for the well-known problem of incompressible flows around a circle will be presented.
Numerical study of SUPG and LPS methods combined with higher order variational time discretization schemes applied to time-dependent convection-diffusion-reaction equations
(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Ahmed, Naveed; Matthies, Gunar
This paper considers the numerical solution of time-dependent convection-diffusion-reaction equations. We shall employ combinations of streamline-upwind Petrov-Galerkin (SUPG) and local projection stabilization (LPS) methods in space with the higher order variational time discretization schemes. In particular, we consider time discretizations by discontinuous Galerkin (dG) methods and continuous Galerkin-Petrov (cGP) methods. Several numerical tests have been performed to assess the accuracy of combinations of spatial and temporal discretization schemes. Furthermore, the dependence of the results on the stabilization parameters of the spatial discretizations are discussed. Finally the long-time behavior of overshoots and undershoots is investigated.
On really locking-free mixed finite element methods for the transient incompressible Stokes equations
(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Ahmed, Naveed; Linke, Alexander; Merdon, Christian
Inf-sup stable mixed methods for the steady incompressible Stokes equations that relax the divergence constraint are often claimed to deliver locking-free discretizations. However, this relaxation leads to a pressure-dependent 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 locking-free, 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. Semi-discrete and fully-discrete a-priori velocity and pressure error estimates are derived, which show beautiful robustness properties. Two numerical examples illustrate the superior accuracy of pressure-robust space discretizations in the case of small viscosities.
On the feasibility of using open source solvers for the simulation of a turbulent air flow in a dairy barn
(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Janke, David; Caiazzo, Alfonso; Ahmed, Naveed; Alia, Najib; Knoth, Oswald; Moreau, Baptiste; Wilbrandt, Ulrich; Willink, Dilya; Amon, Thomas; John, Volker
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 time-averaged stream-wise and vertical-wise 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 root-mean-square 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 time-dependent phenomena in an agricultural context, e.g., like the transport of particulate matter or pathogen-laden aerosols in and around agricultural buildings.
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.
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.
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.
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.

Browse

Browsing WIAS Preprints by Author "Ahmed, Naveed"

Results Per Page

Sort Options