Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук




Скачать 24.11 Kb.
НазваниеУчебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук
страница6/10
Дата03.02.2016
Размер24.11 Kb.
ТипУчебное пособие
1   2   3   4   5   6   7   8   9   10



Выходные данные:

Физический смысл переменной

Название в файле

Фактор эффективности поглощения,

Qa


Фактор эффективности рассеяния,

Qs

Фактор эффективности экстинкции,

Qext

Коэффициент поглощения, , см-1

ba

Коэффициент рассеяния, , см-1

bs

Коэффициент экстинкции, , см-1

bext

Спектральная испускательная способность,

epsl

Спектральная светимость, , Вт/см3

rl

Концентрация, см-3

cont


Файл задания организуется следующим образом. Сначала записываются входные данные в следующем виде [32]:


Тип данных

Запись в файле

Пример

Одно значение

переменная значение

r0 0.02

Набор значений

переменная значение значение …

wavel 0.435 0.465 0.49 0.525

Цикл

переменная

c

начальное_значение

конечное_значение

шаг

rad_change c 0.002 0.2 0.002


Здесь значение может быть как числом (также и в экспоненциальной форме, например 1.e+4, т.е.104), так и именем переменной, объявленной выше. Например:sigma 0.5 sigmabase sigma.


Цикл отличается от набора значений, во-первых, “c” (cycle – цикл) после названия переменной, во-вторых, фиксированным набором чисел – это, соответственно, начальное, конечное значения переменной цикла и шаг.

В файле обязательно должен быть один (и только один) цикл или одна (и только одна) последовательность. Иными словами, только одна из переменных должна изменять свое значение.

После подготовки входных данных нужно поставить заранее предопределенный символ EOD (end of data – конец данных).

Затем нужно записать, какие данные будут выводиться. Названия переменных берутся из приведенной выше таблицы выходных данных.

Оригинальная программа [32] создавалась и тестировалась под операционной системой Unix, где проблем, которые будут описаны ниже, не возникало. При адаптации этой программы под ОС Windows, компилятор Borland С++ 5.02, обнаружились сложности, связанные с вводом из файла. Поэтому при написании файла задания нужно разделять имя переменной и значение(я) или имя переменной и символ «c» не одним, а двумя пробелами, например r0 пробел 0.02 (хотя для первой строки этого можно и не делать). Далее, при перечислении выходных данных нужно размещать список через один пробел от левого края (т. е., например, пробел rl epsl ba bs). Кроме того, после списка выходных данных нужно поставить пустую строку (Enter).

После того, как файл создан, нужно выполнить команду %prog_dir%\w.exe %task_res_dir%\файл_задания %task_res_dir%\файл_результатов, где prog_dir и task_res_dir – пути к директориям, содержащим программные файлы и файлы заданий и результатов, соответственно. Программа может выявить ошибки в построении файла, при этом будет выдан код ошибки и сообщение о ее типе, а также номер содержащей ошибку строки.


Код

ошибки

Диагностика

1

Ошибка при открытии файла …

2

Имя для замены, «…», не определено выше

3

Нельзя начинать строку с числового значения, лишь с имени (выдается также при неправильном построении цикла)




Должен быть лишь один цикл или последовательность, здесь их …, а именно: (имена переменных);




Во входных данных присутствует неизвестная переменная «…» со значением …




В выходных данных присутствует неизвестная переменная …


Если программа все же перестанет в процессе работы выдавать результаты или реагировать на запросы («зависнет»), ее можно прервать последовательным нажатием сначала клавиша Ctrl и, удерживая его, C (Ctrl+C).

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


3.2.3 Излучение индивидуальных частиц микронного и субмикронного размера.


В части 2 (формула 2.21) было показано, что интенсивность источника сплошного спектра связана со спектральной светимостью источника сравнения rλ соотношением:


, (3.27)

