CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Editing RSTM Turbulence Model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 5, 2009, 11:39
Default Editing RSTM Turbulence Model
  #1
Member
 
Sven Winkler
Join Date: May 2009
Posts: 70
Rep Power: 8
sven is on a distinguished road
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 is offline   Reply With Quote

Old   August 5, 2009, 12:09
Default
  #2
Member
 
Sven Winkler
Join Date: May 2009
Posts: 70
Rep Power: 8
sven is on a distinguished road
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
sven is offline   Reply With Quote

Old   September 18, 2009, 05:29
Default
  #3
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
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!
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SimpleFoam case with SpalartAllmaras turbulence model implemented nedved OpenFOAM Running, Solving & CFD 2 November 30, 2014 23:43
Eul-Eul flow, k-e-kp-ep-Theta Turbulence model us FLUENT 5 April 5, 2011 02:29
KOmega Turbulence model from wwwopenFOAMWikinet philippose OpenFOAM Running, Solving & CFD 30 August 4, 2010 10:26
validating turbulence model on flow around a car Pedro CFX 1 February 20, 2008 17:32
SSG Reynolds Turbulence Model Georges CFX 1 February 28, 2007 17:15


All times are GMT -4. The time now is 09:49.