>	restart:#"m11_p87"

Considérese un tanque cilíndrico de propano, de 1 m de diámetro y 3 m de longitud, enterrado bajo 1 m de arena. Determinar el flujo máximo de gas suministrable (o de pérdidas por fisuras).

>	with(RealDomain):read`../therm_eq.m`:read`../therm_const.m`:read`../therm_proc.m`:with(therm_proc):

>	su1:="C3H8":su2:="Arena":dat:=[D=1*m_,L=3*m_,H=1*m_];datP:=get_liq_data(su1);get_sol_data(su2):dat:=k=subs(%,k),op(dat),Const,SI2,SI1:k[arena]=subs(%,k);get_pv_data(su1):

Como el propano se almacena líquido (GLP), al ir saliendo gas se irá reponiendo con la vaporización de líquido, pero como esto requiere un aporte de energía para compensar la entalpía de cambio de fase, el gasto suministrable viene limitado por el flujo de calor entrante, que para un depósito enterrado será por conducción desde la superficie.

Si suponemos que inicialmete está todo atemperado a unos 15 ºC, al empezar a salir gas se irá enfriando (y disminuyendo la presión interior), hasta llegar al equilibrio energético. Si despreciamos la pérdida de presión del gas en el conducto de salida ( o a través de las fisuras), el flujo máximo tendrá lugar cuando el salto térmico sea máximo, l cual implica que el GLP se enfría hasta que la presión interior sea la atmosférica.

Sea 1 las condiciones atmosféricas (100 kPa y 15 ºC), y 2 las condiciones de enfriamiento máximo del depósito (100 kPa, Tb):

>	eqQmax:=T2=T[b];eqQmax:=T2=subs(datP,T[b]);Tini_:=subs(dat,T0);p_ini_:=pv(Tini_);

	(1)

i.e. si la Tinicial=15 ºC, el tanque de propano bifásico estaría a 725 kPa, y cuendo sale gas a flujo máximo su temperatura ha bajado a 231 K (-42 ºC).