code_saturne User's Forum

Hi Antoine

The file cs_user_source_terms.f90 is the new name for ustsns.f90 - if you open cs_user_source_terms.f90 you will see the subroutine called ustsns.

So the modifications you want to make you should do in cs_user_source_terms.f90

Hope that helps,

James

Hi - I think if you've only simulate 0.1s then you have not done enough, even the 1000 iterations at 1D-3 is only 1 second. The flow will take longer to develop. As you're using the code-coupling it will take a bit of time for both domains to synchronise and I don't imagine you will have accurate re...

I don't think it has converged yet, the maximum velocity is oscillating a lot. Though I agree it is large, have you a shot of the flow-field?

You have only performed 10 iterations so you should perform more. And if you are using 5 revs per second then you need to decrease the time-step. But it's hard to know since this listing doesn't correspond to your questions (omegay=0.5 here, and the domain is massive compared to your radius of 0.047...

So looks like you are defining omegay correctly.

It is probably another problem.. does the listing reveal anything?

What is the radius of your turbine blade?

What is your inflow velocity?

Hi - I think you've asked this one before (in any case I have answered it) so try using the forum search feature.

To calculate moment you need the cross product of the force vector and distance vector. Distance vector (assuming you rotate about origin) is stored in cdgfbo array. Just search for it in usproj it is explained well. From this you can calculate the moment at each individual face and then sum the mom...

Yes, for the first time-step.

It is in the listing you provide, line 811 is the first instance.