===================
/MAT/LAW37 (BIPHAS)
===================

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

Формат
------

.. csv-table::
   :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
   :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

   "/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"
   "Пш", "Пш", "", "", "", "", "", "", "", ""
   ":math:`\rho_{0}^{l}`", ":math:`\rho_{0}^{l}`", ":math:`C_{l}`", ":math:`C_{l}`", ":math:`\alpha_{l}`", ":math:`\alpha_{l}`", ":math:`\nu_{l}`", ":math:`\nu_{l}`", ":math:`\nu_{vol​_​l}`", ":math:`\nu_{vol​_​l}`"
   ":math:`\rho_{0}^{g}`", ":math:`\rho_{0}^{g}`", ":math:`\gamma`", ":math:`\gamma`", ":math:`P_{0}`", ":math:`P_{0}`", ":math:`\nu_{g}`", ":math:`\nu_{g}`", ":math:`\nu_{vol​_​g}`", ":math:`\nu_{vol​_​g}`"

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

.. csv-table::
   :header: "Поле", "Содержание", "Единица СИ   Пример"
   :widths: 33, 33, 33

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

Пример
------

.. code-block::

   #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----|

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

1. Для каждого

                        равновесие цикла рассчитывается между жидкостью и газом:
  :math:`\Delta P_{l}=\Delta P_{g}`
  Жидкий ЭОС:
  :math:`P_{l}=P_{0}+C_{l}\mu`
  с объемным модулем 
 :math:`C_{l}=\rho_{0}^{l}(c_{0}^{l})^{2}`
  Где,
  :math:`\rho_{0}^{l}`
  Начальная плотность жидкости
  :math:`c_{0}^{l}`
  Начальная скорость звука
  Газовое ЭОС:
  :math:`Pv^{\gamma}=P_{0}v_{0}^{\gamma}=constant`
  Где,
  P
  Давление
  :math:`v=\frac{V}{V_{0}}`
  Удельный объем
  Стресс Напряжение
  :math:`\sigma=P+\tau`
  :math:`\tau_{ij}=\lambda\dot{\epsilon}_{kk}\delta_{ij}+2\mu\dot{\epsilon}_{ij}`
  Где,
  :math:`\sigma`
  Стрессовое напряжение
  :math:`P`
  Сферическое напряжение напряжения
  :math:`\tau`
  Вязкое стрессовое напряжение
  :math:`\lambda`
  Сдвиговая вязкость (кинематическая)
  :math:`\mu`
  Дилатационная вязкость (кинематическая)
2. Мультиматериал LAW51 (/MAT/MULTIMAT) может имитировать LAW37.Антидиффузионный 
 техника с использованием метода MUSCL находится в стадии реализации. Это уже 
 доступно для LAW51 (/ALE/MUSCL).
3. Если необходим вход/выход, LAW51 должен быть 
 вместо этого используется. LAW51 имеет вход/выход для нескольких материалов. 
 составы.
4. В анимации 
 файлы:/ANIM/ELEM/DENS37 позволяет 
 выходной субматериал 
 плотности./ANIM/BRIC/VFRAC 
 позволяет выводить объемные доли для каждого 
 субматериал. См. временную историю на входе стартера Radioss. 
 Вручную или /ANIM/Eltyp/Retype во входных данных Radioss Engine. 
 Руководство.
5. По умолчанию
  :math:`P_{sh}=−P_{0}`
  . Это означает, что расчетное давление

                        сдвинуто так, чтобы в начальный момент значение выходного давления относительно
  :math:`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}`"
6. Внутренняя энергия не рассчитывается с помощью LAW37. Тем не менее, адиабатический

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

                        уравнение:
  :math:`PV^{\gamma}=const.`
7. LAW37 не совместим 
 с элементарными граничными условиями (/EBCS).