scholarly journals p-Multilevel Preconditioners for HHO Discretizations of the Stokes Equations with Static Condensation

Author(s):  
Lorenzo Botti ◽  
Daniele A. Di Pietro

AbstractWe propose a p-multilevel preconditioner for hybrid high-order (HHO) discretizations of the Stokes equation, numerically assess its performance on two variants of the method, and compare with a classical discontinuous Galerkin scheme. An efficient implementation is proposed where coarse level operators are inherited using $$L^2$$ L 2 -orthogonal projections defined over mesh faces and the restriction of the fine grid operators is performed recursively and matrix-free. Both h- and k-dependency are investigated tackling two- and three-dimensional problems on standard meshes and graded meshes. For the two HHO formulations, featuring discontinuous or hybrid pressure, we study how the combination of p-coarsening and static condensation influences the V-cycle iteration. In particular, two different static condensation procedures are considered for the discontinuous pressure HHO variant, resulting in global linear systems with a different number of unknowns and matrix non-zero entries. Interestingly, we show that the efficiency of the solution strategy might be impacted by static condensation options in the case of graded meshes.

1986 ◽  
Vol 108 (1) ◽  
pp. 68-75 ◽  
Author(s):  
S. Abdallah ◽  
H. G. Smith

The primitive variable formulation originally developed for the incompressible Navier–Stokes equations is applied for the solution of the incompressible Euler equations. The unsteady momentum equation is solved for the velocity field and the continuity equation is satisfied indirectly in a Poisson-type equation for the pressure (divergence of the momentum equation). Solutions for the pressure Poisson equation with derivative boundary conditions exist only if a compatibility condition is satisfied (Green’s theorem). This condition is not automatically satisfied on nonstaggered grids. A new method for the solution of the pressure equation with derivative boundary conditions on a nonstaggered grid [25] is used here for the calculation of the pressure. Three-dimensional solutions for the inviscid rotational flow in a 90 deg curved duct are obtained on a very fine mesh (17 × 17 × 29). The use of a fine grid mesh allows for the accurate prediction of the development of the secondary flow. The computed results are in good agreement with the experimental data of Joy [15].


2020 ◽  
Vol 14 (4) ◽  
pp. 7369-7378
Author(s):  
Ky-Quang Pham ◽  
Xuan-Truong Le ◽  
Cong-Truong Dinh

Splitter blades located between stator blades in a single-stage axial compressor were proposed and investigated in this work to find their effects on aerodynamic performance and operating stability. Aerodynamic performance of the compressor was evaluated using three-dimensional Reynolds-averaged Navier-Stokes equations using the k-e turbulence model with a scalable wall function. The numerical results for the typical performance parameters without stator splitter blades were validated in comparison with experimental data. The numerical results of a parametric study using four geometric parameters (chord length, coverage angle, height and position) of the stator splitter blades showed that the operational stability of the single-stage axial compressor enhances remarkably using the stator splitter blades. The splitters were effective in suppressing flow separation in the stator domain of the compressor at near-stall condition which affects considerably the aerodynamic performance of the compressor.


Author(s):  
Surabhi Rathore ◽  
Tomoki Uda ◽  
Viet Q. H. Huynh ◽  
Hiroshi Suito ◽  
Toshitaka Watanabe ◽  
...  

AbstractHemodialysis procedure is usually advisable for end-stage renal disease patients. This study is aimed at computational investigation of hemodynamical characteristics in three-dimensional arteriovenous shunt for hemodialysis, for which computed tomography scanning and phase-contrast magnetic resonance imaging are used. Several hemodynamical characteristics are presented and discussed depending on the patient-specific morphology and flow conditions including regurgitating flow from the distal artery caused by the construction of the arteriovenous shunt. A simple backflow prevention technique at an outflow boundary is presented, with stabilized finite element approaches for incompressible Navier–Stokes equations.


Processes ◽  
2021 ◽  
Vol 9 (5) ◽  
pp. 841
Author(s):  
Yuzhen Jin ◽  
Huang Zhou ◽  
Linhang Zhu ◽  
Zeqing Li

