Poloidal-toroidal decomposition in a finite cylinder. II: Discretization, regularization and validation

From MaRDI portal
Publication:2465009

DOI10.1016/J.JCP.2007.08.035zbMATH Open1126.76056arXiv0705.0785OpenAlexW2068568713MaRDI QIDQ2465009FDOQ2465009


Authors: Piotr Boronski, Laurette S. Tuckerman Edit this on Wikidata


Publication date: 18 December 2007

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: The Navier-Stokes equations in a finite cylinder are written in terms of poloidal and toroidal potentials in order to impose incompressibility. Regularity of the solutions is ensured in several ways: First, the potentials are represented using a spectral basis which is analytic at the cylindrical axis. Second, the non-physical discontinuous boundary conditions at the cylindrical corners are smoothed using a polynomial approximation to a steep exponential profile. Third, the nonlinear term is evaluated in such a way as to eliminate singularities. The resulting pseudo-spectral code is tested using exact polynomial solutions and the spectral convergence of the coefficients is demonstrated. Our solutions are shown to agree with exact polynomial solutions and with previous axisymmetric calculations of vortex breakdown and of nonaxisymmetric calculations of onset of helical spirals. Parallelization by azimuthal wavenumber is shown to be highly effective.


Full work available at URL: https://arxiv.org/abs/0705.0785




Recommendations




Cites Work


Cited In (6)





This page was built for publication: Poloidal-toroidal decomposition in a finite cylinder. II: Discretization, regularization and validation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2465009)