Skip to content

sscollis/stab

Repository files navigation

STAB: Compressible stability solver

Stab uses a direct approach with either Chebyshev collocation or high-order finite difference. This is complemented by the shoot solver which uses an ODE solver with Conte orthogonalization and shooting to satisfy the boundary conditions. I suggest using shoot to polish the eigensolutions identified by stab.

The stab code has been updated to build with the GCC gfortran compiler on both MacOS and Linux and older builds (SGI, Cray, IBM) have not been recently used so makefiles for those will need updates.

Building

To build use:

    ln -s gcc.mak Makefile
    make USE_NR=1

Notes

stab currently uses a couple of Numerical-Recipes in FORTRAN routines that are commercially licensed and therefore not distributed here. You can easily replace these with public versions and are encouraged to do so.

  1. It would be easy to remove the dependency on NR and I encourage someone to do so.
  2. You may not add the NR software to this repository and you may only use it if you have a valid license to do so.
  3. I have not recently tested the finite-Reynolds number surface roughness solver so please use with caution.

Test Cases

Case Description
thesis Test case from Collis Ph.D. thesis, Chapter 4 (Start here!)
compbl Similar to thesis but using compbl instead of fsc for mean, slight differences
TStest Tollmein-Schlichting test case
CFtest Cross-flow vortex test case
mixl Mixing layer test case
test Test case for experimentation
vk Case for a student

Contact

S. Scott Collis
flow.physics.simulation@gmail.com