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,
https://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