[jiejie]

Hi foamers

I am trying to implement a two-way coupling lagrangian solver by simply modifying the icoUncoupledKinematicParcelFoam.

The solver doesn't solve the flow field at the moment, I am thinking of adding back the flow solving routine from the icoFoam back into it in order for it to solve the flow field.

Regarding to the two-way coupling, it seems the KinematicParcel class has the ability of doing one/two-way coupling with the continuous phase from its description. Since the SU is calculated in KinematicCloudI.H. Can I just do the following to add the momentum source back to the flow solver?

fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
- fvm::laplacian(nu, U)
==
kinematicCloud.SU()/rho
);

Any help will be appreciated.

Jiejie