[shipman]

Hi to all,

I am quite new about implementation and i want to change the kinematic eddy viscosity equation (nut) in k-epsilon turbulence model which is defined:

nut_ = Cmu_*sqr(k_)/epsilon_;

i want to change it as:

nut_ = f(rho)*Cmu_*sqr(k_)/epsilon_;

f(rho) = rho1+ (rho2-rho)/[(rho2-rho1)]^5

here f is function based on the mixture density (rho) which is calculated inside of interphasechangefoam connected with the liquid volume fraction. rho1 and rho2 are densities of phases and constant so there is no problem. However, i just dont know how can i define rho inside the function of f(rho) which can use the calculated mixture density (rho) at every time step in the turbulence equation.

Your help will be appreciated.

Thanks in advance.

B

[alexeym]

Hi,

maybe it'll be easier if you just multiply nut field returned from turbulence model by f(rho) inside the solver?

[Bernhard]

Originally Posted by alexeym

Hi,

maybe it'll be easier if you just multiply nut field returned from turbulence model by f(rho) inside the solver?

I think this is the best advice in terms of the "First get it to work, then make it nice"-spirit

[shipman]

Dear Alexey

Thanks so much for your advice. i have 2 questions.

1. When i look at the U.Eqn. in interphasechangeFoam:

 surfaceScalarField muEff
    (
        "muEff",
        twoPhaseProperties->muf()
      + fvc::interpolate(rho*turbulence->nut())
    );

I really cant understand why the turbulence viscosity multiplied with rho before interpolation. Could you tell me reason?

2. then i inserted my equation here as follows:

surfaceScalarField muEff
    (
        "muEff",
        twoPhaseProperties->muf()
      + fvc::interpolate(double Foam::pow((rho2-rho)/(rho2-rho1))*(rho1-rho2),5.0)+ rho2)*rho*turbulence->nut())
    );

however, it gives following error:

dency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:

UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: no matching function for call to ‘pow(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >)’
/opt/openfoam211/src/OpenFOAM/lnInclude/label.H:257: note: candidates are: Foam::label Foam::pow(Foam::label, Foam::label)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:78: note:                 double Foam::pow(double, double)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:79: note:                 double Foam::pow(double, float)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:80: note:                 double Foam::pow(float, double)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:81: note:                 double Foam::pow(float, float)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:82: note:                 double Foam::pow(double, int)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:83: note:                 double Foam::pow(int, double)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:84: note:                 double Foam::pow(double, long int)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:85: note:                 double Foam::pow(long int, double)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:86: note:                 double Foam::pow(float, int)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:87: note:                 double Foam::pow(int, float)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:88: note:                 double Foam::pow(float, long int)
/opt/openfoam211/src/OpenFOAM/lnInclude/doubleFloat.H:89: note:                 double Foam::pow(long int, float)
/opt/openfoam211/src/OpenFOAM/lnInclude/dimensionSet.H:226: note:                 Foam::dimensionSet Foam::pow(const Foam::dimensionSet&, Foam::scalar)
/opt/openfoam211/src/OpenFOAM/lnInclude/dimensionSet.H:227: note:                 Foam::dimensionSet Foam::pow(const Foam::dimensionSet&, const Foam::dimensionedScalar&)
/opt/openfoam211/src/OpenFOAM/lnInclude/dimensionSet.H:228: note:                 Foam::dimensionSet Foam::pow(const Foam::dimensionedScalar&, const Foam::dimensionSet&)
/opt/openfoam211/src/OpenFOAM/lnInclude/dimensionedScalar.H:58: note:                 Foam::dimensionedScalar Foam::pow(const Foam::dimensionedScalar&, const Foam::dimensionedScalar&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note:                 void Foam::pow(Foam::Field<double>&, const Foam::UList<double>&, const Foam::UList<double>&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::UList<double>&, const Foam::UList<double>&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::UList<double>&, const Foam::tmp<Foam::Field<double> >&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::tmp<Foam::Field<double> >&, const Foam::UList<double>&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::tmp<Foam::Field<double> >&, const Foam::tmp<Foam::Field<double> >&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note:                 void Foam::pow(Foam::Field<double>&, const Foam::scalar&, const Foam::UList<double>&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::scalar&, const Foam::UList<double>&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::scalar&, const Foam::tmp<Foam::Field<double> >&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note:                 void Foam::pow(Foam::Field<double>&, const Foam::UList<double>&, const Foam::scalar&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::UList<double>&, const Foam::scalar&)
/opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note:                 Foam::tmp<Foam::Field<double> > Foam::pow(const Foam::tmp<Foam::Field<double> >&, const Foam::scalar&)
UEqn.H:6: error: expected ‘,’ or ‘;’ before ‘)’ token
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wmake
Making dependency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: expected primary-expression before ‘double’
UEqn.H:5: error: expected ‘)’ before ‘double’
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wmake
Making dependency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: expected primary-expression before ‘double’
UEqn.H:5: error: expected ‘)’ before ‘double’
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wmake
Making dependency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: expected primary-expression before ‘double’
UEqn.H:5: error: expected ‘)’ before ‘double’
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wmake
Making dependency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: expected primary-expression before ‘double’
UEqn.H:5: error: expected ‘,’ or ‘;’ before ‘+’ token
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wmake
Making dependency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: expected primary-expression before ‘double’
UEqn.H:5: error: expected ‘,’ or ‘;’ before ‘*’ token
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1

