CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Rotating channel (http://www.cfd-online.com/Forums/openfoam/83692-rotating-channel.html)

Hagen January 7, 2011 19:24

Rotating channel
 
Hello Everyone,

I want to start a calculation on a rotating channel. In that case the governing equations change due to the Coriolis force.
Has anybody adjusted a solver or turbulence model to capture these effects?

Thanks for your help,
hagen

Curtin University April 6, 2012 00:17

Quote:

Originally Posted by Hagen (Post 289655)
Hello Everyone,

I want to start a calculation on a rotating channel. In that case the governing equations change due to the Coriolis force.
Has anybody adjusted a solver or turbulence model to capture these effects?

Thanks for your help,
hagen

Hi Hagen,

Have you managed to solve this probem? I'm stuck on the same, I've tried using moving reference framce as well as simple UDF to incorporate coriolis and centrifugal forces, however I'm still unable to validate my results with that of establised ASME and other papers....

Any suggestions?

Thanks,

Ivanet April 6, 2012 03:03

Hi Hagen,
MRFSimpleFoam should work for your case.
Greetings
Ivan

Hagen April 10, 2012 03:38

Thanks for your suggestion. However, I think that the MRFSimpleFoam solver alone is not sufficient to account for RANS system rotation thoroughly, since not only the solver, i.e. the momentum equation, but also the turbulence model as well as the scalar flux model (if the energy equation needs to be solved) has to be extended. An additional production term due to Coriolis acceleration appears in both the transport equation for Reynolds stress and the scalar transport equation. More details can be found in the paper by Johnston, Halleen and Lezius (JFM 56, pp. 533-557, 1972). An exemplary extension for an explicit scalar flux model can be found in Younis et al. (Journal Heat Transfer, 132, 2010). Last year I implemented both into OpenFOAM on the basis of SSG and LRR. If anyone is interested in the code, you may send me a private message with your email address so I can provide code + documentation.

Best, Hagen

Curtin University April 10, 2012 04:09

Hi Hagen,

Could I please draw your attention to http://www.cfd-online.com/Forums/fluent/99552-3d-rotating-channel.html#post353897.

I've trying to model a 3D rotating channel using FLUENT. I've been unable to match my Nu plots with that of an expermental analysis conducted by Wagner et al 1991.

I'm currently using a MRF with a constant angular velocity. Are you suggesting that additional source terms for the body forces induced by coriolis acceleration need to be included even if we are using a MRF??

Regards,

Hagen April 10, 2012 05:47

Quote:

Originally Posted by Curtin University (Post 353907)

You mention that you are currently using the standard k-e model for your calculations. Most turbulence models which are based on the eddy viscosity approach are material frame indifferent and therefore "blind" to rotation. In my opinion you might need a rotation sensitized turbulence model for a complete representation of the Coriolis acceleration. I have very little experience with FLUENT and I don't know if such a model is available. However, you'll find the theoretical framework in the papers mentioned overhead.

huangxianbei March 20, 2014 21:52

Hi,Hagen:

I'm just struggling in this problem using DNS. I modified the icoFoam by adding the Coriolics force, and the centrifugal force is considered to be coupling with the pressure gradient term, so no specially term is added with centrifugal force. I ran the rotating case with a satisfied non-rotating case as the initial condition yesterday, until now the velocity is still increasing(my computer is too slow:() , although the tendency of the Reynolds stress is right , the value is not satisfactory and also the U+ distribution along with y+. So can you please share your code with me? I want to see is there anything ignored.My email is:
znhuang@163.com

Huang Xianbei


All times are GMT -4. The time now is 07:05.