
[Sponsors] 
March 11, 2014, 23:16 
Help for the small implementation in turbulence model

#1 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
Hi to all,
I am quite new about implementation and i want to change the kinematic eddy viscosity equation (nut) in kepsilon turbulence model which is defined: Code:
nut_ = Cmu_*sqr(k_)/epsilon_; Code:
nut_ = f(rho)*Cmu_*sqr(k_)/epsilon_; f(rho) = rho1+ (rho2rho)/[(rho2rho1)]^5 Your help will be appreciated. Thanks in advance. B 

March 12, 2014, 05:15 

#2 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
Hi,
maybe it'll be easier if you just multiply nut field returned from turbulence model by f(rho) inside the solver? 

March 12, 2014, 12:08 

#3 
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 14 

March 13, 2014, 01:40 
hi

#4 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
Dear Alexey
Thanks so much for your advice. i have 2 questions. 1. When i look at the U.Eqn. in interphasechangeFoam: Code:
surfaceScalarField muEff ( "muEff", twoPhaseProperties>muf() + fvc::interpolate(rho*turbulence>nut()) ); 2. then i inserted my equation here as follows: Code:
surfaceScalarField muEff ( "muEff", twoPhaseProperties>muf() + fvc::interpolate(double Foam::pow((rho2rho)/(rho2rho1))*(rho1rho2),5.0)+ rho2)*rho*turbulence>nut()) ); Code:
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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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/yusuke2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean yusuke@VENUS:~/OpenFOAM/yusuke2.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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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 primaryexpression 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/yusuke2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean yusuke@VENUS:~/OpenFOAM/yusuke2.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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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 primaryexpression 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/yusuke2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean yusuke@VENUS:~/OpenFOAM/yusuke2.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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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 primaryexpression 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/yusuke2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean yusuke@VENUS:~/OpenFOAM/yusuke2.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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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 primaryexpression 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/yusuke2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean yusuke@VENUS:~/OpenFOAM/yusuke2.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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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 primaryexpression 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 thanks B 

March 13, 2014, 03:32 

#5 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
1. mu is dynamic viscosity, nu is kinematic viscosity, mu = rho*nu. nut is kinematic turbulent viscosity.
2. Code:
fvc::interpolate(double Foam::pow((rho2rho)/(rho2rho1))*(rho1rho2),5.0)+ rho2)*rho*turbulence>nut()) 2.2. You've got problems with brackets. According to your initial post the expressions should be Code:
fvc::interpolate((rho1 + (rho2  rho)*Foam::pow(rho2  rho1, 5))*rho*turbulence>nut()) Code:
(rho1 + (rho2  rho)*Foam::pow(rho2  rho1, 5)) Code:
(rho1.value() + (rho2.value()  rho.field())*Foam::pow((rho2  rho1).value(), 5)) i.e. make your f(rho) nondimensional. 

March 13, 2014, 05:34 
hi

#6 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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): Code:
surfaceScalarField muEff ( "muEff", twoPhaseProperties>muf() + fvc::interpolate((Foam::pow((0.0017rho.field()),5)/Foam::pow((0.0017998),5)*(9980.0017)+ 0.0017)*rho*turbulence>nut()) ); Code:
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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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/yusuke2.1.1/applications/solvers/multiphase/RNGinterPhaseChangeFoam$ wclean yusuke@VENUS:~/OpenFOAM/yusuke2.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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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 Thanks in advance. B 

March 13, 2014, 06:11 

#7 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
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) nondimensional, you can introduce unit density densionedScalar and divide rho by this scalar: Code:
dimensionedScalar rhoUnit("rhoUnit", dimDensity, 1.0); ... surfaceScalarField muEff ( "muEff", twoPhaseProperties>muf() + fvc::interpolate((Foam::pow((0.0017rho/rhoUnit),5)/Foam::pow((0.0017998),5)*(9980.0017)+ 0.0017)*rho*turbulence>nut()) ); Code:
 fvm::laplacian(muEff, U)  (fvc::grad(U) & fvc::grad(muEff)) Code:
volScalarField muEff("muEff", rho*nuEff()); ...  fvm::laplacian(muEff, U)  fvc::div(muEff*dev(T(fvc::grad(U)))) 

