Algorithm for analysis of flows in ribbed annuli

• Published in: International journal for numerical methods in fluids 2012-03, Vol.68 (7), p.805-838
• Main Authors: Moradi, H. V., Floryan, J. M.
• Format: Article
• Language: English
• Subjects: Computational methods in fluid dynamics; differential equations; Exact sciences and technology; Flows in ducts, channels, nozzles, and conduits; Fluid dynamics; Fundamental areas of phenomenology (including applications); Galerkin; immersed boundary; laminar flow; Navier-Stokes; Physics; spectral
• ISSN: 0271-2091; 1097-0363
• DOI: 10.1002/fld.2581

SUMMARY Spectral methods for analyses of steady flows in annuli bounded by walls with either axi‐symmetric or longitudinal ribs are developed. The physical boundary conditions are enforced using the immersed boundary conditions concept. In the former case, the Stokes stream function is used to eliminate pressure and to reduce system of field equations to a single fourth‐order partial differential equation. The ribs are assumed to be periodic in the axial direction and this permits representation of the solution in terms of the Fourier expansion. In the latter case, the problem is reduced to the Laplace equations for the flow modifications that can be expressed in terms of the Fourier expansions. The modal functions, which are functions of the radial coordinate, are discretized using Chebyshev polynomials. The problem formulations are closed using either the fixed volume flow rate constraint or the fixed pressure gradient constraint. Various tests have been carried out in order to demonstrate the spectral accuracy of the discretizations, as well as the spectral accuracy of the enforcement of the flow boundary conditions at the ribbed walls using the immersed boundary conditions concept. Special linear solver that takes advantage of the matrix structure has been implemented in order to reduce computational time and memory requirements. It is shown that the algorithm has superior performance when one is interested in the analysis of a large number of geometries, as part of the coefficient matrix that corresponds to the field equation is always the same and one needs to change only the part of the matrix that corresponds to the boundary relations when changing geometry of the flow domain. Copyright © 2011 John Wiley & Sons, Ltd.