[Sumit07]

Hi
I have combined buoyantPimpleFoam with pimpleDyMFoam to create new solver buoyantPimpleDyMFoam.
In buoyantPimpleFoam, I just modified buoyantPimpleFoam and pEqn to create buoyantPimpleDyMFoam.
In buoyantPimpleDyMFoam.c I added

mesh.update();
// Calculate absolute flux from the mapped surface velocity
phi = mesh.Sf() & Uf;
if (mesh.changing() && correctPhi)
{
#include "correctPhi.H"
}
// Make the flux relative to the mesh motion
fvc::makeRelative(phi, rho, U);
if (mesh.changing() && checkMeshCourantNo)
{
#include "meshCourantNo.H"
}
before Pressure-velocity PIMPLE corrector loop

and in pEqn I added
{
Uf = fvc::interpolate(rho*U);
surfaceVectorField n(mesh.Sf()/mesh.magSf());
Uf += n*(phi/(mesh.magSf()) - (n & Uf));
}
// Make the fluxes relative to the mesh motion
fvc::makeRelative(phi, rho, U);
at the end.

Also I change creatUf.H to compressibleCreateUf.H in which I changed fvc::interpolate(U) to fvc::interpolate(rho*U).

It will show good result until mesh start to move and get diverge. Velocity increase in the portion where mesh get stretched and reduced where mesh get compressed. I have attached result just before the solution get blow up.
Can some one tell me why this is happening ?

Regards

[alexj]

Hi Sumit,

I am currently also working on implementing a buoyantPimpleDyMFoam and got stuck with the same dimensions mismatch that you described in another post (http://www.cfd-online.com/Forums/openfoam-programming-development/154495-buoyantpimpledymfoam-dimension-mismatch-between-uf-phi.html).

Well I did not manage to finish the work yet, say my solver compiles but does not work. Did you progress with your work to the level that the solver was usable. If so would you be willing to share some details either via personal email or here on the forum (surely preferred by the community).

My guess why your solver diverges is that somehow the fluxes do not get corrected during the deforming/moving mesh, did you implement that?

Any discussion and help would be highly appreciated by me.
Kind regards,
Alex