Reply To: Fluid velocity in dpm_drag

[Chris_g]

I was not aware that you could use DPM without an Eulerian phase. I now tried:

Domain *gas_domain;

gas_domain = Get_Domain(2); //Domain ID from multiphase interface

Thread *T0 = Lookup_Thread(gas_domain, 56); //56 from boundary condition Zone-ID for internal Zone

cell_t c = TP_CELL(tp);

real uvel = C_U(c,t0);

This doesn't work, but I am not sure at all if this correct.

 

However in the meantime I found a work-around that seems to work:

cphase_state_t *cell = *(&(tp->cphase));

real uvel = cell->V[1];

In a short test this seems to give plausible values. There is no documentation for this function (I found it in the dpm_types.h file) so it would be nice if you could confirm, that those are the values for the fluid phase. If so I should be able to get all my required values this way.

 

So in the end while it is still interesting what is wrong with my original approach, right now it seems like the problem is solved for me.

Thanks for your help!