Could you help me how can i correct it?

thanks

B

[alexeym]

1. mu is dynamic viscosity, nu is kinematic viscosity, mu = rho*nu. nut is kinematic turbulent viscosity.

2.

fvc::interpolate(double Foam::pow((rho2-rho)/(rho2-rho1))*(rho1-rho2),5.0)+ rho2)*rho*turbulence->nut())

2.1. You don't need that double in front of Foam:ow
2.2. You've got problems with brackets. According to your initial post the expressions should be

fvc::interpolate((rho1 + (rho2 - rho)*Foam::pow(rho2 - rho1, -5))*rho*turbulence->nut())

2.3. If rho1 and rho2 have dimensions (i.e. they are dimensionedScalar), even if the statement compile, you'll get dimensions error during run time. So the expression:

(rho1 + (rho2 - rho)*Foam::pow(rho2 - rho1, -5))

should be rewritten as something like:

(rho1.value() + (rho2.value() - rho.field())*Foam::pow((rho2 - rho1).value(), -5))

(I'm not sure that all calls in the expression are correct, you should check documentation)

i.e. make your f(rho) non-dimensional.

[shipman]

Dear Alexey

Yeah you are right i just mixed with kinematic (nu) . and dynamic (mu) viscosities.

As you know that rho1 and rho2 are constant. In order to avoid to mis write directly i wrote their values and defined the equation as follow( by the way at post #1 i wrote incorrect f(rho) true one is as follows):

surfaceScalarField muEff
    (
        "muEff",
        twoPhaseProperties->muf()
      + fvc::interpolate((Foam::pow((0.0017-rho.field()),5)/Foam::pow((0.0017-998),5)*(998-0.0017)+ 0.0017)*rho*turbulence->nut())
    );

just i tried to write rho as --> rho , rho.fields() , rho() for all the case it gives following error:

Making dependency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: no matching function for call to ‘interpolate(Foam::tmp<Foam::Field<double> >)’
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean
yusuke@VENUS:~/OpenFOAM/yusuke-2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wmake
Making dependency list for source file RNGinterPhaseChangeFoam.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Kunz/Kunz.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Merkle/Merkle.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C
SOURCE=RNGinterPhaseChangeFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam211/src/transportModels -I/opt/openfoam211/src/transportModels/incompressible/lnInclude -I/opt/openfoam211/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam211/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam211/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o
In file included from RNGinterPhaseChangeFoam.C:90:
UEqn.H: In function ‘int main(int, char**)’:
UEqn.H:5: error: no matching function for call to ‘interpolate(Foam::tmp<Foam::Field<double> >)’
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/RNGinterPhaseChangeFoam.o] Error 1

In addition, what you mean as documentation for checking.

Thanks in advance.

B

[alexeym]

Unfortunately I don't have 2.1.1 installation, so these solutions will be quite theoretical (in 2.2.2 and 2.3.0 UEqn.H is slightly different).

As you need to make rho inside your f(rho) non-dimensional, you can introduce unit density densionedScalar and divide rho by this scalar:

    dimensionedScalar rhoUnit("rhoUnit", dimDensity, 1.0);
    ...
    surfaceScalarField muEff
    (
        "muEff",
        twoPhaseProperties->muf()
      + fvc::interpolate((Foam::pow((0.0017-rho/rhoUnit),5)/Foam::pow((0.0017-998),5)*(998-0.0017)+ 0.0017)*rho*turbulence->nut())
    );

Btw in 2.2.2 these two lines of UEqn

      - fvm::laplacian(muEff, U)
      - (fvc::grad(U) & fvc::grad(muEff))

was changed (well, not exactly, there's a call to turbulence->divDevRhoReff(rho, U) which is the code below) to

    volScalarField muEff("muEff", rho*nuEff());

    ...
      - fvm::laplacian(muEff, U)
      - fvc::div(muEff*dev(T(fvc::grad(U))))

By checking documentation I meant going to http://www.openfoam.org/docs/cpp/ and checking if volScalarField indeed has field() method.

[shipman]

Dear Alexey,

Thanks so much for your help. It is running now. Could you give some advice how can i improve myself about this kind of small implementations.

On the other hand, T is transpose in the below equation?if not what is the meaning of it?

fvc::div(muEff*dev(T(fvc::grad(U))))

Thanks in advance.

B.

[alexeym]

Originally Posted by shipman

On the other hand, T is transpose in the below equation?if not what is the meaning of it?

fvc::div(muEff*dev(T(fvc::grad(U))))

Yes, T is transpose (http://foam.sourceforge.net/docs/cpp...a992a1c332bd5c), fvc::grad(U) will return you outer product of operator nabla and velocity (i.e. matrix, http://foam.sourceforge.net/docs/cpp...e60fccb38a70a9) and dev will return deviatoric part (http://foam.sourceforge.net/docs/cpp...362534e3433777) of it.

[shipman]

Dear Alexey,

Thanks so much for your advices. I will have one another problem. I tried to use my same case in the another PC. (which has same version 2.1.1) When i run "wmake", the statement compile however it gives following warning:

/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H: In function ‘int main(int, char**)’:
/opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable]

I checked the readTimeControls.H file but i didnt understand why it is not used??

const bool adjustTimeStep =
    runTime.controlDict().lookupOrDefault("adjustTimeStep", false);

scalar maxCo =
    runTime.controlDict().lookupOrDefault<scalar>("maxCo", 1.0);

scalar maxDeltaT =
    runTime.controlDict().lookupOrDefault<scalar>("maxDeltaT", GREAT);

Also, during run-time this situation results in decrasing the deltaT very small value like 10^-25 and gives error.

Could you tell me what is the mean reason of this problem?

Thanks in advance.

B

[alexeym]

Hi,

Cause you do not use directly maxDeltaT in your solver compiler complains about it. This warning appears during compilation of (almost) every OF solver. So you can ignore it.

Concerning your question about decreasing time step, usually it means that you've got diverging solution. And there're plenty of reasons for it: inconsistent initial conditions, inconsistent boundary conditions, wrong settings in fvSolution, too large initial time step, etc. To tell you something about this problem I need to see the case files.

[shipman]

Dear Alexey,

First of all thank you so much for your instructive informations. Related with time decreasing problems i have attached the my initial conditions. First, i decreased deltaT then changed many times k and epsilon values but still it gives following error:

DILUPBiCG:  Solving for alpha1, Initial residual = 0.00424273, Final residual = 1.84854e-06, No Iterations 1
Liquid phase volume fraction = 0.999545  Min(alpha1) = 0.191732  Max(alpha1) = 1
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00067342, Final residual = 4.82922e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.0577387, Final residual = 3.69088e-07, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 0.175927, Final residual = 0.000876531, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.00159558, Final residual = 8.87093e-05, No Iterations 1
time step continuity errors : sum local = 0.000646625, global = 2.71757e-05, cumulative = -0.0123716
GAMG:  Solving for p_rgh, Initial residual = 0.000603131, Final residual = 3.33279e-05, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 3.97107e-05, Final residual = 3.9419e-06, No Iterations 7
time step continuity errors : sum local = 0.000425533, global = 1.63888e-05, cumulative = -0.0123552
GAMG:  Solving for p_rgh, Initial residual = 0.000117561, Final residual = 9.08278e-06, No Iterations 1
GAMGPCG:  Solving for p_rgh, Initial residual = 1.36361e-05, Final residual = 7.04989e-08, No Iterations 9
time step continuity errors : sum local = 0.000420037, global = 1.71133e-05, cumulative = -0.0123381
ExecutionTime = 10435.8 s  ClockTime = 10636 s

Courant Number mean: 0.0024306 max: 0.935238
deltaT = 7.38217e-16
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 1010
    Increased the timePrecision from 13 to 14 to distinguish between timeNames at time 0.00157894
Time = 0.0015789419679163

DILUPBiCG:  Solving for alpha1, Initial residual = 0.00398756, Final residual = 2.17863e-06, No Iterations 1
Liquid phase volume fraction = 0.999531  Min(alpha1) = 0.191699  Max(alpha1) = 1
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0654812, Final residual = 8.47643e-09, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.489636, Final residual = 7.37747e-07, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 0.125157, Final residual = 0.00217892, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.31948e-06, Final residual = 8.27181e-08, No Iterations 2
time step continuity errors : sum local = 0.000709651, global = 1.68882e-05, cumulative = -0.0123212
GAMG:  Solving for p_rgh, Initial residual = 1.42336e-05, Final residual = 7.77456e-07, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 5.56115e-07, Final residual = 4.50147e-08, No Iterations 2
time step continuity errors : sum local = 0.000680165, global = 4.07245e-05, cumulative = -0.0122805
GAMG:  Solving for p_rgh, Initial residual = 7.13622e-07, Final residual = 3.20422e-08, No Iterations 1
GAMGPCG:  Solving for p_rgh, Initial residual = 6.03872e-08, Final residual = 6.03872e-08, No Iterations 0
time step continuity errors : sum local = 0.000682026, global = 4.08892e-05, cumulative = -0.0122396
ExecutionTime = 10437 s  ClockTime = 10637 s

Courant Number mean: 0.0039758 max: 0.746856
deltaT = 4.94217e-16
Time = 0.0015789419679168

DILUPBiCG:  Solving for alpha1, Initial residual = 0.00188925, Final residual = 1.1283e-06, No Iterations 1
Liquid phase volume fraction = 0.999457  Min(alpha1) = 0.191658  Max(alpha1) = 1
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000332758, Final residual = 6.84687e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.170049, Final residual = 5.61643e-07, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 0.0799552, Final residual = 7.82053e-05, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 4.30074e-12, Final residual = 4.30074e-12, No Iterations 0
time step continuity errors : sum local = 0.256009, global = -0.000611308, cumulative = -0.0128509
GAMG:  Solving for p_rgh, Initial residual = 2.64617e-09, Final residual = 2.64617e-09, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 2.64617e-09, Final residual = 2.64617e-09, No Iterations 0
time step continuity errors : sum local = 0.264605, global = -0.000615829, cumulative = -0.0134667
GAMG:  Solving for p_rgh, Initial residual = 2.6463e-09, Final residual = 2.6463e-09, No Iterations 0
GAMGPCG:  Solving for p_rgh, Initial residual = 2.6463e-09, Final residual = 2.6463e-09, No Iterations 0
time step continuity errors : sum local = 0.268305, global = -0.000638542, cumulative = -0.0141053
ExecutionTime = 10438.2 s  ClockTime = 10639 s

Courant Number mean: 0.411707 max: 39700.8
deltaT = 6.22426e-21
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 1010
    Increased the timePrecision from 14 to 15 to distinguish between timeNames at time 0.00157894
Time = 0.00157894196791677

[3] #0  Foam::error::printStack(Foam::Ostream&)[8] #0  Foam::error::printStack(Foam::Ostream&) in "/op in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int)t/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #2   in "/lib/libc.so.6"
[8] #3   in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2   in "/lib/libc.so.6"
[3] #3  
[3]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[3] #4  

[3]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[3] #5  [8]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[8] #4  
[3]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[3] #6  __libc_start_main
[8]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[8] #5  
 in "/lib/libc.so.6"
[3] #7  
[8]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[8] #6  __libc_start_main in "/lib/libc.so.6"
[8] #7  [3]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[VENUS:16057] *** Process received signal ***
[VENUS:16057] Signal: Floating point exception (8)
[VENUS:16057] Signal code:  (-6)
[VENUS:16057] Failing at address: 0x3e800003eb9
[VENUS:16057] [ 0] /lib/libc.so.6(+0x33af0) [0x7fed744abaf0]
[VENUS:16057] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7fed744aba75]
[VENUS:16057] [ 2] /lib/libc.so.6(+0x33af0) [0x7fed744abaf0]
[VENUS:16057] [ 3] RNGinterPhaseChangeFoam() [0x4a8442]
[VENUS:16057] [ 4] RNGinterPhaseChangeFoam() [0x4b74fd]
[VENUS:16057] [ 5] RNGinterPhaseChangeFoam() [0x435358]
[VENUS:16057] [ 6] /lib/libc.so.6(__libc_start_main+0xfd) [0x7fed74496c4d]
[VENUS:16057] [ 7] RNGinterPhaseChangeFoam() [0x42b849]
[VENUS:16057] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 16057 on node VENUS exited on signal 8 (Floating point exception).

Please let me know what is wrong with this case.

thanks in advance. After checking the simulation i will also send some pictures.

[shipman]

Hello Alexey,

Also you can see the vel, liquid volume fraction and pres. distributions at attached pics. I really dont understand why it gives this kind of distribution. it is very strange that at the inlet pressure is increasing even mass flow (at the inlet) and outlet velocity seem constant???. And also, as you can see that cavitation takes place as attached the wall..

Could you tell me your idea.

thanks in advance

B.

[alexeym]

Hi,

1. Where is inlet, where is outlet?
2. Why do you think that 1 outer corrector step is enough for the simulation to converge?
3. Are you sure about 50 as IC for epsilon?

[shipman]

Dear Alexey,

Thanks for reply.

1. Inlet is upstairs (long width) and down (small width) is outlet.

2. Actaully i dont know whether 1 outCorrecters enough or not? Could you advice what is the good setting? 2 or 3?

3. For epsilon i am not sure. According to initial calculation of epsilon based on k is 0.7 . I thought that it is very small. and before i read one post and one said that you can write 50~100 times higher for epsilon initial condition. So could you tell me your idea?

Thanks in advance..


B

[alexeym]

Well,

1. As usual I'd suggest increasing nOuterCorrectors up to 50 and use residualControl to exit PIMPLE loop, i.e. change current PIMPLE dictionary in fvSolution from

PIMPLE
{
    momentumPredictor          no;
    nOuterCorrectors           1;
    nCorrectors                3;
    nNonOrthogonalCorrectors   1;
        
    cAlpha                     0;
    nAlphaCorr                 1;
    nAlphaSubCycles            1;
}

to

PIMPLE
{
    momentumPredictor          no;
    nOuterCorrectors           50;
    nCorrectors                2;
    nNonOrthogonalCorrectors   1;
        
    cAlpha                     0;
    nAlphaCorr                 1;
    nAlphaSubCycles            1;

    turbOnFinalIterOnly no;

    residualControl
    {
        "(p_rgh|U|k|epsilon)"
        {
            tolerance 1e-2;
            relTol 0;
        }
    }

}

I also added 'turbOnFinalIterOnly no' so turbulent values will be recalculated on every iteration. Maybe it's not necessary though when I run simulations with phase change this addition leads to better convergence.

2. Well, you can keep epsilon=50 if you like but if your simulation is blowing up why not try modifying this value?

3. I assume there's no problem with the mesh. Is it just simple rectangular mesh? Did you estimate y+?

[shipman]

Hi dear Alexey,

Thank so much for your reply.

I started my calculation according to your advice. I will post results here when it finishes.

I have some questions to make sure.

1. To increase the nOuterCorr means that we increase the calculation number for each time step right?

and by adding "turbOnFinalIterOnly no" turbulent values are calculated in the same way for each time step?

2. yeah i changed the epsilon value and put calculated one as 0.7 this time.

3. yeah righ it is rectangular mesh and you can see the checkMesh results:

Mesh stats
    points:           279408
    faces:            669203
    internal faces:   585074
    cells:            195257
    boundary patches: 3
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     167069
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     28188

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
    Patch               Faces    Points   Surface topology                  
    wall                83841    84565    ok (non-closed singly connected)  
    outlet              64       81       ok (non-closed singly connected)  
    inlet               224      261      ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-0.006 -0.008 -2.71051e-20) (0.00194 0.007 0.00194)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-3.85612e-15 -5.87416e-17 3.88763e-16) OK.
    Max cell openness = 2.20581e-16 OK.
    Max aspect ratio = 4.28571 OK.
    Minumum face area = 5.07813e-11. Maximum face area = 3.6e-07.  Face area magnitudes OK.
    Min volume = 6.98242e-16. Max volume = 2.16e-10.  Total volume = 1.37934e-07.  Cell volumes OK.
    Mesh non-orthogonality Max: 43.9569 average: 10.8786
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.683 OK.
    Coupled point location match (average 0) OK.

