Проект установки для дослідження теплопровідності нанофлюідів

Узнать стоимость написания работы

Зміст

Вступ..............................................3

Умовні позначення.....................................6

1. Аналіз літературних даних по теплопровідності нанофлюідів......9

1.1. Огляд і аналіз досліджень теплопровідності нанофлюідов з частинками Al2O3 11

1.2. Вплив концентрації наночастинок, їх розміру і форми..........12

1.3. Вплив температури..................................14

1.4. Вплив методів проведення експерименту...................15

1.5. Вплив конвективної складової..........................18

2. Опис експериментальної установки........................19

2.1. Конструкція вимірювальної комірки......................19

2.2. Опис систем термостатування і вимірювання температури......22

3.Методика проведення досліджень.........................25

4. Технологічні розрахунки елементів установки................32

4.1. Розрахунок вакуумної системи..........................32

4.2. Розрахунок витрати азоту..............................41


5. Охорона праці.......................................42

6. Організаційно – економічна частина.......................56

7. Висновки..........................................61

8. Список використаних джерел............................62

 

 

 

 


Вступ

В останні роки в Україні все більша увага приділяється питанням, пов'язаним з економією енергоресурсів. Питання енергоспоживання є для України надзвичайно складним клубком складних проблем. Тому питання енергозбереження слід вважати найважливішими економічними та екологічними проблемами на даному етапі розвитку.
Існує кілька принципових шляхів енергозбереження:
- Зменшення енергоємності продукції, що випускається на діючих підприємствах за рахунок модернізації технологій;
- Підвищення ефективності та зниження втрат при виробництві і транспортуванні електричної і теплової енергії;
- Впровадження альтернативних екологічно чистих джерел енергії;
- Поступове переведення економіки на так звані «інтелектуальні» технології (комп'ютерні, телекомунікаційні та ін.) Як значно менш енергоємні, більш рентабельні і екологічно чисті;
- Утилізація «енергосомістких» відходів;
- Інтенсифікація теплообмінних процесів в технологічних апаратах і підвищення їх енергетичної ефективності.
Різне енергетичне обладнання, таке як: сверхкритические котли, теплообмінні апарати і ядерні реактори, холодильне обладнання, теплові насоси та теплові труби, потребують інтенсифікації протікають в них теплообмінних процесів.
Аналіз показує, що задовольнити потреби підвищення енергетичної ефективності при тепловому регулювання перерахованих вище пристроїв можна трьома способами.
По-перше, простим збільшенням площі теплопередающей поверхні. Однак такий підхід призведе до збільшення матеріаломісткості, вартості обладнанні, а, отже, і зменшення конкурентоспроможності виробленого обладнання.

По-друге, модифікацією теплообмінних поверхонь апаратів (мікро структуровані поверхні, мікроканали, мікроканальних теплообмінники, мікроструі, нанесення на поверхні теплообмінних апаратів гідрофобних покриттів і т.д.). Однак в роботах [1, 2 EEJET] відзначається, що звичайне збільшення площі поверхні теплообміну і мікроканальних технології не завжди є альтернативним рішенням, наприклад, для теплового регулювання нового покоління напівпровідників і оптоелектронних пристроїв. Варто також враховувати збільшення вартості теплообмінного обладнання використовує зазначені технології, а також можливість засмічення профілю микроканалов і структурованих поверхонь, збільшення гідродинамічного опору теплообмінників.
Третій напрям інтенсифікації теплообміну в апаратах пов'язано з модифікацією вже застосовуваних теплоносіїв або робочих тіл для різного роду обладнання. В рамках цього напрямку пропонується використовувати нанофлюіди (колоїдні системи з твердою дисперсною фазою нано розмірів (до 100 нм) і рідкої дисперсної середовищем). Цей напрямок є досить перспективним. Аналіз літературних даних опублікованих в останні роки показує, що нанофлюіди (традиційні робочі речовини і теплоносії з добавкою частинок від одиниць до декількох десятків нанометрів) можуть істотно змінювати свої властивості і таким чином значно впливати на показники ефективності. Результати досліджень властивостей нанофлюідов, опубліковані в останні роки, показують можливість істотного цілеспрямованої зміни їх теплофізичних і теплообмінних характеристик в порівнянні з традиційними робочими речовинами холодильних систем і систем кондиціонування повітря. У ряді робіт отримано значне збільшення теплопровідності не тільки для нанофлюідов з наночастинками, що володіють високою теплопровідністю (вуглецеві нанотрубки, наночастинки високо теплопровідні металів, як наприклад, золото або мідь), а й для нанофлюідов з наночастинками оксидів металів. Для пояснення причин збільшення теплопровідності в нанорідин до теперішнього часу аналізується кілька основних механізмів: броунівський рух наночастинок, освіта високо теплопровідні рідкого шару (з товщиною кілька молекулярних шарів) на кордоні розділу рідина-наночастинок, балістичний перенос теплової енергії всередині індивідуальної наночастинки і між наночастинками, який відбувається при їх контакті, а також вплив кластеризації наночастинок. У ряді досліджень відзначено значне зростання коефіцієнтів тепловіддачі при різних режимах течії в каналах різної форми і при фазових переходах.
В даний час вивчаються перспективи застосування нанофлюідов для охолодження ядерних реакторів, в автомобільній промисловості для охолодження автоматичної коробки передач і поліпшення роботи системи охолодження двигуна і ін. Вивчаються також можливості поліпшення енергетичних характеристик холодильних систем за рахунок застосування нанофлюідов. Так в ряді робіт показано, що наночастинки оксидів металів можуть бути використані в якості добавки для підвищення розчинності мінерального масла з HFC-хладонами, при цьому поліпшуються повернення масла в компресор і енергетичні характеристики холодильної системи.
Як показують виконані різними авторами дослідження, навіть незначні домішки наночастинок різних металів, їх оксидів, вуглецеві нанотрубки, фулерени і т.п. сприяють значному збільшенню коефіцієнта тепловіддачі при різних процесах теплообміну.

У зв'язку з викладеним, метою дипломної роботи є проект установки для дослідження теплопровідності нанофлюідов.
Для досягнення зазначеної мети були поставлені наступні завдання:
- Розробка і проект експериментальної установки для дослідження теплопровідності методом лінійного джерела (нестаціонарний метод нагрітої нитки), що дозволяє вимірювати теплопровідність нанофлюідов;
- Розробити методику проведення досліджень теплопровідності нанофлюідов при різних температурах і концентраціях наночастинок.
Робота виконана відповідно до Закону України «Про енергозбереження» від 1.07.1994 р .; Постановами Кабінету Міністрів України № 148 від 5.02.1997 р «Про комплексну державну програму енергозбереження України», № 583 від 14.04.1999 р «Про Міжвідомчу комісію Із забезпечення виконан Рамкової Конвенції ООН про зміну клімату», планом науково-дослідних робот МОН Україна «Дослідження, розробка та впровадження нанофлюїдів у холодильну промисловість для Підвищення енергетічної ефектівності устаткування» (тема 12/04).

 


Умовні позначення

Q – тепловий потік;

 - коефіцієнт теплопровідності;

F - площа, поверхня плоского шару досліджуваного речовини;

 - час;

Т - зміна температури;

 - товщина шару досліджуваного речовини;

l - довжина циліндричного шару;

d1 и d2 - відповідно, внутрішній і зовнішній діаметри циліндричного або сферічеcкого шару;

- тиск;

сv - ізохорно теплоємність одиниці маси;

ρ - щільність;

- вектор щільності теплового потоку;

dQ1- кількість теплоти, введене в елементарний об'єм шляхом теплопровідності за час dτ;

dQ2 - кількість теплоти, яке за час dτ виділилося в елементарному обсязі dv за рахунок внутрішніх джерел;

dQ - зміна внутрішньої енергії або ентальпії речовини, що міститься в елементарному обсязі dv, за час dτ;

Ср - ізобарна теплоємність одиниці об'єму;

ср - ізобарна теплоємність одиниці маси;

- Коефіцієнт температуропровідності;

- Оператор Лапласа;

r - радіус-вектор, відстань в радіальному напрямку;

- полярний кут;

z - апліката;

– полярна відстань

- коефіцієнт форми для шару досліджуваного речовини;

R0 - опір термометра при 0 ° С;

А и В - термометричні коефіцієнти;

Т- кінцева температура;

Т0 - початкова температура;

f- коефіцієнт термодифузії;

γ = 0,5772 - константа Ейлера;

н - теплопровідність нанофлюіда;

б - теплопровідність базової рідини;

α - функція величини зазору;

Ra - критерій Релея;

ε – коефіцієнт конвекції;


1.Анализ літературних даних по теплопровідності нанофлюідів

Результати досліджень властивостей нанофлюідов, опубліковані в останні роки, показують можливість істотного цілеспрямованої зміни їх теплофізичних і теплообмінних характеристик в порівнянні з традиційними робочими речовинами холодильних систем і систем кондиціонування повітря. У ряді робіт отримано значне збільшення теплопровідності () не тільки для нанофлюідов з наночастинками, що володіють високою теплопровідністю (вуглецеві нанотрубки, наночастинки високо теплопровідні металів, як наприклад, золото або мідь, але і для нанофлюідов з наночастинками оксидів металів. Для пояснення причин збільшення теплопровідності в нанорідин до теперішнього часу аналізується кілька основних механізмів: броунівський рух наночастинок, освіта високо теплопровідні рідкого шару (з товщиною молекулярного рівня) на кордоні розділу рідина-наночастинок, балістичний перенос теплової енергії всередині індивідуальної наночастинки і між наночастинками, який відбувається при їх контакті, а також вплив кластеризації наночастинок. у ряді досліджень відзначено значне зростання коефіцієнтів тепловіддачі при різних режимах течії в каналах різної форми і при фазових переходах.
 Можливість поліпшення енергетичних характеристик холодильних систем за рахунок цього ефекту в даний час вивчається. Наприклад, в ряді робіт показано, що наночастинки оксидів металів можуть бути використані в якості добавки для підвищення розчинності мінерального масла з HFC-хладонами, при цьому поліпшується повернення масла в компресор і енергетичні характеристики холодильної системи (СОР збільшується до 26%).
Нанофлюіди представляють собою досить складні системи, на властивості яких істотно впливають багато факторів. З огляду на, що більшість публікацій за властивостями нанофлюідов присвячено

 

дослідженню теплопровідності, а отримані в них дані не узгоджуються між собою (див. рис. 2.1), в даній роботі проаналізовано основні фактори, що впливають на результати вимірювань теплопровідності. Як об'єкт експериментального дослідження теплопровідності обрана модельна система ізопропіловий спирт - наночастинки Al2O3.

Мал. 2.1. Вплив концентрації на наведену теплопровідність нанофлюідов нб (де н - теплопровідність нанофлюіда, б - теплопровідність базової рідини)
1 дані [1], 2 - [2], 3 - [3], 4 - [4], 5 - [5], 6 - [6], 7 - [7]

З основних факторів, що впливають на теплопровідність нанофлюідов, найбільш важливими є наступні:
- Спосіб отримання нанофлюіда і методика його підготовки до досліджень;
- Розмір і форма вихідних наночастинок, їх концентрація; здатність наночастинок утворювати кластери, розміри кластерів;
- Температура і тиск;
- Поверхнево-активні речовини (сурфактанти);
- Тип і властивості базових рідин;
- Методика проведення експерименту.

 

При дослідженні теплопровідності крім перерахованих вище факторів важливо враховувати адаптацію методики проведення експерименту до специфіки нанофлюідов, вплив конвективного складової і ефекту «пічної труби», вплив дифузійних ефектів, вплив випромінювання і переизлучения для нанофлюідов на основі базових рідин, напівпрозорих в ІК-області спектра і т . Д.

1.1. Огляд і аналіз досліджень теплопровідності нанофлюідов з частинками Al2O3

В опублікованих роботах досліджувався ефект впливу розміру і концентрації наночастинок з різними базовими рідинами на наведену теплопровідність при об'ємних концентраціях наночастинок, в основному, до 4-5% (див. Табл. 2.1).

Таблиця 2.1 - Огляд основних опублікованих досліджень теплопровідності нанофлюідов з частинками Al2O3

Перший

автор

Макс. концентр. частинок, %

Розмір наночастин,

нм

Макс. еффект

нб, %

Метод

вимірювань

 

Літер. джерело

Базова рідина- вода

Das

4

38

9

МТК

[8]

Putra

4

131

24

МПС

[9]

Masuda

4,4

13

33

НМНН

[2]

Lee

4,3

38

11

НМНН

[6]

Wang

5,5

28

16

МПС

[4]

Chon

1-4

11-150

30

НМНН

[10]

Li

10

36-47

29

МПС

[11]

Zhang

5

11

8

НМНКН

[12]

Timofeeva

5

11-40

13

НМНН

[13]

Xie

5

60

21

НМНН

[3]

Hwang

1

48

4

НМНН

[14]

Murshed 

1

80

12

НМНН

[15]

Oh

4

45

13

3-

[16]

Базова рідина - этиленгліколь

Lee

5

38

19

НМНН

[6]

Wang

8

28

41

МПС

[4]

Eastman

5

35

-

НМНН

[5]

Xie

5

26-300

30

НМНН

[3]

Beck

3

20

14

НМНН

[17]

Murshed 

0,5

150

10

НМНН

[15]

Oh

4

45

10

3-

[16]

Базова рідина – компресорне масло

Wang

7,4

28

30

МПС

[4]

Базова рідина – компресорне масло

Wang

7,1

28

20

МПС

[4]

Xie

5

60

29

НМНН

[3]

Базова рідина - гліцерин

Xie

5

60

38

НМНН

[3]

Базова рідина - трансформаторне масло

Choi

4

13

20

НМНН

[18]

 

*) МТК - метод температурних коливань, МПШ - метод плоского шару, НМНН - нестаціонарний метод нагрітої нитки, НМНКН - нестаціонарний метод короткої нагрітої нитки, 3- - 3 омега метод.

 

1.2. Вплив концентрації наночастинок, їх розміру і форми

 

Максимальна зростання наведеної теплопровідності склав 30-40% для різних базових рідин. Як правило, з ростом концентрації наночастинок наведена теплопровідність збільшувалася, проте в роботі Лі і Петерсона [11] ефект впливу наночастинок на наведену теплопровідність зменшувався при збільшенні концентрації від 6 до 10%. Результати роботи Тимофєєвої та співавт. [13] для цього ж діапазону концентрацій показують протилежний результат.

Аналіз впливу розміру і форми наночастинок на наведену теплопровідність з опублікованих робіт не дозволяє знайти не тільки закономірності, але і навіть тенденції такого впливу. Цей висновок ілюструє рис. 2.2, на якому представлені результати вимірювань з наночастками розміром від 13 до 80 нм (базова рідина - етиленгліколь).

Мал. 2.2. Вплив розмірів наночастинок на наведену теплопровідність нанофлюідов на основі етиленгліколю
1-15 нм [3], 2-35 нм [5], 3-40 нм [13], 4-80 нм [3]

 

Частина опитаних даних показує, що з ростом розміру наночастинок наведена теплопровідність при фіксованій концентрації зростає, інші дані призводять до протилежних висновків. Практично у всіх публікаціях вказана сферична форма наночастинок.

 


1.3. Вплив температури

Результати деяких досліджень впливу температури на наведену теплопровідність наведені на мал. 2.3.

 

Мал. 2.3. Вплив температури на наведену теплопровідність нанофлюідів
1, 2, 3 - дані [2] при х = 1,25%, 2,75% і 4,25%; 4, 5 - дані [8]
при х = 1% і 4%; 6 - дані [11] при х = 2%.

 

Як видно з цього мал., Результати цих вимірювань прямо протилежні - частина з них показує зменшення ефекту впливу наночастинок на теплопровідність базової рідини з ростом температури, інша частина - збільшення цього ефекту. Можливо, це пов'язано зі збільшенням рухливості молекул і розпадом кластерів (конгломератів) наночастинок. До такого ж висновку прийшли Мёрст і ін. [15], вивчаючи вплив добавки сурфактанта на зростання наведеної теплопровідності при різних температурах.

1.4. Вплив методів проведення експерименту.

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

При вимірі теплопровідності нанофлюідов використовувалися як стаціонарні (метод плоского горизонтального шару), так і нестаціонарні (метод нагрітої нитки, метод «короткій» нагрітої нитки, метод нагрітої мікрополоскі, метод температурних коливань, так званий «3ω» метод, метод фототермічна радіометрії, метод теплового компаратора), а також неінвазивні оптичні методи вимірювань. В рамках цієї роботи немає можливості проаналізувати всі зазначені методи вимірювань теплопровідності нанофлюідов з точки зору їх точності, однак слід відзначити їх деякі характерні особливості.

Найбільший обсяг експериментальних даних до теперішнього часу отримано за допомогою нестаціонарного методу нагрітої нитки. Цей метод використовувався багатьма авторами для дослідження теплопровідності різних однофазних рідин. Таке широке застосування нестаціонарного методу нагрітої нитки пояснюється наявністю ряду переваг перед іншими методами досліджень: відносна простота конструкції вимірювального вузла; метод забезпечує достатню точність для однофазних середовищ і вельми експресивний (процес вимірювання займає не більше кількох секунд). Теорія цього методу заснована на одновимірному вирішенні рівняння Фур'є для лінійного джерела теплоти, який являє собою нескінченно довгу і тонку нитку, що розсіює тепло в нескінченну середу. Разом з тим, вплив граничних умов, конвективного, дифузійної та радіаційної складових в дослідах з нестаціонарним методом нагрітої нитки для нанофлюідов не вивчено.
Ряд удосконалень нестаціонарних методів був зроблений в більш пізніх роботах. Наприклад, Xie, Zhang і ін. [3,12] використовували нестаціонарний

метод «короткій» нагрітої нитки, теорія якого заснована на двомірному вирішенні рівняння Фур'є для лінійного джерела теплоти з тим же відношенням довжини до діаметра, що і зразок нанофлюіда. На думку авторів, така методика дозволяє враховувати крайові ефекти.

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

Спробу уникнути агрегації наночастинок і їх осадження зробив Mintsa [7], який помістив пристрій для перемішування нанофлюіда поруч з вимірювальним пристроєм, що реалізує нестаціонарний метод нагрітої нитки. У деяких роботах нагрівач і датчик температури були відокремлені один від одного (в якості нагрівача використовувався лазерний промінь, а в якості датчика температури - термометр опору), однак це не призвело до збільшення точності отриманих результатів.

В останні роки активно розроблялися оптичні методи вимірювання
теплопровідності [19-26], в яких датчик температури і нагрівач відокремлені один від одного, забезпечуючи потенційно більш точні
дані. Вимірювання в такому пристрої відбуваються протягом декількох мікросекунд, тобто набагато швидше, ніж при вимірах за допомогою «класичного» варіанту нестаціонарного методу нагрітої нитки. На думку авторів цих робіт, вплив природної конвекції в таких умовах дуже малий.

 

Порівняння результатів вимірювань теплопровідності нанофлюідов на основі наночастинок Al2O3 різними методами наведено на мал. 2.4.

 

Мал. 2.4. Вплив методів вимірювання на наведену теплопровідність
1-метод плоского шару [4], 2-нестаціонарній метод нагрітої нитки [6], 3-метод температурних коливань [8], 4-3ω метод [16].

Як видно з цього графіка, розбіжності даних різних авторів значно перевершують заявлену ними похибка. Причини таких розбіжностей неочевидні, проте зрозуміло, що найбільшої довіри заслуговують результати, отримані при варіюванні основних факторів, що впливають.

 


1.5. Вплив конвективної складової

При конструюванні вимірювального вузла приладу для визначення теплопровідності рідин його геометричні розміри, конструкція і температурної режим повинні виключити або зменшити ймовірність виникнення конвективного теплообміну. Умови, при яких виникає конвективний теплообмін, вивчені різними дослідниками для однофазних рідин. Зазвичай вплив конвекції враховується введенням відповідної поправки (коефіцієнта конвекції ε=λэф/λ). Для обґрунтованого введення такої поправки необхідно вивчення умов виникнення конвективного переносу тепла в конкретних умовах експерименту.

До перших робіт з вивчення конвекції слід віднести роботи Релея, який досліджував умови виникнення конвективного переносу тепла в шарі рідини між плоскими паралельними пластинами. Їм було встановлено теоретично, що поріг настання конвекції визначається критичною величиною критерію Релея Raкр. Висновком наступних робіт є залежність коефіцієнта конвекції від геометрії шару досліджуваного речовини і його розташування.

У подальших роботах було показано, що для вузьких циліндричних зазорів коефіцієнт конвекції визначається не тільки значеннями Ra, але і величиною вимірювального зазору.

Слід зазначити, що всі дослідження в цій області проводилися для однофазних рідин. У нанофлюідах (двофазних системах) процес виникнення конвективного теплообміну набагато складніше і практично не вивчений. Тому найбільш надійним способом виключити вплив конвекції на результати вимірювання теплопровідності є досліди на декількох значно відрізняються між собою перепадах температур в стаціонарних методах і різних теплових навантаженнях в нестаціонарних методах вимірювань.


2.Опис Експериментальної установки

2.1. Конструкція вимірювальної комірки

Для експериментального дослідження теплопровідності нанофлюідов була розроблена вимірювальна комірка (мал. 2.5).

Мал. 2.5 Конструкція вимірювального пристрою реалізує нестаціонарний метод нагрітої нитки.
1, 10 - струмові дроти; 2 - платинова нитка; 3, 7 - потенційні дроти; 4 пружина; 5 - вузол натягнення нитки; 6, 9, 11 - скляні капіляри; 8 – каркас

Вимірювальна комірка (Мал. 2.5) реалізує варіант нестаціонарного методу нагрітої нитки - метод лінійного джерела. Основним елементом вимірювальної комірки є термометр опору, виготовлений з платинового дроту 2 діаметром 50 мкм і довжиною 61 мм. Робоча ділянка нитки припаяний до струмопровідних проводах 1, 10 діаметром 0,1 мм, верхній з яких запаяний в центрир скляному капілярі 6. Нижній провід, що підводить струм через ізолюючий скляний капіляр 9, приєднаний до вольфрамової пружині 4, яка в свою чергу змонтована в вузлі натягу 5. Натяг пружини здійснювалося вантажем масою 2 г і фіксувалося стопорним гвинтом. Така конструкція дозволяє забезпечити постійне в часі і незалежне від параметрів стану натяг робочої нитки і створити умови для стабільної роботи нитки як термометра опору. Зазначені вище елементи вимірювального центру змонтовані всередині каркаса 8, виготовленого з латуні діаметром 12 мм.

Діаметр дроту визначали двома незалежними методами - виміром діаметра в двох взаємно перпендикулярних площинах через кожні 10 мм по довжині дроту і зважуванням 1 м дроту на аналітичних вагах. Для стабільності показань через напружені режимів роботи осередку платинова дріт попередньо двічі давав жару в муфельній печі при температурі 1100 К протягом 10 годин.

Нижній кінець платинової нитки кріпиться до вольфрамової пружині, що знаходиться в нижньому вузлі комірки. Пружина ізольована від корпусу в гвинті «натягу» за допомогою тонкостінних скляних капілярів. Натяг нитки здійснювалося за допомогою підвішеного вантажу масою 0,02 кг протягом 24 годин з подальшим зняттям вантажу і вивінчіваніем до упору гвинта, який фіксувався контргайкой.

Як струмових і потенційних висновків використовували платинову дріт діаметром 0,1 і 0,02 мм, відповідно. Потенційні висновки, що визначають робочу ділянку нитки, виготовлені з платинового дроту діаметром 0,02 мм. Подпайкі здійснювалася сріблом за допомогою спеціально виготовленої мікрогорелкі. У верхній частині осередку для токового підведення нитки використовується сама нитка. Така конструкція дозволяє значно зменшити торцеві ефекти і спрощує висновок струмових і потенційних проводів в місці їх виведення з робочої зони. Для зниження відтоків тепла по провідникам виведення їх здійснювався по ізотермічної поверхні. Для висновків використовуються свердління діаметром 2 мм, виконані в двох взаємно перпендикулярних площинах в нижньому і верхньому кінцях корпусу вимірювального вузла. Для запобігання вертикально-кільцевої конвекції осередок мала «закриті торці». Це здійснювалося підбором відповідних розмірів скляних капілярів. Висновок платинових струмових і потенційних проводів здійснювався через спеціальні пропили на торцях скляного капіляра.

Для виведення і ізоляції струмових і потенційних проводів поза робочим зазору використовуються скляні тонкостінні капіляри, укладені в спеціально профрезерований канавках по корпусу вимірювальної комірки і фторопластового ізолятора верхнього вузла. Кріплення скляних капілярів до корпусу здійснюється фторопластовой плівкою в проточках корпусу. Величина зазору між скляним капіляром і каркасом осередку становила 0,3 мм, між корпусом вимірювальної комірки і стінкою теплообмінника - 0,2 мм. Таким чином, для усунення можливих конвективних потоків все зазори в корпусі осередки і автоклава зведені до мінімуму.

Центрування вимірювального капіляра щодо нитки вироблялося через свердління в каркасі в двох взаємно перпендикулярних площинах під мікроскопом типу УИМ-21 за допомогою настановних штифтів. Положення «центру» визначалося візуально і за відносним мінімуму ємності системи нитка-капіляр. Висока чутливість вимірювача ємності дозволила виявити і регулювати ексцентриситет до величини 0,001 мм. Щоб уникнути «внутрішнього» ексцентриситету, викликаного можливим викривленням нитки в робочій зоні, вимірювальний ділянку нитки вибирався з умови мінімуму ємності при протягуванні нитки через капіляр. Жорсткі вимоги до центрівці пояснюються тим, що точність визначення теплопровідності істотно залежить від величини ексцентриситету при малих вимірювальних зазорах. Крім того, при дослідженнях теплопровідності нанофлюідов необхідна максимальна симетричність температурного поля.

 


2.2. Опис систем термостатування і вимірювання температури

Загальна схема всієї експериментальної установки представлена ​​на мал.3. Перед заправкою досліджуваної речовини з обсягу заправної камери відкачується повітря за допомогою вакуумної системи. З балона 5 досліджувана рідина заправляється в досліджувану осередок 19 через заправний вентиль 27. Об'єм комірки обмежений притискним фланцем 26 зверху і запірної пробкою 18 знизу. Струмові і потенціометричні дроти, йдуть до вимірювальної комірки через електровводи 29. Також передбачений перетворювач тиску 28 з'єднана з порожниною вимірювальної камери по паровій фазі. Мідний блок 19 служить для вирівнювання поля температур в осередку. Установка призначена для вимірювання теплопровідності при температурах від -60 ° С до 80 ° С. Для забезпечення роботи на низьких температурах (від -60 ° С до -20 ° С) через біфілярного спіральний канал мідного теплообмінника 25 прокачиваются пари азоту, що надходять з системи подачі рідкого азоту. Для температур -20 - +20 - прокачуємо хладоноситель з допоміжного термостата, який охолоджується холодильною машиною. Для температур +20 - 80 - прокачуємо воду з допоміжного термостата. Допоміжні термостати на мал. 3 не вказані.

Для зниження теплопритоков з навколишнього середовища мідний блок поміщений в вакуумну камеру 14, яка відкачується вакуумною системою. Висновок термостата на необхідний температурний рівень здійснюється за допомогою фонового 23 і регулює 22 нагрівачів, які живляться від стабілізованих джерел живлення СП. Датчиком системи регулювання є платиновий термометр опору (R0 = 100 Ом) 20.

Мал. 3. Принципова схема низькотемпературного стенду по дослідженню теплопровідності рідин

1 - сільфонний регулятор тиску; 2 - манометр; 3 - посудина Дьюара; 4 - нагрівач; 5 - заправний балончик; 6, 13 - термопарниє вакуумметри; 7 - іонізаційний вакуумметр; 8, 10 - вакуумні пастки; 9 - дифузний вакуумний насос; 11 - натекатель; 12 - форвакуумний насос; 14 - вакуумна камера; 15 - диференціальна термопара; 16 - платиновий термометр опору 100П; 17 - вимірювальна осередок; 18 - запірна пробка; 19 - мідний блок; 20 - датчик системи регулювання температури; 21, 24 - торцеві нагрівачі; 22 - регулюючий нагрівач; 23 - фоновий нагрівач; 25 - теплообмінник; 26 - фланець; 27 – заправний

вентиль; 28 - перетворювач тиску; 29 - електровводи; 30 - двопозиційний двоканальний перемикач; 31 - трьохпозиційний двоканальний перемикач; 32 - перемикач напряму струму; ІП - джерела живлення; У - підсилювач потужності; ДРТ - блок регулювання температури; МО-62 - міст опору; МСР - магазин опорів; ОКС - зразкова котушка опору; Picotest 1, Picotest 2 - мультиметри.

В термостаті може виникати небажаний градієнт температури по висоті, що призведе до конвективних струмів всередині вимірювальної камери. Для запобігання цьому явищу на торцях термостата встановлена диференційна термопара 15, по сигналу якої включається один з додаткових торцевих нагрівачів 21 або 24. Ці нагрівачі підключені паралельно до джерела живлення фонового нагрівача за допомогою магазину опорів МСР і перемикача 30.

Вимірювання параметрів виробляються в наступному порядку. Після досягнення необхідного температурного рівня за допомогою джерела живлення ІП протягом 10 секунд подається напруга на платинову нитка вимірювального пристрою. Мультиметром Picotest 1 реєструється падіння напруги на зразковою котушці опору ОКС послідовно включеної в ланцюг нитки. Мультиметром Picotest 2 реєструється падіння напруги на самій нитки. У режимі визначення теплопровідності прилади синхронізовані між собою і знімають показання з частотою 50 Гц. У проміжку між вимірами для переходу на інший температурний рівень по приладу Picotest 2 можна стежити за показаннями термометра 16 і термопари 15 використовуючи перемикач 31. Для усунення шумових ефектом пов'язаних з електромагнітними наведеннями в ланцюзі вимірювального пристрою передбачений перемикач напрямку струму 32. До розгляду беруться середні значення параметрів виміряних при прямому і зворотному напрямках струму.

 


3. Методика проведення досліджень

Досліди проводилися на цьому вимірювальному пристрої в нестаціонарному режимі з урахуванням ряду обмежень, пов'язаних з невиконанням граничних умов при вирішенні рівняння Фур'є для лінійного джерела теплоти (діаметр нитки-нагрівача склав 0,1 мм, товщина шару середовища - 0,55 мм).

Методика вимірювання теплопровідності нестаціонарним методом нагрітої нитки заснована на рішенні рівняння Фур'є для лінійного джерела тепла, що представляє собою нескінченно довгу і тонку нитку, що розсіює тепло в нескінченну середу:

 

(2.2)

де, ΔT=T–T0, T - кінцева температура, T0 - початкова температура, r - відстань в радіальному напрямку, τ - час, q - тепловий потік, D - коефіцієнт термодифузії.
Граничні умови можуть бути записані у вигляді:

 

при τ = 0 и r = 0 (2.3)

при τ ≥ 0 и r = ∞ (2.4)

З урахуванням цих умов, зміна температури в радіальному напрямку від джерела тепла може бути представлена у вигляді рівняння

 

 

(2.5)

 


Інтегрування рівняння (7) приводить до формули

 

 

(2.6)

 

а для будь-якого фіксованого в радіальному напрямку відстані зміна температури може бути представлено у вигляді

 

 

(2.7)

звідки випливає, що графік зміни температури як функція логарифма часу являє собою пряму лінію, а нахил цього графіка пропорційний коефіцієнту теплопровідності. З рівняння (2.7) випливає, що теплопровідність може бути розрахована за формулою

 

(2.8)

Таким чином, для розрахунку теплопровідності з даних експерименту необхідно вимірювати потужність теплоти, яка виділяється ниткою, і зміна температури у вре¬мені. На рис. 2.6 і 2.7 показані приклади отриманих в тарувальних експериментах даних (ці вимірювання проводилися з кроком 0,02 сек).


Мал. 2.6. Залежність температури нитки від логарифма часу в перші
7 сек досвіду при різних гріючих токах: 1-400 мА, 2-350 мА, 3-300 мА

Мал. 2.7. Залежність температури нитки від логарифма часу в період часу 0,14-1,0 сек при різних гріючих токах: 1-400 мА, 2-350 мА, 3-300 мА

Як видно з цих графіків, лінійна залежність температури від логарифма часу дотримується в діапазоні -2 ≤ lnτ ≤ 0 при різних теплових навантаженнях. При ln τ ≤ -2і при ln τ ≥ 0 невиконання граничних умов приводить до нелінійності залежності температури від lnτ і, відповідно, до помилок в розрахунку λ.

Приклад результатів вимірювань теплопровідності ізопропілового спирту з наночастинками Al2O3 (об'ємна концентрація наночастинок х = 2,64%) при початковій температурі 8 ° С на вимірювальному пристрої, що реалізує «класичний» варіант нестаціонарного методу нагрітої нитки показаний на рис. 2.8 і 2.9 у вигляді залежності температури нитки від часу нагріву і залежності температури нитки від логарифма часу з "кроком" за часом 0,02 сек при різних значеннях гріє струму. Як видно з цього графіка, практично у всьому діапазоні вимірюваних параметрів кут нахилу лінії, що виражає залежність температури від ln τ, змінюється незначно.

Мал. 2.8 Залежність температури нитки від часу нагріву при різних гріючих токах
1 - 200 мА, 2 - 300 мА, 3 - 400 мА.

Мал. 2.9. Залежність температури нитки від логарифма часу нагріву при різних гріють токах
1-200 мА, 2-300 мА, 3-400 мА

Аналіз цих даних для різних часових діапазонів, представлений в табл. 2.2, показує, що похибка результатів вимірювань теплопровідності не перевищує 1,2%.
Таблиця 2.2. - Результати розрахунку ∂∆T/∂(ln τ), λ і відхилень від середнього (λ-λср)/λср ∙ 100,% для різних часових діапазонів

 

τ, сек

ln τ

∂∆T/∂(ln τ)

λ, Вт/(м∙К)

(λ-λср)/λср∙100, %

0,02 - 0,05

-4…-3

1,251

0,1539

-1,11

0,05 – 0,14

-3…-2

1,232

0,1557

0,04

0,14 – 0,36

-2…-1

1,242

0,1544

-0,79

0,36 – 1

-1…0

1,222

0,1570

0,88

1 – 2,74

0...1

1,234

0,1554

-0,15

2,74 – 7,38

1…2

1,219

0,1574

1,14

3.1. Результати контрольних дослідів

Контрольні досліди були проведені за стандартною рідини - толуолу в діапазоні температур від 20-80оС. Толуол був обраний в якості градуйованної речовини, оскільки для цієї рідини коефіцієнт теплопровідності добре вивчений різними методами (як стаціонарними, так і нестаціонарними) і його теплопровідність визначена з урахуванням напівпрозорої в інфрачервоній області спектра.

Фрагмент результатів вимірювань теплопровідності толуолу представлений в табл. 3.2 і на мал. 3.8.

Таблиця 3.2. Експериментальні значення теплопровідності толуолу

t, 0C

λ, Вт/(мК)

12,63

0,1351

14,35

0,1353

14,98

0,1350

15,53

0,1343

17,54

0,1342

18,56

0,1336

43,04

0,1281

43,05

0,1282

43,06

0,1282

45,12

0,1278

47,13

0,1268

49,13

0,1273

72,90

0,1213

73,06

0,1212

74,07

0,1210

75,91

0,1208

79,91

0,1203

(Табл.3.2.)

(Мал. 3.8)
Мал. 3.8 Результати контрольних дослідів по теплопровідності толуолу

 

Аналіз цих результатів показує, що відхилення, отриманих нами даних по теплопровідності толуолу відрізняються від стандартних довідкових даних в межах похибок експерименту.


4. Технологічні розрахунки елементів установки

4.1. Розрахунок Вакуумної системи

Схема вакуумної системи для розрахунку представлена на рис.4

Рис.4. Схема вакуумної системи.

1, 7, 8, 9, 10, 11 - вентилі вакуумної системи; 4 - вакуумна камера;

5 термостат; 6 - іонізаційний манометр; 12 - трубопровід; 13 - насос для

отримання низького вакууму; 14, 18 - термопарний манометр; 15, 17 –

нітрогену пастка 16 - насос для отримання середнього вакууму.

Основні розміри вакуумної камери наведені нижче.

lтр = 0,43 м - Довжина трубопроводу,

Dтр = 0,32 м - Діаметр трубопроводу,

(2.1)

FTP = πdl = 0,025м2– площа трубопроводу, (2.2)

Розміри вакуумної камери

- обсяг вакуумної камери (2.3)

- площа вакуумної камери, (2.4)

Pраб = 10-3 Па,

Матеріал вакуумної камери - сталь, трубопровід виконаний з нержавіючої сталі Х18Н10Т.

Розрахуємо величину газового потоку

Q= Qгв + Qн , (2.5)

Qгв– потік від газовиділень стінок,

Qн– потік від протікань в систему,

Вк: Qгв1 = qFв = 2·10-4·0,0045=9·10-6 м3Па/с,

qcu = 2·10-4 м3Па/ м2с, [45],

Тр: Qгв2 = qFтр = 4·10-5·0,0025=1·10-6 м3Па/с,

qтр = 4·10-5 м3Па/ м2с, [45],

Qгв = 9·10-6 ·1·10-6 = 10-5 м3Па/с,


Відсутність течі контролюється течошукачем, які мають мінімальний реєстрований потік 10-7 м3Па/с. Тоді величина потоку натекания в систему расчітвается за такою формулою:

 

, (2.6)

кв - ймовірність існування течі, прийнята 0,2;

N = 11 - число з'єднань в системі;

QTИ = 10-7 м3Па/с - мінімальний потік реєстрований течошукачем,

m – 1– число одночасних перевіряються з'єднань,

м3Па/с, (2.7)

Повна величина газового потоку

Q = 10-5 + 2,2·10-7 = 1,0022·10-5 м3Па/с, (2.8)

 

а) Вибираємо насос для роботи в області високого вакууму.

Вибираємо серію паромасляних насосів типу Н, що мають граничний тиск Pпр1 = 3·10-4.. 4·10-5Па. З діапазоном швидкостей дії

SH = 1.5·10-5 .. 30 м3/с, робочий тиск паромасляного Pраб1= 10-3 Па.
Ефективна швидкість відкачування

м3/с, (2.9)

 

Знаходимо значення коефіцієнта використання Ки1= 0,32.

Номінальна швидкість дії пароструйного насоса:

 

м3/с, (2.10)

Найближчий по швидкості дії пароструминний насос Н - 160/170 має наступні характеристики:

Діапазон робочих тисків – 10-4 Па;

Швидкість відкачки в робочому діапазоні тисків – 0,64 м3/с;

Граничне залишковий тиск – 7·10-5 Па;

Найбільший випускний тиск – 33 Па;

б) Вибір насоса для роботи в області низького вакууму.

