[Mat_fr]

Hi foamers,


Everything's in the title.


I was wondering if there's a way to stop an openfoam simulation with a physical criterion (e.g. average pressure on a wall exceeding a threshold), using the dictionaries. And then to write the results, at this last time step.


Otherwise, I suppose I will have to modify the solver.


Thanks,
Mat

[oxadampf]

Hi Mat,


were you able to find a solution to your problem?


Thanks,
Ben

[Mat_fr]

Hello,


No I didn't.

For the moment, I am recording some physical values and looking for which time my physical criterion is achieved. Then if I need I start again the simulation at a previous saved time step untill the criterion time I just found.



It's a bit anoying and takes a bit of time, but it's okay.


Good luck,
Mat

[oxadampf]

Hi,

I think I found a solution. Writing a coded functionObject to the controlDict file can help. Something like this:



T__max__sample
{
enabled true;
type coded;
functionObjectLibs ("libutilityFunctionObjects.so");
region sample00;
redirectType T__max__sample;
writeControl timeStep;
writeInterval 1;
log true;

codeExecute
#{
const volScalarField& T = mesh().lookupObject<volScalarField>("T");
if ( max(T).value() < 270 )
{
bool result=const_cast<Time&>(mesh().time()).writeNow() ;
FatalErrorIn("Stopped")
<< "Bla" << endl
<< exit(FatalError);
}
#};
}


It terminates the simulation, if the max temperature in the region called "sample00" falls below 270 with a fatal error message (I know not pretty, but functional). Then, it writes the last timeStep.



This link helped me a lot:

Mixing #codeStream with coded-function-object


I hope this will help you.
Ben

[Mat_fr]

That's grand! Thanks for your time, Ben.


I haven't made it to the step of conded function objetc yet, and I should.


Regards,
Mat