где γλ , αλ - величины электрических сигналов при регистрации излучения исследуемого источника и источника сравнения, LEFF - эффективная протяженность источника излучения, к методике определения которой мы вернемся ниже. Очевидно, что если источник считать поверхностным со спектральной светимостью rλ и, то:


rλ и = (γλ λ)rλ . (3.28)


Сравнивая выражения (3.27) и (3.28) получаем весьма полезное соотношение:


I0 = 4rλ и /LEFF. (3.29)


Отметим, что приведенные выше соотношения верны для оптически прозрачных источников, когда длина свободного пробега излучения по отношению к поглощению la и рассеянию ls велика по сравнению с протяженностью источника LEFF. Длину свободного пробега для монодисперсных сферических частиц можно оценить согласно соотношению:


la,s = 1/(σ a,sNk)= 1/(Qa,sπr2Nk) . (3.30)


Оценим концентрацию частиц во фронте пламени, для которой можно считать его оптически прозрачным. С учетом сказанного Nk<1/(Qa,sπr2LEFF).

Для характерной величины LEFF ~ 0.1 см [33,34] и Qa,s=1, получаем нижнюю оценку для концентрации частиц:

а) горение газофазное [33,34], r ≈ 5·10-6см , N k< 1011см-3. Полученная величина, можно сказать, является критической. Пламя в коротковолновой части спектра не является оптически прозрачным. Однако, начиная с ближнего ИК диапазона, как показывают оценки с использованием реальных значений факторов эффективности процессов, пламя можно принять оптически тонким.

б) горение гетерогенное, r ≈ 2·10-4см, N k< 108см-3. Полученная величина более чем на два порядка превышает типичное значение Nk во фронте пламени, которое в связи с этим можно считать оптически тонким.

Перейдем к излучательным характеристикам индивидуальных частиц, учитывая то, что под спектральной интенсивностью излучения понимают энергию, излучаемую единицей объема в единицу времени в единичном спектральном интервале длин волн. Полагая, что спектральная испускательная способность частицы равна фактору эффективности поглощения Qa, легко видеть, что для частиц сферической формы c учетом (3.29) :


Iλ=Nkrλb4πr2Qa(r,λ,n,χ)f(r)dr=4rλи/Lэфф, (3.31)


где f(r)- функция распределения частиц по размерам.

Формула (3.31) имеет наиболее простой вид для субмикронных частиц. Принимая во внимание выражение (3.25), можно записать:


Qa(r,λ,n,χ)=а(λ,n,χ)r. (3.32)

И затем:


Nk rλb3>а(λ,n,χ)=4rλи/Lэфф , (3.33)


где 3>=r3f(r)dr.

Выражение (3.33) позволяет по результатам измерений спектральной светимости, определить одну из величин Nk ,3> или а(λ,n,χ) при условии,

что известны какие-то две из них, измеренных не зависимым методом. Например, концентрацию частиц Nk можно определить по рассеянию лазерного излучения, 3> - по данным дисперсного анализа. Протяженность источника излучения для осесимметричных источников Lэфф можно определить из радиальных распределений световых потоков используя преобразования Абеля. К этому вопросу мы вернемся ниже. Особенно просто это можно сделать для оптически тонких источников.

Для относительно крупных частиц, для которых релеевское приближение не годится, выражение (3.31) позволяет определить Nk либо мнимую часть показателя преломления. Принимая во внимание (3.21), а также то, что N(r)=Nkf(r), формулу (3.31) можно привести к виду:


Nk βλ/Lэфф , (3.34)


где ελ=rλи/rλb.