Вибираємо обертальний насос, що мають граничний тиск

Pпр2 = 4·10-1.. 10-2Па. З діапазоном швидкостей дії SH = 10-4 .. 1,5·10-1 м3/с,

Робочий тиск механічного насоса вибираємо по максимальному випускного тиску паромасляного насоса з коефіцієнтом запасу φ = 2, Pраб2= 17 Па.

Ефективна швидкість відкачування:

 

м3/с, (2.11)

Знаходимо значення коефіцієнта використання Ки2= 0,64.

Номінальна швидкість дії пароструйного насоса:

 

м3/с, (2.12)

Найближчий по швидкості дії пароструминний насос ВН – 01 має наступні характеристики:

Номінальна швидкість дії – 10-4 м3/с.

Діаметр вхідного патрубка – 8мм.

граничний тиск – 4 Па.

Граничний тиск насоса з пасткою – 4·10-1 Па.

Максимальна випускний тиск – 105 Па.

Визначимо кордону режимів течії повітря в трубопроводі при Т=293 0К.

Кордон між вязкостним і молекулярно-вязкостним режимами визначаємо за формулою:

 

мм рт. Ст (2.13)

Кордон між і молекулярно-вязкостним і молекулярною режимами визначаємо за формулою:

 

мм рт. ст. (2.14)

