Численный метод расчета нестационарных режимов гидравлических систем (стр. 3 из 4)

Выбирая, для уравнений (15), (19) граничные условия:

, , для (16), (18) , и для (17), (21) , , в результате получаем:

Для вычисления интегралов в значениях коэффициентов

, строим кубический сплайн Эрмита [7] для подынтегральной функции . При этом считается, что кроме значений функции на концах отрезков задаются значения производных и , которые вычислялись при помощи формул (22), (23).

Тогда интегралы в формулах (31), (32) вычисляются с локальной погрешностью

, и глобальная погрешность будет равняться , что и подтверждается вычислительными экспериментами.

Также для сравнения были проведены расчеты уравнений (31), (32) с помощью формулы трапеций.

Отметим, что коэффициенты системы (30) получаются в итоге постоянными величинами.

Исследования свойств схемы

Для исследования устойчивости метода прогонки и монотонности схемы используем принцип максимума, т.е. если в системе уравнений (30)

, , и , то в матрице системы уравнений (30) присутствует диагональное преобладание, матрица является монотонной и выполняется оценка [6].

Рассмотрим коэффициенты системы уравнений (30):

Отсюда следует, что схема является монотонной и метод прогонки устойчив.

Устойчивость по времени нетрудно показать с помощью спектрального признака устойчивости [5].

Перепишем уравнение (9) в следующей форме:

и решение (33) будем искать в виде:

Тогда подставляя это выражение в (33) получим следующее уравнение:

и следовательно,

Очевидно, что

Отсюда следует, что схема абсолютно устойчива.

Для исследования порядка аппроксимации системы (5), (6) схемой (7), (8) разложим точное решение системы по формуле Тейлора в окрестности точки

, предполагая наличие непрерывных вторых производных:

Подставляя эти значения в (7), (8) и, используя, свойства

и системы (5), (6) получим невязку:

Таким образом, учтя погрешность при вычислении коэффициентов (31), (32) , система (5), (6) имеет аппроксимацию

.

Вычислительный эксперимент

Рассматривалась система трубопроводов, состоящая из 3 отрезков и 4 узлов. Длины всех отрезков одинаковы и равны 1, диаметры каналов одинаковы. Расчеты производились приведенным выше методом. Результаты расчетов сравнивались с точным решением системы (5), (6):

где

выбиралась:

1)

.

Начальные условия:

.

Граничные условия: при

при

2)

.

Начальные условия:

.

Граничные условия: при

при

3)

.

Начальные условия:

.

Граничные условия: при

при

Таблицы погрешностей

Расчеты на момент времени

.

1) Для первого случая:

		Сплайн Эрмита		Формула трапеций
tau	h	Общая погрешность	Погрешность по h	Общая погрешность	Погрешность по h
0,1	0,1	0,832851872281	0,000005587977	0,883060158026	0,049862541589
	0,05	0,832846634174	0,000000349870	0,845471064588	0,012273448151
	0,025	0,832846305931	0,000000021627	0,836006949938	0,002809333501
	0,0125	0,832846285402	0,000000001098	0,833636730148	0,000439113711
0,05	0,1	0,411112170340	0,000014415426	0,555115166343	0,142993370234
	0,05	0,411098669010	0,000000914096	0,447750454263	0,035628658154
	0,025	0,411097811602	0,000000056688	0,420301941958	0,008180145849
	0,0125	0,411097757795	0,000000002881	0,413401372060	0,001279575951
0,025	0,1	0,201166720391	0,000030952418	0,635665878379	0,431272044005
	0,05	0,201137819976	0,000002052003	0,316227793596	0,111833959222
	0,025	0,201135896782	0,000000128809	0,230331588708	0,025937754334
	0,0125	0,201135774540	0,000000006567	0,208461455105	0,004067620731

2) Для второго случая: