The linear stability of parallel shear flows of incompressible viscous fluids is classically described by the Orr–Sommerfeld equation in the disturbance streamfunction. This fourth-order equation is obtained by eliminating the pressure from the linearized Navier–Stokes equation. Here we consider retaining the primitive velocity-pressure formulation, as is required for general multidimensional geometries for which the streamfunction is unavailable; this affords a uniform description of one-, two-, and three-dimensional flows and their perturbations. The Orr–Sommerfeld equation is here discretized using Python and scikit- fem, in classical and primitive forms with Hermite and Mini elements, respectively. The solutions for the standard test problem of plane Poiseuille flow show the primitive formulation to be simple, clear, very accurate, and better-conditioned than the classical.
%0 Journal Article
%1 mcbain2022primitive
%A McBain, Geordie Drummond
%D 2021
%J The ANZIAM Journal
%K 65f15-numerical-eigenvalues-eigenvectors 65f35-matrix-norms-conditioning-scaling 76e05-parallel-shear-flows 76m10-finite-element-methods-in-fluid-mechanics scikit-fem
%P C168--C181
%R 10.21914/anziamj.v63.17159
%T The primitive Orr–Sommerfeld equation and its solution by finite elements
%U https://journal.austms.org.au/ojs/index.php/ANZIAMJ/article/view/17159
%V 63
%X The linear stability of parallel shear flows of incompressible viscous fluids is classically described by the Orr–Sommerfeld equation in the disturbance streamfunction. This fourth-order equation is obtained by eliminating the pressure from the linearized Navier–Stokes equation. Here we consider retaining the primitive velocity-pressure formulation, as is required for general multidimensional geometries for which the streamfunction is unavailable; this affords a uniform description of one-, two-, and three-dimensional flows and their perturbations. The Orr–Sommerfeld equation is here discretized using Python and scikit- fem, in classical and primitive forms with Hermite and Mini elements, respectively. The solutions for the standard test problem of plane Poiseuille flow show the primitive formulation to be simple, clear, very accurate, and better-conditioned than the classical.
@article{mcbain2022primitive,
abstract = {The linear stability of parallel shear flows of incompressible viscous fluids is classically described by the Orr–Sommerfeld equation in the disturbance streamfunction. This fourth-order equation is obtained by eliminating the pressure from the linearized Navier–Stokes equation. Here we consider retaining the primitive velocity-pressure formulation, as is required for general multidimensional geometries for which the streamfunction is unavailable; this affords a uniform description of one-, two-, and three-dimensional flows and their perturbations. The Orr–Sommerfeld equation is here discretized using Python and scikit- fem, in classical and primitive forms with Hermite and Mini elements, respectively. The solutions for the standard test problem of plane Poiseuille flow show the primitive formulation to be simple, clear, very accurate, and better-conditioned than the classical.},
added-at = {2022-09-21T08:55:55.000+0200},
author = {McBain, Geordie Drummond},
biburl = {https://www.bibsonomy.org/bibtex/2a97953c4b945068af7f40e8782ec619f/gdmcbain},
doi = {10.21914/anziamj.v63.17159},
id = {17159},
interhash = {9710aae4d3a46b6e7e2760e1e9f1062c},
intrahash = {a97953c4b945068af7f40e8782ec619f},
issn = {1445-8810},
journal = {The ANZIAM Journal},
keywords = {65f15-numerical-eigenvalues-eigenvectors 65f35-matrix-norms-conditioning-scaling 76e05-parallel-shear-flows 76m10-finite-element-methods-in-fluid-mechanics scikit-fem},
pages = {C168--C181},
source = {ANZIAM Journal},
timestamp = {2022-09-21T08:55:55.000+0200},
title = {The primitive Orr–Sommerfeld equation and its solution by finite elements},
uri = {https://journal.austms.org.au/ojs/index.php/ANZIAMJ},
url = {https://journal.austms.org.au/ojs/index.php/ANZIAMJ/article/view/17159},
volume = 63,
year = 2021
}