Відкачування від 650 мм рт.ст. до 1,25 мм рт.ст.

Оскільки провідність трубопроводу більше швидкості відкачки насоса, можна вважати, що в діапазоні тисків від 1,25 мм рт. ст. до 0,0375мм рт.ст., швидкість відкачування дорівнює ефективної швидкості відкачки вакуумної камери.

Отже час відкачування буде розраховуватися за формулою [43]:

 

с (2.15)

 

де V - повний обсяг вакуумованих порожнин

 

V=Vвк + Vтр = Vвк+ π·Rтр2·lтр=0.755·10-3(2.16)

P1/ и P1// - тиску на кордонах режиму (P1// =650 мм рт.ст, P1/= 1,25мм рт.ст).

Відкачування від 1,25 мм рт.ст. до 0,0375 мм рт.ст.

Попередньо розрахувати пропускну здатність трубопроводу:

 

(2.17)

Де Uтр ,Uпов ,Uвент – пропускна здатність трубопроводу, поворотів і вентилів відповідно.

(2.18) л/с

де P1, P2 - тиск кордонів режиму течії [49].

Пропускна здатність вентилів розраховувалася за формулою

(2.19)

 

л/с

 

л/с(20)

 

л/с(2.20)

Пропускна здатність поворотів розраховувалося за формулою

 