Если известен комплексный показатель преломления вещества и проведен дисперсный анализ, то используя теорию Ми можно рассчитать спектральный коэффициент поглощения и затем определить Nk. Можно решить и обратную задачу- найти величину χ(λ,Т), которая имеет исключительную ценность при решении задач о переносе излучения в дисперсных системах. К сожалению, результаты, полученные различными авторами, сильно разнятся как количественно, так и качественно. Сказанное иллюстрирует рисунок 3.9 заимствованный из работы [35]. Ранее мы уже упоминали, что действительная часть показателя преломления слабо зависит от температуры, поэтому считается известной. Для этого расчетным путем находится зависимость β от χ и далее, по экспериментально найденной величине β, графическим путем находим мнимую часть показателя преломления χ(λ,Т).




Рис.3.9 Зависимость мнимой части χ(Т) Al2O3 по данным различных работ суммированных в [35]:а. 1-1473К, 2-1973К, 3-2293К; b.1-2240К, 2-2440К; с- 1-2090К, 2-2550К; d. 2950К; e. 2400-2700К;f. 1-2320К,2-2400К,32600К,4-2900К,5-3000К; g. Данные автора [35] 1-2350К, 2-2500К, 3-2950К; h-расчетный способ.


Обратимся теперь к методике определения эффективной протяженности источника излучения Lэфф.


3.3 Преобразование Абеля. Переход от радиальных распределений световых потоков к локальным интенсивностям. Протяженность оптически тонкого источника излучения.

Кратко остановимся теперь на методике определения ширины осесимметричного источника. Исходной информацией являются радиальные распределения световых потоков, которые можно получить следующим образом. На рисунке 3.10а представлена блок -схема экспериментальной установки, с помощью которой в работе [48] получали радиальные распределения световых потоков зоны конденсации (источник сплошного спектра) горящей одиночной частицы магния.



Рис.3.10а. Блок- схема установки. Рис. 3.10б.

Рис.3.10б. Радиальное распределение излучения (1-частица, 2- зона конденсации).

Излучение от горящей частицы магния 5 после отражения от зеркала 1 фокусировалось объективом 2 в плоскость входной скрещенной щели спектрального прибора 4. Зеркало совершало колебательное движение по пилообразному закону. Благодаря этому изображение объекта смещалось в плоскости входной щели, что и позволяло регистрировать радиальные распределения. Диафрагма 3 использовалась для улучшения пространственного разрешения и увеличения глубины резкости. На рисунке 3.10б. представлено типичное радиальное распределение. Увеличение интенсивности излучения зоны конденсации (II) при смещении от центральной области к периферийной связано с изменение толщины излучающего слоя. В связи с этим становится понятны характер распределения, представленного на рис.3.10б.

Отметим, что изображение источника можно смещать иначе, например, путем поворота плоскопараллельной стеклянной пластинки толщиной не менее 10см. Пластинку следует размещать между фокусирующим объективом и скрещенной щелью спектрального прибора. Можно использовать для этого ПЗС матрицу, на которую фокусируется изображение источника. Излучение при этом следует пропустить через световой фильтр. Обратимся к методике перехода от радиальных распределений световых потоков к локальным интенсивностям с помощью преобразований Абеля. На рисунке 3.11 показан срез осесимметричного источника. Световой поток Ф(y) снимается в направлении оси х для различных значений координаты у. При этом пространственное разрешение должно быть таким, чтобы обеспечить регистрацию пространственной структуры пламени. Например, ширина фронта факела газовзвеси частиц металла составляет приблизительно 1мм. Соответственно пространственное разрешение должно быть не хуже 0.1мм. Функция Ф(y) представляет собой искомое радиальное распределение.




Рис.3.11. Иллюстрация к преобразованиям Абеля.


Пусть I(r) локальная интенсивность. Несложно сообразить, что Ф(y) и I(r) связаны соотношением:

, (3.35 )

где I(r) - локальная интенсивность на расстоянии r от оси симметрии источника z; R - радиус источника, определяемый условием I(r=R) =0. В выражении (3.35) мы учли, что и провели соответствующую замену переменных. Обращение уравнения (3.35) позволяет найти I(r):

. (3.36)

