[Kojote]

Hi,

i want to write a solver which includes LPT and dynamic mesh. And i faced an issue. If there exist a particle which is not in the Domain anymore, the simulation crashes. Is there already a fix for that, even if the particle is just removed i could live with that, i only want to avoid the crashing.

My code looks like this (3.1-ext)

while (runTime.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;
#include "readControls.H"
#include "checkTotalVolume.H"

// Make the fluxes absolute
fvc::makeAbsolute(phi, U);

#include "CourantNo.H"
#include "setDeltaT.H"
{

bool meshChanged = mesh.update();
reduce(meshChanged, orOp<bool>());
# include "volContinuity.H"
if (correctPhi && (mesh.moving() || meshChanged))
{
// Fluxes will be corrected to absolute velocity
// HJ, 6/Feb/2009
# include "correctPhi.H"
}
// Make the fluxes relative to the mesh motion
fvc::makeRelative(phi, U);
if (mesh.moving() && checkMeshCourantNo)
{
# include "meshCourantNo.H"
}

Info<< "Evolving kinematicCloud1" << endl;
kinematicCloud1.evolve();
kinematicCloud1.info();


#include "UEqn.H"

// --- PISO loop
for (int corr=0; corr<nCorr; corr++)
{
#include "pEqn.H"
}
}