Порядок выполнения работы

1. Записать словесную постановку задачи в соответствии с заданием варианта, выданного преподавателем.

2. Составить аналитическую модель описания задачи.

3. Разработать математическую модель решения.

4. Решить задачу с исходными данными, соответствующими вашему варианту. Получить выходные данные.

5. По выходным данным построить графики зависимости выходного параметра во времени в указанном интервале.

Пример расчета

Условие. В проточный аппарат трубчатого вида диаметром d = 0,05 м и длиной = 7 м поступает поток с объемным расходомv = 1,7 м³/ч, содержащий в нулевой момент времени по всей длине аппарата начальную концентрацию c_вых(0) = c_н = 10 моль/м³ (рис. 3).

При этом концентрация вещества в начале аппарата в любой момент времени есть некоторая входная функция от времени, изменяющаяся по следующему закону: . Определить концентрацию вещества в аппаратеc_вых(t) по его длине в момент времени T = 2 ч.

Решение

1. Словесная постановка задачи.

Исходными данными рассматриваемого процесса являются:

- диаметр трубчатого аппарата d = 0,05 м;

- длина аппарата =7 м;

- объемный поток  = 1,7 м³/ч;

- начальная концентрация вещества по всей длине аппарата в момент времени t = 0 была равна начальному значению c (0) = c_н = 10;

- концентрация вещества в начале аппарата в любой момент времени равна концентрации вещества во входном потоке ;

- исследуемый интервал времени T = [0, 2].

2. Синтез аналитической модели описания.

Для рассматриваемого процесса, проходящего в аппарате (см. рис. 3), принимаем следующие допущения:

- диаметр намного меньше длины (d<<);

- в поперечном направлении в каждом слое осуществляется идеальное смешение;

- в продольном направлении перемешивание отсутствует, т. е. вытеснение идеальное.

Концентрация исследуемого вещества в потоке аппарата зависит от двух переменных:

1) концентрация изменяется во времени на каждом бесконечно малом интервале сечения (x);

2) каждый момент времени происходит изменение концентрации вещества по всей длине рассматриваемого аппарата.

Следовательно, исследуемая концентрация есть функция двух переменных - времени (t) и координаты длины аппарата (x), т. е. c_вых = f(t,x). При исследовании динамики изменения c_вых по обеим переменным может быть использовано дифференциальное уравнение 1-го порядка в частных производных.

Исходя из вышесказанного, для описания структуры потока в аппарате (см. рис. 3) применяем следующую математическую модель:

, (5)

где w - линейная скорость передвижения вещества (м/ч),

,

где S - площадь поперечного сечения аппарата (м²).

, тогда .

Для решения дифференциального уравнения в частных производных (5) запишем начальное и граничное условия:

- начальное условие:

с_вых(x, 0) = c_н=10, при t = 0, 0 < x < , (6)

- граничное условие:

с_вых(0, t) = , при x = 0, t > 0. (7)

3. Разработка модели решения.

Математическая модель исследуемого процесса записана в виде дифференциального уравнения 1-го порядка в частных производных (5) с начальным и граничным условиями (6) и (7) соответственно. Для решения воспользуемся методом конечных разностей.

Левая часть уравнения (5) представляет собой производную по времени t (T=[0,2)). В конечных разностях она примет вид:

,

где t - бесконечно малый промежуток интервала времени, , гдеn - число разбиений интервала времени; i - номер интервала разбиения, целые числа, i = 0, 1, 2, …, (n-1), n.

Правая часть уравнения (5) представляет собой производную по координате длины аппарата x (x=[0, )). В конечных разностях она примет вид

,

где x - бесконечно малый участок сечения аппарата, , гдеm - число разбиений аппарата по длине; j - номер интервала разбиения, целые числа, j = 0, 1, 2, …, (m-1), m.

Алгоритм вывода численного решения уравнения (5).

Введем обозначения - безразмерная координата длины;- безразмерная координата концентрации;- безразмерная координата времени. Тогда, уравнение (5), начальное условие (6) и граничное условие (7) примут вид

Заменим в уравнении частные дифференциалы левосторонними конечными разностями относительно:

Таким образом, численное решение уравнения (5), основанное на методе конечных разностей, примет вид

. (8)

Начальное (6) и граничное (7) условия в конечных разностях будут иметь вид:

- начальное условие: _; (9)

- граничное условие: . (10)

Необходимым условием сходимости метода конечных разностей является условие

.

