Эхо Velocimetry Particle Image

Summary

Image Echo частиц велосиметрии (EPIV) система, способная приобретения двумерного поля скорости в оптически непрозрачных жидкостей или через оптически непрозрачных геометрии описаны и проверки измерений в трубе не сообщается.

Abstract

Перенос массы, импульса и энергии в потоках жидкости в конечном итоге определяется пространственно-временным распределением поля скорости жидкости. ¹ Следовательно, предпосылкой для понимания, прогнозирования и управления потоками жидкости является возможность измерения поля скоростей с адекватным пространственным и временным разрешением. ² Для измерения скоростей в оптически непрозрачных жидкостях или с помощью оптически непрозрачной геометрии скорость изображения эхочастиц (EPIV) является привлекательным диагностическим методом для создания «мгновенных» двумерных полей скорости. ^3,4,5,6 В данной статье описан рабочий протокол системы EPIV, построенной путем интеграции коммерческого медицинского ультразвукового аппарата⁷ с ПК, на котором установлено программное обеспечение 8 для измерения скорости изображения частиц (PIV) 8, а также приведены данные о валидационных измерениях в потоке по методу Хагена-Пуазей (т.е. в ламинарной трубе).

Для измерений EPIV преобразователь с фазированной решеткой, подключенный к медицинскому ультразвуковому аппарату, используется для создания двумерного ультразвукового изображения путем импульсного воздействия элементов пьезоэлектрического зонда в разное время. Каждый элемент зонда передает ультразвуковой импульс в жидкость, а частицы индикатора в жидкости (как естественного, так и затравленного) отражают ультразвуковое эхо обратно к зонду, где оно регистрируется. Амплитуда отраженных ультразвуковых волн и их временная задержка относительно передачи используются для создания так называемых двумерных ультразвуковых изображений в B-режиме (режим яркости). В частности, временная задержка используется для определения положения рассеивателя в жидкости, а амплитуда используется для назначения интенсивности рассеивателю. Время, необходимое для получения одного изображения в B-режиме, t, определяется временем, затраченным на импульс всех элементов преобразователя с фазированной решеткой. Для получения нескольких изображений в B-режиме частота кадров системы в кадрах в секунду (fps) = 1/δt. (См. 9 обзор ультразвуковой визуализации.)

Для типичного эксперимента EPIV частота кадров составляет от 20 до 60 кадров в секунду, в зависимости от условий потока, и мы получаем 100-1000 изображений в B-режиме пространственного распределения индикаторных частиц в потоке. После получения ультразвуковые изображения в режиме B передаются через соединение Ethernet на ПК, на котором установлено коммерческое программное обеспечение PIV. С помощью программного обеспечения PIV поля смещения частиц-индикаторов, D(x,y)[пиксели] (где x и y обозначают горизонтальное и вертикальное пространственное положение на ультразвуковом изображении соответственно) получают путем применения алгоритмов взаимной корреляции к последовательным ультразвуковым изображениям в B-режиме. ¹⁰ Поля скорости, u(x,y)[м/с], определяются по полям смещений, зная временной шаг между парами изображений, ΔT[s], и увеличение изображения, M[метр/пиксель], т.е., u(x,y) = MD(x,y)/ΔT. Временной шаг между изображениями ΔT = 1/fps + D(x,y)/B, где B[пикселей/с] — это время, необходимое ультразвуковому датчику для развертки по ширине изображения. В настоящем исследовании M = 77 [мкм/пиксель], fps = 49,5[1/с], а B = 25 047 [пикселей/с]. После получения данные о полях скоростей могут быть проанализированы для вычисления интересующих их объемов потока.

Protocol

1. Создать измеримый поток

1. Измерения EPIV проверки будет показано в трубопроводе раствора глицерина, воды (50% глицерина - 50% воды). Схема экспериментальной установки показана на рисунке 1.
2. Полые стеклянные сферы с номинальным диаметром 10 мкм добавляют в жидкость при концентрации примерно 17 весовых частей на миллион. Полые стеклянные сферы служить в качестве ультразвуковых контрастных средств, а также их размеров и плотности выбраны так, что они пассивно следуют за потоком жидкости ^10.
3. Фиксированное напряжение подается на насос ввести известно расхода. Расход выбрана такой, что U << Ах / t, где U ì средняя скорость в трубе, Ах является линейной длины объем EPIV измерения и δt является шагом по времени между изображениями, т. е. поток необходимости "медленные" по сравнению с кадров ультразвукового сишток ^3.

2. Калибровка Ультразвуковое

