This code was written by Anand Radhakrishnan
This code solves for an incompressible Couette flow between two cylinders, with the inner one rotating at a constant angular speed.
The theoretical formulation is in the JCP_Paper.pdf
paper in the root directory.
Mesh files .msh
in conc_circle_geoms/
were generated with gmsh
.
The code has 4 main components:
- Operator discretization for pressure Poisson with all Neumann boundary conditions in
Grid.cpp
- Multigrid preconditioned GMRES/ ILU-GMRES solvers in
FracStepMultirid.cpp
- Fractional Step procedure in
fractionalStepGrid.cpp
- Main file which generates the simulation in
FractionalStepSim.cpp
Build using make all