(2.21)

 

где lрас=l+1.33·d·=·2.00+1.33·0.004·8=2.042 м

P1, P2 - тиск кордонів режиму течії;

d- діаметр трубопроводу;

i- кількість поворотів.

 

л/с (2.22)

Ефективність відкачування в молекулярно - вязкостному режимі розраховувалася за формулою:

л/с (2.23)

 

Так як , то час відкачування в молекулярно-вязкостной режимі знайдемо за формулою:

с(2.24)

Відкачування від 0,0375 мм рт.ст. до 0,01 мм рт.ст:

Пропускна здатність трубопроводу при молекулярному режимі:

л/с(2.25)

Пропускна здатність поворотів при молекулярному режимі:

л/с(28)

Де л/с

Общая пропускная способность при молекулярном режиме:

л/с (2.26)

Ефективність відкачування в молекулярному режимі розраховувалася за формулою:

 

л/с (2.27)

час відкачування в молекулярному режимі знайдемо за формулою:

 

с

Загальний час відкачування для трьох режимів:

 

с, (2.29)

 

, (2.30)

м3/с,

(36)

м3/с, (2.31)

Час відкачування

Отриманий час задовольняє вимогам реалізації методики проведення досвіду.

4.2.Розрахунок витрати азоту

В ході досліджень експериментально було встановлено, що азот випаровувався повністю протягом 50 хв = 3000 с. Так як заливаємий обсяг азоту дорівнює:

То витрата азоту дорівнює:


5. Охорона праці

Охорона праці - це система правових соціально-економічних, організаційно-технічних, санітарно-гігієнічних і лікувально-профілактичних заходів та засобів, спрямованих на збереження здоров'я і працездатності людини в процесі праці. В даний час охорона праці підтримується кодексом законів про працю України.

Небезпечні і шкідливі речовини і фактори, що застосовуються в лабораторії. Їх вплив на людину і методи захисту

Основними шкідливими факторами на проектованому об'єкті являють-ся: електричний струм і R600a (при витоках в робочу зону).

Проходячи через організм, електричний струм надає термічне, електролітичне і біологічне дії. Це різноманіття дій електричного струму призводить до електротравма і електричним ударам.

Розрахунковим опором тіла людини є для непошкоджений-ної шкіри 1000 Ом. Величина струму, що протікає через організм, є головним фактором, від якого залежить результат поразки (чим більше струм, тим небезпечніше його дія).

Основними заходами захисту від ураження струмом є: забезпе-ня недоступності струмоведучих частин, що знаходяться під напругою, для випадкового дотику; електричне поділ мережі; усунення небезпеки поразки у разі напруги на корпусах, кожухах та інших частинах електроустаткування, що досягається застосуванням малих напруг, використанням подвійної ізоляції, вирівнюванням потенціалу, захисним заземленням, занулением, відключенням.

При атмосферному тиску холодоагент R600a - безбарвний газ, легший за повітря, практично без запаху.

Найбільш небезпечною властивістю R600a є його вибухонебезпечність.

Крім того, газоподібний R600a подразнює слизові оболонки очей і верхніх дихальних шляхів, на спітнілі ділянки шкіри. При важких отруєннях R600a спостерігаються утруднене дихання, сильний кашель, задуха, спазми голосової щілини.

Гранично допустима концентрація (ГДК) R600a в повітрі робочої зони складає 20 мг/м3.

Проникнення R600a в повітря робочої зони можливо при розгерметизації лабораторної експериментальної установки.

У лабораторії обов'язкова установка сигналізаторів витоку і аварійної концентрації парів R600a в повітрі приміщення. Сигналізатори витоку повинні давати попереджувальний сигнал і включати припливно-витяжну вентиляцію при концентрації R600a вище 0.5 - 1 мг / л. При досягненні концентрації 1.5 мг / л сигналізатори аварійної концентрації повинні вимикати електроживлення всієї лабораторної експериментальної установки і одночасно включати аварійну вентиляцію і світлозвукову сигналізацію.

 

Захист від ураження електричним струмом

 

Для захисту від впливу електричного струму рекомендується ис-користувати захисне заземлення. Захисне заземлення - це навмисне електричне з'єднання з землею або її еквівалентом металевих неструмоведучих частин електрообладнання. В електричних мережах напругою до 1000 В опір заземлювальної системи не повинно перевищувати 4 Ом. Це в 250 разів менше, ніж розрахунковий опір тіла людини. При випадковому зіткненні людини до такого заземленого корпусу струм, який піде через нього буде відповідно менше.

