Incompact3d
  • Home
  • About
  • Contact
  • Download
  • Impact
Picture

About Incompact3d

Incompact3d is a powerful high-order flow solver for academic research. Dedicated to Direct and Large Eddy Simulations (DNS/LES), it can combine the versatility of industrial codes with the accuracy of spectral codes. It scales with up to one million cores. 
The incompressible Navier-Stokes equations are discretized with finite-difference sixth-order schemes on a Cartesian mesh.  Explicit or semi-implicit temporal schemes can be used for the time advancement depending on the flow configuration. To treat the incompressibility condition, a fractional step method requires to solve a Poisson equation. This equation is fully solved in spectral space via the use of relevant 3D Fast Fourier transforms(FFTs), allowing any kind of boundary conditions for the velocity field in each spatial direction. Using the concept of the modified wavenumber, the divergence free condition is ensured up to machine accuracy.  The pressure field is staggered from the velocity field by half a mesh to avoid spurious oscillations.  
The modelling of a solid body inside the computational domain is performed with a customised Immersed Boundary Method. It is based on a direct forcing to ensure a no-slip boundary condition at the wall of the solid body while imposing non-zero velocities inside the solid body to avoid discontinuities on the velocity field. This customised IBM, fully compatible with the 2D domain decomposition and with a possible mesh refinement at the wall, is based on a 1D expansion of the velocity field from fluid regions into solid regions using Lagrange polynomials. 
To reach realistic Reynolds numbers, an implicit LES strategy can be implemented to solve the Navier-Stokes equations without any extra explicit modelling. In order to mimic a subgrid-scale model, artificial dissipation can be added via the viscous term thanks to the artificial dissipative features of the high-order compact schemes.

Acknowledgement
Sylvain Laizet and Eric Lamballais would like to thank Imperial College London and The University of Poitiers for agreeing to make the solver Incompact3d available for the scientific community. Christos Vassilicos, Sylvain Lardeau, Ning Li, Charles Moulinec, Véronique Fortuné, Jorge Silvestrini and Thibault Dairay are acknowledged for helping in the development of Incompact3d. 

About the numerical methods

Laizet S. & Lamballais E., High-order compact schemes for incompressible flows: a simple and efficient method with the quasi-spectral accuracy, J. Comp. Phys., vol 228-15, pp 5989-6015, 2009

Lamballais E., Fortune V. & Laizet S., Straightforward high-order numerical dissipation via the viscous term for Direct and Large Eddy Simulation, J. Comp. Phys., Vol 230-9, pp 3270-3275, 2011

Gautier R., Laizet S. & Lamballais E.,  A DNS  study of jet control with microjets using an alternating direction  forcing strategy, Int. J. of Computational Fluid  Dynamics, 28, pp 393-410, 2014

About the parallel strategy

Laizet S. & Li N., Incompact3d, a powerful tool to tackle turbulence problems with up to 0(10^5) computational cores, Int. J. of Numerical Methods in Fluids, Vol 67-11, pp 1735-1757, 2011

Li N. & Laizet S., 2DECOMP&FFT - a highly scalable 2D decomposi on library and FFT interface, Cray User Group meeting: Simulation comes of age, Edinburgh, Scotland -- 24/05-27/05, 2010
​
Laizet S., Lamballais E. & Vassilicos J.C., A numerical strategy to combine high-order schemes, complex geometry and parallel computing for high resolution DNS of fractal generated turbulence, Computers & Fluids, vol 39-3, pp 471-484, 2010

Incompact3d is using the open-source library 2D Decomp&FFT, available at http://www.2decomp.org

Proudly powered by Weebly