1. Установите ультразвуковой зонд к внешней стенке трубы. На водной основе актуальных гель наносится на ультразвуковой датчик, чтобы свести к минимуму потери передачи ультразвукового луча между торцом датчика и стенкой трубы.
2. Включите аппарат УЗИ. В прямом эфире поток ультразвуковых изображений начинается автоматически, как только все системы нагрузки.
3. Установите глубину изображения использовании Глубина ручку управления на панели управления ультразвукового аппарата.
4. Установить общий коэффициент усиления изображения с помощью 2D ручки Gain на панели управления ультразвукового аппарата.
5. Отрегулируйте время получить компенсацию (ТГК) ползунки для ослабления разброс от стенок труб и для компенсации глубины связанных затухания ультразвукового сигнала.
6. Ширина изображения, фокусировка, частота зонда операционной, и частота кадров настраиваются с помощью назначаемых ручки управления. Этичетыре ручки, расположенной в левом верхнем углу панели управления, меняется в зависимости от режима, в котором работает система. В 2D-режиме (как в настоящее время используется), слева-направо ручки соответствуют ширине, фокус, частоту и частоту кадров, соответственно. Обратите внимание, что в связи с основополагающими принципами ультразвукового ^9, этих четырех параметров по своей природе связаны. Следовательно, для данного сканирования изображений ультразвука (т. е. эксперимент EPIV) существует компромисс между пространственным и временным разрешением.
7. На рисунке 2 представителя изображений УЗИ трубе высевают с 10 мкм полые стеклянные сферы. Обратите внимание, что в связи с ограниченным пространственным разрешением, стеклянные шары размываются в поперечном направлении и появляются в виде эллипсоида в изображении.

3. Сбор данных

1. Нажмите кнопку Создать экзамен на панели управления ультразвуковым, чтобы начать новый эксперимент.
2. Создатьновая "пациента", введя в трубе Фамилия и сегодняшней даты в имя и номер теста у пациента ID.
3. После создания «пациент», ультразвуковое сканирование начинается до заданного максимума между 1000-1500 изображение будет достигнута, после чего новый цикл начинается сканирование. Нажатие на кнопку Стоп на панели управления дважды ультразвуковых будет перезапустить сканирование в любое время до достижения максимального заданного количества изображений.
4. После хорошего набора ультразвуковых изображений были приобретены (например, резкое семян изображения частиц и достаточной плотности частиц семян), нажмите кнопку Стоп на панели управления ультразвук, чтобы остановить получения изображений.
5. Нажмите кнопку кинопетли на панели управления ультразвука. Выберите набор ультразвуковых изображений быть проанализированы с помощью первого регулятора цикла на панели управления ультразвук, чтобы выбрать первое изображение в набор, и последний цикл ручку, чтобы выбратьПоследнее изображение в набор.
6. Нажмите кнопку Image Store на панели управления ультразвук, чтобы сохранить выбранный набор ультразвуковых изображений.
7. Нажмите кнопку архива на панели управления УЗИ и с помощью мыши курсор, чтобы выбрать конец экзамена. Это будет предлагать пользователю для выбора изображения или cineloops для сохранения на локальном жестком диске. Выберите кинопетли (ы), представляющие интерес затем выйти на экзамене.
8. Нажмите кнопку архива на панели управления ультразвука и использовать курсор мыши сначала выбрать все больше и затем выберите Управление дисками. Управление дисками будет передавать сохраненные кинопетли (ы) с ПК под управлением программного обеспечения PIV.

4. Преобразование Filetype

1. Ультразвуковое изображение хранится в виде цифровых коммуникаций изображений в медицине (DICOM) типа файлов на аппарат УЗИ. Для того, чтобы открыть и прочитать с помощью программного обеспечения PIV, DICOM файлы должны быть преобразованы в графические файлы. В настоящее времяСценарий Matlab работает DICOM2JPG.m используется для преобразования файлов в формате DICOM Joint Photographic Experts Group (JPEG) тип файла.
2. Изображений JPEG ультразвуковой затем анализируются с помощью программного обеспечения с Дэвис LaVision.

5. Вычисление поля смещения, D (х, у), используя Дэвиса