March 13, 2014, 22:01 
Hi

#8 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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. 

March 14, 2014, 02:09 

#9  
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
Quote:


March 15, 2014, 03:09 
hi

#10 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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: Code:
/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’ [Wunusedvariable] HTML Code:
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); Could you tell me what is the mean reason of this problem? Thanks in advance. B 

March 15, 2014, 09:57 

#11 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
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. 

March 16, 2014, 20:59 
hi

#12 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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: HTML Code:
DILUPBiCG: Solving for alpha1, Initial residual = 0.00424273, Final residual = 1.84854e06, 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.82922e07, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.0577387, Final residual = 3.69088e07, 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.87093e05, No Iterations 1 time step continuity errors : sum local = 0.000646625, global = 2.71757e05, cumulative = 0.0123716 GAMG: Solving for p_rgh, Initial residual = 0.000603131, Final residual = 3.33279e05, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 3.97107e05, Final residual = 3.9419e06, No Iterations 7 time step continuity errors : sum local = 0.000425533, global = 1.63888e05, cumulative = 0.0123552 GAMG: Solving for p_rgh, Initial residual = 0.000117561, Final residual = 9.08278e06, No Iterations 1 GAMGPCG: Solving for p_rgh, Initial residual = 1.36361e05, Final residual = 7.04989e08, No Iterations 9 time step continuity errors : sum local = 0.000420037, global = 1.71133e05, cumulative = 0.0123381 ExecutionTime = 10435.8 s ClockTime = 10636 s Courant Number mean: 0.0024306 max: 0.935238 deltaT = 7.38217e16 > 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.17863e06, 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.47643e09, No Iterations 2 DILUPBiCG: Solving for k, Initial residual = 0.489636, Final residual = 7.37747e07, 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.31948e06, Final residual = 8.27181e08, No Iterations 2 time step continuity errors : sum local = 0.000709651, global = 1.68882e05, cumulative = 0.0123212 GAMG: Solving for p_rgh, Initial residual = 1.42336e05, Final residual = 7.77456e07, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 5.56115e07, Final residual = 4.50147e08, No Iterations 2 time step continuity errors : sum local = 0.000680165, global = 4.07245e05, cumulative = 0.0122805 GAMG: Solving for p_rgh, Initial residual = 7.13622e07, Final residual = 3.20422e08, No Iterations 1 GAMGPCG: Solving for p_rgh, Initial residual = 6.03872e08, Final residual = 6.03872e08, No Iterations 0 time step continuity errors : sum local = 0.000682026, global = 4.08892e05, cumulative = 0.0122396 ExecutionTime = 10437 s ClockTime = 10637 s Courant Number mean: 0.0039758 max: 0.746856 deltaT = 4.94217e16 Time = 0.0015789419679168 DILUPBiCG: Solving for alpha1, Initial residual = 0.00188925, Final residual = 1.1283e06, 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.84687e07, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.170049, Final residual = 5.61643e07, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 0.0799552, Final residual = 7.82053e05, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.30074e12, Final residual = 4.30074e12, 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.64617e09, Final residual = 2.64617e09, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.64617e09, Final residual = 2.64617e09, 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.6463e09, Final residual = 2.6463e09, No Iterations 0 GAMGPCG: Solving for p_rgh, Initial residual = 2.6463e09, Final residual = 2.6463e09, 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.22426e21 > 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/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [3] #4 [3] in "/home/yusuke/OpenFOAM/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [3] #5 [8] in "/home/yusuke/OpenFOAM/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [8] #4 [3] in "/home/yusuke/OpenFOAM/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [3] #6 __libc_start_main [8] in "/home/yusuke/OpenFOAM/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [8] #5 in "/lib/libc.so.6" [3] #7 [8] in "/home/yusuke/OpenFOAM/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [8] #6 __libc_start_main in "/lib/libc.so.6" [8] #7 [3] in "/home/yusuke/OpenFOAM/yusuke2.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). thanks in advance. After checking the simulation i will also send some pictures. 

March 17, 2014, 01:27 
hi

#13 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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. Last edited by shipman; March 17, 2014 at 03:44. 

March 17, 2014, 03:48 

