Dear Everyone,
I am writing a new solver that couples interfoam with cfdemsolverIB solver.
I have followed this paper "https://pubs.acs.org/doi/abs/10.1021/acs.iecr.0c00108" to build the solver.
In the interFoam solver, the pressure is not divided to rho, so I have defined the p_rgh with the dimension of kg/m.s2.
in the cfdemCloudIB.C there is a line to correct pressure as follows:
"p=p+phiIB/U.mesh().time().deltaT(); "
Because p=(p/rho) in the original code, i have changed it to "p=p+rho*phiIB/U.mesh().time().deltaT();" to account for rho
I complied with no error .
but when running the program shows that the dimension of both side is not compatible as follows:
LHS and RHS of - have different dimensions
dimensions : [0 1 -2 0 0 0 0] - [1 -2 -2 0 0 0 0]
Any help will be much appreciated.