- **OpenFOAM Programming & Development**
(*https://www.cfd-online.com/Forums/openfoam-programming-development/*)

- - **Implementation of system rotation**
(*https://www.cfd-online.com/Forums/openfoam-programming-development/84085-implementation-system-rotation.html*)

Implementation of system rotationDear 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: Omega_ ( IOobject ( "Omega", runTime_.timeName(), mesh_, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh_ ) later: 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! |

All times are GMT -4. The time now is 03:28. |