A Stokes-Brinkman solver in Julia

View the Project on GitHub bmwilly/brinkman-stokes

Intro

High-order, parallel, matrix-free finite element solver for Stokes and Brinkman flow written in Julia.

Details and benchmarks can be found here.

Basic Usage

There are two implementations of the software, a parallel Julia implementation and a serial Matlab implementation. For each implementation, there are three versions that correspond to various finite element assembly methods, using assembled matrices, matrix-free methods using matrix-vector products, and matrix-free methods using matrix-matrix products.

From the top-level directory, cd julia-parallel/efficient-operators. Start julia and run include(sss.jl) to start an interactive session that will prompt you for a type of test problem to solve (lid-driven cavity or brinkman flow through obstacles), mesh size, preconditioning steps, etc.