========= /BEM/FLOW ========= Ключевое слово формата блока Описывает поток несжимаемой жидкости по границе. метод элементов. Формат ------ .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID", "/BEM/FLOW/flow_ID/unit_ID" "flow_title", "flow_title", "flow_title", "flow_title", "flow_title", "flow_title", "flow_title", "flow_title", "flow_title", "flow_title" "surf_IDex", "Нио", "внутри", "fct_IDfsp", "Фскалефсп", "Фскалефсп", "Аскалефсп", "Аскалефсп", "", "" "grn_IDaux", "Итест", "Толе", "Толе", "", "", "", "", "", "" "Ро", "Ро", "Ивинф", "", "", "", "", "", "", "" If N io > 0, вставить N io временная линия .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "surf_IDio", "fct_IDnv", "fct_IDp", "", "Фскаленв", "Фскаленв", "Фскалеп", "Фскалеп", "Аскалет", "Аскалет" Флаг формулировки .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "Яформа", "Ипри", "Dtflow", "Dtflow", "", "", "", "", "", "" "fct_IDv", "Фскалев", "Фскалев", "Аскалёв", "Аскалёв", "", "", "", "", "" "Диркс", "Диркс", "Дири", "Дири", "Дирз", "Дирз", "", "", "", "" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "flow_ID", "Несжимаемый проточный блок идентификатор.(Целое число, максимум 10 цифр)", "" "unit_ID", "(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)", "" "flow_title", "Несжимаемый проточный блок заголовок.(Символ, максимум 100 символов)", "" "surf_IDex", "Поток внешней поверхности идентификатор.(Целое число)", "" "Нио", "Количество притоков-оттоков поверхности.(Целое число)", "" "внутри", "Флаг внутреннего или внешнего потока. = 1 (по умолчанию) Поток вычисляется внутри поверхности, определяемой формулой surf_IDex. Нормали элемента поверхности должны быть ориентированы наружу. = 2 Поток вычисляется за пределами поверхности, определяемой формулой surf_IDex. Нормали элемента поверхности должны быть ориентированы внутрь. (Целое число)", "" "fct_IDfsp", "Кривая давления застоя идентификатор.(Целое число)", "" "Фскалефсп", "Шкала давления застоя коэффициент.По умолчанию = 1,0 (Реальный)", ":math:`[Pa]`" "Аскалефсп", "Масштабный коэффициент по оси абсцисс для кривая давления застоя. По умолчанию = 1,0 (реальное)", ":math:`[s]`" "grn_IDaux", "Группа вспомогательных узлов идентификатор. 2 (целое число)", "" "Итест", "Флаг проверки вспомогательных узлов. 2(Целое число > 0)", "" "Толе", "Размерный допуск. 2По умолчанию = 1.e-5 (Реал)", "" "Ро", "жидкость плотность.(Реальная)", ":math:`[\frac{kg}{m^{3}}]`" "Ивинф", "Дополнительное поле скоростей флаг. 3(Целое число > 0)", "" "surf_IDio", "Поверхность притока-оттока идентификатор. 4 (целое число)", "" "fct_IDnv", "Нормальная кривая скорости. 4 (целое число)", "" "fct_IDp", "Наложенная кривая давления. 5 (целое число)", "" "Фскаленв", "Нормальная шкала скоростей коэффициент.По умолчанию = 1,0 (Реальный)", ":math:`[\frac{m}{s}]`" "Фскалеп", "Шкала наложенного давления коэффициент.По умолчанию = 1,0 (Реальный)", ":math:`[Pa]`" "Аскалет", "Масштабный коэффициент по оси абсцисс для кривая нормальной скорости и кривая вынужденного давления. По умолчанию = 1,0. (Реал)", ":math:`[s]`" "Яформа", "Флаг формулировки. 6 = 1 (по умолчанию) Поток жидкости рассчитывается с использованием БЭМ с коллокацией подход к решению интегрального уравнения. = 2 Поток жидкости рассчитывается с использованием БЭМ с Галеркиным. подход к решению интегрального уравнения. (Целое число > 1)", "" "Ипри", "Выходной уровень.(Целое число > 1)", "" "Dtflow", "Шаг по времени для матриц БЭМ сборка. 7 = 0 (по умолчанию) :math:`\Delta t` использовался для обновления до BEM матрицы. ≠ 0 :math:`max(Dt_{flow},\Delta t)` использовался для обновления до BEM матрицы. (Настоящий)", ":math:`[s]`" "fct_IDv", "Кривая скорости идентификатор.(Целое число)", "" "Фскалев", "Шкала скорости коэффициент.По умолчанию = 1,0 (Реальный)", ":math:`[\frac{m}{s}]`" "Аскалёв", "Масштабный коэффициент по оси абсцисс для кривая скорости. По умолчанию = 1,0 (Реальное)", ":math:`[s]`" "Диркс", "X-компонент вектор направления дополнительного поля.(Реальный)", "" "Дири", "Y-компонент вектор направления дополнительного поля.(Реальный)", "" "Дирз", "Z-компонент вектор направления дополнительного поля.(Реальный)", "" Пример ------ В этом примере установите нормальную скорость на поверхности. surf_ID io = 2 и давление на поверхность surf_ID io = 1 . Вспомогательные узлы находятся внутри замкнутой поверхности. surf_ID ex = 3 . .. code-block:: #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /BEM/FLOW/1 Flow 1 #surf_IDex Nio Iinside fct_IDfsp Fscale_fsp Ascale_fsp 3 2 1 0 0. 0. #grn_IDaux Itest Tole 1 0 1e-5 # Rho Ivinf 1.0 0 #surf_IDio fct_IDnv fct_IDp Fscale_nv Fscale_p Ascale_t 2 1 0 10.0 0.0 1.0 #surf_IDio fct_IDnv fct_IDp Fscale_nv Fscale_p Ascale_t 1 0 1 0.0 101325.0 1.0 # Iform Ipri Dt_flow 1 1 1e-3 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/1 Function 1 # X Y 0 1 100 1 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #enddata .. image:: images/bem_flow_starter_obsolete_r_bem_flow_example.jpg *(Рисунок 1.)* Комментарии ----------- 1. surf_IDex должен определить закрытый поверхность. 2. Используя БЭМ, поток потенциал, скорость и давление вычисляются для узлов, принадлежащих поверхности определяется surf_IDex.Для визуальные проблемы и проблемы после обработки, можно вычислить характеристики потока для множества узлов внутри потока, принадлежащего grn_IDaux.Если Итест = 1, действительно ли расположены вспомогательные узлы внутри (если Iinside =1) или снаружи (если внутри =2), поверхность, определяемая формулой surf_IDex на каждом шаг по времени проверяется. Неправильные узлы затем отменяются на текущий момент. Step.Tolerance Tole используется для выполнения Тест «точка внутри закрытой поверхности». 3. Флаг Ивинф эффективен только для расчета потока в неограниченной области вне поверхности, определяемой формулой surf_IDex (Внутри =2).Если Ивинф = 1, условие притока определяется дополнительным однородный поток, определенный в свободном пространстве. Рассчитанный поток будет идентичен дополнительному потоку на бесконечном расстоянии от поверхности, определяемому формулой surf_IDex. 4. Если внутри = 0: поток вычисляется внутри поверхности, определяемой surf_IDex. Должно быть в хотя бы одна поверхность, на которой задана нормальная скорость, и только одна поверхность где нормальную скорость можно оставить свободной. Скорость на свободном поверхность будет рассчитываться благодаря равновесию потоков на глобальной поверхности определяется surf_IDex.Если внутри = 2 и Ивинф = 0: внутри = 0, но поток вычисляется за пределами поверхности, определяемой формулой surf_IDex.Если внутри = 2 и Ивинф = 1: номера поверхностей могут быть свободными и нормальными Скорость должна быть задана всем из них. 5. Чтобы снизить давление из поля скоростей необходимо приложить одно и только одно давление для весь расчет потока: это может быть либо глобальное давление застоя, либо давление на одной из поверхностей притока-оттока. 6. Коллокационный подход – это быстрее, но может оказаться недостаточно надежным для обработки очень сложной геометрии. Подход Галеркина работает в любой ситуации, но существенно медленнее. 7. Матрицы БЭМ зависят только от геометрия поверхности. Если Dtпоток = 0 (по умолчанию), они собираются в каждом цикле моделирование (шаг по времени классически задается условием устойчивости конечных элементов). Если Dtflow ≠ 0: макс(Dtflow, Dt) используется для обновления матриц БЭМ; где Dt — шаг по времени для конечных элементов.