I have finally completed the implementation and the testing phases of a two-fluid code for the simulation of gas particle flows based on OpenFOAM®. The code is derived from the existing twoPhaseEulerFoam, with significant differences in the solution algorithm in order to stabilize it, and with the addition of a variety of models. This is the evolution of the twoPhaseEulerPimpleFoam code discussed previously here and here. I decided to change name to the code, and call it fluidParticleTwoPhase4Foam, since the implementation has become more and more specific for particulate flows, and future developments will move further in this direction. A summary of the improvements is reported below: Solution algorithm Adoption…

The development of twoPhaseEulerPimpleFoam proceeds well, and a good number of milestones, including Implementation of the conservative form of the momentum equation Implementation of the partial elimination algorithm Stabilization of the algorithm in the dense limit Boundary conditions for the particle phase Electrostatic model have been completed, and are under testing (they are not released yet to the public). Further information is available on the Extend project (You need an account to see this page). At the moment the only released development is represented by the implementation of the Johnson and Jackson boundary conditions, developed by Juho Peltola and me. They are available in my git repository. These boundary conditions…

OpenFOAM® 1.7.x (and previous versions) provides a wide range of LES SGS models, however the implementation of the models based on the dynamic procedure, originally proposed by Germano, is not suitable for flows in arbitrary geometries, since the coefficient, determined dynamically, is averaged over the whole computational domain. I implemented a modified version of the dynamic model, based on the work of Lilly (1992), and I made it available through my git repository for those interested. Please, refer to the README file for instructions. This offering is not approved or endorsed by OpenCFD Limited, the producer of the OpenFOAM software and owner of the OPENFOAM®  and OpenCFD®  trade marks.

I have been working for some time on the development of an improved version of the Euler-Euler two-fluid solver in OpenFOAM®, in order to improve its stability when a dense granular phases is present. An initial version of the solver, derived from twoPhaseEulerFoam, with the following modifications Adoption of the PIMPLE (PISO/SIMPLE) algorithm instead than PISO, to allow equation under-relaxation and ensure the convergence of all the equations at each time-step Solution of the momentum predictor, accounting for the strongly varying terms using flux-reconstruction is available on github. You can download it by typing git clone git://github.com/AlbertoPa/twoPhaseEulerPimpleFoam.git and keep it up to date typing git pull in a terminal opened…