TEXT( NACA2414 Airfoil at RE=60200, Alpha=2.8
TITLE     
mesg( this is a test case of drag force on airfoil) 
  DISPLAY
  This is a test case of drag force on airfoil with PARSOL.           
  The geometry, Cd and Cl are all from NASG,
  http://www.nasg.com/.  
  This case used imbalance patch printed in result file
  to get the drag and lift forces like:
  
  x-wise force on IMBL     is     1.922913E-03
  y-wise force on IMBL     is     2.185478E-02

  Other parameters for calculations:
  Velocity 0.929 [m/s]; density 1.189 [kg/m^3];
  Planform area 0.1 [m^2]

  Cd= (x-wise force)/(Planform area*0.5*density*Velocity^2) 
  Cl= (y-wise force)/(Planform area*0.5*density*Velocity^2)
  Comparision to NASG's experiemtal data was made in table:
               NASG-exp         Simulation 
   ------------------------------------------------------- 
       Cl       0.315              0.426
       Cd       0.0316             0.0375
   --------------------------------------------------------
  
   Click to see pressure contour 
  ENDDIS

 TALK=T;RUN( 1, 1) 
 ************************************************************
   Q1 created by VDI menu, Version 3.5, Date 07/03/03
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =         1 ;LIBREF =        0
 ************************************************************
  Group 1. Run Title
 TEXT(NACA wing Cd test ( angle 2.8)          )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,180,140,1)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
 RHO1    = 1.189000E+00
 ENUL    = 1.544000E-05
 CP1     = 1.005000E+03
 ENUT    = 0.000000E+00
 DVO1DT  = 3.410000E-03
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
   No PATCHes used for this Group
 
 
 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
PATCH(IMBL,VOLUME,3,NX-3,3,NY-3,1,1,1,1)
COVAL(IMBL,U1,0.0,0.0)
COVAL(IMBL,V1,0.0,0.0)
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      1000
 RESFAC  = 1.000000E-04
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 RELAX(U1  ,FALSDT, 1.000000E-02)
 RELAX(V1  ,FALSDT, 1.000000E-02)
 SARAH=10.0
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 ASAP    =    T
 PARSOL  =    T
 RG(  1) = 6.830000E+01
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =        13 ;IYMON  =        10 ;IZMON  =         1
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        1.000000E+01, 6.000000E+00, 1.000000E-01
> DOM,    MONIT,       1.499490E+00, 1.156200E+00, 5.000000E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,     60,-1.040000E+00,G
> GRID,   RSET_X_2,     60, 1.000000E+00,G
> GRID,   RSET_X_3,     60, 1.050000E+00,G
> GRID,   RSET_Y_1,     60,-1.055000E+00,G
> GRID,   RSET_Y_2,     20, 1.000000E+00
> GRID,   RSET_Y_3,     60, 1.055000E+00,G
 
> OBJ,    NAME,        IN
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 6.000000E+00, 1.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    VELOCITY,      9.290000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,   0.000000E+00
 
> OBJ,    NAME,        OUT
> OBJ,    POSITION,    1.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 6.000000E+00, 1.000000E-01
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   SAME
> OBJ,    COEFFICIENT,   1.000000E+03
 
> OBJ,    NAME,        B3
> OBJ,    POSITION,    3.500000E+00, 2.800000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.401000E-01, 1.000000E-01
> OBJ,    ROT-ANGLE,   3.572000E+02, 0.000000E+00, 0.000000E+00
> OBJ,    CLIPART,     naca2414
> OBJ,    ROTATION24,        0
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
STOP
distil=t
 EX(U1  )=   8.946E-01 
 EX(PRPS)=   6.668E+00 
 EX(P1  )=   3.121E-02 
 EX(V1  )=   2.416E-02