Python Kappa.Coefficientsの例

コード例 #1

ファイル: kappa_p.py プロジェクト: govtmirror/proteus-mprans

from proteus import *
from proteus.default_p import *
from jet_test2D import *
from proteus.mprans import Kappa

LevelModelType = Kappa.LevelModel

coefficients = Kappa.Coefficients(
    V_model=0,
    ME_model=3,
    LS_model=1,
    RD_model=None,
    dissipation_model=4,
    dissipation_model_flag=dissipation_model_flag,
    useMetrics=useMetrics,
    rho_0=rho_0,
    nu_0=nu_0,
    rho_1=rho_1,
    nu_1=nu_1,
    g=g,
    c_mu=0.09,
    sigma_k=1.0,
    sc_uref=kappa_sc_uref,
    sc_beta=kappa_sc_beta)


def getDBC_k(x, flag):
    if flag == boundaryTags['inflow']:
        return lambda x, t: kInflow
    if flag == boundaryTags['wall']:
        return lambda x, t: 0.0

コード例 #2

    LS_model = 2
    ME_model = 5
    dissipation_model = 6
#
dissipation_model_flag = 1
if useRANS == 2:
    dissipation_model_flag = 2
coefficients = Kappa.Coefficients(
    V_model=0,
    ME_model=ME_model,
    LS_model=LS_model,
    RD_model=RD_model,
    dissipation_model=dissipation_model,
    dissipation_model_flag=dissipation_model_flag,  #1 -- K-epsilon, 2 -- K-omega
    useMetrics=useMetrics,
    rho_0=rho_0,
    nu_0=nu_0,
    rho_1=rho_1,
    nu_1=nu_1,
    g=g,
    c_mu=0.09,
    sigma_k=1.0,
    sc_uref=kappa_sc_uref,
    sc_beta=kappa_sc_beta)


def getDBC_k(x, flag):
    if flag == boundaryTags['left']:
        return lambda x, t: kInflow
    if flag == boundaryTags['right']:
        return lambda x, t: 0.0

コード例 #3

ファイル: kappa_p.py プロジェクト: hhc2tech/air-water-vv

dissipation_model_flag = 1
if ct.useRANS == 2:
    dissipation_model_flag = 2
elif ct.useRANS == 3:
    dissipation_model_flag = 3

coefficients = Kappa.Coefficients(
    V_model=int(ct.movingDomain) + 0,
    ME_model=ME_model,
    LS_model=LS_model,
    RD_model=RD_model,
    dissipation_model=dissipation_model,
    dissipation_model_flag=
    dissipation_model_flag,  #1 -- K-epsilon, 2 -- K-omega 1998, 3 -- K-omega 1988
    useMetrics=ct.useMetrics,
    rho_0=ct.rho_0,
    nu_0=ct.nu_0,
    rho_1=ct.rho_1,
    nu_1=ct.nu_1,
    g=ct.g,
    nd=nd,
    c_mu=0.09,
    sigma_k=1.0,
    sc_uref=ct.kappa_sc_uref,
    sc_beta=ct.kappa_sc_beta)

