scholarly journals Scalable Block Preconditioners for Linearized Navier-Stokes Equations at High Reynolds Number

Algorithms ◽  
2020 ◽  
Vol 13 (8) ◽  
pp. 199
Author(s):  
Filippo Zanetti ◽  
Luca Bergamaschi

We review a number of preconditioners for the advection-diffusion operator and for the Schur complement matrix, which, in turn, constitute the building blocks for Constraint and Triangular Preconditioners to accelerate the iterative solution of the discretized and linearized Navier-Stokes equations. An intensive numerical testing is performed onto the driven cavity problem with low values of the viscosity coefficient. We devise an efficient multigrid preconditioner for the advection-diffusion matrix, which, combined with the commuted BFBt Schur complement approximation, and inserted in a 2×2 block preconditioner, provides convergence of the Generalized Minimal Residual (GMRES) method in a number of iteration independent of the meshsize for the lowest values of the viscosity parameter. The low-rank acceleration of such preconditioner is also investigated, showing its great potential.

1958 ◽  
Vol 8 ◽  
pp. 966-974
Author(s):  
H. E. Petschek

Analyses of aerodynamic dissipation in ordinary un-ionized gases are all based upon the Navier-Stokes equations. These equations relate the rate of dissipation to the local gradients in velocity and temperature through the viscosity and heat conduction coefficients. Although it is true that in many flow situations the magnitude of the total dissipation in the gas does not depend on the magnitude of the viscosity coefficient, this coefficient does determine the minimum scale of variations observed in the gas and the form of the Navier-Stokes equations determines the type of phenomena which are observed on a small scale. In order to discuss dissipation in an ionized gas in the presence of a magnetic field, it is therefore necessary to re-examine the derivation of the basic flow equations. This paper attempts to do this for a case of a completely ionized gas and demonstrates that the basic microscopic dissipation mechanism is appreciably different. For example, it is shown that the minimum length in which the properties of the flow field can change noticeably is appreciably less than one mean free path.


Author(s):  
Iraj Saeedpanah ◽  
M. Shayanfar ◽  
E. Jabbari ◽  
Mohammad Haji Mohammadi

Free surface flows are frequently encountered in hydraulic engineering problems including water jets, weirs and around gates. An iterative solution to the incompressible two-dimensional vertical steady Navier-Stokes equations, comprising momentum and continuity equations, is used to solve for the priori unknown free surface, the velocity and the pressure fields. The entire water body is covered by a unstructured finite element grid which is locally refined. The dynamic boundary condition is imposed for the free surface where the pressure vanishes. This procedure is done continuously until the normal velocities components vanish. To overcome numerical errors and oscillations encountering in convection terms, the SUPG (streamline upwinding Petrov-Galerkin) method is applied. The solution method is tested for different discharges onto a standard spillway geometries. The results shows good agreement with available experimental data.


2013 ◽  
Vol 13 (5) ◽  
pp. 1309-1329 ◽  
Author(s):  
Laura Lazar ◽  
Richard Pasquetti ◽  
Francesca Rapetti

AbstractSpectral element methods on simplicial meshes, say TSEM, show both the advantages of spectral and finite element methods, i.e., spectral accuracy and geometrical flexibility. We present aTSEM solver of the two-dimensional (2D) incompressible Navier-Stokes equations, with possible extension to the 3D case. It uses a projection method in time and piecewise polynomial basis functions of arbitrary degree in space. The so-called Fekete-Gauss TSEM is employed,i.e., Fekete (resp. Gauss) points of the triangle are used as interpolation (resp. quadrature) points. For the sake of consistency, isoparametric elements are used to approximate curved geometries. The resolution algorithm is based on an efficient Schur complement method, so that one only solves for the element boundary nodes. Moreover, the algebraic system is never assembled, therefore the number of degrees of freedom is not limiting. An accuracy study is carried out and results are provided for classical benchmarks: the driven cavity flow, the flow between eccentric cylinders and the flow past a cylinder.


2013 ◽  
Vol 734 ◽  
pp. 275-316 ◽  
Author(s):  
Rashad Moarref ◽  
Ati S. Sharma ◽  
Joel A. Tropp ◽  
Beverley J. McKeon

AbstractWe study the Reynolds-number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier–Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behaviour with Reynolds number in the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it is shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. Under the practical assumption that the mean velocity consists of a logarithmic region, we identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_{\tau } \approx 1{0}^{3} {\unicode{x2013}} 1{0}^{10} $). Results from this low-rank model of the Navier–Stokes equations compare favourably with experimental results in the literature.


Sign in / Sign up

Export Citation Format

Share Document