Расчетная сетка криволинейная

Сильное локальное изменение кривизны профиля в окрестности сопряжения сферы и обратного конуса приводит к необходимости выделять при численных расчетах три различные области течения газа — на сфере, в области сильного локального изменения кривизны профиля, на поверхности обратного конуса — и в соответствии с этим строить криволинейную расчетную сетку на поверхности тела. [c.301]

В случае расчета ОСН в сварных узлах при наличии криволинейных границ наиболее удобен МКЭ, что обусловлено отсутствием недостатков, присущих МКР (основные из которых трудность аппроксимации криволинейной области прямоугольной сеткой и равномерность шага сетки), иначе очень усложняется расчетная схема и теряется основное достоинство метода — простота. [c.278]

Использование ONDU T в случае сложной геометрической формы области возможно с применением особой технологии, описанной в этом параграфе. Эта технология является всего лишь вычислительной уловкой, разработанной в целях использования программы в случае областей сложной формы. За счет использования традиционных систем координат вычислительный метод и программа очень просты, удобны и эффективны. В то же время программу довольно непросто применить для областей сложной геометрической формы. Если бы в центре внимания находились области произвольной геометрической формы, то мы могли бы построить программу, использующую криволинейные неортогональные расчетные сетки или метод конечных элементов. Но тогда программа стала бы намного сложнее по структуре, а также трудна для понимания и использования. В представленной книге, сфокусированной на физическом понимании процессов, выбор простых систем координат кажется вполне оправданным. [c.116]

В 1961 году А. Ф. Сидоров возглавил один из ведущих математических отделов — отдел по расчету критических параметров и энерговыделения ядерных сборок, ему тогда было всего 28 лет. К этому времени относится начало его работы по созданию эффективных вариационных методов построения оптимальных криволинейных сеток, к которой он неоднократно возвращался впоследствии при решении самых различных прикладных задач. Оценивая результаты, полученные на этом направлении в самом начале его научной карьеры, достаточно лишь заметить, что методика и соответствующая программа, автоматизирующие процесс выбора расчетной сетки, разработанные А.Ф. Сидоровым тогда (в конце 50-х годов), до сих пор являются частью комплексов программ при расчете задач энерговыделения как в РФЯЦ-ВНИИТФ (г. Снежинск), так и в РФЯЦ-ВНИИЭФ (г. Саров). [c.6]

А.Ф. Сидоров уделял много внимания исследованиям, связанным с разработкой эффективных вариационных методов построения оптимальных криволинейных адаптивных сеток в двумерных и трехмерных областях сложных конфигураций, использующихся для решения задач механики сплошных сред (это было его хобби). Исследования были начаты А.Ф. Сидоровым в конце 50-х годов, когда он работал во ВНИИТФ. Им были предложены одномерный функционал, отвечающий за близость сетки к равномерной, и алгоритм построения сетки, обладающей достаточно хорошими аппроксимационными свойствами, созданы методика и программа, автоматизирующие процесс выбора одномерной расчетной сетки. [c.11]

Для проверки возможности переноса результатов, полученных для одномерного модельного уравнения (Б.1), на двумерные уравнения гидродинамики был проведен численный экспе-)имент с использованием программы Моретти (см. Моретти и Злейх [1968]) расчета обтекания затупленного тела невязким газом. Рассматривалось обтекание сферически затупленного конуса с полууглом раствора 6° совершенным газом с показателем адиабаты 7 = 1.4 при числе Маха невозмущенного потока, равном 10. Программа осуществляет выделение ударной волны на криволинейной расчетной сетке, перестраивающейся по мере изменения решения во времени. Поскольку ударная волна в процессе расчета все время сохраняется как разрыв, представленные результаты не искажаются послескачковыми всплесками, характерными для методов сквозного счета, или размазывания скачка. Для усиления влияния величины aes была выбрана чрезвычайно грубая сетка она содержала только три узла (две ячейки) между поверхностью тела и ударной волной и только пять узлов вдоль тела. Целью эксперимента являлось доказательство того, что стационарное решение, полученное по схеме Моретти, зависит от выбранной величины At, как это следует нз стационарного анализа величины е. (В этом состоит отличие схемы Моретти от схемы конечных разностей против потока, обсуждавшейся ранее, а также от ряда других конечно-разностных схем.) [c.524]

