Skip to main content

A Scipy-Fenics interface for incompressible Navier-Stokes

Project description

# dolfin_navier_scipy


This python module dns provides an interface between the FEM toolbox [FEniCS]( and [SciPy]( in view of simulation and control of incompressible flows. Basically, FEniCS is used to discretize the incompressible Navier-Stokes equations in space. Then dns makes the discretized operators available in SciPy for use in model reduction, simulation, or control and optimization.

dns also contains a solver for the steady state and time dependent problems.

## Quick Start

To get started, create the needed subdirectories and run one of the tests/ files, e.g.

` cd tests mkdir data mkdir results # export PYTHONPATH="$PYTHONPATH:path/to/repo/" # add the repo to the path python3 `

Then, to examine the results, launch ` paraview results/vel_TH__timestep.pvd `

## Test Cases and Examples

A selection:

  • tests/ a minimal setup for a steady-state simulation

  • tests/ the 2D steady-state cylinder wake benchmark by Schäfer/Turek

  • tests/ the 2D cylinder wake with a freely rotating cylinder as benchmarked in Richter et al.

  • tests/ time integration with Picard and Newton linearization

  • tests/ time integration with explicit treatment of the nonlinearity

  • tests/ time integration of the cylinder wake with boundary controls

  • tests/ time integration with iterative solves of the state equations via [krypy](

## Dependencies

The latter is my home-brew module that includes the submodule lin_alg_utils with routines for solving the saddle point problem as it arises in the (v,p) formulation of the NSE.

Note: the branch lau-included already contains the module sadptprj_riclyap_adi

## Documentation

Documentation of the code goes [here](

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

dolfin_navier_scipy-1.0.0.tar.gz (36.7 kB view hashes)

Uploaded Source

Built Distribution

dolfin_navier_scipy-1.0.0-py3-none-any.whl (51.0 kB view hashes)

Uploaded Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page