[asfik]

Dear All,

I would like to add dynamic mesh refinement with MPPICInterFoam Solver(v2206) and did the following task:

a. Modified the .C file according to InterFoam/InterDyMFoam
b. Added dynamicRefine2DFvMesh @dynamicMeshDict

Solver Complied successfully, tried to simulate the dam brake 2D refinement tutorial along with dynamicRefine2DFvMesh.
the solver refined the mesh for a few timesteps and then got the error which is coming from PEqn.H file of the solver.


*************************************************
fvScalarMatrix p_rghEqn
(
fvm::laplacian(Dp, p_rgh)
==
fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA)
);

p_rghEqn.setReference(pRefCell, getRefCellValue(p_rgh, pRefCell));

p_rghEqn.solve(mesh.solver(p_rgh.select(pimple.fin alInnerIter())));

*** Works until this stage then got singularity from following part***********


if (pimple.finalNonOrthogonalIter())
{
phi = phiHbyA - p_rghEqn.flux()/alphacf;

p_rgh.relax();

U =
HbyA
+ rAUc*fvc::reconstruct((phig - p_rghEqn.flux()/alphacf)/rAUcf);

U.correctBoundaryConditions();

fvOptions.correct(U);
}
}


**************************************************

Please show me some light of hope.
Thanks