Computational results and python files for the work "Divergence-conforming velocity and vorticity approximations for incompressible fluids obtained with minimal facet coupling"

This repository contains data accompanying the paper "Divergence-conforming velocity and vorticity approximations for incompressible fluids obtained with minimal facet coupling".

The implementation is based on the python-interface of the NGSolve open source Finite Element library (ngsolve.org).

The file solve_problem_allione.py represents a minimum working example where the proposed MCS/HDG (set the use_MCS flag) method is used to solve the problem from the numerics section of the paper.

The files FlowTemplates.py and krylovspace_extension.py contain a somewhat larger and more modular implementation of the proposed method that also features preconditioned iterative solvers, including support for the NgsAMG NGSolve extension library as well as the NGSolve-PETSc interface.

The files errors_hdg.pickle, errors_mcs.pickle and kappas.pickle contain the raw data the tables and pictures in the paper were generated from.

This data was generated with the scripts conv3d_hdg.py, conv3d_mcs.py and calc_kappas.py which use the FlowTemplates.py infrastructure.