source.for

关于二维射流计算源程序

   
   
        SUBROUTINE Wsource
	  USE VOF2D        
        IMPLICIT  REAL*8(A-H,O-Z),INTEGER(I-N)	  

	   fa=1.0
	   TJ=T-TMAKE
	   if(TJ.lt.twave)fa=0.5*(1.0-dcos(pi*TJ/twave))

	   IF(IWAVE.EQ.1)THEN
	    CWAVE=2*PI/(WNK*TWAVE)
	   	WSC=(THETA0*CWAVE/SAREA)*DSIN(WOMG*TJ)
	   ELSE
	   AS  = 0.5*THETA0
	
	   SHKH   = DSINH(WNK*FLHT0)
	   CHKH   = DCOSH(WNK*FLHT0)
	   CH2KH  = DCOSH(2*WNK*FLHT0)

	   BS=THETA0**2*WNK*CHKH*(2+CH2KH)/(16*SHKH**3)

	   PS=DACOS((-AS+DSQRT(AS**2+8*BS**2))/(4*BS))

         SETA=0.5*PI-WOMG*TJ-PS
	   CS  =DCOS(SETA)
	   CS2 =DCOS(2*SETA)
	   CW  =2*PI/(WNK*TWAVE)
	   WSC  =CW*THETA0*CS/SAREA+
     *       CW*THETA0**2*WNK*CHKH*(2+CH2KH)*CS2/(8*SAREA*SHKH**3)
	  ENDIF
	   WSC=WSC*FA
	   SCALL=SCALL+WSC

	END