/MAT/K-EPS

Ключевое слово формата блока описывает: \(k−\epsilon\) вязкий турбулентный материал для жидкости.

Формат


/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID	/MAT/K-EPS/mat_ID/unit_ID
mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title
\(\rho_{i}\)	\(\rho_{i}\)	\(\rho_{0}\)	\(\rho_{0}\)
\(\nu\)	\(\nu\)	Пмин	Пмин
\(\rho_{0}k_{0}\)	\(\rho_{0}k_{0}\)	SSL	SSL
\(c_{\mu}\)	\(c_{\mu}\)	\(\sigma_{k}\)	\(\sigma_{k}\)	\(\sigma_{\epsilon}\)	\(\sigma_{\epsilon}\)	\(P_{r}/P_{rt}\)	\(P_{r}/P_{rt}\)
\(C_{1\epsilon}\)	\(C_{1\epsilon}\)	\(C_{2\epsilon}\)	\(C_{2\epsilon}\)	\(C_{3\epsilon}\)	\(C_{3\epsilon}\)
\(\kappa\)	\(\kappa\)	E	E	\(\alpha\)	\(\alpha\)	\(\chi_{t}\)	\(\chi_{t}\)

Определение

Поле	Содержание	СИ Пример устройства
mat_ID	Идентификатор материала.(Целое число, максимум 10 цифры)
unit_ID	Идентификатор объекта.(Целое число, максимум 10 цифр)
mat_title	Название материала.(Символ, максимум 100 персонажи)
\(\rho_{i}\)	Начальная плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
\(\rho_{0}\)	Эталонная плотность, используемая в E.O.S (уравнение состояние).По умолчанию = \(\rho_{i}\) (Реал)	\([\frac{kg}{m^{3}}]\)
\(\nu\)	Кинематическая вязкость.(Реальная)	\([\frac{m^{2}}{s}]\)
Пмин	Отключение давления.(Реальное)	\([Pa]\)
\(\rho_{0}k_{0}\)	Начальная турбулентная энергия (первая часть).(Реальная)	\([J]\)
SSL	Длина шкалы подсетки (первая часть). По умолчанию = 1e+10 (Реал)	\([m]\)
\(c_{\mu}\)	Коэффициент турбулентной вязкости (второй часть).По умолчанию = 0,09 (Реальное)
\(\sigma_{k}\)	к коэффициент диффузии (вторая часть).По умолчанию = 1,00 (Реал)
\(\sigma_{\epsilon}\)	Число диссипации Прандтля (второе часть).По умолчанию = 1,30 (Реальное)
\(P_{r}/P_{rt}\)	Ламинарное/турбулентное соотношение Прандтля (второе часть).По умолчанию = 0,7/0,9 (реальный)
\(C_{1\epsilon}\)	:math:`epsilon ` коэффициент уравнения 1 (третий часть).По умолчанию = 1,440 (Реальное)
\(C_{2\epsilon}\)	:math:`epsilon ` коэффициент уравнения 2 (третий часть).По умолчанию = 1,920 (Реальное)
\(C_{3\epsilon}\)	:math:`epsilon ` коэффициент уравнения 3 (третий часть).По умолчанию = -0,375 (Реальный)
\(\kappa\)	Постоянная стены Каппа (четвертая часть). По умолчанию = 0,4187 (Реал)
E	Постоянная стенки E (четвертая часть).По умолчанию = 9,7930 (Реальное)
\(\alpha\)	\(\kappa\) , \(\epsilon \)tau` эксцентрация (четвертая часть).По умолчанию = 0,5000 (Реальное)
\(\chi_{t}\)	Фактор исходного термина (четвертый часть).(Реальная)

Пример (Газ)

#RADIOSS STARTER

/UNIT/1

unit for mat

                  kg                   m                   s

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

/MAT/K-EPS/4/1

GAS

#              RHO_I               RHO_0

               .3828                   0

#                KNU                Pmin

             1.05E-4                   0

#            RHO0_K0                 SSL

                  20                   0

#               C_MU               SIG_k             SIG_EPS         P_R_ON_P_RT

                   0                   0                   0

#             C_1eps              C_2eps              C_3eps

                   0                   0                   0

#              KAPPA                   E              ALPHA                GSI_T

                   0                   0                   0                   0

/EOS/POLYNOMIAL/4/1

GAS

#                 C0                  C1                  C2                  C3

                   0                   0                   0                   0

#                 C4                  C5                  E0                Pmin               RHO_0

                 0.4                 0.4              253300                   0                1.22

/ALE/MAT/4

#     Modif. factor.

                   0

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

#enddata

/END

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Комментарии

\(S_{ij}=2\rho\nu_{eq}\dot{e}_{ij}\)

Где, \(S_{ij}\) Девиаторный тензор напряжений \(e_{ij}\) Тензор девиаторной деформации - Если элемент соединен с границей

В этом случае используется модель турбулентного пограничного слоя:

\(v_{eq}=max⁡(v,v⋅\kappa(\frac{y^{+}}{ln⁡(E⋅y^{+})}))\) \(y^{+}=\frac{c_{\mu}k^{2}}{\epsilon\alpha\kappa\nu}\) \(\chi=(1−\chi_{t})+\chi_{t}\alphaln(E⋅y^{+})\) Где,

\(\kappa\) это неспокойный кинетическая энергия.

Если соотношение между ламинарным и

турбулентные числа Прантля выше, чем

\(P_{r}/P_{rt}\) , затем: - Для ламинарного потока:

\(\nu_{eq}=\nu\)

Для турбулентного потока: \(\nu_{eq}=\nu+\frac{c_{\mu}k^{2}}{\epsilon}\)

Где,

\(\dot{\epsilon}\) это турбулентная диссипация и рассчитывается по следующей формуле уравнения:

\(\frac{d}{dt}\underset{V}{\int}\rho\epsilondV=\underset{S}{\int}\rho\epsilon(v−w)ndS+\underset{S}{\int}\frac{\mu_{t}}{\sigma_{\epsilon}}grad(\epsilon)ndS+\underset{V}{\int}S_{\epsilon}dV\) С, \(\mu_{t}=\frac{C_{\mu}⋅\kappa^{2}}{\epsilon}\)

(бурный вязкость)

\(G=\frac{\partialv_{i}}{\partialx_{j}}[\mu_{t}(\frac{\partialv_{i}}{\partialx_{j}}+\frac{\partialv_{j}}{\partialx_{i}}−\frac{2}{3}\rho\kappa\delta_{ij})−\frac{2}{3}\rho\kappa\delta_{ij}]\) \(S_{\epsilon}=\frac{\epsilon}{\kappa}(C_{1\epsilon}G−C_{2\epsilon}\rho\epsilon+C_{3\epsilon}\rho\kappa\frac{\partialv_{j}}{\partialx_{j}})\) Где, \(v\) Скорость материала \(w\) Скорость сетки

Уравнение состояния гидродинамического

давление необходимо задавать через

/EOS карта. Если

\(P=cst=0\)

, тогда

\(C_{1}\mu+\alpha_{\nu}T=0\)

, так что

\(\mu=−\frac{\alpha_{\nu}T}{C_{1}}\): Где, \(\mu\) Коэффициент дилатации \(\mu<0\) Расширение В этом случае параметры C 2 и C 3 не будет учитываться.

При использовании LAW6 в сочетании с

/MAT/LAW37 (BIPHAS) для жидкости

фаза (без газовой фазы), совместимость

жидкий EOS:

\(\Delta P_{1}=C_{1}\mu\)

для /MAT/LAW37 (BIPHAS)
\(p=C_{0}+C_{1}\mu+C_{2}\mu^{2}+C_{3}\mu^{3}+(C_{4}+C_{5}\mu)E\)

для LAW6, через полиномиальное EOS, определенное в примере выше,

тогда,

\(p=C_{1}\mu\)

При использовании LAW6 в сочетании с

/MAT/LAW37 (BIPHAS) для газа

фаза (без жидкой фазы), совместимость

EOS газа:

\(PV\gamma=const.\)

для /MAT/LAW37 (BIPHAS)
\(p=(\gamma−1)(\mu+1)E\)

, для LAW6, через

/EOS/IDEAL-GAS уравнение государство.

Где,

\(E\) это энергия на единицу объем.

Все тепловые данные (

\(\rho_{0}C_{p},T_{0},A,andB\) ) можно определить с помощью

ключевое слово

/HEAT/MAT .