A three-dimensional numerical study of a single droplet splashing vertically on a liquid film is presented. The numerical method is based on the finite volume method (FVM) of Navier–Stokes equations coupled with the volume of fluid (VOF) method, and the adaptive local mesh refinement technology is adopted. It enables the liquid–gas interface to be tracked more accurately, and to be less computationally expensive. The relationship between the diameter of the free rim, the height of the crown with different numbers of collision Weber, and the thickness of the liquid film is explored. The results indicate that the crown height increases as the Weber number increases, and the diameter of the crown rim is inversely proportional to the collision Weber number. It can also be concluded that the dimensionless height of the crown decreases with the increase in the thickness of the dimensionless liquid film, which has little effect on the diameter of the crown rim during its growth.


Author(s):  
Eiman B Saheby ◽  
Xing Shen ◽  
Anthony P Hays ◽  
Zhang Jun

This study describes the aerodynamic efficiency of a forebody–inlet configuration and computational investigation of a drone system, capable of sustainable supersonic cruising at Mach 1.60. Because the whole drone configuration is formed around the induction system and the design is highly interrelated to the flow structure of forebody and inlet efficiency, analysis of this section and understanding its flow pattern is necessary before any progress in design phases. The compression surface is designed analytically using oblique shock patterns, which results in a low drag forebody. To study the concept, two inlet–forebody geometries are considered for Computational Fluid Dynamic simulation using ANSYS Fluent code. The supersonic and subsonic performance, effects of angle of attack, sideslip, and duct geometries on the propulsive efficiency of the concept are studied by solving the three-dimensional Navier–Stokes equations in structured cell domains. Comparing the results with the available data from other sources indicates that the aerodynamic efficiency of the concept is acceptable at supersonic and transonic regimes.


Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 756
Author(s):  
Federico Lluesma-Rodríguez ◽  
Francisco Álcantara-Ávila ◽  
María Jezabel Pérez-Quiles ◽  
Sergio Hoyas

One numerical method was designed to solve the time-dependent, three-dimensional, incompressible Navier–Stokes equations in turbulent thermal channel flows. Its originality lies in the use of several well-known methods to discretize the problem and its parallel nature. Vorticy-Laplacian of velocity formulation has been used, so pressure has been removed from the system. Heat is modeled as a passive scalar. Any other quantity modeled as passive scalar can be very easily studied, including several of them at the same time. These methods have been successfully used for extensive direct numerical simulations of passive thermal flow for several boundary conditions.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Niklas Ericsson

Abstract We develop a framework for solving the stationary, incompressible Stokes equations in an axisymmetric domain. By means of Fourier expansion with respect to the angular variable, the three-dimensional Stokes problem is reduced to an equivalent, countable family of decoupled two-dimensional problems. By using decomposition of three-dimensional Sobolev norms, we derive natural variational spaces for the two-dimensional problems, and show that the variational formulations are well-posed. We analyze the error due to Fourier truncation and conclude that, for data that are sufficiently regular, it suffices to solve a small number of two-dimensional problems.


2000 ◽  
Vol 411 ◽  
pp. 1-38 ◽  
Author(s):  
C. ROSS ETHIER ◽  
SUJATA PRAKASH ◽  
DAVID A. STEINMAN ◽  
RICHARD L. LEASK ◽  
GREGORY G. COUCH ◽  
...  

Numerical and experimental techniques were used to study the physics of flow separation for steady internal flow in a 45° junction geometry, such as that observed between two pipes or between the downstream end of a bypass graft and an artery. The three-dimensional Navier–Stokes equations were solved using a validated finite element code, and complementary experiments were performed using the photochromic dye tracer technique. Inlet Reynolds numbers in the range 250 to 1650 were considered. An adaptive mesh refinement approach was adopted to ensure grid-independent solutions. Good agreement was observed between the numerical results and the experimentally measured velocity fields; however, the wall shear stress agreement was less satisfactory. Just distal to the ‘toe’ of the junction, axial flow separation was observed for all Reynolds numbers greater than 250. Further downstream (approximately 1.3 diameters from the toe), the axial flow again separated for Re [ges ] 450. The location and structure of axial flow separation in this geometry is controlled by secondary flows, which at sufficiently high Re create free stagnation points on the model symmetry plane. In fact, separation in this flow is best explained by a secondary flow boundary layer collision model, analogous to that proposed for flow in the entry region of a curved tube. Novel features of this flow include axial flow separation at modest Re (as compared to flow in a curved tube, where separation occurs only at much higher Re), and the existence and interaction of two distinct three-dimensional separation zones.


Sign in / Sign up

Export Citation Format

Share Document