AGA3 1990 Algorithm
AGA3 1990 Algorithm
org/
AGA 3 takes meter settings, compressibilities, a differential pressure, static pressure and temperature, and calculates base volume, base volume flow rate or volume correction factor. It requires flowing and base pressure, temperature and compressibility.
AGA3: Enter a flowing temperature, pressure and compressibility, a base temperature, pressure and compressibility, a flowing volume and a meter factor, then press Calculate.
Value Units ? Diff Press = Flowing Pressure = Flowing Temp = Downstream Pressure Tap? Orifice diameter = Pipe diameter = Flowing Compressibility = Base Compressibility = Ideal specific gravity = Base Pressure = Base Temp = Orifice diameter temp = Orifice therm exp. coeff. = Pipe diameter temp = Pipe therm exp. coeff. = Input
SI
10. 0 4101. 325 10. 0
US
1. 665e-5 20 1. 116e-5
mm in mm in Deg C Deg F kPaA PSIA Deg C Deg F Deg C Deg F mm/(mm Deg C) in/(in Deg F) Deg C Deg F mm/(mm Deg C) in/(in Deg F)
Results Table
Set unit conversion constants for the unit system chosen Calculate flowing density: Dflow = Pf * Mrair * Gi / (Zf * R * (Tf + N5)) Calculate base density: Dbase = Pb * Mrair * Gi / (Zb * R * (Tb + N5)) Calculate size of orifice at flowing conditions: d' = d * [1 + alpha1 * (Tf - Tdm)] Calculate size of meter run at flowing conditions: D' = D * [1 + alpha2 * (Tf TDm)] Calculate beta = d' / D' Calculate velocity of approach factor: Ev = 1 / sqrt( 1 - beta 4) Calculate factors required for coefficient of discharge / Reynolds number iteration loop. o Upstream tap position L1 = N4/D o Downstream tap position L2 = N4/D ( -8.5 L ) (-6.0 L ) o Upstream tap correction factor Tu = [ S2 + S3 e 1 + S4 e 1 ] * beta
1.1
o o o o o o o o o o o o o
Dimensionless downstream dam height M2 = (2 L2 ) / (1 - beta) Downstream tap correction factor Td = S6 [ M2 + S7 M2 1.3 ] beta 1.1 Small pipe correction factor: if D <= (A4 N4) then Ts = A3 (1 - beta) (A4 D/N4) else Ts = 0 Cd0 = A0 + A1 beta 2 + A2 B 8 + Tu + Td + Ts Cd1 = A5 * beta 0.7 (250) 0.7 Cd2 = A6 * beta 4 (250) 0.35 Cd3 = S1 * beta 4 * beta 0.8 (4.75) 0.8 (250) 0.35 Cd4 = (S5Tu + S8Td ) beta 0.8 (4.75) 0.8 Orifice differential pressure to flowing pressure ratio x = dP / (N3 * Pf) Yp = (0.41 + 0.35 beta 4) / k Expansion factor Y = 1 - Yp x Calculate Ftmp = Ev Y d 2 Calculate FIc = 4000 NIc D * mu / Ftmp
Calculate FIp = sqrt ( 2 * Dflow * dP) Calculate FI = FIc / FIp Limit FI to 1000 maximum. Set up coefficient of discharge loop by setting Cd(FT) to Cd0 Perform loop: o X = Fl / Cd(FT) o If X < Xc then 0.35 Fc = Cd0 + ( Cd1 X + Cd2 + Cd3 X 0.8 ) X 0.35 + Cd4 X 0.8 0.35 Dc = ( 0.7 Cd1 X + 0.35 Cd2 + 1.15 Cd3 X 0.8 ) X 0.35 + 0.8 Cd4 X
0.8
o o o
else
Fc = Cd0 + Cd1 X 0.7 + ( Cd2 + Cd3 X 0.8 ) ( A - B / X ) + Cd4 X 0.8 Dc = 0.7 Cd1 X 0.7 + ( Cd2 + Cd3 X 0.8 ) * B / X + 0.8 Cd3 ( A - B / X) X 0.8 + 0.8 Cd4 X 0.8 o Calculate the amount to vary Cd(FT) by: dCd = [ Cd (FT) - Fc] / [ 1 + Dc / Cd(FT) ] o Change Cd(FT) : Cd(FT) = Cd(FT) - dCd Until the absolute value of dCd is less than 0.000005. Reynolds number is Re = 4000 / X Calculate mass flow factor: F mass = ( pi / 4) * Nc Ev d2 Calculate mass flow rate: qm = F mass Cd(FT) Y * sqrt (2 Dflow * dP ) Calculate standard flow rate: qv = qm / D base