![]() |
Главная Случайная страница Контакты | Мы поможем в написании вашей работы! | |
|
Следующее логическое уточнение должно аппроксимировать параболой переменный профиль между узлами Ρ и Ε, а не прямой. Чтобы создать параболу, мы должны использовать данные в еще одной точке; в соответствии с природой конвекции, третья точка взята в расположенном выше по потоку месте, то есть W, если поток от Ρ до Ε (то есть ux> 0) или EE, если ux < 0, см. рис. 1.10. Таким образом, получаем:
, (1.15)
где D, U, и UU обозначают основное направление, первый расположенный выше по потоку узел, и второй расположенный выше по потоку узел, соответственно (E, P, и W или Ρ, Ε, и EE, в зависимости от направления потока). Коэффициенты g2 и g1 могут быть выражены в терминах узловых координат:
;
.
Для регулярных сеток берутся коэффициенты трех узловых значений, включенных в интерполяцию: 3/8 для точки основного направления, 6/8 для первого расположенного выше по потоку узла, 1/8 для второго расположенного выше по потоку узла. Эта схема является несколько более сложной, чем схема CD: это расширяет вычислительную область на еще один узел в каждом направлении (в двухмерном случае, включены узлы EE, WW, NN и SS), а на неортогональных и/или неравномерных сетках, выражения для коэффициентов gi не просты. Леонард (1979) сделал эту схему популярной и дал ей название QUICK (Квадратичная противоточная интерполяция для конвективной кинематики).
У этой квадратичной схемы интерполяции есть погрешность метода третьего порядка и на равномерных и на неравномерных сетках. Это можно показать, устраняя вторую производную от уравнения (1.15) использование fw, который, на равномерной сетке в декартовой СК с ux> 0, приводит:
.
Первые три члена на правой стороне представляют аппроксимацию QUICK, в то время как последний член – основную погрешность метода. Когда используется эта схема интерполяции вместе с аппроксимацией правила середины поверхностного интеграла, полная аппроксимация второго порядка точности (точность аппроксимации квадратуры). Аппроксимация QUICK немного более точна, чем CD.
Дата публикования: 2015-02-22; Прочитано: 400 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!