Расчеты можно провести с помощью компьютера с использованием какого либо пакета программы, например," Origin", "Maple" или "Mathcad". За протяженность источника можно принять полуширину распределения I(r) или определить ее согласно выражению: LEFF =(см. часть 2),где ζ(r)=I(r)/I0. Типичные зависимости Ф(y) и I(r) вблизи вершины ламинарного двухфазного диффузионного (ЛДД) факела Al показаны на рисунке 3.12.



Рис.3.12 Радиальные распределения светового потока.


Задание 1.

В таблице 3.2 представлены сглаженные результаты обработки радиальных распределений световых потоков ЛДД факела Al (y0=9мм). Используя преобразования Абеля найти:

  1. зависимость интенсивности от радиальной координаты факела I(r);

  2. эффективную ширину фронта факела L'эфф.


Таблица 3.2 Радиальное распределение светового потока.

y/y0

Ф(y/y0) отн

y/y0

Ф(y/y0) отн

y/y0

Ф(y/y0) отн

0

0,47

0,51

0,52

0,79

0,68

0,051

0,47

0,54

0,56

0,82

0,63

0,1

0,47

0,56

0,65

0,85

0,56

0,15

0,47

0,59

0,81

0,87

0,5

0,21

0,47

0,62

0,87

0,9

0,41

0,26

0,47

0,64

1

0,92

0,3

0,31

0,47

0,67

0,94

0,95

0,2

0,36

0,47

0,69

0,88

0,97

0,15

0,41

0,47

0,72

0,87

1

0,07

0,46

0,47

0,74

0,8







0,49

0,5

0,77

0,73









Задание 2.

В таблице 3.3 представлены результаты измерений спектральной светимости ЛДД факела частиц алюминия с наиболее вероятным диаметром 14 мкм. В результате газофазного горения образуются частица Al2O3 сферической формы, функция распределения которых по размерам хорошо описывается логарифмически нормальным распределением:


f(d)= ((2πσ2)0.5d)-1exp(-(ln(d/d0))2/2σ2), где: σ=0.4; d0=0,08 мкм.


Таблица 3.3 Спектральная светимость ЛДД факела частиц алюминия.

λ, мкм

rλ, Вт/см3

λ, мкм

rλ, Вт/см3

0.45

2800

0.8

27800

0.5

7350

0.85

29500

0.55

11400

0.9

28350

0.6

14650

0.95

30000

0.65

17850

1.0

28000

0.75

25500

1.05

30000


Принимая концентрацию частиц равной Nk=1011см-3 и Lэфф=0.1см [33], вычислите коэффициент пропорциональности между факторами эффективности поглощения и радиусом частицы Al2O3 для различных длин волн при температуре горения Т=3150±100К. Определите величину Qa(λ) для частицы с наиболее вероятным радиусом.

Используя приближение Релея, рассчитайте согласно соотношению (3.24) факторы эффективности рассеяния Qs(λ) для заданного интервала длин волн. Действительную часть показателя преломления примите равной n=1.76 [36].

Найдите длину свободного пробега излучения по отношению к поглощению и рассеянию. Можно ли считать рассмотренную систему как оптически тонкую?

Сравните эффективность процессов рассеяния и поглощения излучения.


3.4 Измерение экстинкции излучения.


Экстинкция излучения определяется по отношению интенсивности падающего на среду и прошедшего сквозь нее излучения. Если многократным рассеянием можно пренебречь, то интенсивность светового пучка в среде с частицами затухает в соответствии с законом Бугера:


(3.37)

где: b - коэффициент экстинкции. Экстинкция зависит как от поглощения, так и от рассеяния и для монодисперсных частиц, равен:


(3.38)


