# Help for the small implementation in turbulence model

 Register Blogs Members List Search Today's Posts Mark Forums Read

 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 k-epsilon turbulence model which is defined: Code: `nut_ = Cmu_*sqr(k_)/epsilon_;` i want to change it as: Code: ```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

 March 12, 2014, 05:15 #2 Senior Member   Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,423 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? shipman likes this.

March 12, 2014, 12:08
#3
Senior Member

Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 14
Quote:
 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

 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()) );``` 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: Code: ```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: 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 -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 >)’ /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&, const Foam::UList&, const Foam::UList&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note: Foam::tmp > Foam::pow(const Foam::UList&, const Foam::UList&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note: Foam::tmp > Foam::pow(const Foam::UList&, const Foam::tmp >&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note: Foam::tmp > Foam::pow(const Foam::tmp >&, const Foam::UList&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:85: note: Foam::tmp > Foam::pow(const Foam::tmp >&, const Foam::tmp >&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note: void Foam::pow(Foam::Field&, const Foam::scalar&, const Foam::UList&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note: Foam::tmp > Foam::pow(const Foam::scalar&, const Foam::UList&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note: Foam::tmp > Foam::pow(const Foam::scalar&, const Foam::tmp >&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note: void Foam::pow(Foam::Field&, const Foam::UList&, const Foam::scalar&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note: Foam::tmp > Foam::pow(const Foam::UList&, const Foam::scalar&) /opt/openfoam211/src/OpenFOAM/lnInclude/scalarField.H:86: note: Foam::tmp > Foam::pow(const Foam::tmp >&, 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

 March 13, 2014, 03:32 #5 Senior Member   Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,423 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((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 Code: `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: Code: `(rho1 + (rho2 - rho)*Foam::pow(rho2 - rho1, -5))` should be rewritten as something like: Code: `(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.

 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.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: 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 -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 >)’ /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 >)’ /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

 March 13, 2014, 06:11 #7 Senior Member   Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,423 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) non-dimensional, 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.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 Code: ``` - 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 Code: ``` 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.

 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,423
Rep Power: 25
Quote:
 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.

 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’ [-Wunused-variable]``` I checked the readTimeControls.H file but i didnt understand why it is not used?? HTML Code: ```const bool adjustTimeStep = runTime.controlDict().lookupOrDefault("adjustTimeStep", false); scalar maxCo = runTime.controlDict().lookupOrDefault("maxCo", 1.0); scalar maxDeltaT = runTime.controlDict().lookupOrDefault("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

 March 15, 2014, 09:57 #11 Senior Member   Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,423 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.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] [ 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.
Attached Files
 0.tar.gz (1.6 KB, 5 views) system.tar.gz (1.9 KB, 4 views)

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.

B.
Attached Images
 Vel.jpg (24.8 KB, 13 views) vel2.jpg (11.0 KB, 13 views) p.jpg (14.8 KB, 13 views) alpha.jpg (12.6 KB, 11 views)

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,423 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,423 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; }``` to Code: ```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 likes this.

 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 (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

 March 17, 2014, 05:54 #18 Senior Member   Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,423 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.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.

 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.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&, Foam::Field const&, Foam::Field const&, Foam::Field const&) const in "/lib/libc.so.6" [10] #3 Foam::GAMGSolver::scalingFactor(Foam::Field&, Foam::Field const&, Foam::Field const&, Foam::Field const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [11] #4 Foam::GAMGSolver::scalingFactor(Foam::Field&, Foam::lduMatrix const&, Foam::Field&, Foam::FieldField const&, Foam::UPtrList const&, Foam::Field const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [10] #4 Foam::GAMGSolver::scalingFactor(Foam::Field&, Foam::lduMatrix const&, Foam::Field&, Foam::FieldField const&, Foam::UPtrList const&, Foam::Field const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [11] #5 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList >&, Foam::PtrList >&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [10] #5 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList >&, Foam::PtrList >&, unsigned char) const in "/lib/libc.so.6" [8] #3 Foam::GAMGSolver::scalingFactor(Foam::Field&, Foam::Field const&, Foam::Field const&, Foam::Field const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [11] #6 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [10] #6 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [11] #7 Foam::fvMatrix::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [10] #7 Foam::fvMatrix::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [8] #4 Foam::GAMGSolver::scalingFactor(Foam::Field&, Foam::lduMatrix const&, Foam::Field&, Foam::FieldField const&, Foam::UPtrList const&, Foam::Field const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [8] #5 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList >&, Foam::PtrList >&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [8] #6 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [8] #7 Foam::fvMatrix::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 ***```

 March 18, 2014, 03:14 #20 Senior Member   Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,423 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 Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post zhengjg Main CFD Forum 0 July 24, 2013 03:43 idefix OpenFOAM Programming & Development 5 March 27, 2013 02:22 ganesh Main CFD Forum 0 March 6, 2006 13:23 suneesh Main CFD Forum 4 November 23, 2005 18:35 Mohammad Kermani Main CFD Forum 2 December 26, 1999 03:56

All times are GMT -4. The time now is 06:20.