[cute]

Hi All,

I want to stop my simulation when C < tolerance_Value or C < 0 at any grid point. I don't know, how to write the condition in while loop? I need your help, code is given below.

while (????)    
   {        
         Info<< "Time = " << runTime.timeName() << nl << endl;               

         forAll(CX, cellI) 
            {   
                   CX[cellI] = alpha.value()*C[cellI]*C[cellI];            
            }                          

         solve(fvm::ddt(C) == fvm::laplacian(DT, C) - k*CX*C);   

         runTime.write();  
    }

[anishtain4]

You can use panicDump during the run rather than putting it in the solver.

[Jerryfan]

Hi cute,


You can simply add a line that gives the maximum of C, and then determine based on this maximum value.

while (maxC > tolerance_value)
{
Info<< "Time = " << runTime.timeName() << nl << endl;

forAll(CX, cellI)
{
CX[cellI] = alpha.value()*C[cellI]*C[cellI];
}

dimensioned<Type> maxC = max(C); //Type depends on the Type of variable in field C.

solve(fvm::ddt(C) == fvm::laplacian(DT, C) - k*CX*C);

runTime.write();
}

[cute]

Thanks for the input.
I defined C as volScalarField, is my Type equals volScalarField?

dimensioned<volScalarField> maxC = max(C);

But all the values in C are floating point numbers, I am confused what to write?

[Jerryfan]

What I mean is Type of variable in field C, not the Type of field C itself. In this case, you should use

dimensioned<scalar>

.

[cute]

Thanks for a quick response.

[cute]

I am using the following code as suggested, but I am getting error:

CODE:

dimensioned<scalar> maxC = max(C);
dimensioned<scalar> minC = min(C);
    
    while (maxC > 1E-6 && minC > 0.0)
    { -------------------- }

ERROR:

 error: no match for ‘operator>’ (operand types are ‘Foam::dimensioned<double>’ and ‘double’)
     while (maxC > 1E-6 && minC > 0.0)

[Jerryfan]

I am using the following code as suggested, but I am getting error:

CODE:
Code:
dimensioned<scalar> maxC = max(C);
dimensioned<scalar> minC = min(C);

while (maxC > 1E-6 && minC > 0.0)
{ -------------------- }
ERROR:

Code:
error: no match for ‘operator>’ (operand types are ‘Foam::dimensioned<double>’ and ‘double’)
while (maxC > 1E-6 && minC > 0.0)

You can try to get the value out from the dimensionedType variable by:


while (maxC.value() > 1E-6 && minC.value() > 0.0)

[cute]

Thanks, it works.