У приміщенні лабораторії заземлення виконують наступним обра-зом: всередині приміщення, по периметру, на відстані 20-30 см від підлоги прокладають металлопроводнік, який з'єднаний з землею. У землі знаходяться електроди довжиною від 2 до 3 м. Металеві корпуси обладнання, що знаходяться в лабораторії з'єднані провідниками з заземлюючим контуром за допомогою зварювання. З корпусом провідник з'єднується різьбових з'єднань з пружиною і шайбою для запобігання самовідгвинчування.

Класифікація приміщення за ступенем небезпеки ураження електричним струмом

 

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

 

Особливо небезпечні приміщення за ступенем небезпеки ураження електричним-струмом характеризуються підвищеною вогкістю, з відносною вологістю повітря, близькою до 100%, наявністю хімічно активної середовища.

 

розрахунок заземлення

 

Розрахунок заземлюючого пристрою для електрообладнання напругою до 1000 В при грунті - глина. Як штучного заземлювача використовуються сталеві труби:

– діаметр и довжиною ,

– відстань між трубами ,

– ширина полоси ,

– опір грунту ,

– кліматичний коефіцієнт .

 

Для проведення заземлення в лабораторії використовувалися матеріали з наступними геометричними і іншими параметрами: м

м, м, м, , .

Визначаємо розрахункове значення питомої опору грунту:

(5.1)

Ом∙м

Вибираємо систему розподілу вертикальних заземлювачів. Для малих беремо в ряд.

Задаємо в межах . Нехай м.

 

(6.2) м (5.2)

Опір одного вертикального заземлювача становить:

 

(5.3)

 

Ом.

Визначаємо кількість одиночних заземлювачів:

 

,(5.4)

 

В електричних мережах напругою до 1000 В відповідно до вимог ПУЕ опір заземлювальної системи не повинно перевищувати 4 Ом. Це в 250 разів менше, ніж розрахунковий опір тіла людини. При випадковому дотику людини до такого заземленого корпусу струм, який піде через нього буде відповідно менше. Виходячи з того, що розраховане опір системи складає 2,37 Ом, можна зробити висновок про правильність розрахунку, так як розрахункове значення задовольняє умові ≤ 4 Ом.

 

шт.

Приймаємо =4; вибираємо =4 шт.

Визначаємо коефіцієнт використання вертикальних стрижневих заземлювачів.


(5.5)

При значення коефіцієнту Ом.

Розраховуємо довжину горизонтальної електрода:

(5.6)

м.

Визначаємо опір заземлювача:

,(5.7)

Де - коефіцієнт використання горизонтального заземлювача. Визначаємо, що = = 0,89[ ]тоді

 

Ом.

Розрахунок загального опору системи:

 

(5.8)

Ом.

 

Згідно з вимогами ПУЕ - в електричних мережах напругою до 1000 В опір заземлювальної системи не повинно перевищувати 4 Ом. Виходячи з того, що розраховане опір системи становить Ом, можна зробити висновок про правильність розрахунку, так як рас-парне значення задовольняє умові.


Виробнича санітарія

 

Освітлення

 

У приміщенні лабораторії використовується штучний тип освітлений-ня. Освітлювальна система повинна забезпечувати високу якість осве-щенности, яке визначається найменшим розміром розглянутого предмета і контрастом його з фоном; рівномірний розподіл яскравості по робочій поверхні; відсутність різких тіней на робочій поверхні; відсутність прямої і відбитої блескости; відсутність пульсації коливання і освітленості в часі.

 

Розрахунок штучного освітлення

 

Завданням розрахунку є визначення потрібної потужності електричної освітлювальної установки для створення в приміщенні лабораторії заданої освітленості. Розрахунок проводиться за такою методикою.

Вибираємо тип джерела світла - газорозрядна лампа.

Система освітлення - загальна.

Тип світильника - ВВПЗ (влаговзривопилезахищення). У цих све-тільніках встановлено дві і більше ламп, що дає можливість зменшити пульсацію сумарного світлового потоку світильника і виключити стробо-скопически ефект. У приміщенні, де знаходиться установка, встановлено штучне освітлення.

Необхідно розподілити світильники і визначити їх кількість. Забезпечення рівномірного розподілу джерела досягається в тому випадку, якщо відношення відстаней між центрами світильників (L) до висоти їх підвісу над робочою поверхнею складає 1,5 для світильника типу ПВЛ.

Нормована освітленість для зорової роботи високої точності становить Ен = 400 люкс. Висота приміщення Н = 3 м, ширина В = 5 м, довжина А = 6м.

Нехай відстань від центру люстри до робочого місця hр = 2 м. Величина відстані між центрами світильників становить

 

м.(5.9)

Кількість світильників становить

 

(5.10)

 

шт

Визначити індекс приміщення

 

(5.11)

 

 

При коефіцієнтах відображення стелі і стіни . Визначаємо коефіцієнт використання світлового потоку ламп: [ ]

Світловий потік виражається

 

.(5.12)

Коефіцієнт мінімальної освітленості для люмінесцентних ламп дорівнює Z = 1.1, коефіцієнт запасу дорівнює K = 1.5.
Світловий потік становить

 

лм.

 

Кількість ламп в світильнику складе 4 шт., Підбираємо необхідну лампу: лампа ЛДЦ 80.

потужність електричної освітлювальної установки становить

 

Вт.

Таким чином, в приміщення лабораторії рекомендується оборудо-вать електричної освітлювальної установкою з використанням чотирьох світильників ПВЛ з чотирма лампами ЛДЦ 80.

Порахуємо відхилення:

 

(5.13)

Розраховане відхилення знаходиться в межах норми, Δ має бути в межах від -10% до + 20%.

Догляд за освітлювальними приладами

 

Необхідно стежити за цілісністю плафонів, періодично (раз в шість місяців) витирати пил з освітлювальних приладів. Потрібно стежити за цілісністю проводки, а також змінювати розбиті плафони.

 

Розрахунок вентиляції приміщення

 

Розрахувати продуктивність системи вентиляції в приміщенні з заданими параметрами, з урахуванням перебування в приміщенні персоналу і експериментального устаткування

Початкові дані:

а = 6 м – длина приміщення;

b = 5м - ширина;

h = 3 м – висота.

Розрахунок:

Розрахуємо обсяг приміщення по формулі:

 

;(5.14)

м3

Витрата повітря при надходженні надлишкового тепла визначається за формулою:

 

,(5.15)

Де - продуктивність системи вентиляції ;

- Питома теплоємність повітря при постійному тиску;

– щільність припливного повітря;

0C – темпера тураповітря, що видаляється.

Для теплого періоду року і робіт «легка 1б» (робота, виконувана сидячи, якимось чином пов'язана з невеликими переміщеннями, що не вимагає великих фізичних витрат, тепловиділення порядку 185 Вт) температура повітря повинна бути від 0C [ ];

 0Cробоча різниця температур,

 - температура припливного повітря, 0C.Обчислюється за формулою:
 

;(5.16)

0C.

- Загальна кількість тепла, що надходить в приміщення, Вт;

 

,(5.17)

 

Де - Теплопритоки від навколишнього середовища, Вт;

- Теплопритоки від експериментального обладнання, Вт;

- Теплопритоки від персоналу, Вт;

- Теплопритоки від освітлення, Вт;

- кількість тепла, що виділяється установкою;

- Тепловиділення людини;

 

Qoc=Vnom∙30 = 320∙30 = 9600;(5.18)

Qob= qnk∙n = 1200∙1 = 1200;(5.19)

Qч = qч∙n = 170∙4 = 680;(5.20)

Qосв = 1280– з відповідних розрахунків.

∑Q = 9600+1200+680+1280 = 12760 ;

L = 12380 / (1,006∙1,2∙ (24 - 20))=2642,5 .

Визначити встановлену потужність електродвигуна вентилятора за формулою:

,(5.21)

де – коефіцієнт запасу ();

- Аеродинамічний опір вентиляційної мережі повітропроводів, ;

- ККД вентилятора (приймається за влучним висловом вентилятора );

- ККД приводу, який при клинопасовій передачі дорівнює

=0,95:

Висновок - приймаємо для вентиляції вентиляційний агрегат Повітрявиробництво L = 2642,5 вентилятор потужністю N = 0,35 .

Система вентиляції не повинна викликати переохолодження або перегріву робітників. Система вентиляції не повинна створювати шум на робочих місцях, що перевищує гранично допустимий рівень. Система вентиляції повинна бути електро-, пожежо-і вибухобезпечної, проста в пристрої, надійна в експлуатації і ефективна.

 

Противопожарная профилактика

 

Матеріали, що застосовуються для огороджувальних конструкцій та оздоблення приміщень, повинні бути вогнестійкими. У будівлі повинно бути передбачено не менше двох евакуаційних виходів. У приміщеннях з установками необхідно обмежити застосування дерева. Проходи, коридори і робочі місця не слід захаращувати архівними матеріалами, приладами і т.д. В системі вентиляції повинні бути передбачені протипожежні клапани в місцях перетину вогнестримувальних перекриттів і стін з нормованим межею вогнестійкості. Повітроводи, вентиляційні камери та регулюючі пристрої систем виконуються з негорючих матеріалів.

Протипожежна автоматична сигналізація

Пожежна автоматична сигналізація служить для швидкого сповіщення служби пожежної охорони про яка виникла пожежі в будь-якому приміщенні або спорудженні підприємства. У системи автоматичного пожежогасіння включається також і пожежна сигналізація. При необхідності пожежна сигналізація може бути поєднана з охоронною сигналізацією.

Пожежні сповіщувачі перетворюють неелектричні фізичні величини (випромінювання теплової і світлової енергії, рух частинок диму) в електричні, які у вигляді сигналу певної форми направляються по проводах на приймальну станцію.

Система пожежогасіння

В практика гасіння пожеж найбільшого поширення полоучілі наступні Принцип припинення горіння: 1) ізоляція вогнища горіння від повітря негорючими газами; 2) охолодження вогнища горіння нижче певних температур; 3) інтенсивне гальмування (інгібування) швидкості хімічної реакції в полум'ї; 4) механічний зрив полум'я в результаті впливу на нього сильного струменя газу; 5) створення умов огнепрегражденія, тобто таких умов при яких полум'я поширюється через вузькі канали.