где Nк – концентрация частиц; sa , ss – сечения поглощения и рассеяния. Некоторые формулы для удобства чтения мы повторяем намерено. Оба процесса происходят одновременно, но нередко один из них является определяющим. Для сажистых частиц решающую роль играет поглощение излучения [30,37]. И наоборот, для частиц оксидов металлов роль поглощения в экстинкции излучения в видимой и ультрафиолетовой спектральных областях несущественна [38]. Более того, в зависимости от длины волны вклад упомянутых процессов различен. В реальных системах частицы имеют некоторое распределение по размерам. В этом случае удобнее работать не с сечениями процессов, а с факторами их эффективности Qa и Qs. Указанные характеристики можно рассчитать на основе теории Ми. Для малых частиц (в сравнении с длиной волны) справедливо приближение Релея, в рамках которого Qa и Qs зависят от длины волны и размера частиц по-разному : Qa~ r/l, Qs~ (r/l)4 (см. соотношения (3.24, 3.25)). Поэтому для малых частиц можно выявить, какой из процессов – рассеяние или поглощение – играет в экстинкции решающую роль , если построить график зависимости lnbl от lnl. Тангенс угла наклона в первом случае должен быть близок к 1, во втором – к 4.

Для полидисперсной системы следует учитывать распределение частиц по размерам f(r). Для сферических частиц:


(3.39)


Для малых частиц величины bа и bs будут зависеть соответственно от третьего (поскольку Qa~r) и шестого (поскольку Qa~r4) моментов функции распределения частиц по размерам. Эта информация, как видим, является чрезвычайно важной и может быть получена с использованием электронного микроскопа или другим методом.

Обратимся теперь методике эксперимента. На рис.3.13 представлена принципиальная схема установки для измерения экстинкции излучения. Свет от "точечного" источника (1) колимируется линзой (2), проходит сквозь Рис.3.13 Принципиальная схема установки.


исследуемый объект (3) и фокусируется на скрещенную входную щель спектрального прибора (5). Щель в данном случае дополнительно играет роль весьма узкой диафрагмы, так что свет, рассеянный на малый угол, вообще не попадет на фотоприемник (6). Спектральный прибор позволяет проводить исследования на различных длинах волн. Благодаря диафрагмам (7) и (8) можно сформировать узкий параксиальный световой пучок, обеспечивающий, в частности, необходимое пространственное разрешение установки. Описанная установка позволяет исследовать частицы, размеры которых меньше длины волны λ [27]. Для частиц, больших по сравнению с λ, необходимо использовать устройства для устранения рассеянного вперед излучения (за счет дифракции), которое может давать основной вклад в экстинкцию излучения.

В качестве примера, иллюстрирующего особенности поглощения и рассеяния излучения дисперсной средой, приведем результаты исследования работы [38]. Измерения экстинкции проводились согласно описанной выше методике. Рассеивающими центрами являлись частицы окиси магния с характерным размером около 0,1 мкм с концентрацией Nk≈5·1011см-3. В качестве источника излучения использовались газоразрядные лампы высокого давления. Расчеты проводились согласно теории Ми (кривые1и2, рис.3.14) и в приближении Релея (кривая 3). Как видно, теория Ми хорошо описывает результаты измерений. Причем для длин волн λ< 0.7 мкм определяющий вклад в экстинкцию излучения вносит рассеяние. Вклад поглощения не превышает в этой области длин волн 10%. Отметим, что для наглядности на рисунке представлены зависимости 10·βа·l (намерен увеличен масштаб) и βs·l, так что βа·l на порядок величины меньше и βs·l. Для λ~1мкм роль обоих процессов сравнивается. Понятна и дальнейшая тенденция указанных зависимостей. Из рис.3.14 видно, что релеевское приближение хорошо "работает" для λ>0.7 мкм. При этом дифракционный параметр 2πr/λ для частиц с размером 0.1 мкм оказывается < 0.45, но не <<1, как этого требует приближение Релея.




Рис.3.14. Оптическая плотность зоны горения, обусловленная поглощением (10bаl) - кривая 1, и рассеянием излучения (bsl) – кривые 2, 3. Кривые 1, 2 – расчет по теории Ми, 3 – расчет в релеевском приближении. Точки – результаты измерений(bl): (bl) » (bsl).

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