Так как ,и, получаем условие выбора значений интерваловn и m:

. (11)

Задаваясь одним из параметров (n или m), определяется другой. Например, если примем n = 100, тогда из условия (12) имеем:

.

Аналогично, задаваясь параметром m=const, определим параметр n.

Таким образом, модель решения уравнения (5) с начальным (6) и граничным (7) условиями будет состоять из (8), (9), (10), (11).

4. Решение задачи в среде MathСad.

Введем исходные данные задачи (рис. 4).

Рис. 4. Ввод исходных данных

Рассчитаем значение концентрации исходного вещества на входе в аппарат в каждой точке разбиения интервала времени. Для определения выражений, переменных и функций в несколько строк (при возвращении различных значений, как в нашем случае) в системе MathCAD используются программные модули.

Для вставки программного кода в документы используется специальная панель инструментов Programming (Программирование, рис. 5), которую можно вызвать на экран нажатием кнопки Programming Toolbar на панели Math (Математика).

Рис. 5. Панель инструментов Programming

Присваивание в пределах программ, в отличие от остальных документов MathСAD, производится с помощью оператора Local Definition (Локальное присваивание), который вставляется нажатием кнопки с изображением стрелки на панели Programming. Ни оператор присваивания «:=», ни оператор вывода «=» в пределах программ не применяются.

Определим с помощью программы функцию с_вх(t) (рис. 6).

Рис. 6. Вычисление значений функции с_вх(t)

Программный модуль обозначается в MathСAD вертикальной чертой, справа от которой последовательно записываются операторы языка программирования.

Чтобы создать программный модуль (см. рис. 6):

1. Введите часть выражения, которая будет находиться слева от знака присваивания и сам знак присваивания. В нашем примере это имя функции cv.

2. При необходимости вызовите на экран панель инструментов Programming (Программирование) (см. рис. 5).

3. Нажмите на этой панели кнопку Add Line (Добавить линию).

4. Если приблизительно известно, сколько строк кода будет содержать программа, можно создать нужное количество линий повторным нажатием кнопки Add Line (Добавить линию) соответствующее число раз (рис.).

Рис. 7. Начало создания программного модуля

5. В появившиеся местозаполнители введите желаемый программный код, используя программные операторы. В нашем примере в верхнем местозаполнителе вводим оператор For (при этом автоматически добавляется местозаполнитель для записи тела цикла), имя параметра цикла и диапазон его значений. В появившейся строке с помощью оператора Local Definition, вводим тело цикла, т. е. выражение, которое должно выполняться циклически. В нижнем местозаполнителе записано имя заполняемого массива. Заметим, что переменная s существует только внутри программы, выделенной вертикальной чертой. Из других мест документа получить ее значение невозможно.

После того как программный модуль полностью определен, и ни один местозаполнитель не остался пустым, функция может использоваться обычным образом, как в численных, так и в символьных расчетах. Для того чтобы вывести посчитанные значения после программного кода, введем выражение “cv=” и нажмем клавишу Enter, на экран будет выведена таблица решений (см. рис. 6).

Далее для решения дифференциального уравнения в частных производных (5) запишем начальное и граничное условия (6), (7) и численное решение уравнения (5) – формулу (10). Для этого необходимо записать следующий программный код (рис. 8).

Рис. 8. Заполнение массива концентрации исходного

вещества

Для вывода посчитанных значений, запишем выражение “cw=” и нажмем клавишу Enter (рис. 9).

5. Построение графиков.

Теперь по полученным данным построим графики зависимости выходного параметра во времени в указанном интервале. Сначала построим график изменения концентрации по длине аппарата в конечный момент времени. Для этого на макете для построения графиков в позиции оси абсцисс вводим параметр i, пределы изменения которого задаются перед построением, а в положении оси ординат – обозначение последнего столбца матрицы значений функции с_вых(t) (рис. 10).

Рис. 9. Массив концентрации исходного вещества

Рис. 10. Зависимость распределения концентрации вещества с_вых(t)

по длине аппарата в момент времени t= 2 ч

Аналогично строим график изменения концентрации исходного вещества в конце аппарата в каждый момент времени, при этом на оси абсцисс отображается изменение параметра j (начальное и конечное значение которого предварительно задаются), на оси ординат – интересующие нас значения выходного параметра, что соответствует последней строке матрицы значений функции с_вых(t) (рис. 11).

Рис. 11. Зависимость распределения концентрации вещества с_вых(t)

в каждый момент времени tв конце аппарата