У лабораторіях гасіння пожеж краще вести першим способом.

Розрахунок кількості балонів з вуглекислотою для внутрішнього пожежогасіння

Кількість вогнегасної газового складу (кг):

(5.22)

Де - вогнегасна концентрація газового складу для вуглекислоти;

- обсяг приміщення, м3;

- Коефіцієнт, який враховує особливості газообміну, витікання вуглекислого газу через нещільності приміщення (приймаємо )

кг.

Необхідна кількість робочих балонів (шт.) З вуглекислотою:

(5.23)

Де = 25 л - обсяг балона (при 25 л в балоні міститься 15,6 кг вуглекислоти);

 

= 0.625 кг/л - щільність засоб гасіння;

= 1 - коефіцієнт наповнення.

шт.

Кількість резервних балонів (приймаємо рівним числу робочих балонів):

шт.

Загальна кількість балонів:

шт.

Для можливості ліквідації джерела пожежі установка газового пожежогасіння всередині приміщення повинна мати 14 балонів з вуглекислотою.

 

Висновок:

У цій частині дипломної роботи були викладені вимоги до робо-чого місця інженера. Створені умови повинні забезпечити комфортну роботу. Був проведений розрахунок вентиляційної системи, розрахунок оптимального освітлення виробничого приміщення, а також розрахована система штучного заземлення та кількість балонів з вуглекислотою для пожежогасіння.

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

 

 

 

 

 


6. Організаційно – економічна частина

Проект має на меті, сфера використання та можливі користувачі проектного продукту (підприємства, організації).

Метою проекту є вивчення коефіцієнта тепловіддачі при кипінні у вільному обсязі хладагентов, а так само їх сумішей з маслами або нанофлюідамі з метою поліпшення енергетичної ефективності холодильних машин і зменшення екологічного впливу на навколишнє середовище.

Сферою використання є побутові та комерційні компресорні холодильні установки.

Користувачами є виробники побутових і комерційних холодильних установок.

Основні технічні, соціальні, екологічні, економічні показники проекту і їх зіставлення з аналогічними (альтернативними варіантами).

Збільшення коефіцієнта тепловіддачі за рахунок исспользование хлада-гент нового покоління, а так само вибору оптимальної концентація масла в реальному робочому тілі холодильної установки покращує її енергетичні характеристики та зменшує антропогенний вплив на навколишнє середовище в порівнянні з існуючими робочими тілами цих холодильних машин.

Можливі ризики і конкурентоспроможність проектного продукту на внутрішньому та зовнішньому ринках.

Можливим ринком збуту результату дослідження є холодильна промисловість, в тому числі провідні фірми як нашої країни, так і за кордоном, які виробляють холодильні компресора і інше обладнання для дование (випарники і конденсотрои) для холодильних машин.

Учасники проекту, фактори зовнішнього середовища. Можливі джерела фінансування і фактори ризику.

Учасниками проектіов виступають: розробник, дипломних керівник, ОГАХ, фірма NORD. Фактори зовнішнього середовища такі як економічна і політична ситуація на проведення досліджень не впливають.

Можливості та умови реалізації проекту.

В ході розробки дипломного проекту кінцевим продуктом розробки ми отримаємо інформацію про коефіцієнт тепловіддачі при кипінні хла-дагентов їх розчинів з мастилами при різних кончентраціях, що в свою чергу дозволить точно провести розрахунок холодильного обладнання та вибрати оптимальні речовини для використання. З точки зору економіки це дозволить з економити значні кошти при роботі холодильного обладнання.

Організація виконання проекту

Календарний план розробки проекту (основні етапи, терміни їх виконання)

Побудова і розрахунок мережного плану-графіка.

При виконанні технічних виробів, будівництва різних со-оружений, підготовки і проведення наукових і практичних експериментів доводиться виконувати велику кількість робіт. Виконання їх пов'язане в часі: почати одні роботи можна тільки тоді, коли виконано одну або декілька інших робіт, підготовчих до них робіт. Інформацію про взаємозв'язок і порядок виконання всіх робіт в комплексі можна уявити за допомогою орієнтованого графіка, який доповнюють певними числовими даними, такими, наприклад, як кількість часу, необхідного для завершення конкретних робіт, час початку певної роботи, номери вершин графіка і т.д. Такий графік називають сітковим графіком або гратчастої моделлю. У таблиці 1 представлений перелік подій і робіт сіткового графіка.

На основі даних таблиці 1 проведено розрахунок параметрів мережного графіка, представлений в таблиці 1 і мал. 1

Таблиця 1 - Перелік подій і робіт сіткового графіка

№ кода роботи

Назва роботи

Т (дні)

№ коду попередньої роботи

1

Збір даних і аналіз існуючого стану

7

0-1

2

Вивчення предметної області

7

0-2

 

фіктивна робота

0

1-3

 

фіктивна робота

0

2-3

3

Планування термінів необхідних робіт

7

3-4

4.

Затвердження зведеного плану

2

4-6

5.

Організаційне виконання робіт

3

5-6

6.

Проектування дослідницької установки

21

6-7

7.

Проектування дослідницької частини робочої установки

14

6-8

 

фіктивна робота

0

7-9

 

фіктивна робота

0

8-9

8.

Установка ділового контакту

2

9-10

9.

Контроль за виконанням робіт

2

10-11

10.

Розробка схеми осередку

14

11-12

11.

Розробка схеми електровимірювань

7

12-13

12.

Аналіз поправок до вимірювання коефіцієнта тепловіддачі

7

13-14

13.

Проведення розрахунків за отриманими дослідними даними

7

14-15

14.

Аналіз результатів

7

15-16

15.

Підтвердження завершення робіт

2

16-17

16.

Порівняння отриманих даних з існуючими моделями

2

17-18

17.

Підготовка до здачі

2

18-19

18.

Оцінка результатів проекту та підтвердження підсумків

2

19-20

Всього

 

115

 

 

Таблиця 2 - Результати розрахунків параметрів мережного графіка

Кількість робіт що входять в подію

Работи

Ранній термін виконання події

ТPI

тривалість роботи

Тij

 

Ранній термін виконання події ТPI

Пізній термін виконання події Тin

тривалість роботиТij

тривалість події

повний резерв Rj

i

j

1

2

3

4

5

6

7

8

9

10

-

0

1

0

7

7

0

7

7

0

1

0

2

0

7

7

0

7

7

0

1

1

3

7

0

7

7

0

7

0

1

2

3

7

0

7

7

0

7

0

2

3

4

7

7

14

7

7

14

0

1

4

5

14

2

16

14

2

16

0

1

5

6

16

3

19

16

3

19

0

1

6

7

19

21

40

19

21

40

0

1

6

8

19

14

33

19

14

40

7

1

7

9

40

0

40

40

0

40

0

1

8

9

33

0

33

33

0

40

7

2

9

10

40

2

42

40

2

42

2

1

10

11

42

2

44

42

2

44

2

1

11

12

44

14

58

44

14

58

14

1

12

13

58

7

65

58

7

65

0

1

13

14

65

7

72

65

7

72

0

1

14

15

72

7

79

72

7

79

0

1

15

16

79

7

86

79

7

86

0

1

16

17

86

2

88

86

2

88

0

1

17

18

88

2

90

88

2

90

0

1

18

19

90

2

92

90

2

92

0

1

19

20

92

2

94

92

2

94

0


7. Висновки

В ході дослідження були вирішені такі завдання:

- Розроблена та проектована експериментальна установка для дослідження теплопровідності методом лінійного джерела (нестаціонарний метод нагрітої нитки), що дозволяє вимірювати теплопровідність нанофлюідів;
- Розроблена методика проведення досліджень теплопровідності нанофлюідів при різних температурах і концентраціях наночастинок.

- Проаналізовано літературні дані по теплопровідності нанофлюідів.

- Досліджені та Розраховані технологічні елементи установки.

8. Список використаних джерел

 

 

1.                 Chiesa M, Simonsen AJ: The importance of suspension stability for hotwire measurements of thermal conductivity of colloidal suspensions. 16th Australasian Fluid Mechanics Conference, Gold Coast, Australia 2010.

2.                 Masuda H., Ebata A., Teramae K., Hishinuma N., Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles (dispersion of c-Al2 O3, SiO2 and TiO2 ultra-fine particles), Netsu Bussei 4 (1993) 227–233.

3.                 Xie H, Wang J, Xi T, Liu Y, Ai F: Thermal conductivity enhancement of suspensions containing nanosized alumna particles. J Appl Phys 2002, 91:4568-4572.