Last question is that how can i estimate y+?

Thanks in advance.

B

[alexeym]

If you take a look at the interPhaseChangeFoam source code, you'll find there:

        while (pimple.loop())
        {
            #include "UEqn.H"

            // --- Pressure corrector loop
            while (pimple.correct())
            {
                #include "pEqn.H"
            }

            if (pimple.turbCorr())
            {
                turbulence->correct();
            }
        }

- pimple.loop() will return false when either number of iterations exceed nOuterCorrectors or residuals in the beginning of the iteration will go below the values you set in residualControl. So if I get your question right, yes, when you increase number of outer correctors, you increase amount of computation during every time step.

- pimple.turbCorr() will return true on every iteration if turbOnFinalIterOnly is no, otherwise turbulence->correct() will be executed only once per time step.

About y+: there are plenty of calculators. One of them in Tools menu - http://www.cfd-online.com/Tools/yplus.php.

[shipman]

Dear Alexey,

Finally, program gave again following floating point error. Mass flow, velocity inside of nozzle seem OK, but at the inlet pressure are increasing so excessively.

Also inside of error i recognized that time step cont. error became crazy?? what is the meaning of this error and how can i correct it?

Thanks in advance.

B

time step continuity errors : sum local = 4.64985e+125, global = -7.72928e+121, cumulative = -1.28369e+122

