(function(){ var PipeID; PipeID=fieldname125; var CSArea; CSArea=Math.PI/4*POW(PipeID/12,2); var LiqFlow; if(fieldname16=="lb/hr") LiqFlow=fieldname14; if(fieldname16=="Kg/hr") LiqFlow=fieldname14*2.20462; var GasFlow; if(fieldname59=="lb/hr") GasFlow=fieldname58; if(fieldname59=="Kg/hr") GasFlow=fieldname58*2.20462; var LiquidDensity; if(fieldname27=="Kg/m3") LiquidDensity=fieldname24/16.01846; if(fieldname27=="lb/ft3") LiquidDensity=fieldname24; var GasDensity; if(fieldname76=="Kg/m3") GasDensity=fieldname75/16.01846; if(fieldname76=="lb/ft3") GasDensity=fieldname75; /* Convert Temp */ if(fieldname17=="°C") Temp=fieldname2; if(fieldname17=="°F") Temp=(fieldname2-32)/1.8; if(fieldname17=="K") Temp=(fieldname2-273.15); /* Convert Gas Viscoisty */ var GasViscosity; if(fieldname19=="cP") GasViscosity=fieldname4; if(fieldname19=="mPa-s") GasViscosity=fieldname4; if(fieldname19=="Poise") GasViscosity=fieldname4*100; /* Convert Liquid Viscoisty */ var LiqViscosity; if(fieldname127=="cP") LiqViscosity=fieldname126; if(fieldname127=="mPa-s") LiqViscosity=fieldname126; if(fieldname127=="Poise") LiqViscosity=fieldname126*100; /* Convert Pressure */ if(fieldname18=="barg") GasPressure=fieldname15*1.01972; if(fieldname18=="psig") GasPressure=fieldname15/14.2233; if(fieldname18=="Kg/cm2g") GasPressure=fieldname15; /* Convert the barometric pressure */ var Barometric; if(fieldname132=="bar a") Barometric=fieldname131*1.01972; if(fieldname132=="psi a") Barometric=fieldname131/14.2233; if(fieldname132=="Kg/cm2 a") Barometric=fieldname131; /* Roughness */ var Roughness; if(fieldname130=="mm") Roughness=fieldname129/25.4; if(fieldname130=="inch") Roughness=fieldname129; /* Vapor Density calculation */ Compressibility=fieldname47; var GasDensityCalculated; GasDensityCalculated=(348.54*(Barometric+GasPressure)*fieldname13/28.97/(Temp+273.15)/Compressibility)/16.01846; /* Logic for using the Vapor Density */ var GasDensityUse; if(fieldname77=="Vapor Density Available") GasDensityUse=GasDensity; if(fieldname77=="Calculate Vapor Density") GasDensityUse=GasDensityCalculated; var LiquidVelocity; LiquidVelocity=LiqFlow/(LiquidDensity*3600)/CSArea; /* Liquid Pressure Drop alone */ var NreLiq; NreLiq=19.83*LiqFlow/(Math.PI*PipeID*LiqViscosity); var A,B,FFTurbLiq,FFUseLiq; A=POW(2.457*LOG(1/(POW(7/NreLiq,0.9)+0.27*Roughness/PipeID)),16); B=POW(37530/NreLiq,16); FFTurbLiq=8*POW(POW(8/NreLiq,12)+1/POW(A+B,1.5),1/12); if(NreLiq<=2100) FFUseLiq=64/NreLiq; if(NreLiq>2100) FFUseLiq=FFTurbLiq; var DelPLiq; DelPLiq=4.167*FFUseLiq*POW(LiquidVelocity,2)*LiquidDensity/(32.1741*PipeID); /* Gas Pressure Drop alone */ var GasVelocity; GasVelocity=GasFlow/(GasDensityUse*3600)/CSArea; var NreGas; var AA,BB; NreGas=19.83*GasFlow/(Math.PI*PipeID*GasViscosity); AA=POW(2.457*LOG(1/(POW(7/NreGas,0.9)+0.27*Roughness/PipeID)),16); BB=POW(37530/NreGas,16); FFTurbGas=8*POW(POW(8/NreGas,12)+1/POW(AA+BB,1.5),1/12); var FFUseGas; if(NreGas<=2100) FFUseGas=64/NreGas; if(NreGas>2100) FFUseGas=FFTurbGas; var DelPGas; DelPGas=4.167*FFUseGas*POW(GasVelocity,2)*GasDensityUse/(32.1741*PipeID); /* Gas viscosity ratio */ var GasVisRatio; GasVisRatio=LiqViscosity/GasViscosity; var TotalMassFlow; TotalMassFlow=GasFlow+LiqFlow; var AvgDensity; AvgDensity=TotalMassFlow/((GasFlow/GasDensityUse)+(LiqFlow/LiquidDensity)); var AvgViscosity; AvgViscosity=TotalMassFlow/((GasFlow/GasViscosity)+(LiqFlow/LiqViscosity)); var TotalMassFlux; TotalMassFlux=TotalMassFlow/(3600*CSArea); var NreTotal; NreTotal=19.83*TotalMassFlow/(Math.PI*PipeID*AvgViscosity); /* 2 phase method based on Viscosity Ratio */ var TwoPhMethod; if(AND(GasVisRatio>=1000,TotalMassFlux>20.5)) TwoPhMethod="Chisholm-Baroczy (C-B)"; if(AND(GasVisRatio>=1000,TotalMassFlux<20.5)) TwoPhMethod="Lockhart-Martinelli (L-M)"; if(GasVisRatio<1000) TwoPhMethod="Friedel"; var AvgVelocity; AvgVelocity=TotalMassFlow/(AvgDensity*3600*CSArea); var NreLiqTotal; NreLiqTotal=19.83*TotalMassFlow/(Math.PI*PipeID*LiqViscosity); var AAA,BBB,FFTotalLiqTurb,FFTotalLiqUse; AAA=POW(2.457*LOG(1/(POW(7/NreLiqTotal,0.9)+0.27*Roughness/PipeID)),16); BBB=POW(37530/NreLiqTotal,16); FFTotalLiqTurb=8*POW(POW(8/NreLiqTotal,12)+1/POW(AAA+BBB,1.5),1/12); if(NreLiqTotal<=2100) FFTotalLiqUse=64/NreLiqTotal; if(NreLiqTotal>2100) FFTotalLiqUse=FFTotalLiqTurb; /* Liquid Friction drop based on total flow*/ var DelPLiqTotal; DelPLiqTotal=4.167*FFTotalLiqUse*POW(TotalMassFlux,2)/(32.1741*PipeID*LiquidDensity); /* Gas Friction drop based on total flow*/ var NreGasTotal; var AAAA,BBBB; var FFTurbGasTotal; NreGasTotal=19.83*TotalMassFlow/(Math.PI*PipeID*GasViscosity); AAAA=POW(2.457*LOG(1/(POW(7/NreGasTotal,0.9)+0.27*Roughness/PipeID)),16); BBBB=POW(37530/NreGasTotal,16); FFTurbGasTotal=8*POW(POW(8/NreGasTotal,12)+1/POW(AAAA+BBBB,1.5),1/12); /* Friction factor Gas to be used*/ var FFUseGasTotal; if(NreGasTotal<=2100) FFUseGasTotal=64/NreGasTotal; if(NreGasTotal>2100) FFUseGasTotal=FFTurbGasTotal; /* Gas Friction drop*/ var DelPGasTotal; DelPGasTotal=4.167*FFUseGasTotal*POW(TotalMassFlux,2)/(32.1741*PipeID*GasDensityUse); /* C-B Method in detail*/ var PrFactor; PrFactor=POW(DelPGasTotal/DelPLiqTotal,0.5); var CBFactor; if(PrFactor<9.5) CBFactor=24.9/POW(TotalMassFlux,0.5); if(AND(PrFactor>=9.5,PrFactor<=28)) CBFactor=235.3/(PrFactor*POW(TotalMassFlux,0.5)); if(PrFactor>28) CBFactor=6788.5/(PrFactor*PrFactor*POW(TotalMassFlux,0.5)); var GasMassFraction; var Nfactor; GasMassFraction=GasFlow/(LiqFlow+GasFlow); Nfactor=(2-0.25)/2; var FiFactor; FiFactor=1+(PrFactor*PrFactor-1)*(CBFactor*POW(GasMassFraction,Nfactor)*POW(1-GasMassFraction,Nfactor)+POW(GasMassFraction,1.75)); var CBPressureDrop; CBPressureDrop=4.167*FiFactor*FFTotalLiqUse*POW(TotalMassFlux,2)/(32.1741*PipeID*LiquidDensity); /* L-M Method */ var PRLM; PRLM=LOG(POW(DelPLiq/DelPGas,0.5)); var FiLiqLM,FiGasLM; if(AND(NreLiq>2100,NreGas>2100)) FiLiqLM=1.44-0.508*PRLM+0.0579*PRLM*PRLM-0.000376*PRLM*PRLM*PRLM-0.000444*PRLM*PRLM*PRLM*PRLM; if(AND(NreLiq>2100,NreGas>2100)) FiGasLM=1.44+0.409*PRLM+0.0577*PRLM*PRLM-0.000352*PRLM*PRLM*PRLM-0.000432*PRLM*PRLM*PRLM*PRLM; if(AND(NreLiq>2100,NreGas<2100)) FiLiqLM=1.25-0.458*PRLM+0.067*PRLM*PRLM-0.00213*PRLM*PRLM*PRLM-0.000585*PRLM*PRLM*PRLM*PRLM; if(AND(NreLiq>2100,NreGas<2100)) FiGasLM=1.25+0.542*PRLM+0.067*PRLM*PRLM-0.00212*PRLM*PRLM*PRLM-0.000583*PRLM*PRLM*PRLM*PRLM; if(AND(NreLiq<2100,NreGas>2100)) FiLiqLM=1.24-0.484*PRLM+0.072*PRLM*PRLM-0.00127*PRLM*PRLM*PRLM-0.00071*PRLM*PRLM*PRLM*PRLM; if(AND(NreLiq<2100,NreGas>2100)) FiGasLM=1.24+0.516*PRLM+0.072*PRLM*PRLM-0.00126*PRLM*PRLM*PRLM-0.000706*PRLM*PRLM*PRLM*PRLM; if(AND(NreLiq<2100,NreGas<2100)) FiLiqLM=0.979-0.444*PRLM+0.096*PRLM*PRLM-0.00245*PRLM*PRLM*PRLM-0.00144*PRLM*PRLM*PRLM*PRLM; if(AND(NreLiq<2100,NreGas<2100)) FiGasLM=0.979+0.555*PRLM+0.096*PRLM*PRLM-0.00244*PRLM*PRLM*PRLM-0.00144*PRLM*PRLM*PRLM*PRLM; var DelLMLiqSep; var DelLMGasSep; DelLMLiqSep=POW(Math.exp(FiLiqLM),2)*DelPLiq; DelLMGasSep=POW(Math.exp(FiGasLM),2)*DelPGas; var DeltaPLMMethod; DeltaPLMMethod=Math.max(DelLMGasSep,DelLMLiqSep); /* Friedel Method */ var NreFriedelL; var NreFriedelG; NreFriedelL=19.83*TotalMassFlow/(Math.PI*PipeID*LiqViscosity); NreFriedelG=19.83*TotalMassFlow/(Math.PI*PipeID*GasViscosity); var FroudeN; FroudeN=12*TotalMassFlux*TotalMassFlux/(32.1741*AvgDensity*AvgDensity*PipeID); var SurfaceT; SurfaceT=fieldname128; var WeberN; WeberN=37.8*PipeID*TotalMassFlux*TotalMassFlux/(AvgDensity*SurfaceT); var FriedelCo1Liq,FriedelCo1Gas; if(NreFriedelL<1055) FriedelCo1Liq=64/NreFriedelL; if(NreFriedelL>1055) FriedelCo1Liq=POW((0.86859*LOG(NreFriedelL/(1.964*LOG(NreFriedelL)-3.8215))),-2); if(NreFriedelG<1055) FriedelCo1Gas=64/NreFriedelG; if(NreFriedelG>1055) FriedelCo1Gas=POW((0.86859*LOG(NreFriedelG/(1.964*LOG(NreFriedelG)-3.8215))),-2); var FriedelCo2; FriedelCo2=(1-GasMassFraction)*(1-GasMassFraction)+GasMassFraction*GasMassFraction*(LiquidDensity*FriedelCo1Gas/(GasDensityUse*FriedelCo1Liq)); var FiFriedleHorVerUp; FiFriedleHorVerUp=FriedelCo2+3.24*POW(GasMassFraction,0.78)*POW((1-GasMassFraction),0.24)*POW(LiquidDensity/GasDensityUse,0.91)*POW(GasViscosity/LiqViscosity,0.19)*POW((1-GasViscosity/LiqViscosity),0.7)*POW(FroudeN,-0.045)*POW(WeberN,-0.035); var FiFriedleVerDown; FiFriedleVerDown=FriedelCo2+38.5*POW(GasMassFraction,0.75)*POW((1-GasMassFraction),0.314)*POW(LiquidDensity/GasDensityUse,0.86)*POW(GasViscosity/LiqViscosity,0.73)*POW((1-GasViscosity/LiqViscosity),6.84)*POW(FroudeN,-0.0001)*POW(WeberN,-0.037); var DelPFriedelHorVerUp; DelPFriedelHorVerUp=4.167*FiFriedleHorVerUp*FriedelCo1Liq*TotalMassFlux*TotalMassFlux/(32.1741*LiquidDensity*PipeID); var DelPFriedelVerDown; DelPFriedelVerDown=4.167*FiFriedleVerDown*FriedelCo1Liq*TotalMassFlux*TotalMassFlux/(32.1741*LiquidDensity*PipeID); /* Pressure Drop to be shown */ var PressureDropShown; if(AND(GasVisRatio>=1000,TotalMassFlux>20.5)) PressureDropShown=CBPressureDrop; if(AND(GasVisRatio>=1000,TotalMassFlux<20.5)) PressureDropShown=DeltaPLMMethod; if(GasVisRatio<1000) PressureDropShown=DelPFriedelHorVerUp; /* Display the results */ getField(133).setVal(PREC(LiquidVelocity,3)); getField(134).setVal(PREC(NreLiq,3)); getField(135).setVal(PREC(FFUseLiq,3)); getField(138).setVal(PREC(DelPLiq,3)); getField(43).setVal(PREC(GasDensityUse,3)); getField(136).setVal(PREC(GasVelocity,3)); getField(137).setVal(PREC(NreGas,3)); getField(139).setVal(PREC(FFUseGas,3)); getField(140).setVal(PREC(DelPGas,3)); getField(142).setVal(PREC(FFTotalLiqUse,3)); getField(143).setVal(TwoPhMethod); getField(144).setVal(DelPLiqTotal); getField(145).setVal(PREC(TotalMassFlux,3)); getField(146).setVal(DelPGasTotal); getField(148).setVal(CBFactor); getField(147).setVal(FiFactor); getField(149).setVal(GasMassFraction); getField(150).setVal(PREC(GasVisRatio,3)); getField(151).setVal(GasDensityUse,3); getField(152).setVal(PRLM,3); getField(153).setVal(FiLiqLM,3); getField(154).setVal(FiGasLM,3); getField(155).setVal(DelLMLiqSep); getField(156).setVal(DelLMGasSep); getField(157).setVal(PREC(DeltaPLMMethod,3)); getField(158).setVal(PREC(FriedelCo1Liq,4)); getField(159).setVal(PREC(FriedelCo1Gas,4)); getField(160).setVal(PREC(FriedelCo2,4)); getField(161).setVal(PREC(FiFriedleHorVerUp,4)); getField(162).setVal(PREC(FiFriedleVerDown,4)); getField(163).setVal(PREC(DelPFriedelHorVerUp,3)); getField(164).setVal(PREC(DelPFriedelVerDown,3)); getField(165).setVal(PREC(AvgDensity,3)); getField(166).setVal(PREC(AvgViscosity,3)); getField(167).setVal(PREC(AvgVelocity,3)); getField(168).setVal(PREC(NreTotal,0)); getField(172).setVal(PREC(PressureDropShown,2)); return PREC(CBPressureDrop,3); })()