dirichletConditions = {0: lambda x, flag: domain.bc[flag].k_dirichlet}
advectiveFluxBoundaryConditions = {
    0: lambda x, flag: domain.bc[flag].k_advective
}
diffusiveFluxBoundaryConditions = {

コード例 #4

ファイル: kappa_p.py プロジェクト: burgreen/proteus_cases

RD_model = None
LS_model = None
ME_model = 1
dissipation_model = 2

coefficients = Kappa.Coefficients(
    V_model=int(movingDomain) + 0,
    ME_model=ME_model,
    LS_model=LS_model,
    RD_model=RD_model,
    dissipation_model=dissipation_model,
    #1 -- K-epsilon, 2 -- K-omega 1998, 3 -- K-omega 1988
    dissipation_model_flag=dissipation_model_flag,
    useMetrics=user_param.useMetrics,
    rho_0=user_param.rho_water,
    nu_0=user_param.nu_water,
    rho_1=user_param.rho_air,
    nu_1=user_param.nu_air,
    #g=user_param.gravity,
    g=numpy.array(
        [user_param.gravity[0], user_param.gravity[1], user_param.gravity[2]],
        dtype='d'),
    nd=user_param.nd,
    c_mu=0.09,
    sigma_k=1.0,
    sc_uref=user_param.kappa_sc_uref,
    sc_beta=user_param.kappa_sc_beta)

dirichletConditions = {
    0: lambda x, flag: domain.bc[flag].k_dirichlet.init_cython()
}
advectiveFluxBoundaryConditions = {

コード例 #5

ファイル: turb_k_p.py プロジェクト: burgreen/proteus-msu-cases

RD_model = None
LS_model = None
ME_model = 1
dissipation_model = 2

coefficients = Kappa.Coefficients(
    V_model=0 + int(movingDomain),
    ME_model=ME_model,
    LS_model=LS_model,
    RD_model=RD_model,
    dissipation_model=dissipation_model,
    dissipation_model_flag=
    dissipation_model_flag,  # 1=K-epsilon, 2=K-omega 1998, 3=K-omega 1988
    useMetrics=useMetrics,  # main_param.useMetrics,
    rho_0=rho_0,  # main_param
    nu_0=nu_0,  # main_param
    rho_1=rho_1,  # main_param
    nu_1=nu_1,  # main_param
    g=gravity,
    nd=nd,  #main_param.nd,
    c_mu=0.09,
    sigma_k=1.0,
    sc_uref=kappa_sc_uref,  # main_param
    sc_beta=kappa_sc_beta  # main_param
)
'''
dirichletConditions = {
  0: lambda x, flag: domain.bc[flag].k_dirichlet.init_cython()
}

advectiveFluxBoundaryConditions = {

コード例 #6

ファイル: kappa_p.py プロジェクト: hhc2tech/air-water-vv

    ME_model = 5
    dissipation_model = 6
#
dissipation_model_flag = 1
if ct.useRANS >= 2:
    dissipation_model_flag = 2

coefficients = Kappa.Coefficients(
    V_model=0 + int(ct.movingDomain),
    ME_model=ME_model + int(ct.movingDomain),
    LS_model=LS_model + int(ct.movingDomain),
    RD_model=RD_model + int(ct.movingDomain),
    dissipation_model=dissipation_model + int(ct.movingDomain),
    dissipation_model_flag=dissipation_model_flag +
    int(ct.movingDomain),  #1 -- K-epsilon, 2 -- K-omega
    useMetrics=useMetrics,
    rho_0=rho_0,
    nu_0=nu_0,
    rho_1=rho_1,
    nu_1=nu_1,
    g=g,
    c_mu=ct.opts.c_mu,
    sigma_k=ct.opts.sigma_k,
    sc_uref=kappa_sc_uref,
    sc_beta=kappa_sc_beta)

kInflow = ct.kInflow

dirichletConditions = {
    0: lambda x, flag: domain.bc[flag].k_dirichlet.init_cython()
}

コード例 #7

dissipation_model_flag = 1
if ct.useRANS >= 2:
    dissipation_model_flag = 2

coefficients = Kappa.Coefficients(
    V_model=ct.V_model,
    ME_model=ct.K_model,
    LS_model=ct.LS_model,
    RD_model=ct.RD_model,
    dissipation_model=ct.EPS_model,
    SED_model=ct.SED_model,
    dissipation_model_flag=dissipation_model_flag +
    int(ct.movingDomain),  #1 -- K-epsilon, 2 -- K-omega
    useMetrics=useMetrics,
    rho_0=rho_0,
    nu_0=nu_0,
    rho_1=rho_1,
    nu_1=nu_1,
    g=g,
    nd=ct.nd,
    c_mu=ct.opts.Cmu,
    sigma_k=ct.opts.sigma_k,
    sc_uref=kappa_sc_uref,
    sc_beta=kappa_sc_beta,
    closure=ct.sedClosure)

kInflow = ct.kInflow

dirichletConditions = {
    0: lambda x, flag: domain.bc[flag].k_dirichlet.init_cython()