コード例 #1
0
 def __init__(self, wrm=None, coef=1, dT=1.):
     if wrm is None:
         wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
     self.smooth_coef = coef
     self.dp = capillary_pressure.pc_ice(273.15-dT, 101325., 1)
     print self.dp
     self.pm = permafrost_model_explicit_fpd.PermafrostModel(wrm)
コード例 #2
0
 def __init__(self, wrm=None, coef=1, dT=1.):
     if wrm is None:
         wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
     self.smooth_coef = coef
     self.dp = capillary_pressure.pc_ice(273.15 - dT, 101325., 1)
     print self.dp
     self.pm = permafrost_model_explicit_fpd.PermafrostModel(wrm)
コード例 #3
0
    def wc(self, T, p):
        pcliq = capillary_pressure.pc_liq(T, p)
        pcice = capillary_pressure.pc_ice(T, p)
        sg, sl, si = self.pm.saturations_Tp(T, p)

        ni = eos.n_ice(T, p)
        nl = eos.n_water(T, p)
        ng = eos.n_gas(T, p)
        om = eos.mol_frac_gas(T)

        return sg * om * ng + sl * nl + si * ni
コード例 #4
0
ファイル: wc_T.py プロジェクト: amanzi/ats-dev 
    def wc(self, T, p):
        pcliq = capillary_pressure.pc_liq(T,p)
        pcice = capillary_pressure.pc_ice(T,p)
        sg,sl,si = self.pm.saturations_Tp(T,p)

        ni = eos.n_ice(T,p)
        nl = eos.n_water(T,p)
        ng = eos.n_gas(T,p)
        om = eos.mol_frac_gas(T)

        return sg*om*ng + sl*nl + si*ni
コード例 #5
0
 def __init__(self, wrm=None, coef=1., dT=1.):
     if wrm is None:
         wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
     self.wrm = wrm
     self.coef = coef
     self.dp = capillary_pressure.pc_ice(273.15 - dT, 101325., coef)
コード例 #6
0
 def saturations_Tp(self, T, p):
     pcliq = capillary_pressure.pc_liq(T, p)
     pcice = capillary_pressure.pc_ice(T, p, self.coef)
     return self.saturations_pc(pcliq, pcice)
コード例 #7
0
 def dsaturations_dpcice_Tp(self, T, p):
     pcliq = capillary_pressure.pc_liq(T,p)
     pcice = capillary_pressure.pc_ice(T,p)
     return self.dsaturations_dpcice_pc(pcliq, pcice)
コード例 #8
0
ファイル: permafrost_model.py プロジェクト: amanzi/ats-dev 
 def saturations_Tp(self, T, p):
     pcliq = capillary_pressure.pc_liq(T,p)
     pcice = capillary_pressure.pc_ice(T,p)
     return self.saturations_pc(pcliq,pcice)
コード例 #9
0
 def __init__(self, wrm=None, coef=1., dT=1.):
     if wrm is None:
         wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
     self.wrm = wrm
     self.coef = coef
     self.dp = capillary_pressure.pc_ice(273.15-dT, 101325., coef)