Приложение 4. "Методические указания по составлению энергетической характеристики для систем транспорта тепловой энергии по показателю "Удельный расход сетевой воды". СО 153-34.20.523(2)-2003" (утв. Приказом Минэнерго РФ от 30.06.2003 N 278)

    Mixed scheme with circulation, with PT, direct scheme CO.
    Tn=(or<)Tnc,Ktp>0.            gt=?
    SetAtribute[{Tv,Tnp,T1p,T2p,T3p,gop,Kot,Tgp,Tx,Tc,Tn,T1,DTtp1,
    y,Ktp,a,m1,m2,n1,n2},Constant]
    Tv=
    Tnp=
    T1p=
    T2p=
    T3p=
    gop=
    Kot=
    Tgp=
    Tx=
    Tc=
    Tn=
    T1=
    DTtp1=
    y=
    Ktp=
    a=
    m1=
    m2=
    n1=
    n2=
    FindRoot[{T1k==T3k+(T1p-Т3р)*q,
    T2k==T3k-(T3p-T2p)*q,
                                                /\
    T3k==Tv+0.5*(T3p-T2p)*q+(0.5*(T3p+t2p)-Tv)*q  0.8,
    q==(Tv-Tn)/(Tv-Tnp),
    x==(T1tp-Tn)/(0.5*(T3k+T2k)+(T1k-0.5*(T3k+T2k))/(y*Kot)-Tn),
    T1tp==T1-DTtp1,
    T2o==T1tp-x*(T1k-T2k)/(y*Kot),
    c==0.8*Ktp*(0.5*(Tgp+Tc)-Tx)/(Tgp-Tc),
    P1==(1000/(0.5*(Tgp+Tc)-Tx))*a/(y*Kot*gop+a*gt),
    P2==(1000/(0.5*(Tgp+Tc)-Tx))*(1+c)/gt,
    A1==0.49*m1*n1*(1-P1)/Sqrt[P1],
    A2==0.63*m2*n2*(1-P2)/Sqrt[P2],
    Ts2==T1tp-(T1tp-Tgp)*(Exp[A2]-1)*P2/(1-P2),
    Th2==Tgp-(T1tp-Ts2)/P2,
    Th1==Tx+(Ts1-Tx)*(Exp[A1]-1)/(Exp[A1]-P1),
    Ts1==(a*gt*Ts2+y*Kot*gop*T2o)/(a*gt+y*Kot*gop),
    Th2==(Th1+Tc*c)/(1+c),
    T2==Ts1-(Th1-Tx)*P1},
    {T1k,120},{T2k,60},{T3k,80},{q,0.7},{x,1},{T1tp,115},
{T2o,60},{c,1},{P1,0.3},{P2,2.5},
    {A1,3},{A2,1.5},{Ts2,50},{Ts1,50},{Th2,50},{Th1,50},
{T2,40},{gt,15},MaxIterations->100]
    ClearAtribute[{Tv,Tnp,T1p,T2p,T3p,gop,Kot,Tgp,Tx,Tc,Tn,T1,
DTtp1,y,Ktp,a,m1,m2,n1,n2}]

    Mixed scheme
    with (without) circulation, with (without) PT, direct scheme
