CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Implementation of system rotation

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

LinkBack Thread Tools Display Modes
Old   January 19, 2011, 19:24
Default Implementation of system rotation
Hagen Müller
Join Date: Nov 2010
Posts: 34
Rep Power: 9
Hagen is on a distinguished road
Dear Foamers,

I need to extend an existing RAS turbulence model, to capture effects of system rotation. Therefore a new production term needs to be defined, which goes into the differential equation for ddt(R_).
This term is: G_ = 2.0 * (Omega_ & R_ - R_ & Omega_)
where R_ is the Re stress tensor and Omega_ is a new tensor that depends on the rate of system rotation. I want to make the turb.model read Omega_ from a file in the caseSetup. runTime or transportProperties.
So far I added these lines to the existing code:

In the beginning, where k,R,epsilon are initialized, I added:



volTensorField OR_ = (Omega_ & R_);

volTensorField RO_ = (R_ & Omega_);

volSymmTensorField G_rot = 2.0*symm(OR_- RO_);

Unfortunately, that doesn't work. The error text is to long for my terminal window. I can't read the beginning.
Can anybody help me with this problem? Basically it's just reading a tensor and calculate G_rot. But that can be hard for someone like me with little C++ knowledge.

Thank you!
Hagen is offline   Reply With Quote


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
CFX11 + Fortran compiler ? Mohan CFX 20 March 30, 2011 18:56
Need ideas-fuel discharge system Jan FLUENT 1 October 10, 2006 23:05
need ideas-fuel discharge system Jan Siemens 0 October 9, 2006 04:31
Need ideas-fuel delivery system Jan NUMECA 0 October 9, 2006 04:30
Need ideas-fuel discharge system Jan Main CFD Forum 0 October 9, 2006 04:27

All times are GMT -4. The time now is 13:25.