Oasis: A high-level/high-performance open source Navier–Stokes solver

Oasis is a high-level/high-performance finite element Navier–Stokes solver written from scratch in Python using building blocks from the FEniCS project (fenicsproject.org). The solver is unstructured and targets large-scale applications in complex geometries on massively parallel clusters. Oasis utilizes MPI and interfaces, through FEniCS, to the linear algebra backend PETSc. Oasis advocates a high-level, programmable user interface through the creation of highly flexible Python modules for new problems. Through the high-level Python interface the user is placed in complete control of every aspect of the solver. A version of the solver, that is using piecewise linear elements for both velocity and pressure, is shown to reproduce very well the classical, spectral, turbulent channel simulations of Moser et al. (1999). The computational speed is strongly dominated by the iterative solvers provided by the linear algebra backend, which is arguably the best performance any similar implicit solver using PETSc may hope for. Higher order accuracy is also demonstrated and new solvers may be easily added within the same framework.Program summaryProgram title: OasisCatalogue identifier: AEUW_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEUW_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: GNU Lesser GPL version 3 or any later versionNo. of lines in distributed program, including test data, etc.: 3491No. of bytes in distributed program, including test data, etc.: 266924Distribution format: tar.gzProgramming language: Python/C++.Computer: Any single laptop computer or cluster.Operating system: Any (Linux, OSX, Windows).RAM: a few Megabytes to several hundred GigabytesClassification: 12.External routines: FEniCS 1.3.0 (www.fenicsproject.org, that in turn depends on a number of external libraries like MPI, PETSc, Epetra, Boost and ParMetis)Nature of problem:Incompressible, Newtonian fluid flow.Solution method:The finite element method.Unusual features:FEniCS automatically generates and compiles low-level C++ code based on high-level Python code.Running time:The example provided takes a couple of minutes on a single processor.