[JM27]

Hi there,


I am running some simulations in compressibleInterFoam and I want to make the solver more robust by preventing the temperature from becoming negative.



After some research on how this can be done, I have tried using the elegant fvOptions function to limit the temperature in the domain. I am using OpenFOAM v5.



The fvOptions file I set is as follows (values are arbitrary, just want to test functionality with compressibleInterFoam):

 limitT
    {
        type                  limitTemperature;
        active               yes;     
        selectionMode   all;   
        min                   200;            
        max                  20000;   
    }

Upon trying to run this I get the following error:

Selecting finite volume options model type limitTemperature
    Source: limitT
    - selecting all cells
    - selected 571536 cell(s) with volume 9e-07
[0] 
[0] 
[0] --> FOAM FATAL ERROR: 
[0] Not implemented
[0] 
[0]     From function const Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> &Foam::twoPhaseMixtureThermo::he() const
[0]     in file twoPhaseMixtureThermo.H at line 138.
[0] 
FOAM parallel run aborting
[0] 
[2] 
[2] 
[2] --> FOAM FATAL ERROR: 
[2] Not implemented
[2] 
[2]     From function const Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> &Foam::twoPhaseMixtureThermo::he() const
[2]     in file twoPhaseMixtureThermo.H at line 138.
[2] 
FOAM parallel run aborting
[2] 
[1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] Not implemented
[1] 
[1]     From function [0] #0  Foam::error::printStack(Foam::Ostream&)const Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> &Foam::twoPhaseMixtureThermo::he() const
[1]     in file twoPhaseMixtureThermo.H at line 138.
[1] 
FOAM parallel run aborting
[1] 
[1] #0  Foam::error::printStack(Foam::Ostream&)[2] #0  Foam::error::printStack(Foam::Ostream&)[3] 
[3] 
[3] --> FOAM FATAL ERROR: 
[3] Not implemented
[3] 
[3]     From function const Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> &Foam::twoPhaseMixtureThermo::he() const
[3]     in file twoPhaseMixtureThermo.H at line 138.
[3] 
FOAM parallel run aborting
[3] 
[3] #0  Foam::error::printStack(Foam::Ostream&) at ??:?

I understand that the error is in the twoPhaseMixtureThermo however I cannot seem identify the cause of the problem yet. The phases involved are air and water. I use perfectGas (compressible) to model air and rhoConst for water (incompressible)...



Could this be the cause of the problem ? Does limitTemperature function not work with incompressible phases?


I will keep on looking into this and post the solution if found, however any help will be greatly appreciated!



Thanks

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

Note:

1. I tried running the case with fvOptions specified in /constant or in /system, however this does not seem to make a difference (as expected).



2. I have been through some tutorials to understand the fvOptions functionality however there are none available on limitTemperature.


3. I have double-checked the pEqn.H, UEqn.H and TEqn.H files in compressibleInterFoam allow for the fvOptions functionality.