Splitting of Taylor Hood Discretizations and Minimal Extension for Stable Time Integration of Nonviscous Flow
For the theory see our Preprint.
To check invertibility of the local FEM divergence matrices of regular n-gons use the script detOfRegTayHooCluster.py
.
The simulations are started from through file
python ct_main.py
There all parameters are set. The default values of internal parameters of the first commit are the ones used in the simulations presented in the preprint.
You will need to create the directories data
, mats
, results
.
To compute and plot the errors for the cylinder wake setup use
python check_cyl_errs
You will need a Dolfin/FeNiCs installation. You may want to install Krypy, but you can solve the linear systems with a different solver.
Other features:
- Paraview integration for visualization
- Convergence plots in the ipython notebook
- Storing of simulation data in json files