cd_pk_ms.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cd_pk_ms (rad, d, xl, reynolds, p2, p1, beta, kappa, cd, u, T1, R)

Function/Subroutine Documentation

cd_pk_ms()

subroutine cd_pk_ms	(	real*8	rad,
		real*8	d,
		real*8	xl,
		real*8	reynolds,
		real*8	p2,
		real*8	p1,
		real*8	beta,
		real*8	kappa,
		real*8	cd,
		real*8	u,
		real*8	T1,
		real*8	R
	)

!
!     This subroutines enable to calculate the compressible discharge 
!     coefficient for thin and long orifices with corner radiusing;
!
!     author: Yannick Muller
!     
      implicit none
!
      real*8 rad,d,xl,lqd,rqd,reynolds,p2,p1,p2p1,beta,beta_cor,kappa,
     &     cd,cdc_cl1,cdc_cl3,rldb,r,u,t1,c1,c2,
     &     c3,ms_cdr,rv,vid
!
      p2p1=p2/p1
      rqd=rad/d
      lqd=xl/d
      rldb=max(lqd,0.d0)
!
!     the method of cd calculation for a sharp edged aperture is only valid 
!     for beta comprised between 0 and 0.7 
!      
      if (beta.gt.0.7d0) then
         beta_cor=0.7d0
      else
         beta_cor=beta
      endif
!
!     differences between class1 or class2 or class3
!
      if (lqd.eq.rqd) then
!     
!     class1
!
         call pk_cdc_cl1(lqd,reynolds,p2p1,beta_cor,kappa,cdc_cl1)
         cd=cdc_cl1
      else
!     
!     class2 or class3 (clas2 is a sub class of class3 )
!     
         call pk_cdc_cl3(lqd,rqd,reynolds,p2p1,beta_cor,kappa,cdc_cl3)
         cd=cdc_cl3
      endif
!      
!     if rotating orifice with Mac Greehan & Scotch
!     The decription of the method can be found in :
!     "Flow characteristics of long orifices with rotation and 
!     corner radiusing" ASME 87-GT-16
!
!     rotating case eq 17
      
      if (u.ne.0) then
         vid=dsqrt(2.d0*kappa/(kappa-1.d0)*r*t1*
     &        (1.d0-(p2/p1)**((kappa-1.d0)/kappa)))
         rv=u/vid*(cd/0.6)**(-3)
         c1=exp(-rv**1.2d0)
         c2=0.5d0*rv**0.6d0*dsqrt(0.6/cd)
         c3=exp(-0.5*rv**0.9d0)
         ms_cdr=cd*(c1+c2*c3)
         cd=ms_cdr
         cd=min(max(cd,0.d0),1.d0)

!         vid=dsqrt(2.d0*kappa/(kappa-1.d0)*R*T1*
!     &        (1.d0-p2/p1**((kappa-1.d0)/kappa)))
!         rv=1000*u/vid*(cd/0.6)**(-3)
!         c1=exp(-rv**1.2d0)
!         c2=0.5*rv**(0.6d0)*(cd/0.6)**(-0.5d0)
!         c3=exp(-0.5d0*rv**0.9d0)
!         ms_cdr=cd*(c1+c2*c3)
!         cd=ms_cdr
      endif
!
      return