/MAT/LAW37 (BIPHAS)

Ключевое слово формата блока Описывает гидродинамический двухкомпонентный жидкий газовый материал. Не рекомендуется: использовать законы мультиматериалов с одинарной точностью Radioss двигатель.

Формат


/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID	/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/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_{0}^{l}\)	\(\rho_{0}^{l}\)	\(C_{l}\)	\(C_{l}\)	\(\alpha_{l}\)	\(\alpha_{l}\)	\(\nu_{l}\)	\(\nu_{l}\)	\(\nu_{vol_l}\)	\(\nu_{vol_l}\)
\(\rho_{0}^{g}\)	\(\rho_{0}^{g}\)	\(\gamma\)	\(\gamma\)	\(P_{0}\)	\(P_{0}\)	\(\nu_{g}\)	\(\nu_{g}\)	\(\nu_{vol_g}\)	\(\nu_{vol_g}\)

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

Поле	Содержание	Единица СИ Пример
mat_ID	Материал идентификатор.(Целое число, максимум 10 цифры)
unit_ID	Идентификатор объекта.(Целое число, максимум 10 цифр)
mat_title	Материал титул.(Персонаж, максимум 100 персонажи)
Пш	Давление смена. 5По умолчанию: \(P_{sh}=−P_{0}\) (Реал)	\([Pa]\)
\(\rho_{0}^{l}\)	Жидкость эталонная плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
\(C_{l}\)	Жидкий объем модуль.(Реальный)	\([Pa]\)
\(\alpha_{l}\)	Начальная масса фракция жидкая доля. = 0,0 Газ = 1,0 Жидкость (Реал)
\(\nu_{l}\)	сдвиг кинематическая вязкость ( \(=\mu/\rho_{0}^{l}\) ).(Реал)	\([\frac{m^{2}}{s}]\)
\(\nu_{vol_l}\)	Объемный вязкость (кинематическая) для жидкости: \(=\frac{3\lambda+2\mu}{\rho}\) По умолчанию = 0: Стокса Гипотеза ( \(3\lambda+2\mu=0\) )(Реал)	\([\frac{m^{2}}{s}]\)
\(\rho_{0}^{g}\)	Ссылка плотность газа.(Реальная)	\([\frac{kg}{m^{3}}]\)
\(\gamma\)	Идеальный газ константа.(Реальная)
\(P_{0}\)	Ссылка Начальное давление газа.(Реальное)	\([Pa]\)
\(\nu_{g}\)	сдвиг кинематическая вязкость = \(=\mu/\rho_{0}^{g}\) .(Реал)	\([\frac{m^{2}}{s}]\)
\(\nu_{vol_g}\)	Объемный вязкость (кинематическая) для газа: \(=\frac{3\lambda+2\mu}{\rho}\) По умолчанию = 0: Стокса Гипотеза ( \(3\lambda+2\mu=0\) ).(Реал)	\([\frac{m^{2}}{s}]\)

Пример

#RADIOSS STARTER

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

#-  1. MATERIALS:

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

/UNIT/1

unit for mat

                   g                  mm                  ms

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

/MAT/LAW37/1/1

air (unit g_mm_ms)

#                Psh

                   0

#             RHO_l0                 C_l             ALPHA_l                NU_l            NU_VOL_l

                .001                2089                   0            1.787e-3                   0

#             RHO_G0               GAMMA                  P0                NU_g            NU_VOL_g

             1.22E-6                 1.4                  .1             1.48e-2

/ALE/MAT/1

#     Modif. factor.

                   0

/MAT/LAW37/2/1

water (unit g_mm_ms)

#                Psh

                   0

#             RHO_l0                 C_l             ALPHA_l                NU_l            NU_VOL_l

                .001                2089                   1            1.787e-3                   0

#             RHO_G0               GAMMA                  P0                NU_g            NU_VOL_g

             1.22E-6                 1.4                  .1             1.48e-2

/ALE/MAT/2

#     Modif. factor.

                   0

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

#enddata

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

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

Для каждого

равновесие цикла рассчитывается между жидкостью и газом:

\(\Delta P_{l}=\Delta P_{g}\) Жидкий ЭОС: \(P_{l}=P_{0}+C_{l}\mu\) с объемным модулем

\(C_{l}=\rho_{0}^{l}(c_{0}^{l})^{2}\)
Где, \(\rho_{0}^{l}\) Начальная плотность жидкости \(c_{0}^{l}\) Начальная скорость звука Газовое ЭОС: \(Pv^{\gamma}=P_{0}v_{0}^{\gamma}=constant\) Где, P Давление \(v=\frac{V}{V_{0}}\) Удельный объем Стресс Напряжение \(\sigma=P+\tau\) \(\tau_{ij}=\lambda\dot{\epsilon}_{kk}\delta_{ij}+2\mu\dot{\epsilon}_{ij}\) Где, \(\sigma\) Стрессовое напряжение \(P\) Сферическое напряжение напряжения \(\tau\) Вязкое стрессовое напряжение \(\lambda\) Сдвиговая вязкость (кинематическая) \(\mu\) Дилатационная вязкость (кинематическая)

Мультиматериал LAW51 (/MAT/MULTIMAT) может имитировать LAW37.Антидиффузионный

техника с использованием метода MUSCL находится в стадии реализации. Это уже доступно для LAW51 (/ALE/MUSCL).

Если необходим вход/выход, LAW51 должен быть

вместо этого используется. LAW51 имеет вход/выход для нескольких материалов. составы.

В анимации

файлы:/ANIM/ELEM/DENS37 позволяет выходной субматериал плотности./ANIM/BRIC/VFRAC позволяет выводить объемные доли для каждого субматериал. См. временную историю на входе стартера Radioss. Вручную или /ANIM/Eltyp/Retype во входных данных Radioss Engine. Руководство.

По умолчанию

\(P_{sh}=−P_{0}\) . Это означает, что расчетное давление

сдвинуто так, чтобы в начальный момент значение выходного давления относительно

\(P_{0}\) значение 0,0. Чтобы использовать полное давление, P sh должно быть

обнулить, используя небольшое значение, например -1e-20 Давление

P можно экспортировать в анимацию с помощью /ANIM/ELEM/P . .. csv-table:
:header: "", "Относительное давление   Пш=-P0   (по умолчанию)", "Общее давление   Пш   = -1e-20"
:widths: 33, 33, 33
"Начальное состояние:", ":math:`P=\Delta P=0`", ":math:`P=P_{0}`"
"Расширение:", ":math:`P=\Delta P<0`", ":math:`P<P_{0}`"
"Сжатие:", ":math:`P=\Delta P>0`", ":math:`P>P_{0}`"

Внутренняя энергия не рассчитывается с помощью LAW37. Тем не менее, адиабатический

состояние

\(dE=−P⋅dV\) используется посредством политропного процесса

уравнение:

\(PV^{\gamma}=const.\)

LAW37 не совместим

с элементарными граничными условиями (/EBCS).