ERROR:

Courant Number mean: 1153.13 max: 617400
deltaT = 2.1961e-18
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 1010
    Increased the timePrecision from 8 to 9 to distinguish between timeNames at time 0.000458572
Time = 0.000458572308

DILUPBiCG:  Solving for alpha1, Initial residual = 0.186596, Final residual = 0.000478897, No Iterations 1
Liquid phase volume fraction = 0.999853  Min(alpha1) = 0.352503  Max(alpha1) = 1
DILUPBiCG:  Solving for epsilon, Initial residual = 1, Final residual = 1.78371e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 1.90975e-07, No Iterations 2
PIMPLE: iteration 1
GAMG:  Solving for p_rgh, Initial residual = 0.996129, Final residual = 0.00234635, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.3641e-05, Final residual = 5.29413e-10, No Iterations 1
time step continuity errors : sum local = 0.213357, global = 1.01812e-06, cumulative = 8.80345
GAMG:  Solving for p_rgh, Initial residual = 0.000155482, Final residual = 7.27544e-08, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.00180428, Final residual = 3.9623e-06, No Iterations 1
time step continuity errors : sum local = 10447.6, global = 8.63117e-07, cumulative = 8.80345
PIMPLE: iteration 2
GAMG:  Solving for p_rgh, Initial residual = 0.00207729, Final residual = 6.19547e-05, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.00286512, Final residual = 0.000244628, No Iterations 1
time step continuity errors : sum local = 3.31121e+09, global = -5.84454, cumulative = 2.95891
GAMG:  Solving for p_rgh, Initial residual = 0.116709, Final residual = 0.0062631, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.104965, Final residual = 0.00992775, No Iterations 1
time step continuity errors : sum local = 1.13935e+09, global = -2.36284e+06, cumulative = -2.36284e+06
PIMPLE: iteration 3
GAMG:  Solving for p_rgh, Initial residual = 0.159014, Final residual = 0.00972769, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.939243, Final residual = 0.0262975, No Iterations 2
time step continuity errors : sum local = 1.37437e+15, global = -1.1007e+14, cumulative = -1.1007e+14
GAMG:  Solving for p_rgh, Initial residual = 0.139423, Final residual = 0.00878953, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.831211, Final residual = 0.0467336, No Iterations 1
time step continuity errors : sum local = 3.92446e+16, global = -8.8333e+13, cumulative = -1.98403e+14
PIMPLE: iteration 4
GAMG:  Solving for p_rgh, Initial residual = 0.242785, Final residual = 0.0136361, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.973402, Final residual = 0.054481, No Iterations 1
time step continuity errors : sum local = 4.48943e+23, global = -5.0509e+19, cumulative = -5.05092e+19
GAMG:  Solving for p_rgh, Initial residual = 0.157715, Final residual = 0.00695594, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.346942, Final residual = 0.0191784, No Iterations 1
time step continuity errors : sum local = 4.63282e+23, global = -7.63235e+19, cumulative = -1.26833e+20
PIMPLE: iteration 5
GAMG:  Solving for p_rgh, Initial residual = 0.181717, Final residual = 0.0104671, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.97619, Final residual = 0.0542755, No Iterations 1
time step continuity errors : sum local = 8.67454e+30, global = -9.74491e+26, cumulative = -9.74491e+26
GAMG:  Solving for p_rgh, Initial residual = 0.170678, Final residual = 0.00688831, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324195, Final residual = 0.0186419, No Iterations 1
time step continuity errors : sum local = 8.8706e+30, global = -1.47446e+27, cumulative = -2.44895e+27
PIMPLE: iteration 6
GAMG:  Solving for p_rgh, Initial residual = 0.183186, Final residual = 0.0103863, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976201, Final residual = 0.0542631, No Iterations 1
time step continuity errors : sum local = 1.6764e+38, global = -1.88316e+34, cumulative = -1.88316e+34
GAMG:  Solving for p_rgh, Initial residual = 0.170606, Final residual = 0.0068871, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186373, No Iterations 1
time step continuity errors : sum local = 1.71433e+38, global = -2.84964e+34, cumulative = -4.7328e+34
PIMPLE: iteration 7
GAMG:  Solving for p_rgh, Initial residual = 0.183132, Final residual = 0.0103844, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976203, Final residual = 0.0542605, No Iterations 1
time step continuity errors : sum local = 3.23962e+45, global = -3.63913e+41, cumulative = -3.63913e+41
GAMG:  Solving for p_rgh, Initial residual = 0.17059, Final residual = 0.00688685, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186363, No Iterations 1
time step continuity errors : sum local = 3.31293e+45, global = -5.50695e+41, cumulative = -9.14608e+41
PIMPLE: iteration 8
GAMG:  Solving for p_rgh, Initial residual = 0.183121, Final residual = 0.010384, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.05426, No Iterations 1
time step continuity errors : sum local = 6.26046e+52, global = -7.03248e+48, cumulative = -7.03248e+48
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.0068868, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 6.40213e+52, global = -1.0642e+49, cumulative = -1.76745e+49
PIMPLE: iteration 9
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542599, No Iterations 1
time step continuity errors : sum local = 1.20981e+60, global = -1.359e+56, cumulative = -1.359e+56
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 1.23719e+60, global = -2.05654e+56, cumulative = -3.41554e+56
PIMPLE: iteration 10
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542599, No Iterations 1
time step continuity errors : sum local = 2.33792e+67, global = -2.62622e+63, cumulative = -2.62622e+63
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 2.39082e+67, global = -3.97418e+63, cumulative = -6.6004e+63
PIMPLE: iteration 11
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 4.51794e+74, global = -5.07508e+70, cumulative = -5.07508e+70
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 4.62018e+74, global = -7.67996e+70, cumulative = -1.2755e+71
PIMPLE: iteration 12
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 8.73076e+81, global = -9.80741e+77, cumulative = -9.80741e+77
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 8.92833e+81, global = -1.48413e+78, cumulative = -2.46487e+78
PIMPLE: iteration 13
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 1.68719e+89, global = -1.89525e+85, cumulative = -1.89525e+85
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 1.72537e+89, global = -2.86802e+85, cumulative = -4.76326e+85
PIMPLE: iteration 14
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 3.26043e+96, global = -3.6625e+92, cumulative = -3.6625e+92
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 3.33421e+96, global = -5.54234e+92, cumulative = -9.20484e+92
PIMPLE: iteration 15
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 6.30066e+103, global = -7.07764e+99, cumulative = -7.07764e+99
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 6.44324e+103, global = -1.07104e+100, cumulative = -1.7788e+100
PIMPLE: iteration 16
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 1.21758e+111, global = -1.36773e+107, cumulative = -1.36773e+107
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 1.24513e+111, global = -2.06974e+107, cumulative = -3.43747e+107
PIMPLE: iteration 17
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 2.35293e+118, global = -2.64309e+114, cumulative = -2.64309e+114
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 2.40618e+118, global = -3.9997e+114, cumulative = -6.64279e+114
PIMPLE: iteration 18
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.976204, Final residual = 0.0542598, No Iterations 1
time step continuity errors : sum local = 4.54695e+125, global = -5.10767e+121, cumulative = -5.10767e+121
GAMG:  Solving for p_rgh, Initial residual = 0.170587, Final residual = 0.00688679, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.324194, Final residual = 0.0186361, No Iterations 1
time step continuity errors : sum local = 4.64985e+125, global = -7.72928e+121, cumulative = -1.28369e+122
PIMPLE: iteration 19
GAMG:  Solving for p_rgh, Initial residual = 0.183118, Final residual = 0.0103839, No Iterations 1
[11] #0  Foam::error::printStack(Foam::Ostream&)[8] #0[10] #0    Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/li in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAb/libOpenFOAM.so"
[8] #1M.so"
[10] #1  Foam::sigFpe::sigHandler(int)  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[11] #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[11] #2   in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[10] #2   in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #2   in "/lib/libc.so.6"
[11] #3  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/lib/libc.so.6"
[10] #3  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[11] #4  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[10] #4  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[11] #5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[10] #5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/lib/libc.so.6"
[8] #3  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[11] #6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[10] #6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[11] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[10] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #4  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[11] #8   in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[10] #8   in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[8] #8  