#14 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
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? 

March 17, 2014, 04:15 

#15 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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 

March 17, 2014, 05:12 

#16 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
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 Code:
PIMPLE { momentumPredictor no; nOuterCorrectors 1; nCorrectors 3; nNonOrthogonalCorrectors 1; cAlpha 0; nAlphaCorr 1; nAlphaSubCycles 1; } Code:
PIMPLE { momentumPredictor no; nOuterCorrectors 50; nCorrectors 2; nNonOrthogonalCorrectors 1; cAlpha 0; nAlphaCorr 1; nAlphaSubCycles 1; turbOnFinalIterOnly no; residualControl { "(p_rghUkepsilon)" { tolerance 1e2; relTol 0; } } } 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+? 

March 17, 2014, 05:40 
hi

#17 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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: HTML Code:
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 (nonclosed singly connected) outlet 64 81 ok (nonclosed singly connected) inlet 224 261 ok (nonclosed singly connected) Checking geometry... Overall domain bounding box (0.006 0.008 2.71051e20) (0.00194 0.007 0.00194) Mesh (nonempty, nonwedge) directions (1 1 1) Mesh (nonempty) directions (1 1 1) Boundary openness (3.85612e15 5.87416e17 3.88763e16) OK. Max cell openness = 2.20581e16 OK. Max aspect ratio = 4.28571 OK. Minumum face area = 5.07813e11. Maximum face area = 3.6e07. Face area magnitudes OK. Min volume = 6.98242e16. Max volume = 2.16e10. Total volume = 1.37934e07. Cell volumes OK. Mesh nonorthogonality Max: 43.9569 average: 10.8786 Nonorthogonality check OK. Face pyramids OK. Max skewness = 1.683 OK. Coupled point location match (average 0) OK. Thanks in advance. B 

March 17, 2014, 05:54 

#18 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
If you take a look at the interPhaseChangeFoam source code, you'll find there:
Code:
while (pimple.loop()) { #include "UEqn.H" //  Pressure corrector loop while (pimple.correct()) { #include "pEqn.H" } if (pimple.turbCorr()) { turbulence>correct(); } }  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.cfdonline.com/Tools/yplus.php. 

March 17, 2014, 21:39 
hi

#19 
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 5 
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: HTML Code:
Courant Number mean: 1153.13 max: 617400 deltaT = 2.1961e18 > 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.78371e10, No Iterations 2 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 1.90975e07, 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.3641e05, Final residual = 5.29413e10, No Iterations 1 time step continuity errors : sum local = 0.213357, global = 1.01812e06, cumulative = 8.80345 GAMG: Solving for p_rgh, Initial residual = 0.000155482, Final residual = 7.27544e08, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 0.00180428, Final residual = 3.9623e06, No Iterations 1 time step continuity errors : sum local = 10447.6, global = 8.63117e07, cumulative = 8.80345 PIMPLE: iteration 2 GAMG: Solving for p_rgh, Initial residual = 0.00207729, Final residual = 6.19547e05, 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/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [8] #9 __libc_start_main[11] in "/home/yusuke/OpenFOAM/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [11] #9 __libc_start_main in "/lib/libc.so.6" [11] #10 [10] in "/home/yusuke/OpenFOAM/yusuke2.1.1/platforms/linux64GccDPOpt/bin/RNGinterPhaseChangeFoam" [10] #9 __libc_start_main in "/lib/libc.so.6" [8] #10 [11] in "/home/yusuke/OpenFOAM/yusuke2.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 *** 

March 18, 2014, 03:14 

#20 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,419
Rep Power: 25 
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. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
implementation of SpalartAllmaras Turbulence Model  zhengjg  Main CFD Forum  0  July 24, 2013 03:43 
Implementation of turbulence model  idefix  OpenFOAM Programming & Development  5  March 27, 2013 02:22 
SA Turbulence model implementation  ganesh  Main CFD Forum  0  March 6, 2006 13:23 
kw Turbulence model implementation  suneesh  Main CFD Forum  4  November 23, 2005 18:35 
A reference on implementation of SpalartAlam. Turbulence Model?  Mohammad Kermani  Main CFD Forum  2  December 26, 1999 03:56 