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 inncompressible fro thin orifices with corner radiusing (eq 5) 20! 21! author: Yannick Muller 22! 23 subroutine pk_cdi_r (rqd,reynolds,beta,cdi_r) 24! 25 implicit none 26! 27 real*8 rqd,reynolds,beta,cdi_r,frqd,cdi_se,cdi_noz 28! 29 call pk_cdi_noz(reynolds,cdi_noz) 30 call pk_cdi_se(reynolds,beta,cdi_se) 31 32 frqd=0.008d0+0.992d0*exp(-5.5d0*rqd-3.5d0*rqd**2.d0) 33! 34 cdi_r=cdi_noz-frqd*(cdi_noz-cdi_se) 35! 36 return 37 end 38