1! 2! CalculiX - A 3-dimensional finite element program 3! Copyright (C) 1998-2021 Guido Dhondt 4! 5! This program is free software; you can redistribute it and/or 6! modify it under the terms of the GNU General Public License as 7! published by the Free Software Foundation(version 2); 8! 9! 10! This program is distributed in the hope that it will be useful, 11! but WITHOUT ANY WARRANTY; without even the implied warranty of 12! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 13! GNU General Public License for more details. 14! 15! You should have received a copy of the GNU General Public License 16! along with this program; if not, write to the Free Software 17! Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 18! 19! cd compressible for class 3 orifice where, l/d>0 and r/d>0 20!typ a) with 0 <= l/d<=0.28 (eq. 25 modified) 21! 22 subroutine pk_cdc_cl3a(lqd,rqd,reynolds,p2p1,beta,kappa,cdc_cl3a) 23! 24 implicit none 25! 26 real*8 lqd,rqd,reynolds,p2p1,beta,kappa,cdc_cl3a,cdi_noz, 27 & cdi_rl,cdi_se,y0,yg,cdqcv_noz,cdqcv_rl 28! 29! cd incompressible nozlle eq 4a 4b 30 call pk_cdi_noz(reynolds,cdi_noz) 31! cd incompresible eq.6 32 call pk_cdi_rl(lqd,rqd,reynolds,beta,cdi_rl) 33! cd incompressible sharp edge eq.3 34 call pk_cdi_se(reynolds,beta,cdi_se) 35! y0,yg ,eq. 15-17, eq.18 36 call pk_y0_yg(p2p1,beta,kappa,y0,yg) 37! 38 cdqcv_noz=cdi_noz/(0.0718d0*cdi_noz+0.9282d0) 39 cdqcv_rl=cdi_rl/(0.0718d0*cdi_rl+0.9282d0) 40! 41! eq.26 modified for class 3a 42! 43 cdc_cl3a=cdi_rl*((cdqcv_noz-cdqcv_rl)/(cdqcv_noz-cdi_se/0.971d0) 44 & *(y0/yg-1.d0)+1.d0) 45! 46 return 47! 48 end 49