proteus.TwophaseDarcyCoefficients module

An extension of the TransportCoefficients module for two-phase flow in porous media

Inheritance diagram of proteus.TwophaseDarcyCoefficients

class proteus.TwophaseDarcyCoefficients.TwophaseDarcyFlow_base(g=9.8, rhon=1.0, rhow=0.0, mun=1.0, muw=1.0, Ksw=1.0, psk_model='VGM', vg_alpha=5.0, vg_m=0.75, bc_pd=0.18281535648994515, bc_lambda=0.5, omega=1.0, Sw_max=1.0, Sw_min=0.0, density_w_model='Constant', density_n_model='Constant')[source]

Bases: proteus.TransportCoefficients.TC_base

setMaterialTypes(Ksw_types=[1.0], omega_types=[0.4], Sw_max_types=[1.0], Sw_min_types=[0.0], bc_lambda_types=None, bc_pd_types=None, vg_alpha_types=None, vg_m_types=None)[source]
setMaterialFunction(setParams)[source]
initializeMesh(mesh)[source]
initializeElementQuadrature(t, cq)[source]
initializeElementBoundaryQuadrature(t, cebq, cebq_global)[source]
initializeGlobalExteriorElementBoundaryQuadrature(t, cebqe)[source]
initializeGeneralizedInterpolationPointQuadrature(t, cip)[source]
class proteus.TwophaseDarcyCoefficients.TwophaseDarcy_fc(g=9.8, rhon=1.0, rhow=0.0, mun=1.0, muw=1.0, Ksw=1.0, psk_model='VGM', vg_alpha=5.0, vg_m=0.75, bc_pd=0.18281535648994515, bc_lambda=0.5, omega=1.0, Sw_max=1.0, Sw_min=0.0, density_w_parameters=None, density_n_parameters=None, diagonalHet=False, sparseDiffusionTensors={}, sd=True)[source]

Bases: proteus.TwophaseDarcyCoefficients.TwophaseDarcyFlow_base

evaluate(t, c)[source]
evaluate_sd_diagonal(t, c)[source]
twophaseDarcy_fc()[source]

Evaluate the coefficients of the fully coupled formulation of incompressible, two-phase Darcy flow

twophaseDarcy_fc_het_matType()[source]

Evaluate the coefficients of the fully coupled formulation of incompressible, two-phase Darcy flow for a heterogeneous medium

twophaseDarcy_fc_sd_diag()[source]

Evaluate the coefficients of the fully coupled formulation of incompressible, two-phase Darcy flow, sparse diffusion rep. assuming isotropic Ks

twophaseDarcy_fc_sd_diag_het_matType()[source]

Evaluate the coefficients of the fully coupled formulation of incompressible, two-phase Darcy flow for a heterogeneous medium, sparse diffusion rep. assuming isotropic Ks

class proteus.TwophaseDarcyCoefficients.TwophaseDarcy_fc_ff(g=9.8, rhon=1.0, rhow=0.0, mun=1.0, muw=1.0, Ksw=1.0, psk_model='VGM', vg_alpha=5.0, vg_m=0.75, bc_pd=0.18281535648994515, bc_lambda=0.5, omega=1.0, Sw_max=1.0, Sw_min=0.0, diagonalHet=False, sparseDiffusionTensors={}, sd=True)[source]

Bases: proteus.TwophaseDarcyCoefficients.TwophaseDarcyFlow_base

evaluate(t, c)[source]
evaluate_sd_diagonal(t, c)[source]
twophaseDarcy_fc_ff()[source]

Evaluate the coefficients of the fully coupled fractional flow formulation of incompressible, two-phase Darcy flow

twophaseDarcy_fc_ff_het_matType()[source]

Evaluate the coefficients of the fully coupled fractional flow formulation of incompressible, two-phase Darcy flow for a heterogeneous medium

twophaseDarcy_fc_ff_sd_diag()[source]

Evaluate the coefficients of the fully coupled fractional flow formulation of incompressible, two-phase Darcy flow, sparse diff. diagonal het

twophaseDarcy_fc_ff_sd_diag_het_matType()[source]

Evaluate the coefficients of the fully coupled fractional flow formulation of incompressible, two-phase Darcy flow for a heterogeneous medium, sparse diff. diagonal het

class proteus.TwophaseDarcyCoefficients.TwophaseDarcy_split_pressure(g=9.8, rhon=1.0, rhow=0.0, mun=1.0, muw=1.0, Ksw=1.0, psk_model='VGM', vg_alpha=5.0, vg_m=0.75, bc_pd=0.18281535648994515, bc_lambda=0.5, omega=1.0, Sw_max=1.0, Sw_min=0.0, swConstant=0.5, capillaryDiffusionScaling=1.0, nModel=1, diagonalHet=False, sparseDiffusionTensors={}, sd=True)[source]

Bases: proteus.TwophaseDarcyCoefficients.TwophaseDarcyFlow_base

attachModels(modelList)[source]
initializeElementQuadrature(t, cq)[source]
initializeElementBoundaryQuadrature(t, cebq, cebq_global)[source]
initializeGlobalExteriorElementBoundaryQuadrature(t, cebqe)[source]
initializeGeneralizedInterpolationPointQuadrature(t, cip)[source]
evaluate(t, c)[source]
evaluate_sd_diagonal(t, c)[source]
twophaseDarcy_split_pressure()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow

twophaseDarcy_split_pressure_het_matType()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow for a heterogeneous medium

twophaseDarcy_split_sd_diag_pressure()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow, sparse diffusion diagonal het

twophaseDarcy_split_sd_diag_pressure_het_matType()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow for a heterogeneous medium, sparse diffusion diagonal het

class proteus.TwophaseDarcyCoefficients.TwophaseDarcy_split_saturation(g=[9.8], rhon=1.0, rhow=1.0, mun=1.0, muw=1.0, Ksw=1.0, psk_model='VGM', vg_alpha=5.0, vg_m=0.75, bc_pd=0.18281535648994515, bc_lambda=0.5, omega=1.0, Sw_max=1.0, Sw_min=0.0, qScalarConstant=1.0, capillaryDiffusionScaling=1.0, nModel=0, diagonalHet=False, sparseDiffusionTensors={}, sd=True)[source]

Bases: proteus.TwophaseDarcyCoefficients.TwophaseDarcyFlow_base

twophaseDarcy_split_saturation()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow

twophaseDarcy_split_saturation_het_matType()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow for a heterogeneous medium

twophaseDarcy_split_sd_diag_saturation()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow, sparse diff. diagonal het

twophaseDarcy_split_sd_diag_saturation_het_matType()[source]

Evaluate the coefficients of the split fractional flow formulation of incompressible, two-phase Darcy flow for a heterogeneous medium, sparse diffusion diagonal

attachModels(modelList)[source]
initializeElementQuadrature(t, cq)[source]
initializeElementBoundaryQuadrature(t, cebq, cebq_global)[source]
initializeGlobalExteriorElementBoundaryQuadrature(t, cebqe)[source]
initializeGeneralizedInterpolationPointQuadrature(t, cip)[source]
evaluate(t, c)[source]
evaluate_sd_diagonal(t, c)[source]