По мере прохождения технологического маршрута в результате выполнения локального окисления кремниевая поверхность перестанет быть планарной. При моделировании непланарной поверхности кремния или границы раздела Si - Si02 на прямоугольной сетке возникает криволинейная граница, произвольным образом пересекающая дискретные ячейки, на которые разделена область моделирования расчетной сеткой. В случае инертной окружающей среды на этой линии ставятся граничные условия отражения. Эти условия также можно легко представить с помощью фиктивных точек. Область моделирования ограничена при этом узловыми точками, находящимися в кремнии и являющимися соседними с криволинейной границей, а соседние с границей узловые точки, находящиеся вне кремния, являются фиктивными. Фиктивные точки применяются также при моделировании процессов предварительного осаждения примесей, однако в этом случае значение концентрации в них соответствует значению приповерхностной концентрации. [c.285]

Уравнения интегрировались по времени конечно-разностным методом второго порядка типа ЕУО в варианте, близком к схеме [12]. Течение в сопле с внезапным сужением имеет довольно сложную структуру. Для адекватного разрегнения его деталей, критичных в плане учета вязких эффектов, был развит достаточно общий подход [13], допускающий разбиение расчетной области на блоки четырех- или треугольной формы с криволинейными границами. Внутри блока сетка строилась посредством интерполяции. Вдоль каждой из границ блока возможно заданное сгущение сетки, что обеспечивало необходимую гибкость при описании областей сложной формы. [c.335]

При расчете трехмерных течений определяюгцая система уравнений записывалась в консервативной форме в произвольной неортогональной системе координат. Это позволяло использовать расчетную область с криволинейными границами и сгугцать сетки в областях с болыпими градиентами параметров. Параметры потока рассчитывались в центрах ячеек, а потоки — на их гранях. Конвективные потоки вычислялись с использованием противопоточной схемы с третьим порядком аппроксимации, диффузионные потоки на гранях определялись при помогци центральных разностей второго порядка точности [22]. [c.588]

Заметим, что при условии М>6В( =6АЬ практически отсутствует влияние водоотбора на поток противоположного берега, так что при выполнении этого условия водоток представляется внешней границей потока. Эффективным при этом может быть использование вблизи водотока криволинейной ортогональной сетки, переходящей иа удалении от водотока в прямолинейную сетку. Такую криволинейную сетку можно трансформировать в прямоугольную, исходя из равенства фильтрационных сопротивлений. При этом ближайшие к границе узловые точки следует задавать на линии водозаборных скважин, относя к ним водоотбор, равный суммарному дебиту попадающих в блок скважин. Расчетное (модельное) значение параметра сопротивления ложа водотока ALм при переходе к прямоугольной сетке с шириной блока Ах будет А щ=АЬАх 1в. [c.220]

Метод расчета. Примененный расчетный алгоритм основан на обобщенной процедуре глобальных итераций, предназначенной для решения конечно-объемным факторизованным методом уравнений переноса на многоблочных пересекающихся сетках О- и Н-типа. Система исходных уравнений записьшается в дельта-форме в криволинейных, согласованных с границами расчетной области координатах относительно приращений зависимых переменных, включающих декартовые составляющие скорости. После линеаризации система исходных уравнений решается с помощью согласованной неявной конечно-объемной процедуры коррекции давления [1], основанной на концепции расщепления по физическим процессам и записанной в -факторной формулировке. При этом для дискретизации временных производных используется схема второго порядка аппроксимации [10]. Для уменьшения влияния численной диффузии в расчетах течений с организованным отрывом потока, весьма чувствительных к ошибкам аппроксимации конвективных членов, в явной части уравнений переноса используется одномерный аналог противопоточной схемы с квадратичной интерполяцией [11]. Одновременно, чтобы избежать ложных осцилляций при воспроизводстве течений с тонкими сдвиговыми слоями, в неявной части уравнений использован механизм искусственной диффузии в сочетании с применением односторонних противопоточных схем для представления конвективных членов. В свою очередь, для устранения немонотонностей в распределении давления при дискретизации градиента давления по схеме с центральными разностями на согласованном (с совмещенными узлами для скалярных переменных и декартовых составляющих скорости) шаблоне в блок коррекции давления введен монотонизатор с эмпирическим сомножителем. Его величина 0.1 определена в ходе численных экспериментов на задаче обтекания цилиндра и шара потоком вязкой несжимаемой жидкости. Высокая эффективность вычислительной процедуры для решения дискретных алгебраических уравнений обеспечена применением метода неполной матричной факторизации. Более подробно детали описанной процедуры расчета течения на моноблочных сетках изложены в [11]. [c.46]

Энциклопедия по машиностроению XXL

Оборудование, материаловедение, механика и ...