[Jinzhihui]

Hi
I'm trying to add a body force to momentum source through CFX.
I have a constant body force f(fx fy 0) which Located in a specific area（such as x1<x<x2,y1<y<y2）. I want to add them to the momentum equation.I have added only one direction body force(X or Y,but not all)to Subdomain Source form with User CEL function and Routine in CFX successfully.

PHP Code:

#include "cfx5ext.h"
dllexport(user_source)
      SUBROUTINE USER_SOURCE (
     & NLOC,NRET,NARG,RET,ARGS,CRESLT,CZ,DZ,IZ,LZ,RZ)

      INTEGER NLOC, NRET, NARG
      CHARACTER*(4) CRESLT
      REAL    RET(1:NLOC,1:NRET), ARGS(1:NLOC,1:NARG)

      INTEGER IZ(*)
      CHARACTER CZ(*)*(1)
      DOUBLE PRECISION DZ(*)
      LOGICAL LZ(*)
      REAL RZ(*)

      CALL USER_SOURCE_SUB (NLOC,RET(1,1),ARGS(1,1),ARGS(1,2))
C
      CRESLT = 'GOOD'
      END

      SUBROUTINE USER_SOURCE_SUB (NLOC,SOURCE,X,Y)

      DO ILOC=1,NLOC
         SOURCE(ILOC) = 0.0
         IF (X(ILOC).GE.1.5  .AND. X(ILOC).LE.2.5  .AND.
     &    Y(ILOC).GE.1.5 .AND. Y(ILOC).LE.2.5) THEN
            SOURCE(ILOC) =  1000000.0
         END IF
C---------------------
      END DO
C
      END 


Now I need to add Fx and Fy simultaneously ,The NRET of CEL Routine is always 1 in CFX according to User CEL Function and Routines in the CFX Modeling Guide.So I tried to use two CEL Function from two different Routine,and define of the Functions as below:

(1)this work done successfully,these Functions from the same Routine.




(2)this work can't run




Unfortunately,I got error message.Notice that the Function 1 and 2 from two Routines,because i need two kinds of return Values.





1.If there are any mistakes in what I have finished(add one dirction body force)?

2.What is the reason for this mistake?

3.How to add different direction momentum source at same time?
Appreciate your response....

[Opaque]

You have shown the UserFortran code, but not the context on how it is used.

Please show how it is accessed in the setup, i.e.

MOMENTUM SOURCE:
Option = Cartesian Components
X Component =
Y Component =
Z Component =
END

[Jinzhihui]

Originally Posted by Opaque

You have shown the UserFortran code, but not the context on how it is used.

Please show how it is accessed in the setup, i.e.

MOMENTUM SOURCE:
Option = Cartesian Components
X Component =
Y Component =
Z Component =
END

Thinks for your reply.
this UserFortran code is used to define one User Routine whose Option is User CEL Function, in CFX Pre.The Function which has three Arguments and onlu one return variavle is defined as below,for this routine structure of User CEL Functions only have one renturn variable :




So do you have some suggestions....

[Gert-Jan]

Are you convinced that function 2 is working properly? Did you for example try running only function 2? Or using functions 2 twice, for x and y component?

[Jinzhihui]

Originally Posted by Gert-Jan

Are you convinced that function 2 is working properly? Did you for example try running only function 2? Or using functions 2 twice, for x and y component?

Yes,I have already run only function 2 for X ,for X and Y component both,It's OK.

I guess the reason for error is the two User Function for X and Y component come from two different User Routines, and CFX need to call two different User CEL Function Routines(Both of these two Routines can run independently,and the Structure of Routines same as the first Fig) .This means CFX call the Function using the same Routine(only one) for x and y component, It can run.

[Gert-Jan]

I would like know what ANSYS says........

[Jinzhihui]

Originally Posted by Gert-Jan

I would like know what ANSYS says........

ANSYS prompt this error message


I don't get what you mean....

[Gert-Jan]

I mean the ANSYS support help desk