CO.
    Tn=(or>)Tnu.              gt=?
    SetAtribute[{Tv,Tnp,T1p,T2p,T3p,gop,Kot,Tgp,Tx,Tc,Tn,T1,DTtp1,
    y,Ktp,a,m1,m2,n1,n2},Constant]
    Tv=
    Tnp=
    T1p=
    T2p=
    T3p=
    gop=
    Kot=
    Tgp=
    Tx=
    Tc=
    Tn=
    T1=
    DTtp1=
    y=
    Ktp=
    a=
    m1=
    m2=
    n1=
    n2=

    FindRoot[{T1k==T3k+(T1p-Т3р)*q,
    T2k==T3k-(T3p-T2p)*q,
                                                /\
    T3k==Tv+0.5*(T3p-T2p)*q+(0.5*(T3p+T2p)-Tv)*q  0.8,
    q==(Tv-Tn)/(Tv-Tnp),
    x==(T1tp-Tn)/(0.5*(T3k+T2k)+(T1k-0.5*(T3k+T2k))/(y*Kot)-Tn),
    T1tp==T1-DTtp1,
    T2o==T1tp-x*(T1k-T2k)/(y*Kot),
    c==0.8*Ktp*(0.5*(Tgp+Tc)-Tx)/(Tgp-Tc),
    P1==(1000/(0.5*(Tgp+Tc)-Tx))*a/(y*Kot*gop+a*qt),
    P2==(1000/(0.5*(Tgp+Tc)-Tx))*(1+c)/gt,
    A1==0.46*m1*n1*(1-P1)/Sqrt[P1],
    A2==0.58*m2*n2*(1-P2)/Sqrt[P2],
    Ts2==T1tp-(T1tp-Tgp)*(Exp[A2]-1)*P2/(1-P2),
    Th2==Tgp-(T1tp-Ts2)/P2,
    Th1==Tx+(Ts1-Tx)*(Exp[A1]-1)/(Exp[A1]-P1),
    Ts1==(a*gt*Ts2+y*Kot*gop*T2o)/(a*gt+y*Kot*gop),
    Th2==(Th1+Tc*c)/(1+c),
    T2==Ts1-(Th1-Tx)*P1},
    {T1k,70},{T2k,40},{T3k,50},{q,0.2},{x,1.2},{T1tp,70},{T2o,50},
{c,1},{P1,0.3},{P2,0.9},
    {A1,2},{A2,0.4},{Ts2,50},{Ts1,50},{Th2,50},{Th1,50},{T2,30},
{gt,25},MaxIterations->100]
    ClearAtribute[{Tv,Tnp,T1p,T2p,T3p,gop,Kot,Tgp,Tx,Tc,Tn,T1,
DTtp1,y,Ktp,a,m1,m2,n1,n2}]

    3)  Алгоритм  определения  эксплуатационного удельного расхода
сетевой воды на ГВС при смешанной схеме включения подогревателей с
РТ,  с  циркуляцией  воды  и  при ее отсутствии и непосредственной
схеме  присоединения  систем  отопления  в  диапазонах  температур
наружного воздуха, прилегающих к t   ~= -5 °С
                                  нв
    Mixed scheme with (without) circulation, with PT, direct
scheme CO.
    Tn=-5(or><-5).           gt=?

    SetAtribute[{Tv,Tnp,T1p,T2p,T3p,gop,Kot,Tgp,Tx,Tc,Tn,T1,DTtp1,
    y,Ktp,a,m1,m2,n1,n2},Constant]
    Tv=
    Tnp=
    T1p=
    T2p=
    T3p=
    gop=
    Kot=
    Tgp=
    Tx=
    Tc=
    Tn=
    T1=
    DTtp1=
    y=
    Ktp=
    a=
    m1=
    m2=
    n1=
    n2=
    FindRoot[{T1k==T3k+(T1p-Т3р)*q,
    T2k==T3k-(T3p-T2p)*q,
                                               /\
    T3k==Tv+0.5*(T3p-T2p)*q+(0.5(T3p+T2p)-Tv)*q  0.8,
    q==(Tv-Tn)/(Tv-Tnp),
    x==(T1tp-Tn)/(0.5*(T3k+T2k)+(T1k-0.5*(T3k+T2k))/(y*Kot)-Tn),
    T1tp==T1-DTtp1,
    T2o==T1tp-x*(T1k-T2k)/(y*Kot),
    c==0.8*Ktp*(0.5*(Tgp+Tc)-Tx)/(Tgp-Tc),
    P1==(1000/(0.5*(Tgp+Tc)-Tx))*a/(y*Kot*gop+a*gt),
    P2==(1000/(0.5*(Tgp+Tc)-Tx))*(1+c)/gt,
    A1==0.49*m1*n1*(1-P1)/Sqrt[P1],
    A2==0.63*m2*n2*(1-P2)/Sqrt[P2],
    Ts2==T1tp-(T1tp-Tgp)*(Exp[A2]-1)*P2/(1-P2),
    Th2==Tgp-(T1tp-Ts2)/P2,
    Th1==Tx+(Ts1-Tx)*(Exp[A1]-1)/(Exp[A1]-P1),
    Ts1==(a*gt*Ts2+y*Kot*gop*T2o)/(a*gt+y*Kot*gop),
    Th2==(Th1+Tc*c)/(1+c),
    T2==Ts1-(Th1-Tx)*P1},
    {T1k,90},{T2k,50},{T3k,60},{q,0.5},{x,1},{T1tp,85},{T2o,50},
{c,1},(P1,0.3},{P2,4},{A1,2},{A2,-5},{Ts2,50},{Ts1,50},{Th2,50},
{Th1,50},{T2,40},{gt,10},MaxIterations->100]
    ClearAtribute[{Tv,Tnp,T1p,T2p,T3p,gop,Kot,Tgp,Tx,Tc,Tn,T1,
DTtp1,y,Ktp,a,m1,m2,n1,n2}]