[8]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[8] #9  __libc_start_main[11]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[11] #9  __libc_start_main in "/lib/libc.so.6"
[11] #10  
[10]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[10] #9  __libc_start_main in "/lib/libc.so.6"
[8] #10  [11]  in "/home/yusuke/OpenFOAM/yusuke-2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam"
[VENUS:21120] *** Process received signal ***
[VENUS:21120] Signal: Floating point exception (8)
[VENUS:21120] Signal code:  (-6)
[VENUS:21120] Failing at address: 0x3e800005280
[VENUS:21120] [ 0] /lib/libc.so.6(+0x33af0) [0x7fd7d9a32af0]
[VENUS:21120] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7fd7d9a32a75]
[VENUS:21120] [ 2] /lib/libc.so.6(+0x33af0) [0x7fd7d9a32af0]
[VENUS:21120] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver13scalingFactorERNS_5FieldIdEERKS2_S5_S5_+0x74) [0x7fd7da9c8664]
[VENUS:21120] [ 4] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver13scalingFactorERNS_5FieldIdEERKNS_9lduMatrixES3_RKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xa6) [0x7fd7da9c87d6]
[VENUS:21120] [ 5] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_RNS1_IS8_EESD_h+0x126c) [0x7fd7da9caa6c]
[VENUS:21120] [ 6] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x492) [0x7fd7da9cc182]
[VENUS:21120] [ 7] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x157) [0x7fd7db548bf7]
[VENUS:21120] [ 8] RNGinterPhaseChangeFoam() [0x432aa8]
[VENUS:21120] [ 9] /lib/libc.so.6(__libc_start_main+0xfd) [0x7fd7d9a1dc4d]
[VENUS:21120] [10] RNGinterPhaseChangeFoam() [0x42b849]
[VENUS:21120] *** End of error message ***

[alexeym]

Hi,

well, now we arrived to the point where you can describe the problem you are trying to simulate Cause 1e5 for p_rgh seems to be quite large. Also you can try to change relaxation settings: 0.7 for U, 0.3 for p_rgh, k and epsilon can be left as they are.