[chenbo]

Hi guys,
I'm trying to merge the solver DPMFoam with reactingTwoPhaseEulerFoam.
However, I don't know how to access basicKinematicCollidingCloud.theta() in twoPhaseSystem::solve().
In the other thread

Access particle data from functionObject

They suggested using

[CODE] const kinematicCloud& myCloud = obr_lookupObject<kinematicCloud> ("kinematicCloud"); [CODE]

However, when I include "#include basicKinematicCollidingCloud.H" in twoPhaseSystem/twoPhaseSystem.C and add that code, I get errors:

[CODE] error: 'obr_' was not declared in this scope [CODE]

So, can anyone give me some hints about how to change the code to let me get the value of "kinematicCloud.theta()" in twoPhaseSystem::solve() function, just as the theta() could be calculated in DPMFoam?
Thanks!

[raumpolizei]

Hey chenbo,
well, it looks like the line you cited is missing a point (literally)

  const kinematicCloud& myCloud = obr_.lookupObject<kinematicCloud> ("kinematicCloud");

However, I am not sure if it is going to compile as obr_ has to be known within the class. What could be worth trying is this (mesh_ is a protected member of phaseSystem, which is the base class of twoPhaseSystem so it should be possible to access it):

#include "kinematicCloud.H"
/* ... */
const kinematicCloud& lagrangeCloud = this->mesh_.lookupObject<kinematicCloud> ("kinematicCloud");
auto lagrangeTheta = lagrangeCloud.theta();

Good luck with your task.
Cheers
RP

[chenbo]

Originally Posted by raumpolizei

Hey chenbo,
well, it looks like the line you cited is missing a point (literally)

  const kinematicCloud& myCloud = obr_.lookupObject<kinematicCloud> ("kinematicCloud");

However, I am not sure if it is going to compile as obr_ has to be known within the class. What could be worth trying is this (mesh_ is a protected member of phaseSystem, which is the base class of twoPhaseSystem so it should be possible to access it):

#include "kinematicCloud.H"
/* ... */
const kinematicCloud& lagrangeCloud = this->mesh_.lookupObject<kinematicCloud> ("kinematicCloud");
auto lagrangeTheta = lagrangeCloud.theta();

Good luck with your task.
Cheers
RP

Dear raumpolizei,
Thanks for your kindly reply.
1. I did miss a "." in the original post, but I input it right in the code, there's still several errors while compiling, because, as you said, "obr_" is not callable (the compiler suggested using "db_", which was not working either).

2. I did your code, now it compiles well, no errors occur, thanks a lot, really appreciate your help. My research time may be halved by your help.

The forum is really a nice and warm community. I wish I could know OpenFOAM as much as you guys do, so I can post replies to help others.

Now I'll go make some experiments based on the code.
Thanks again. Have a nice day.

[Jasper Z]

Dear chenbo,

I know this thread is a long time ago, but I am also coupling these two solvers these days. So I wonder did you couple these two solvers to simulate three-phase flow? where two phases used Euler and the other used lagrangian method? If so, how did you calculate the alpha1 and alpha2 in twoPhaseSystem.C? Are they related to 1 - kinematicCloud.theta()? Thanks for your help, I am really struggled with this recently

Kind regards,
Jasper