|
[Sponsors] |
April 28, 2010, 21:13 |
User sub routine for Stokes number
|
#1 |
New Member
Mathieu Etcheparre
Join Date: Aug 2009
Posts: 14
Rep Power: 16 |
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 |
|
April 29, 2010, 19:42 |
|
#2 |
New Member
Mathieu Etcheparre
Join Date: Aug 2009
Posts: 14
Rep Power: 16 |
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 20:58. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
RPM in Wind Turbine | Pankaj | CFX | 9 | November 23, 2009 04:05 |
Problems with Courant number (LaunderGibsonTurbulence Model) | sven | OpenFOAM | 3 | August 10, 2009 03:12 |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 07:36 |
user defined function | cfduser | CFX | 0 | April 29, 2006 10:58 |
CFX 10 User Sub Routine | Claudia | CFX | 6 | February 15, 2006 08:32 |