1. Двухместный мыши на значке Дэвиса на ПК. Выберите New Project. Выберите PIV.
2. Выберите Импорт изображения в панели инструментов и выберите пункт Импорт с помощью пронумерованных файлов. В выпадающем меню, найдите папку, в которой изображений JPEG ультразвуковой хранятся, и дважды щелкните на первое изображение из набора. Это будет импортировать все ультразвуковых изображений в этой номером набора.
3. Обычно изображение маски будут определены выделить область интереса (ROI) в ультразвуковом изображении должны быть обработаны. Для трубе, маска используется, чтобы определить ROI между стенок труб (например, жидкости).
4. К главной панели управления в Дэвисе, выберите вкладку расположенную под Текущий проект содержащие импортированные изображения и выберите вкладку пакетной обработки. Это позволяет окну векторной обработки Дэвиса для пакетной обработки импортированных изображений ультразвука.
5. Из списка операций, с помощью PIV-временных рядов деревьев, выберите вектора параметров расчета, а также выбрать параметры, которые будут использоваться для векторной обработки. Если маска используется, проверьте Range окно Data = использовать масках области в меню параметров расчета вектора. Обратите внимание, что оптимальный выбор параметров расчета вектора зависит от потока геометрии, свойств потока, разрешение изображения, индикатор плотности частиц, и желаемый количественный анализ потока ^10.
   Для измерения трубе, параметры, которые, как правило, дала лучшие результаты многопроходной с уменьшением допроса размер 32 х 32 пикселя ² до 8 х8 пиксельных ^2, с перекрытием 50%. Относительная ограниченность диапазона вектор был установлен на ± все (размер окна / 2) и абсолютное ограничение диапазона вектор был установлен в пределах ± 5 пикселей. Наконец, 3 х ² 3pixel медианного фильтра была использована для подавления шума и сглаживания векторных полей.
6. В левой части экрана, пакетная обработка, выделите общего количества изображений для обработки и выбрать начальную обработку. Это позволит вычислить смещение поля, D (х, у), между последовательными ультразвуковых изображений с использованием кросс-корреляцию алгоритмов.

6. Анализ векторных полей

1. Для анализа пост-обработки и передачи данных, поля EPIV вектор которые экспортируются из Дэвис. Файлы TXT. Это достигается путем выбора филиала вектор смещения под филиал JPEG изображений в проекте экране. В панели инструментов, выберите вкладку Экспорт, выберите тип файла ASCII. TXT, выберите / создайте папку экспорт,D Выберите экспорт.
2. Экспортировать векторные поля называются Bxxxxx.txt, где 00001 ≤ ≤ ххххх 99999, с B обозначает буфер. Каждый файл содержит четыре столбца данных: (1) X-положение вектора в образ, (2) у-положение вектора в изображение, (3) х-компонента смещения (т.е. продольного смещения), (4) Y-компоненту смещения (то есть, стены нормальные смещения). Bxxxxx.txt файлы открываются и обрабатываются в MATLAB, чтобы сначала вычислить поле скоростей, зная шаг по времени между изображением пары, ΔT [с], а увеличение изображения, M [м / пиксель], т. е. и (х, у ) = MD (х, у) / T,. Шаг по времени между изображениями T = 1/fps + D (х, у) / B, где B [пикселов / с] является время, необходимое для ультразвукового зонда прокатиться по всей ширине изображения. В настоящем исследовании, M = 77 [мкм / пиксель], FPS = 49,5 [1 / с], а B = 25047 [пикселов / с]. Далее, ансамблейBLE средняя скорость векторные поля, нормальной к стенке профилей средней скорости, среди других величин потока интерес вычисляются. (См. раздел Результаты представитель.)

Results

Мгновенным эхом скорости Изображения Частиц (EPIV) векторное поле показано на рисунке 3. Вектор график показывает векторы скорости каждого четвертого столбца, а цвет фона контурная карта соответствует скорости величина. Вектор среднее по ансамблю участок среднем более 1000 мгновенное вектор участков EPIV показано на рисунке 4. В соответствии с трубе, векторы скорости в основном в продольном направлении, крупнейший скоростей происходят в центральной трубе, и скорость уменьшится до нуля н...

Discussion

Операционная протокол Image Echo частиц велосиметрии (EPIV) система, способная приобретения двумерного поля скорости в оптически непрозрачных жидкостей или через оптически непрозрачных геометрии было описано. Практическое применение EPIV хорошо подходит для изучения промышленных и биологических проточных системах, где поток непрозрачные жидкости происходит в очень многих приложений. Особая система представленные здесь был построен специально для изучения реологических свойств сжиженных жидкости биомассы, используемой в ...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Название реагента	Компания	Номер в каталоге	Комментарии (по желанию)
Ультразвуковая машина	GE	Vivid 7 Pro
Линейный массив Ультразвуковое	GE	10 L
DC водяного насоса	KNF	NF 10 KPDC
Векторной обработки программного обеспечения	Lavision	Дэвис 7,2
Сообщение обработка программного обеспечения	Mathworks	MATLAB 7,12
Акриловые трубы	МакМастер-Carг	8486K531
Ультразвуковой гель	Паркер	Aquasonic 100