Cтраница 2
Основные особенности алгоритма построения расчетных сеток и постановки граничных условий на контактной поверхности в настоящей методике определяются тем, что вычислительная сетка строится в зависимости от геометрии препятствия. Большая формоизменяемость мягкой оболочки во время ее взаимодействия с газом обуславливает необходимость классификации различных положений оболочки с целью проведения однообразных вычислений. [16]
Расположение границы препятствия по отношению к семейству координатных линий и номера зон, разбивающих расчетную область на подобласти. [17] |
В этом случае применяется перестроечная фаза ПЛЭ метода. Она основана на рассмотрении движения среды относительно движущейся системы координат, скорость которой в конечно-разностных соотношениях интерпретируется как скорость движения узлов вычислительной сетки. Такой подход довольно прост и экономичен. [18]
Прогресс в развитии вычислительной техники и создание многопроцессорных вычислительных систем позволяют в приемлемые сроки получить решение рассмотренных задач с помощью алгоритмов интегрирования уравнений Эйлера модифицированным методом С. К. Годунова на подвижных сетках. Координаты узлов вычислительной сетки на нижней границе ( поверхности обтекаемого тела) изменяются в соответствии с законом его движения, а положение верхней границы в абсолютной системе координат определяется размером возмущенной области. Вследствие подвижности расчетной области вычислительная сетка перестраивается на каждом шаге интегрирования системы уравнений движения газа. [19]
Четыре типа форм оболочки. [20] |
Алгоритм метода в значительной степени зависит и от способа постановки граничных условий на поверхности контакта оболочки с газом. Непосредственное применение способов, используемых для объемных тел и заключающихся во введении в тело слоя фиктивных ячеек, в данном случае затруднено из-за того, что тонкая оболочка моделируется поверхностью, не имеющей толщины. Для того, чтобы можно было проводить расчеты как для внутренних, так и граничных ячеек на основе единого алгоритма, в вычислительной сетке оболочка моделируется в виде совокупности вырожденных ячеек нулевой толщины. Расчеты показали, что такая неравномерность площадей вносит заметную погрешность. Эта погрешность является следствием ошибки, допускаемой при вычислении массы в узле, прилегающем к треугольной ячейке. [21]
Если желательно резко уменьшить объем вычислений, то можно использовать даже модель со ступенчатой функцией из разд. Для электростатических линз можно, например, положить W ( z) Ez ( z), но, как мы увидим в разд. W ( z) производной наивысшего порядка, которая появляется в интеграле аберраций. Таким образом, задача сводится к поиску N ( 2M l) точек пересечения вычислительной сетки, которые будут задавать линейные отрезки оптимизированной функции. [22]