Preconditioning Techniques for Nonsymmetric Linear Systems in the Computation of Incompressible Flows

2009 ◽  
Vol 76 (2) ◽  
Author(s):  
Murat Manguoglu ◽  
Ahmed H. Sameh ◽  
Faisal Saied ◽  
Tayfun E. Tezduyar ◽  
Sunil Sathe

In this paper we present effective preconditioning techniques for solving the nonsymmetric systems that arise from the discretization of the Navier–Stokes equations. These linear systems are solved using either Krylov subspace methods or the Richardson scheme. We demonstrate the effectiveness of our techniques in handling time-accurate as well as steady-state solutions. We also compare our solvers with those published previously.

AIAA Journal ◽  
2000 ◽  
Vol 38 ◽  
pp. 1603-1614
Author(s):  
Martin Scholtysik ◽  
Bernhard Mueller ◽  
Torstein K. Fannelop

2014 ◽  
Author(s):  
P. Bigay ◽  
A. Bardin ◽  
G. Oger ◽  
D. Le Touzé

In order to efficiently address complex problems in hydrodynamics, the advances in the development of a new method are presented here. This method aims at finding a good compromise between computational efficiency, accuracy, and easy handling of complex geometries. The chosen method is an Explicit Cartesian Finite Volume method for Hydrodynamics (ECFVH) based on a compressible (hyperbolic) solver, with a ghost-cell method for geometry handling and a Level-set method for the treatment of biphase-flows. The explicit nature of the solver is obtained through a weakly-compressible approach chosen to simulate nearly-incompressible flows. The explicit cell-centered resolution allows for an efficient solving of very large simulations together with a straightforward handling of multi-physics. A characteristic flux method for solving the hyperbolic part of the Navier-Stokes equations is used. The treatment of arbitrary geometries is addressed in the hyperbolic and viscous framework. Viscous effects are computed via a finite difference computation of viscous fluxes and turbulent effects are addressed via a Large-Eddy Simulation method (LES). The Level-Set solver used to handle biphase flows is also presented. The solver is validated on 2-D test cases (flow past a cylinder, 2-D dam break) and future improvements are discussed.


1987 ◽  
Vol 109 (4) ◽  
pp. 345-352 ◽  
Author(s):  
M. Reggio ◽  
R. Camarero

A numerical procedure to solve three-dimensional incompressible flows in arbitrary shapes is presented. The conservative form of the primitive-variable formulation of the time-dependent Navier-Stokes equations written for a general curvilinear coordiante system is adopted. The numerical scheme is based on an overlapping grid combined with opposed differencing for mass and pressure gradients. The pressure and the velocity components are stored at the same location: the center of the computational cell which is used for both mass and the momentum balance. The resulting scheme is stable and no oscillations in the velocity or pressure fields are detected. The method is applied to test cases of ducting and the results are compared with experimental and numerical data.


2005 ◽  
Vol 9 (1) ◽  
pp. 67-78 ◽  
Author(s):  
J. Socolowsky

iscous two‐fluid channel flows arise in different kinds of coating technologies. The corresponding mathematical models represent two‐dimensional free boundary value problems for the Navier‐Stokes equations. In this paper the solvability of the related stationary problems is discussed and computational results are presented. Furthermore, it is shown that depending on the flow parameters like viscosity or density ratios and on the fluxes there can happen nonexistence of steady‐state solutions. For other parameter sets the solution is even unique. Dvieju, tekančiu kanale, klampiu skysčiu srauto uždavinys iškyla taikant ivairias skirtingu rušiu paviršiu padengimo technologijas. Atitinkamas matematinis modelis išreiškiamas dvimačiu kraštiniu uždaviniu su laisvu paviršiumi Navje-Stokso lygtims. Straipsnyje nagrinejamas santykinai stacionaraus uždavinio išsprendžiamumas ir pateikiami skaičiavimo rezultatai. Be to parodoma, kad priklausomai nuo sroves parametru kaip ir nuo klampumo ir tankio santykio stacionarus sprendiniai gali neegzistuoti. Su kitais parametrais egzistuoja tiksliai vienas sprendinys.


1969 ◽  
Vol 37 (4) ◽  
pp. 727-750 ◽  
Author(s):  
Gareth P. Williams

A method of numerically integrating the Navier-Stokes equations for certain three-dimensional incompressible flows is described. The technique is presented through application to the particular problem of describing thermal convection in a rotating annulus. The equations, in cylindrical polar co-ordinate form, are integrated with respect to time by a marching process, together with the solving of a Poisson equation for the pressure. A suitable form of the finite difference equations gives a computationally-stable long-term integration with reasonably faithful representation of the spatial and temporal characteristics of the flow.Trigonometric interpolation techniques provide accurate (discretely exact) solutions to the Poisson equation. By using an auxiliary algorithm for rapid evaluation of trigonometric transforms, the proportion of computation needed to solve the Poisson equation can be reduced to less than 25% of the total time needed to’ advance one time step. Computing on a UNIVAC 1108 machine, the flow can be advanced one time-step in 2 sec for a 14 × 14 × 14 grid upward to 96 sec for a 60 × 34 × 34 grid.As an example of the method, some features of a solution for steady wave flow in annulus convection are presented. The resemblance of this flow to the classical Eady wave is noted.


2018 ◽  
Vol 855 ◽  
pp. 43-66 ◽  
Author(s):  
Ke Wu ◽  
Bruno D. Welfert ◽  
Juan M. Lopez

The dynamic response to shear of a fluid-filled square cavity with stable temperature stratification is investigated numerically. The shear is imposed by the constant translation of the top lid, and is quantified by the associated Reynolds number. The stratification, quantified by a Richardson number, is imposed by maintaining the temperature of the top lid at a higher constant temperature than that of the bottom, and the side walls are insulating. The Navier–Stokes equations under the Boussinesq approximation are solved, using a pseudospectral approximation, over a wide range of Reynolds and Richardson numbers. Particular attention is paid to the dynamical mechanisms associated with the onset of instability of steady state solutions, and to the complex and rich dynamics occurring beyond.


Sign in / Sign up

Export Citation Format

Share Document