"""turbulence model variable entity (Turb class definition)"""
from typing import Any
import os
if os.environ.get('__USE_CUDA', '0') == '1':
import cupy as np
from cupy import ndarray
else:
import numpy as np
from numpy import ndarray
from g3d.defs import G3DEntity
from g3d.defs import FunctionByPass
from g3d.flow.vars import Auxv
from g3d.gas import Gas
[docs]
class Turb(G3DEntity):
def __init__(self, *args: Any, **kwargs: Any) -> None:
self._visct = 0.
super().__init__(*args, **kwargs)
def _set_default(self) -> None:
super()._set_default()
self._tag = 'g3d_turb'
def _set_from_args(self, *args: Any, **kwargs: Any) -> ndarray:
bps = self._get_function_evaluate_status(**kwargs)
if bps is FunctionByPass.SELF:
args_stat = np.ones(len(args))
else:
args_stat = super()._set_from_args(*args)
if bps in [FunctionByPass.EVALUATE, FunctionByPass.SELF]:
if np.sum(args_stat) > 0:
flow_data = []
for i in range(len(args)):
if args_stat[i] == 1:
arg = args[i]
if type(arg) is Turb:
args_stat[i] = 0
self.copy(arg)
elif type(arg) is list:
args_stat[i] = 0
if self._list_is_all_floats(arg) and len(arg) == 1:
self._visct = arg[0]
elif type(arg) in [float, np.float64, np.float32, int, np.int64, np.int32]:
args_stat[i] = 0
flow_data.append(np.float64(arg))
if len(flow_data) == 1:
self._set_from_args(flow_data, bps=FunctionByPass.SELF)
return args_stat
[docs]
def set_from_key_value(self, key: str, value: Any) -> bool:
if super().set_from_key_value(key, value):
return True
elif key == 'visct':
self._visct = value
return True
return False
@property
def visct(self) -> float:
"""
:return:
"""
return self._visct
@visct.setter
def visct(self, val: float) -> None:
self._visct = val
[docs]
def copy(self, turb: 'Turb') -> None:
super().copy(turb)
self._visct = turb.visct
[docs]
def update(self, auxv: Auxv, gas: Gas):
"""
:param auxv:
:param gas:
"""
pass
[docs]
def print(self) -> None:
print(' turbulence viscosity %11.5E [] ' % self._visct)