[FBenz]

Dear OpenFoam Users,


I am trying to implement a subroutine inside the main time loop which looks something like:



while(runTime.loop())
{
solve first equation
while (some convergence criteria)
{
solve second equation
}
}


The solver compiles and the case is working, but it seems to be that openfoam does not like it if I call solve over and over again in the same time step.



If I run the two pieces of code:

1:


while (runTime.loop())
{
dU *= (0.0,0.0,0.0);
solve
(
fvm::laplacian(mu,dU) + fvc::div(sigma) //mu, dU and sigma are just dummy fields
);
prevTime = time;
time = runTime.elapsedCpuTime();
Info << "Time needed : " << time - prevTime << nl;
}


and
2:


while(true)
{
dU *= (0.0,0.0,0.0);
solve
(
fvm::laplacian(mu,dU) + fvc::div(sigma)
);
prevTime = time;
time = runTime.elapsedCpuTime() ;
Info << "Time needed : " << time - prevTime << nl;
}


the first code will give me a constant "Time needed" as expected, while the second will increase with every iteration.
I guess runTime.loop() will kind of reset the field or the objectRegistry?


Is there a way to circumvent this problem or to access the registry manually?


Thanks in advance
Felix