Главная страница
Навигация по странице:

  • КУРСОВАЯ РАБОТА по дисциплине «

  • Задача.

  • ∆σ=σ

  • Г=(В/

  • Решение задачи по расчёту удельной адсорбции на границе раздела жидкость-газ. Курсовая работа по дисциплине Информационные технологии в защите окружающей среды


    Скачать 423.5 Kb.
    НазваниеКурсовая работа по дисциплине Информационные технологии в защите окружающей среды
    АнкорРешение задачи по расчёту удельной адсорбции на границе раздела жидкость-газ
    Дата03.12.2021
    Размер423.5 Kb.
    Формат файлаdoc
    Имя файлаPEChATAJ.doc
    ТипКурсовая
    #289753

    МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ
    КАЗАНСКИЙ ГОСУДАРСТВЕННЫЙ АРХИТЕКТУРНО-СТРОИТЕЛЬНЫЙ УНИВЕРСИТЕТ
    Кафедра химии и инженерной экологии в строительстве

    КУРСОВАЯ РАБОТА
    по дисциплине «Информационные технологии в защите окружающей среды»
    на тему: "Решение задачи по расчёту удельной адсорбции

    на границе раздела жидкость-газ"

    Выполнил: студент гр. 0ИЗ01

    Мингазов Д. Р.

    Преподаватель: Громаков Н.С.

    Защитил: _________ ________

    (оценка) (дата)


    Казань

    2020

    Задача. Рассчитать удельную адсорбцию гептановой кислоты и построить изотерму адсорбции по зависимости σ = f(с) для её водных растворов при Т = 293 К.

    Решение. 1) Работа проводится с применением средств информационных технологий в виде компьютерной программы Excel. На основании экспериментальных данных, приведенных в таблице 1 строю график зависимости (рис.1) поверхностного натяжения раствора от концентрации гептановой кислоты σ=f(c).

    Таблица 1. Поверхностное натяжение водных растворов гептановой кислоты
    при Т=293 К.


    C*10,моль/мᶾ

    σ·10ᶾ , Н/м

    0

    72,75

    0,100

    72,40

    0,398

    72,00

    1,000

    71,20

    1,259

    70,60

    1,585

    69,70

    1,995

    68,40

    2,512

    67,00

    3,162

    65,30

    3,981

    63,60

    5,012

    61,80

    6,310

    60,10

    7,943

    58,40

    10,000

    56,70

    

    Таблица 2. Расчётные значения σ растворов гептановой кислоты при равноотстоящих значениях концентрации

    c,моль/мᶾ

    σ·10ᶾ , Н/м

    0

    73,25

    1,25

    69,96

    2,5

    67,01

    3,75

    64,42

    5

    62,17

    6,25

    60,26

    7,5

    58,71

    8,75

    57,50

    10

    56,63

    2) По точкам на рис.1 с помощью редактора электронных таблиц MicrosoftExcel провожу линию тренда в виде полиноминальной кривой второй степени как наиболее точно описывающей экспериментальную кривую, об этом свидетельствует величина достоверности аппроксимации близкая к единице (R2=0,9944). Уравнение линии тренда (1) приведено на рис.1. С помощью полученного уравнения зависимости σ=f(с) рассчитываю значения σ при равноотстоящих значениях концентрации с (шаг l=0,125). Полученные расчетные данные представлены в таблице 2. Ниже приведены алгоритмы для расчета этих величин:

    1) 69,9563=0,1109*1,252-2,7704*1,25+73,246

    2) 67,0131=0,1109*2,52-2,7704*2,5+73,246

    3) 64,4165=0,1109*3,752-2,7704*3,75+73,246

    4) 62,1665=0,1109*52-2,7704*5+73,246

    5) 60,263=0,1109*6,252-2,7704*6,25+73,246

    6) 58,7061=0,1109*7,52-2,7704*7,5+73,246

    7) 57,4958=0,1109*8,752-2,7704*8,75+73,246

    8) 56,632=0,1109*102-2,7704*10+73,246

    По рассчитанным значениям на этом же рисунке (рис.1) строю график для равноотстоящих точек. Эти данные потребуются при расчете производной dσ/dc.



    3) Для расчёта и построения изотермы адсорбции по адсорбционному уравнению Гиббса требуется знание зависимости dσ /dc от концентрации гептановой кислоты.

    Величина производной рассчитывается не менее чем по пяти равноотстоящим точкам. Для этого первым четырём точкам приписываем значения х=0, l, 2l, 3l. Производные в точках х=0 и х=l рассчитываем по формулам:

    (dy/dx)x=0 = (-21f(x0)+13f(x0+l)+17f(x0+2l)-9f(x0+3l))/20l, (2)

    (dy/dx)x=l = (-11f(x0)+3f(x0+l)+7f(x0+2l)-f(x0+3l))/20l, (3)

    где x0- первое значение абсциссы (х=0); l=1,25.

    Для определения производной в двух последних точках производим те же операции, начиная с последней (n-ой) точкой, которой приписываем значение х=0, а величины f(x0), f(x0+l), f(x0+2l), f(x0+3l) заменяем, соответственно, на f(xn),f(xn-l), f(xn-2l), f(xn-3l). Кроме того, знак производной изменяем на противоположный.

    Производные в остальных точках определяем следующим образом. Пяти выбранным точкам условно приписываем значения х, соответственно:

    -2l, -l, 0, l, 2l, а производную находим в точке х=0.

    Значение производной в точке х=0 рассчитываем по формуле:

    (dy/dx)x=0 = (-2f(x0-2l)-f(x0-l)+f(x0+l)+2f(x0+2l))/10l, (4)

    Алгоритм для расчета dσ /dc приведём ниже:

    1) 0,00249= -(-21*69,96 +13*67,01 +17*64,42-9*62,17)/(20*1,25*10^3)

    2) 0,00222=-(-11*69,96 +3*67,01 +7*64,42 +62,17)/(20*1,25*10^3)

    3) 0,00194=-(-2*69,96 -67,01 +62,17 +2*60,26)/(10*1,25*10^3)

    4) 0,00166=-(-2*67,01 -64,42 +60,26 +2*58,71)/(10*1,25*10^3)

    5) 0,00138=-(-2*64,42 -62,17 +58,71 +2*57,5)/(10*1,25*10^3)

    6) 0,00111=-(-2*62,17 -60,26 +57,5 +2*56,63)/(10*1,25*10^3)

    7) 0,00083=(-11*56,63 +3*57,5 +7*58,71 +60,26)/(20*1,25*10^3)

    8) 0,00055=(-21*56,63 +13*57,5 +17*58,71 -9*60,26)/(20*1,25*10^3)

    4) Полученные знания производной dσ /dc подставляем в уравнение Гиббса и рассчитываем величину адсорбции Г для разных концентраций с:

    (5)

    Рассчитываем Г следующим образом:

    1) 1,2799=-(-1,25 /(8,31*293))* 0,0025 *10^6

    2) 2,2752 =-(-2, 5 /(8,31*293))* 0,0022*10^6

    3) 2,9858=-(-3,75/(8,31*293))* 0,0019 *10^6

    4) 3,4117=-(-5/(8,31*293))*0,0017 *10^6

    5) 3,553 =-(-6,25/(8,31*293))*0,0014*10^6

    6) 3,4096 =-(-7, 5/(8,31*293))*0,0011*10^6

    7) 2,9815 =-(-8,75/(8,31*293))*0,0008*10^6

    8) 2,2687 =-(-10/(8,31*293))*0,0006*10^6

    Все рассчитанные значения сведены в таблицу 3.


    c,моль/мᶾ

    dσ/dc, Н/м

    Г*106,моль/м²

    0

    0

    0

    1,25

    0,002493

    1,280

    2,5

    0,002216

    2,275

    3,75

    0,001939

    2,986

    5

    0,001661

    3,412

    6,25

    0,001384

    3,553

    7,5

    0,001107

    3,410

    8,75

    0,00083

    2,981

    10

    0,000552

    2,269
    Таблица 3. Результаты расчёта адсорбции Г с примечанием численного дифференцирования

    По данным вычислений, приведенным в таблице 3, строим график изотермы адсорбции Г=f(c) (рис. 2).



    5) Для проверки применимости экспериментальных данных к теории Ленгмюра строим график в линейной форме уравнения Ленгмюра (с/Г от с). Данные для расчета с/Г берем из таблицы 3 и считаем следующим образом:

    1) 0,997=1,25/1,2799

    2) 1,099=2,5/2,2752

    3) 1,256=3,75/2,9858

    4) 1,466=5/3,4117

    5) 1,759=6,25/3,553

    6) 2,2=7, 5/3,4096

    7) 2,935=8,75/2,9815

    8) 4,408=10/2,2687

    Таблица 4. Результаты расчета адсорбции Г
    по уравнению Ленгмюра

    c, моль/мᶾ

    c/Г

    0

    0

    1,25

    0,977

    2,5

    1,099

    3,75

    1,256

    5

    1,466

    6,25

    1,759

    7,5

    2,200

    8,75

    2,935

    10

    4,408

    По полученным данным строим график зависимости с/Г от с:



    Видно, что через полученные точки можно провести прямую только при малых концентрациях, сама прямая тренда и ее уравнение даны на рис. 3.

    6)Далее проводим проверку применимости уравнения Шишковского

    ∆σ=σ0-σ=Bln(l+c/A) (7)

    Для расчета удельной адсорбции на основе экспериментальных данных по измерению поверхностного натяжения на границе раздела: раствор-воздух. Вначале оцениваются его применимость в области относительно высоких концентраций, когда С>>1 его можно упростить до уравнения прямой в координатах ∆ σ=f(lnc):

    ∆σ=σ0-σ=Bln(с/А) (8)

    Это позволяет определить константы уравнения А и В.

    Рассчитываем lnc:

    1) 0,921=LN(2,512)

    2) 1,151=LN(3,162)

    3) 1,382=LN(3,981)

    4) 1,612=LN(5,012)

    5) 1,842=LN(6,31)

    6) 2,042=LN(7,943)

    7) 2,303=LN(10)

    Рассчитаем ∆ σ по формуле ∆σ=σ0-σ:

    1) 5,75=72,75-67

    2) 7,45=72,75-65,3

    3) 9,15=72,75-63,6

    4) 10,95=72,75-61,8

    5) 12,65=72,75-60,1

    6) 14,35=72,75-58,4

    7) 16,05=72,75-56,7

    Все рассчитанные значения сведены в таблицу 4.

    Таблица 4. Данные для проверки применимости уравнения Шишковского

    c,моль/мᶾ

    lnc

    Δσ=σₒ-σ

    σ·10ᶾ , Н/м

    0

    0

    0

    72,75

    0,1

    -2,303

    0,35

    72,4

    0,398

    -0,921

    0,75

    72

    1

    0

    1,55

    71,2

    1,259

    0,230

    2,15

    70,6

    1,585

    0,461

    3,05

    69,7

    1,995

    0,691

    4,35

    68,4

    2,512

    0,921

    5,75

    67

    3,162

    1,151

    7,45

    65,3

    3,981

    1,382

    9,15

    63,6

    5,012

    1,612

    10,95

    61,8

    6,31

    1,842

    12,65

    60,1

    7,943

    2,072

    14,35

    58,4

    10

    2,303

    16,05

    56,7

    По приведенным данным строю график зависимости ∆σ от lnc (рис.4):



    Из полученных данных видно, что линейная зависимость справедлива только для относительно высоких концентраций адсорбата. Для этих точек с помощью Excel провожу прямую линию тренда, уравнение которой представлено на полях рис.4. Это свидетельствует, с одной стороны, о возможности применения уравнения Шишковского в области относительно высоких концентраций, а с другой, позволяет упростить его до вида

    σ=σ0-σ=Bln(Ас)=Blnс+BlnА

    и рассчитать константы А и В, после чего по ним проверить применимость данного уравнения для малых концентраций гептановой кислоты. Исходя из уравнения прямой у(σ)=7,476x+1,1428  нахожу, что В=7,476  и 

     А=EXP(1,1428/7,476)=1,16

    Проверим применимость полученных значений А и В для малых концентраций, например, при с=0,1. Поскольку ∆σ=σ0, то σ=σ0-∆σ. Отсюда σ=72,75 -(7,476*ln(1+0,1*EXP(1,1428/7,476)))=71,93.

     

    При σ  (экспериментальное) =72,4, ошибка в определении σ(расчётное) составляет:

    (72,4-71,93)*100/72,4=0,65%

    Как видно различие между экспериментальными и рассчитанными значениями σ незначительно (равно 0,65%). Это говорит о применимости уравнения Шишковского (7) во всем диапазоне концентраций и возможности расчета удельной адсорбции по уравнению Ленгмюра и константам Аи В (10).

    Г=(В/RT)*((c/A)/(1+c/A)) (10)

    Для этого надо произвести соответствующие расчеты, заполнить таблицу 5 и на основе ее данных построить рис.5.

    Таблица 5. Результаты расчета Гиббсовкской адсорции по уравнению 10

    c,моль/мᶾ

    Г*103,моль/м²

    0

    0

    0,1

    0,244

    0,398

    0,784

    1

    1,422

    1,259

    1,598

    1,585

    1,773

    1,995

    1,942

    2,512

    2,100

    3,162

    2,246

    3,981

    2,378

    5,012

    2,493

    6,31

    2,594

    7,943

    2,679

    10

    2,751



    Видно, что в этом случае график имеет теоретически правильную форму. Считается, что уравнение Шишковского позволяет определить поверхностную активность веществ более строго, чем уравнение Ленгмюра, из-за отсутствия концентрационных ограничений.

    7) В заключение проведем расчеты размеров молекулы адсорбата. На рис.3 по углу наклона линейной зависимости с/Г=f(с) равному 0,1117 рассчитываю предельную адсорбцию Г=1/0,1117=8,95*10-6 моль/м2. Если предположить, что при Г=Гобразуется насыщенный монослой из молекул гептановой кислоты, то площадь, занимаемая одной молекулой, составит:

    S0=(1/Г*Na)=1/(8,95*10-6*6,03*1023)=18,5*10-20 м2.

    Этот результат соответствует литературным данным, согласно которым площадь, занимаемая одной группой OH, равна (20-22)* 10-20 м2.

    Таким образом, в области малых концентраций к поверхностному слою изученного раствора в первом приближении применима теория мономолекулярной адсорбции.





    написать администратору сайта