Search found 72 matches

by JamesMcNaughton
Wed Sep 05, 2012 5:08 pm
Forum: code_saturne usage
Topic: effortts calculation (F=PxS)
Replies: 37
Views: 15443

Re: effortts calculation (F=PxS)

Without knowing more about the set-up it's just a guess to the problem.. it's possible you are running in parallel and Face_201 is not on process 0 (the process that writes to the listing) and so you will never enter the loop. In your example xfor is the sum of forces, the parrsm sums across all pro...
by JamesMcNaughton
Wed Sep 05, 2012 3:48 pm
Forum: code_saturne usage
Topic: effortts calculation (F=PxS)
Replies: 37
Views: 15443

Re: effortts calculation (F=PxS)

Yeah like Yvan says you have what you need to calculate the moment so you just need to print your result to a file, if it's just a number then you can put it in the listing, for example: write(nfecra,*)xfor(1),xfor(2),xfor(3) Or could put this in a specific file, look at the "profile" exam...
by JamesMcNaughton
Tue Sep 04, 2012 3:06 pm
Forum: code_saturne usage
Topic: effortts calculation (F=PxS)
Replies: 37
Views: 15443

Re: effortts calculation (F=PxS)

Change this line:

Code: Select all

xfor(ii) = xfor(ii) + forbr(ii, ifac)
To:

Code: Select all

xfor(ii) = xfor(ii) + ra(iforbr + (ifac-1)*ndim + ii-1)
That should work as per the v2.0.1 example of usproj..

Cheers
by JamesMcNaughton
Tue Sep 04, 2012 11:12 am
Forum: code_saturne usage
Topic: effortts calculation (F=PxS)
Replies: 37
Views: 15443

Re: effortts calculation (F=PxS)

You're using saturne 2.0 so you should use the usproj of that version.. check your previous post here http://code-saturne.org/forum/viewtopic ... 6077#p6077 (use the line you commented out).
by JamesMcNaughton
Tue Sep 04, 2012 9:19 am
Forum: code_saturne usage
Topic: effortts calculation (F=PxS)
Replies: 37
Views: 15443

Re: effortts calculation (F=PxS)

The forces are displayed in the post-processing by default (efforts).

If you want to print them to a file or something then you can use usproj as you suggest though delete the lines:

Code: Select all

iutile = 0

if (iutile.eq.0) return
Or change the value of iutile...
by JamesMcNaughton
Wed Aug 29, 2012 9:07 am
Forum: code_saturne usage
Topic: yplus output
Replies: 10
Views: 6748

Re: yplus output

Hi, to interupt a calculation make a file called ficstp in the tmp_Saturne/CASE.12345678 folder. Leave first line of ficstp blank then on the second line enter the new time-step number to stop at, if you enter a value less than the value already reached it will stop at the next time-step. This ends ...
by JamesMcNaughton
Wed Jul 25, 2012 5:06 pm
Forum: code_saturne usage
Topic: Refined mesh and scalable wall function
Replies: 5
Views: 4047

Re: Refined mesh and scalable wall function

Just so you know you can change the wall-function in usini1. Type the command

Code: Select all

 ideuch = 
in usini1 in the usipsu part of the subroutine.

Thanks,

James
by JamesMcNaughton
Thu Jul 19, 2012 8:07 am
Forum: code_saturne usage
Topic: turbine modelling
Replies: 123
Views: 51281

Re: turbine modelling

Hi Stephane - if you search for efforts on this forum there's quite a few discussions and examples for calculating forces. There's an example in usproj too (again search for efforts).

Cheers
by JamesMcNaughton
Tue Jul 17, 2012 10:50 am
Forum: code_saturne usage
Topic: turbine modelling
Replies: 123
Views: 51281

Re: turbine modelling

Glad it is working now..

To print the post-processing (ensight / paraview) files change ntchr=X with X the frequency of time-steps to print at.
by JamesMcNaughton
Tue Jul 17, 2012 8:54 am
Forum: code_saturne usage
Topic: turbine modelling
Replies: 123
Views: 51281

Re: turbine modelling

It looks like it is looking for the folder: DATA.1/RESTART/suiava when it should be looking for the folder DATA.1/RESTART (this contains the files suiava and suiavx). Have hard coded something in runcase that might be an issue? Does it say suiava anywhere in the file? Otherwise you'll have to ask so...