Обратим внимание еще на одно обстоятельство. Измерение экстинкции относительно- простая задача. Гораздо труднее выделить вклад рассеяние и поглощения по отдельности. Для этого необходимо проводить измерения экстинкции и полного рассеяния излучения независимым образом. Тогда разность полученных результатов позволит найти, например, оптическую плотность по отношению к поглощению. Понятно, что, если рассеяние вносит основной вклад в экстинкцию ( βext·l) , то ошибки при определении βа·l могут оказаться неконтролируемо большими. Действительно, пусть вклад рассеяния в экстинкцию составляет > 90%. Причем погрешности измерений составляют 5% . Понятно, в этом случае нельзя ничего сказать о величине βа·l, кроме того, что она < 0.1· βext·l. Для определения полного рассеяния необходимо провести измерения интенсивности рассеянного света по всем направлениям, например, с помощью интегрирующей сферы. Возможно измерение поглощения непосредственно, так называемым фотоакустическим методом [27]. Суть его сводится к следующему. Объект освещается промодулированным светом (например, с помощью обтюратора) на частоте звукового диапазона. При этом частицы благодаря поглощению излучения периодически будут нагреваться, что приведет к колебаниям давления (к звуку). Последний можно зарегистрировать с помощью микрофона и при его надлежащей калибровке найти поглощение излучения.


1   2   3   4   5   6   7   8   9   10

Похожие:

Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconУчебное пособие Чита 2012 министерство образования и науки российской федерации федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Учебное пособие предназначено для студентов 1 2 курсов направлений 081100 Государственное и муниципальное управление, 080100 Экономика...
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconУчебное пособие разработано для подготовки студентов по специиальностям очной, очно-заочной форм обучения и экстерната по специальности «документоведение в управленческой деятельности»
...
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconІi международная научно-практическая конференция “ полупроводниковые материалы, информационные технологии и фотовольтаика
Украины «Киевский политехнический институт», Национальный технический университет «Харьковский политехнический институт», Харьковский...
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconУчебное пособие Санкт-Петербург 2010 удк 616-058: 93/99+61(06) isbn 57645-0189
В. П. Романюк, В. С. Лучкевич, И. Л. Самодова. История мировой и отечественной медицины: учебное пособие. – Спб.: Спбгма им. И. И....
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconУчебное пособие министерство образования российской федерации гоу впо «уральский государственный педагогический университет»
Учебное пособие предназначено для курса «История музыкального образования», который входит в федеральный компонент учебного плана...
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconМинистерство здравоохранения украины национальный медицинский университет имени А. А. Богомольца
Международная научно-практическая конференция ко всемирному дню здоровья, который
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconВера Савченко «В глухом переулке»
Дерибасовская – Ришельевская: Одесский альманах (сб.) Книга 38.– Одесса, 2009. С. 225-232
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconМинистерство образования и науки украины национальная академия наук украины
Б-63 Биоразнообразие и роль животных в экосистемах: Материалы V международной научной конференции. – Днепропетровск: Лира, 2009....
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconМинистерство образования и науки украины национальная академия наук украины
Б-63 Биоразнообразие и роль животных в экосистемах: Материалы IV международной научной конференции. – Днепропетровск: Изд-во дну,...
Учебное пособие Одесса 2006 министерство образования и науки украины одесский национальный университет им. И. И. Мечникова Кафедра общей и химической физики. А. В. Флорко, В. Г. Шевчук iconМинистерство образования и науки украины национальная академия наук украины
Биоразнообразие и роль зооценоза в естественных и антропогенных экосистемах: Материалы III международной научной конференции. – Днепропетровск:...
Разместите кнопку на своём сайте:
Библиотека


База данных защищена авторским правом ©lib2.znate.ru 2012
обратиться к администрации
Библиотека
Главная страница