Beispiel #1
0
from proteus import *
from proteus.default_p import *
from dambreak_Ubbink_medium import *
from proteus.mprans import MCorr

LevelModelType = MCorr.LevelModel

coefficients = MCorr.Coefficients(LSModel_index=2,
                                  V_model=0,
                                  me_model=4,
                                  VOFModel_index=1,
                                  applyCorrection=applyCorrection,
                                  nd=nd,
                                  checkMass=False,
                                  useMetrics=useMetrics,
                                  epsFactHeaviside=epsFact_consrv_heaviside,
                                  epsFactDirac=epsFact_consrv_dirac,
                                  epsFactDiffusion=epsFact_consrv_diffusion)


class zero_phi:
    def __init__(self):
        pass

    def uOfX(self, X):
        return 0.0

    def uOfXT(self, X, t):
        return 0.0

from proteus import *
from proteus.default_p import *
from vortex import *
name = soname + "_phicor"
from proteus.mprans import MCorr

LevelModelType = MCorr.LevelModel
coefficients = MCorr.Coefficients(applyCorrection=applyCorrection,
                                  epsFactHeaviside=epsFactHeaviside,
                                  epsFactDirac=epsFactDirac,
                                  epsFactDiffusion=epsFactDiffusion,
                                  LSModel_index=0,
                                  V_model=0,
                                  me_model=3,
                                  VOFModel_index=2,
                                  checkMass=checkMass,
                                  nd=nd)


class zero_phi:
    def __init__(self):
        pass

    def uOfX(self, X):
        return 0.0

    def uOfXT(self, X, t):
        return self.uOfX(X)


analyticalSolutions = None
Beispiel #3
0
domain = ct.domain
nd = ct.domain.nd
mesh = domain.MeshOptions

genMesh = mesh.genMesh
movingDomain = ct.movingDomain
T = ct.T

LevelModelType = MCorr.LevelModel

coefficients = MCorr.Coefficients(LSModel_index=int(ct.movingDomain) + 2,
                                  V_model=int(ct.movingDomain) + 0,
                                  me_model=int(ct.movingDomain) + 4,
                                  VOFModel_index=int(ct.movingDomain) + 1,
                                  applyCorrection=ct.applyCorrection,
                                  nd=nd,
                                  checkMass=True,
                                  useMetrics=ct.useMetrics,
                                  epsFactHeaviside=ct.epsFact_consrv_heaviside,
                                  epsFactDirac=ct.epsFact_consrv_dirac,
                                  epsFactDiffusion=ct.epsFact_consrv_diffusion)


class zero_phi:
    def __init__(self):
        pass

    def uOfX(self, X):
        return 0.0

    def uOfXT(self, X, t):
Beispiel #4
0
from __future__ import absolute_import
from builtins import object
from proteus import *
from proteus.default_p import *
from .thelper_cons_ls import *
name = soname + "_phicor"
from proteus.mprans import MCorr

LevelModelType = MCorr.LevelModel
coefficients = MCorr.Coefficients(
    applyCorrection=applyCorrection,
    epsFactHeaviside=epsFactHeaviside,
    epsFactDirac=epsFactDirac,
    epsFactDiffusion=epsFactDiffusion,
    LSModel_index=LS_model,  #0
    V_model=V_model,  #0
    VOFModel_index=VOF_model,
    me_model=MCorr_model,
    checkMass=checkMass,
    nd=nd,
    useQuadraticRegularization=ct.STABILIZATION_TYPE_vof > 0,
    edgeBasedStabilizationMethods=ct.STABILIZATION_TYPE_vof > 0)


class zero_phi(object):
    def __init__(self):
        pass

    def uOfX(self, X):
        return 0.0
Beispiel #5
0
from __future__ import absolute_import
from builtins import object
from proteus.default_p import *
from proteus.mprans import MCorr
from .multiphase import *

LevelModelType = MCorr.LevelModel
coefficients = MCorr.Coefficients(LSModel_index=int(movingDomain) + 2,
                                  V_model=int(movingDomain) + 0,
                                  me_model=int(movingDomain) + 4,
                                  VOFModel_index=int(movingDomain) + 1,
                                  applyCorrection=applyCorrection,
                                  nd=nd,
                                  checkMass=False,
                                  useMetrics=useMetrics,
                                  epsFactHeaviside=ecH,
                                  epsFactDirac=epsFact_consrv_dirac,
                                  epsFactDiffusion=epsFact_consrv_diffusion)


class zero_phi(object):
    def __init__(self):
        pass

    def uOfX(self, X):
        return 0.0

    def uOfXT(self, X, t):
        return 0.0

domain = ct.domain
nd = domain.nd
mesh = ct.domain.MeshOptions

genMesh = mesh.genMesh
movingDomain = ct.movingDomain
T = ct.T

LevelModelType = MCorr.LevelModel

coefficients = MCorr.Coefficients(LSModel_index=2,
                                  V_model=0,
                                  me_model=4,
                                  VOFModel_index=1,
                                  applyCorrection=ct.applyCorrection,
                                  nd=ct.domain.nd,
                                  checkMass=False,
                                  useMetrics=ct.useMetrics,
                                  epsFactHeaviside=ct.ecH,
                                  epsFactDirac=ct.epsFact_consrv_dirac,
                                  epsFactDiffusion=ct.epsFact_consrv_diffusion)


class zero_phi:
    def __init__(self):
        pass

    def uOfX(self, X):
        return 0.0

    def uOfXT(self, X, t):