pyCALC-LES
pyCALC-LES is written in Python (3.12). The code solves the incompressible momentum equations, the continuity equation and transport equations for modeled turbulent quantities such as k, eps and omega. The density is assumed to be constant and equal to one, i.e. rho = 1. The transport equations are solved in 3D and the grid may be curvi-linear in the x-y plane. In the z direction the grid is Cartesian but dz may vary. The code is suitable for DNS, LES or DES (hybrid LES-RANS). For LES, the Smagorinsky model and the WALE model are implemented. For DES, a k-omega DES model and a PANS k-eps and IDDES model are implemented. pyCALC-LES is a finite volume code. It is fully vectorized (i.e. no for loops). The solution procedure is based on fractional step. Second-order central differencing is used in space and the Crank-Nicolson scheme in time. The discretized equations are solved with Pythons sparse matrix solvers. For users who have an Nvidia graphics card, the entire code runs on the GPU. Using IDDES on a 2.3M mesh the speed-up is a factor of 70 on the GPU compared to the CPU.
Department of Mechanics and Maritime Sciences Division of Fluid Dynamice |