Hi all, I am working on multi-phase euler foam problem. I have implemented population balance module in the multiPhaseEulerFoam model. Successfully compiled the new solver, but the case file is not running. In fact it is not showing any error. It is showing... Build : 2.2.2-9240f8b967db Exec : multiPhaseEulerFoam Date : Jan 17 2017 Time : 23:10:13 Host : "sarojdip-OptiPlex-7010" PID : 9546 Case : /home/sarojdip/OpenFOAM/sarojdip-2.2.2/run/testCase nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading g Reading transportProperties Calculating face flux field phia Calculating face flux field phib Calculating face flux field phia1 Calculating face flux field phia2 Calculating face flux field phia3 Calculating face flux field phia4 Reading field alpha Reading field alpha1 Reading field alpha2 Reading field alpha3 Reading field alpha4 Reading field beta --> FOAM Warning : From function GeometricField<Type, PatchField, GeoMesh>::readIfPresent() in file /opt/openfoam222/src/OpenFOAM/lnInclude/GeometricField.C at line 108 read option IOobject::MUST_READ or MUST_READ_IF_MODIFIEDsuggests that a read constructor for field beta would be more appropriate. Reading field p Reading field k Reading field epsilon Calculating field nutb Calculating field nuEffa Calculating field nuEffa2 Calculating field nuEffa3 Calculating field nuEffa4 Calculating field nuEffb Calculating field DDtUa DDTUb DDtUa1 DDtUa2 DDtUa3 DDtUa4 Calculating field g.h Selecting dragModel for phase a: Tsuchiya Selecting dragModel for phase b: Tsuchiya Selecting dragModel for phase a1: Tsuchiya Selecting dragModel for phase a2: Tsuchiya Selecting dragModel for phase a3: Tsuchiya Selecting dragModel for phase a4: Tsuchiya dragPhasea is a dragPhasea1 is a1 dragPhasea2 is a2 dragPhasea3 is a3 dragPhasea4 is a4 liftModel is constant liftPhase is b Selecting viscosityModel Syamlal Selecting conductivityModel HrenyaSinclair Selecting radialModel Gidaspow Selecting granularPressureModel Lun Selecting frictionalStressModel JohnsonJackson Courant Number mean: 0 max: 0 PIMPLE: Operating solver in PISO mode Starting time loop fieldAverage fieldAverage1: Starting averaging at time 0 Courant Number mean: 0 max: 0 Max Ura Courant Number = 0 Max Ura1 Courant Number = 0 Max Ura2 Courant Number = 0 Max Ura3 Courant Number = 0 Max Ura4 Courant Number = 0 deltaT = 1.19999e-05 Time = 1.19999e-05 ............The simulation stopped immediately here.................. If any body has Population balance code implemented in OpenFOAM, then help me by sharing. I'll provide my email-ID in my next thread. Any help would help me in a great way... Thank you in advance Saroj..