CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Editing RSTM Turbulence Model (http://www.cfd-online.com/Forums/openfoam/67147-editing-rstm-turbulence-model.html)

sven August 5, 2009 11:39

Editing RSTM Turbulence Model
 
Hi everyone,

I want to edit the LaunderGibsonRSTM turbulence model which comes with OpenFOAM. More exactly I want to change the equation for the Reynolds Stresses in this model, so that it becomes a new modelling approach. I copied the Folder with the model and renamed it with the name of the new model and I also renamed the files within this copied folder. Then I began to edit the .C and .H file and it seemed to work. But unfortunately I get an error while compiling the edited file. The eror tells me that LHS and RHS (Left- and right Handside) of the equation I edited no longer have the same Dimensions. I think this is because I deleted some Coefficients in this equation and added some new ones. Sadly I cant find where the dimensions of these constants are implemented. Did anyone encounter this problem too or does anyone know how to solve it?
Thanls a lot

sven August 5, 2009 12:09

I thought the posted error message might probably help:

Code:


LHS and RHS of - have different dimensions
    dimensions : [0 0 0 0 0 0 0] - [0 0 -1 0 0 0 0]
#0  Foam::error::printStack(Foam::Ostream&) in "/home/sven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1  Foam::error::abort() in "/home/sven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2  Foam::operator-(Foam::dimensionSet const&, Foam::dimensionSet const&) in "/home/sven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#3  Foam::tmp<Foam::GeometricField<Foam::typeOfSum<Foam::SymmTensor<double>, Foam::SphericalTensor<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::operator-<Foam::SymmTensor<double>, Foam::SphericalTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<Foam::SphericalTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/sven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#4  Foam::incompressible::RASModels::LaunderGibsonModified::correct() in "/home/sven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#5  main in "/home/sven/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/turbFoam"
#6  __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7  _start in "/home/sven/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/turbFoam"


    From function operator-(const dimensionSet& ds1, const dimensionSet& ds2)
    in file dimensionSet/dimensionSet.C at line 423.

FOAM aborting


makaveli_lcf September 18, 2009 05:29

Sven,
from your error out put you have some mismatch of dimensions in substraction of two values in Foam::incompressible::RASModels::LaunderGibsonModi fied::correct() method.
Analyze it if all parts have the same dimensions.
As a hint: it is better first to analyze equations on paper if they are correctly stated before transitioning them to source code.
Good luck!


All times are GMT -4. The time now is 12:33.