4.                 Wang X., Xu X., Choi S.U.S., Thermal conductivity of nanoparticles–fluid mixture, J. Thermophys. Heat Transfer 13 (1999) 474–480.

5.                 Eastman JA, Choi SUS, Li S, Yu W, Thomson LJ: Anomalously increased effective thermal conductivities of ethylene glycol-based nanofluids containing copper nanoparticles. Appl Phys Lett 2001, 78:718-720.

6.                 Lee S., Choi S.U.S., Li S., Eastman J.A., Measuring thermal conductivity of

fluids containing oxide nanoparticles, ASME J. Heat Transfer 121 (1999)

280–289.

7.                 Mintsa HA, Roy G, Nguyen CT, Doucet D: New temperature dependent thermal conductivity data for water-based nanofluids. International Journal of Thermal Sciences 2009, 48:363-371.

8. Das S, Putra N, P. Thiesen, W. Roetzel, Temperature dependence of thermal conductivity enhancement for nanofluids, J. Heat Transfer, 125, 2003, p. 567–574.

9. Putra N., Roetzel W., Das S.K., Natural convection of nanofluids, Heat Mass Transfer 39 (2003) 775–784.

10. Chon C.H., Kihm K.D., Lee S.P., Choi S.U.S., Empirical correlation finding the role of temperature and particle size for nanofluid (Al2O3) thermal conductivity enhancement, Appl. Phys. Lett. 87 (2005) 153107.

11. Li C.H., Peterson G.P., The effect of particle size on the effective thermal

conductivity of Al2O3–water nanofluids, J. Appl. Phys. 101 (2007) 044312.

12. Zhang X., Gu H., Fujii M., Effective thermal conductivity and thermal

diffusivity of nanofluids containing spherical and cylindrical nanoparticles,

J. Appl. Phys. 100 (2006) 1–5.

13. Timofeeva E.V., Gavrilov A.N., McCloskey J.M., Tolmachev Y.V., Thermal

conductivity and particle agglomeration in alumina nanofluids: experiment

and theory, Phys. Rev. E 76 (2007) 061203.

14. Hwang D, Hong KS, Yang HS: Study of thermal conductivity nanofluids for the application of heat transfer fluids. Thermochim Acta 2007, 455:66-69.

15. Murshed SMS, Leong KC, Yang C: Invesitions of thermal conductivity and viscosity of nanofluids. Int J Therm Sci 2008, 47:560-568.

16. Oh DW, Jain A, Eaton JK, Goodson KE, Lee JS: Thermal conductivity measurement and sedimentation detection of aluminum oxide nanofluids by using 3ω method. Int J Heat Fluid Fl 2008, 29:1456-1461.

17. Beck MP, Sun T, Teja AS: The thermal conductivity of alumina nanoparticles dispersed in ethylene glycol. Fluid Phase Equilib 2007, 260:275-278.

18. Choi C, Yoo HS, Oh JM: Preparation and heat transfer properties of nanoparticles-in-transformer oil dispersions as advanced energy-efficient coolants. Curr Appl Phys 2008, 8:710-712.

19. Rusconi R, Rodari E, Piazza R: Optical measurements of the thermal

properties of nanofluids. Applied Physics Letters 2006, 89, 261916-1-2619163.

20. Putnam SA, Cahill DG, Braun PV: Thermal conductivity of nanoparticle

suspensions. Journal of Applied Physics 2006, 99, 084308-1-084308-6.

21. Venerus DC, Kabadi MS, Lee S, Perez-Luna V: Study of thermal transport in

nanoparticle suspensions using forced Rayleigh scattering. Journal of

Applied Physics 2006, 100, 094310-1-094310-5.

22. Buongiorno J, Venerus DC, Prabhat N: A benchmark study on the thermal

conductivity of nanofluids. Journal of Applied Physics 2009, 106, 094312-1-

094312-14.

23. Bazan JAN: Thermal conductivity of poly-aelpha-olefin (PAO)-based nanofluids. Ph.D. Thesis, University of Dayton, Dayton, OH, USA 2010.

24. Schmidt AJ, Chiesa M, Torchinsky DH, Johnson JA, Nelson KA, Chen G:

Thermal conductivity of nanoparticle suspension in insulating media measured with a transient optical grating and a hotwire. Journal of Applied Physics 2008, 103, 083529-1-083529-5.

25. Ali FM, Yunus WMM, Moksin MM, Talib ZA: The effect of volume fraction

concentration on the thermal conductivity and thermal diffusivity of

nanofluids: numerical and experimental. Review of Scientific Instruments 2010, 81, 074901-1-074901-9.

26. Kolade B, Goodson KE, Eaton JK: Convective performance of nanofluids in

a laminar thermally developing tube flow. Journal of Heat Transfer 2009, 131, 052402-1-052402-8.

27. Geller, V., Nemzer, B., Cheremnykh, U. Thermal Conductivity of the Refrigerant Mixtures R404A, R407C, R410A, and R507A. International Journal of Thermophysisc, Vol. 32, N. 5, pp. 48-56, 2002.

28. Li C.H., Peterson G.P. Mixing effect on the enhancement of the effective thermal conductivity of nanoparticle suspensions (nanofluids). // Inter. J. of Heat and Mass Transfer. 2007. Vol. 50. P. 4668−4677.

29. Keblinski P, Phillpot S.R., Choi S.U.S., Eastman J.A. Mechanisms of heat flow in suspensions of nano-sized particles (nanofluids) // Inter. J. of Heat and Mass Transfer. 2002. Vol. 45.Р. 855−863.

30. Ding Y.L, Chen H., Wang L., Yang C-Y, He Y., Yang W., Lee W.P., Zhang L. and Huo R. Heat Transfer Intensificaton Using Nanofluids // Powder and Particle. 2007. No. 25. P. 23−36.

31. Yu W., Choi S.U.S. The role of interfacial layers in the enhanced thermal conductivity of nanofluids: a renovated Maxwell model // J. of Nanoparticle Research.2003. No. 5. P. 167−171.

32. Wang B.-X., Zhou L.-P., X.-F. Peng X.-F.A fractal model for predicting the effective thermal conductivity of liquid with suspension of nanoparticles // Inter. J. of Heat and Mass Transfer. 2003. Vol. 46. P. 2665−2672.

33. Xue L., Keblinski P., Phillpot S.R., Choi S.U.S., Eastman J.A. Effect of liquid layering at the liquid - solid interface on thermal transport // Inter. J. of Heat and Mass Transfer. 2004. Vol. 47, No. 19-20. P. 4277−4284.

34. Gao L., Zhou X., Ding Y.L. Effective thermal and electrical conductivity of carbon nanotube composites //Chem. Phys. Lett. 2007. Vol. 434. P. 297−300.

35. Елецкий А.В. Транспортные свойства углеродных нанотрубок// УФН. 2009. Т. 179, Nо. 3. С. 225−242.

36. Maxwell J.C., A Treatise on Electricity and Magnetism, second ed., Clarendon Press, Oxford, UK, 1881.

37. Bruggeman D.A.G., Berechnung verschiedener physikalischer Konstanten von heterogenen Substanzen, I. Dielektrizitatskonstanten und Leitfahigkeiten der Mischkorper aus isotropen Substanzen, Annalen der Physik, Leipzig 24 (1935) 636–679.

38. Hamilton R.L., Crosser O.K., Thermal conductivity of heterogeneous two component systems, I&EC Fundam. 1 (1962) 182–191.

39. Xue Q.-Z., Model for effective thermal conductivity of nanofluids, Physics Letters A307 (2003) 313–317.

40. Xue Q., Xu W.-M., A model of thermal conductivity of nanofluids with interfacial shells, Materials Chemistry and Physics 90 (2–3) (2005) 298–301.

41. Patel H.E., Das S.K., Sundararagan T., Nair A.S., Geoge B., Pradeep T., Thermal conductivities of naked and monolayer protected metal nanoparticle based nanofluids: Manifestation of anomalous enhancement and chemical effects, Applied Physics Letters 83 (2003) 2931–2933.

42. Xuan Y., Li Q., Hu W., Aggregation structure and thermal conductivity of nanofluids, AIChE Journal 49 (4) (2003) 1038–1043.

43. Kumar D.H., Patel H.E., Kumar V.R.R., Sundararajan T., Pradeep T., Das S.K., Model for heat conduction in nanofluids, Physical Review Letters 93 (14) (2004) 144301.

44. Jang S.P., Choi S.U.S., Role of Brownian motion in the enhanced thermal conductivity of nanofluids, Applied Physics Letters 84 (2004) 4316–4318.

45. Prasher R., Bhattacharya P., Phelan P.E., Thermal conductivity of nanoscale colloidal solutions (nanofluids), Physical Review Letters 94 (2) (2005) 025901.

46. Koo J., Kleinstreuer C., A new thermal conductivity model for nanofluids, Journal of Nanoparticle Research 6 (6) (2004) 577–588.

47. Koo J., Kleinstreuer C., Laminar nanofluid flow in micro-heat sinks, International Journal of Heat and Mass Transfer 48 (13) (2005) 2652–2661.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Источник: портал www.KazEdu.kz

Каталог учебных материалов

Свежие работы в разделе

Наша кнопка

Разместить ссылку на наш сайт можно воспользовавшись следующим кодом:

Контакты

Если у вас возникли какие либо вопросы, обращайтесь на email администратора: admin@kazreferat.info