Python FEM and Multiphysics Simulations with FEniCS and FEATool
FEniCS is a flexible and comprehensive finite element FEM and partial differential equation PDE modeling and simulation toolkit with Python and C++ interface...
In computational fluid dynamics (CFD) simulations accuracy is often very important, especially for drag and lift values. References [1-3] describe a CFD benchmark problem for time-dependent flow around a cylinder where the Reynolds number varies between 0<=Re<=100. The references contain very accurate reference values for the drag and lift coefficients, and pressure difference between the front and rear of the cylinder.
This benchmark problem has been simulated with FEATool with several grid and time step sizes. The results are shown in the table below.
As can be seen the computed values on the finest grid and time step are very accurate and close to the reference solution. The time evolution of the lift coefficient is the most difficult quantity to capture correctly, but from the figure below the computed solution is again very close to the reference.
The FEATool model script file to run the time-dependent cylinder benchmark can be found in the FEATool examples directory or from the link below. Note that it might take quite some time for the simulation to finish depending on the choice of grid and time step sizes.
FEATool Instationary Flow Around Cylinder CFD Model Tutorial
 Schaefer M, Turek S. Benchmark computations of laminar flow around a cylinder. In Flow Simulation with High-Performance Computers II, Notes on Numerical Fluid Dynamics, vol. 52. Vieweg, 1996; 547-566.
 John V. Reference values for drag and lift of a two-dimensional time-dependent flow around a cylinder. International Journal for Numerical Methods in Fluids 2004; 44:777-788.
 John V, Rang J. Adaptive time step control for the incompressible Navier–Stokes equations. Comput. Methods Appl. Mech. Engrg. 199 (2010) 514–524.