CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

User sub routine for Stokes number

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 28, 2010, 22:13
Default User sub routine for Stokes number
  #1
New Member
 
Mathieu Etcheparre
Join Date: Aug 2009
Posts: 14
Rep Power: 7
cheche is on a distinguished road
Hi,

I'm am modeling an airblast atomizer and would like to look at the Stokes number of my particles amongst other things in the subsequent solved flow. According to the CFX help file I have to wrtie a user sub routine in Fortran to do so.

Can anyone point me as to where I could find similar examples or a nice user guide on how to implement this.

Thanks in advacne for any time spent on this.

Cheers

Mat
cheche is offline   Reply With Quote

Old   April 29, 2010, 20:42
Default
  #2
New Member
 
Mathieu Etcheparre
Join Date: Aug 2009
Posts: 14
Rep Power: 7
cheche is on a distinguished road
Ok after doing some research over the web on this issue I found that I have to use the USER_GETVAR function in Fortran.

Can anyone point me to somewhere where an example of a subroutine with this could be found or help me with the issue?

Cheers

Edit-- this is what I have so far:

#include "cfx5ext.h"

dllexport(stokes)

SUBROUTINE STOKES (NLOC, NRET, NARG, RET, ARGS, CRESLT,

& CZ,DZ,IZ,LZ,RZ)

IMPLICIT NONE

C Calculate magnitude of Stokes number

#include "stack_point.h"

C Result flag

CHARACTER CRESLT*(*)

INTEGER NLOC

INTEGER NARG

INTEGER NRET

REAL ARGS(NLOC,NARG)

REAL RET(NLOC,NRET)

INTEGER IZ(*)

CHARACTER CZ(*)*(1)

DOUBLE PRECISION DZ(*)

LOGICAL LZ(*)

REAL RZ(*)

__stack_point__ pST

C Obtain Stokes number

CALL USER_GETVAR('ptstt',CRESLT,
& pST,CZ,DZ,IZ,LZ,RZ)

IF (CRESLT.NE. 'GOOD') STOP


END

Last edited by cheche; April 29, 2010 at 21:58.
cheche is offline   Reply With Quote

Reply

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
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Problems with Courant number (LaunderGibsonTurbulence Model) sven OpenFOAM 3 August 10, 2009 04:12
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 08:36
user defined function cfduser CFX 0 April 29, 2006 11:58
CFX 10 User Sub Routine Claudia CFX 6 February 15, 2006 09:32


All times are GMT -4. The time now is 12:17.