Summary
네 도면 스플리터 고속 카메라 기반 입체 입자 추적 속도계 (3D-PTV) 시스템은 여기에 설명된다. 기술 하류 ≈ 7000 다시 레이놀즈 수에서 열 지름 근방 원형 파이프에서 제트류에인가된다.
Abstract
3D-PTV는 이미지 시퀀스의 스테레오 촬영에 입체적으로 입자의 집합의 라그랑 경로를 추적 할 목적으로 정량 유량 측정 기술이다. 구성 요소, 기능, 제약 및 네보기 스플리터와 고속 카메라로 구성된 3D-PTV 토폴로지 최적화 팁 기본 설명이 문서에서 설명합니다. 이 기술은 다시 ≈ 7,000에서 원형 제트의 중간 유동장 (5 <x / d의 <25)에 적용된다. 오일러 프레임의 라그랑 플로우 기능과 난류 수량은 제트 기원의 다운 스트림 열 직경 주위와 제트 코어에서 다양한 반경 거리에서 추정된다. 라그랑 특성 궤적, 속도 및 가속도 선택된 입자뿐만 아니라 Serret-Frenet 식으로부터 얻어지는 흐름 경로의 곡률을 포함한다. 10시에 위치한 크로스면에서 제트 핵심 축을 중심으로 3D 속도와 난류 필드의 추정제트의 하류 직경은 문헌과 비교되고, 대규모 유선 속도 움직임의 파워 스펙트럼은 제트 코어로부터 다양한 방사상 거리에서 얻어진다.
Introduction
난류 제트 흐름은 엔지니어링 분야에서 유비쿼터스 있습니다. 이러한 흐름의 상세한 특성화는 전자 마이크로 스케일 디바이스에 대규모 환경 배출 시스템에서 실질적인 문제 걸치는 넓은 스펙트럼에서 매우 중요하다. 4 - 때문에 다양한 응용 프로그램의 수에 미치는 영향의 제트 흐름은 깊이 1에서 연구되고있다. 열선 풍속계 등 여러 실험 기법 4 - 8, 레이저 도플러 속도계 (LDV) (4), 9-12 및 입자상 유속계 (PIV) 12-16은 제트 레이놀즈 수와 경계 다양한 흐름 특성화하는데 사용되어왔다 정황. 최근 몇몇 연구는 제트 (17) 흐름의 18 난류 / 비 난류 인터페이스를 연구하기 위해 3D-PTV를 사용되었습니다. 3D-PTV는 복잡한 난류 파이를 설명하는 데 적합하는 기술이다다른 관점에서 필드들. 그것은 멀티 뷰 입체를 사용하여 참조의 라그랑 프레임에 볼륨 내의 입자 궤적의 재구성을 할 수 있습니다. 첫번째 기술은 창 (19)에 의해 도입 된 추가 Racca 듀이 (20)에 의해 개발되었다. 24 - 그 이후로 많은 개선이 3D-PTV 알고리즘과 실험 장치 (21)에 이루어졌다. 이러한 성과 이전 작품으로, 시스템이 성공적으로 4 MX 2 MX 2m (25), 실내 공기 흐름 필드 (26)의 영역에서 대규모의 유체 운동 등 다양한 유체 현상을 연구하는 데 사용되었습니다, 박동성 (27)과 대동맥 혈액 흐름 (28) 흐름 .
3D-PTV 측정의 작동 원리는 데이터 수집 시스템 설정, 기록 / 재생의 전처리, 교정, 3D 대응, 시간 추적 및 후 처리로 구성된다. 정확한 보정은 입자 위치의 정확한 검출이 가능에스. 세 개 이상의 이미지보기 검출 입자의 대응 에피 폴라 기하에 기초하여 3 차원 입자 위치의 재구성을 허용한다. 연속적인 이미지 프레임의 결합은, 입자 궤도의 S (t)를 정의 시간 추적 결과. 3D-PTV 시스템의 최적화 다중 입자 추적 가능성을 최대화하기 위해 중요하다.
최적화의 첫 번째 단계는 고속 카메라, 조명 원 시드 입자의 특징을 포함 적절한 데이터 수집 시스템을 획득하는 것이다. 질의 볼륨의 크기와 함께 카메라 해상도가 있으므로, 화소 크기를 정의하고, 단일 픽셀보다 커야 필요한 시드 입자 크기. 검출 된 입자의 무게 중심은 명도 (21)에 의해 가중 된 입자 화소의 평균 위치를 취함으로써 서브 픽셀 정확도로 추정된다. 카메라의 프레임 비율이 밀접되어 associat레이놀즈 수 검출 입자를 연결할 수있는 능력 에디션. 높은 프레임 레이트는 이미지 사이의 평균 변위가 입자의 평균 분리를 넘으면 트랙킹이 어려워지기 때문에 빠른 흐름 또는 입자들보다 많은 수의 해결을 허용한다.
셔터 속도, 조리개 및 감도는 이미지 캡처에서 고려해야 할 세 가지 요소입니다. 셔터 속도가 빠른 충분해야 입자 중심 위치의 불확실성을 감소시키는 입자 주위 흐리게 최소화합니다. 카메라 조리개 볼륨 외부 입자를 검출 할 확률을 감소시키는 질의 볼륨의 심도로 조정되어야한다. 카메라의 최대 감도가 고정되어 있기 때문에, 프레임 레이트가 증가함에 따라, 필요에 따라 빛이 증가한다 입자를 조명하는 데 필요한. PIV 달리 복잡한 광학 설정 및 고출력 레이저는 엄격하게 한 광원이 충분히 배설물 같이, 3D-PTV 필요하지카메라에 추적 입자에서 거행 한. 연속 LED 또는 할로겐 조명 동기화 (21)의 필요성을 우회 좋은 비용 효율적인 옵션입니다.
3D-PTV에 다른 광학 유동 측정 기술처럼 추적 입자 속도는 로컬 순간 유속 29으로 가정한다. 그러나,이 널 직경 관성 이상적인 트레이서 경우만이다 추적 입자는 충분히 큰 카메라로 포착 할 수 있어야합니다. 유한 입자의 충실도 스톡스 번호 S t에 의해 결정될 수 있고, 입자의 완화 시간 스케일 비율 관심 난류 구조의 시간 척도 즉. 일반적으로, S t는 S의 t은 1보다 훨씬 작아야한다 ≤0.1 흐름 추적 오차는 1 % 30 이하. 깊이있는 논의는 메이 등의 알에서 찾아 볼 수있다 29 -. 31 (예를 들면 50 ~ 200 μm의) 작은 입자 (32) 반면, (예를 들어 1-50 μm의) 33, 34 (고출력 레이저를 사용할 수있다 예를 들어 80-100 와트 사용 CW 레이저). 특정 파장의 빛에 대한 높은 반사율을 가진 입자는, 할로겐 조명 아래 코팅은 같은 이미지로 자신의 마크를 증폭 할 수있다. 파종 밀도 성공적인 3D-PTV 측정하는 또 다른 중요한 변수이다. 입자의 과도한 수의 대응을 수립하고 추적에서 모호성을 야기하면서 몇 가지 입자는 궤도의 낮은 숫자 초래한다. 대응을 수립 모호성은 중복 입자와 정의 에피 폴라 라인을 따라 여러 후보를 검출을 포함한다. 추적 과정에서의 모호성 때문에 높은 seedin 행 g 밀도 때문에 입자의 평균 비교적 짧은 간격으로 발생한다.
두 번째 단계는 이미지 품질을 향상시키기 위해 기록 / 전처리에 최적의 설정이다. 이러한 흑색 레벨 이득 (G 및 B)와 같은 촬영 설정, 이미지 품질을 최적화하는 중요한 역할을한다. 이득은 화상의 밝기를 증폭하는 반면 흑색 레벨은 화상의 가장 어두운 부분에서 휘도 레벨을 정의한다. G & B 수준의 약간의 변화가 현저하게 추적 가능성에 영향을 미칠 수 있습니다. 사실, 높은 G & B 이미지를 과도하게 밝게 결국 카메라 센서가 손상 될 수 있습니다. 이를 설명하기 위해, 흐름 재건에 G & B 수준의 영향은이 문서에서 검사된다. 전처리 단계에서, 이미지로부터의 광 분산 입자를 강조하는 고역 통과 필터로 필터링된다. 픽셀 크기는 계조 질의 볼륨 내의 입자 검출을 최대화하도록 조정된다.
t "> 최적화의 세 번째 단계는 에피 폴라 기하학, 카메라 파라미터 (초점 거리 원칙적 가리킨 왜곡 계수) 및 굴절율의 변화에 기반 입체 영상의 정확한 보정된다.이 과정은 3 차원을 최소화하기 위해 중요 기점 목표 지점의 재구성 에러. 에피 폴라 기하 타겟 이미지로부터 경사각 (카메라 질의 볼륨 사이의) 상대 거리를 사용한다. 질의 볼륨을 통해 카메라 뷰를 따라 굴절률 변화가 질량 절차에 기초하여 고려 될 수있다 등. (21).이 실험에서, 규칙적으로 목표 지점으로 3D 계단 형 구조를 대상으로 사용된다.두 이미지는 3D 입자 위치를 결정하는 데 필요한 있지만 3D-PTV 실험에서 전형적 이상의 카메라는 21 모호성을 줄이기 위해 사용된다. 여러 고속 카메라와 비싼 설정에 대한 대안은 VI입니다으으 3D-PTV의 사용에 대한 호이어 등의 알에 의해 제안 된 스플리터. (35) 최근 생물 의학 응용 프로그램에 대한 Gulean 등. (28)에 의해인가. 보기 스플리터는 피라미드 모양의 거울 (이에 관해서 차 거울)과 네 개의 조정 가능한 미러 (보조 미러 이에 관해서)로 구성되어 있습니다. 본 연구에서는 네보기 스플리터와 하나의 카메라는 4 대의 카메라에서 입체 영상을 모방하는 데 사용되었다. 시스템은 직경 (D)의 H = 1cm 및 R6 ≈ 7000 라그랑에서 함께 토출 관의 중간 유동장의 특성을 사용하고, 오일러 제트 기원의 하류 직경 약 14.5-18.5에서 프레임.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. 실험실 안전
- 선택된 조명 원의 안전 가이드 라인 (예를 들어 레이저, LED 산업, 할로겐)을 검토한다.
주 :이 실험에서, 다섯 250 와트 할로겐 스포트라이트 조명 세트로 사용된다. 다음과 같이 광원의 기본 안전 및 추천 측면이 설명되어 있습니다.- 높은 온도 (~ 3000 K 색 온도)에서 작동 할로겐 조명과 직접 접촉을 피하십시오.
- 고려 흐름을 가열 방지하기 위해 데이터를 수집하는 경우에만 ON 빛을 유지합니다.
- 어떤 종류의 용지를 포함하여, 광원 근처에 떨어져 모든 가연성 물질을 보관하십시오.
2. 실험 셋업
- 적절한 렌즈를 선택
- 교정 문제를 방지하기 위해 낮은 수차 렌즈를 선택합니다. 추천 렌즈는 망원 타입 또는 마이크로 렌즈이다.
- 렌즈가 OBJ 내에서 원하는 시야 (FOV)을 커버 있는지 확인요법 거리, O, 필요한 확대, M을 추정하여.
주 : 배율 FOV 카메라 칩의 길이의 비이며, 물체 거리 f는 렌즈의 초점 거리 O = F (1 / M + 1),로서 계산 될 수있다. 이 실험에서는, 카메라 칩의 길이 20.34 mm이고, 대응 FOV 또는 주경이 제한된 목적 거리 O ≤250 mm 50 mm이다. (피사체 거리 인해 카메라 뷰 스플리터 장착되는 슬라이더의 한정된 길이로 억제된다.) 배율은 M = 20.34 / 50 = 0.41, 및 피사체 거리의 소정 범위의 대략적인 초점 길이 f를이다 ≤72.7 mm. 따라서, 60mm의 초점 거리의 마이크로 렌즈는 F / 2.8D의 초점 비율로 사용된다.
- 마운트보기 분배기에서 카메라를 조정합니다.
- 미르을 밀어 심문 볼륨의 그와 기본 거울의 중심을 수평수직 장착 후 함께 ROR 및 포스트 홀더 미러를 고정. 이 단계는 보조 미러를 설치하기 전에 수행합니다.
- 카메라를 탑재하고 기본 거울의 중심과 이미지 일치의 중심을 설정합니다.
- 단지 기록 소프트웨어 설정이자 (ROI)의 영역을 제어함으로써, 주 미러를 커버하도록 카메라 뷰의 폭과 높이를 조정한다. 이 프로세스는 이미지 사이즈와 소음을 감소시키는 것이다. 주 :이 실험에서, 주 미러의 크기는 5 × 5cm 2 (1728 X 1728 픽셀)이다.
- (정의) 3D 교정 대상을 디자인합니다. 그것은 전체 조사 볼륨을 묶어야한다. 스플리터의 각보기가 균일 한 교정을 할 수 있도록 모든 대상 마크를 캡처 있는지 확인합니다.
참고 :이 데모에서, 대상 베로 다시 플라스틱을 사용하여 3D 프린터로했다. 은 1 mm 직경 화이트 목표 지점의 해부와 사이즈 35 X 35 X 30mm 3 계단 형상을 갖는다2.5 mm, 5mm, 수직, 유선과 스팬 방향 10mm 에드. 이 수조에 대하여 캘리브레이션 모델 및 카메라의 위치로 돌출로 타겟의 정확한 구조는 중요하다.
- 심문 볼륨으로 교정 대상을 놓습니다.
- 카메라를 향해 직면 높이 조절 플랫폼에서 보정 대상을 놓습니다.
- 타겟 플랫폼의 높이를 조정함으로써 질의 볼륨의 중심과 보정 대상의 중심 높이와 일치.
주 :이 예에서, 보정 대상의 중심 마크는 분출 노즐의 중심 20cm 높이로 평평하게된다. 버블 레벨 게이지 목표 수준으로 사용될 수있다.
- 마운트 및 네보기 분배기의 보조 미러를 조정합니다.
- 그 전체 캡처를 보장 심문으로부터의 거리에서 차 거울을 찾습니다. 그것은 (Figur이 데모 0.2 m이다E 1).
- 기본 미러의 각면에서 카메라 뷰가 약 보조 미러와 정렬 대략적인 위치에서 보조 미러를 탑재합니다. 보조 미러의 수직 장착 포스트에 고정하여 보조 미러를 고정합니다.
- 다른 세 개의 거울이 과정을 반복합니다. 차 거울에 대한 모든 보조 거울의 기하학적 대칭을 확인합니다.
- 네 개의 각보기는 전체 보정 대상을 둘러싸 수 있도록 보조 미러의 미러 마운트를 조정하여 최종 조정을합니다. 거울 '위치와 각도를 확인하는 하나의 효과적인 방법은 각각의 뷰의 이미지 경로를 시각화 레이저 포인터를 사용하는 것이다.
- 하나의 미러를 이동하여 서브 이미지의 중첩을 확인합니다.
참고 : 하나의 뷰가 변경되면, 중복 영역을 무시할 수있다. 그렇지 않으면, 반복 중복 영역이 최소화 될 때까지 2.4.5에 2.4.2 단계를 반복합니다.
- (광원을 배치(S)에 직접 질의 볼륨 접함). 카메라 센서에 손상을 방지하기 위해 빛을 조정하면 카메라가 커버 한해야합니다.
- 광원이 균일하게 전체 조사 볼륨에 분산되어 있는지 확인합니다.
- 필요한 경우, 광원 바로 아래 확대 렌즈를 배치함으로써, 광 강도를 향상시킨다. 참고 :이 실험에서, 초점 거리 f를 0 = 450mm의 평 볼록 확대 렌즈가 조명을 강화하는 데 사용됩니다.
3. 업 최적화
- 켜기 및 이미지 품질을 높이기 위해 카메라에 설정을 조정합니다.
- 기본 거울을 통해 반사가 똑같이 보조 미러의 네보기에 초점을 맞추고 될 때까지 렌즈의 배율을 조정합니다.
- 보기 스플리터에서 이미지가 대칭인지 확인 4 개의보기에서 보정 대상 이미지의 대칭성을 관찰하여 심문 볼륨을 캡처.
- 카메라로부터 가장 가깝고 먼 교정 목표 지점을 포착 할 F 값을 조정한다.
주 :이 카메라에만 질의 볼륨의 깊이 추적 입자를 포착 할 수있다. 이러한 예에서, F 수치는 11이다. - 550 Hz에서 (이 소개를 참조하십시오, 특정 응용 프로그램에 따라 다릅니다) 등의 원하는 프레임 속도를 설정하고 레코딩 소프트웨어에 따라 빛 감도를 극대화 할 수 있습니다.
- 라이브 카메라 뷰를 통해 스플리터의 각보기의 입자 밀도의 차이를 관찰하여 기본 거울의 각각의 관점에서 조명을 확인합니다.
주 : 여러 광원 질의 볼륨을 조명하는 데 사용되는 경우에는 분할 각 뷰의 차이를 가질 가능성이 높다. 조명 상부로부터 오기 때문에 본 실험에서는 상위 두 보조 미러 적은 광 받았다. 수조의 바닥에 평면 미러의 이용 광 변동을 줄이는 데 도움보기에서. - 3D-PTV 광원을 사용하기 전에 방에 배경 조명을 끄십시오.
- 더 나은 입자의 빛 산란을 캡처하는 카메라의 G & B 레벨을 조정합니다. 다양한 G & B 수준의 몇 가지 짧은 시퀀스를 기록하고 입자 궤적의 분포와 밀도를 관찰하여 최적의 하나를 찾을 수 있습니다.
주 :이 실험에서, G 및 B 레벨의 범위는 0-500이고, 블랙 (B) 레벨이 적당히 이미지를 증폭하는 반면, 게인 (G)이 설정되어 매체 (300) 디머 광 산란을 밝게하기 위해 500으로 설정했다 신호 않도록 이미지를 오버 밝게.
4. 교정
- 추적 입자를 추가하기 전에 조사 볼륨 보정 대상을 놓고 약간의 보정 이미지를 가져 가라. 대상을 조명 디머 광원 (예 : LED 플래시 라이트)를 사용합니다.
- 네 개의 독립 서브 이미지로 캘리브레이션 화상 나누고 참조 포함한 텍스트 파일을줬어 대상 마크의 위치를 좌표입니다. OpenPTV 소프트웨어 (http://www.openptv.net)이 목적을 위해 여기에 사용됩니다.
주 : 이후 처리는 다중 카메라 셋업을 이용하는 사용자와 동일하다. - 이미지 및 소프트웨어의 '칼'폴더의 단계 4.2에서 얻어진 텍스트 파일을 저장 한 후 교정 프로세스를 시작하는 '교정 만들기'탭을 클릭합니다.
- '편집 교정 매개 변수'탭을 클릭하고 각 분할 뷰의 중심과 보정 대상의 기원 사이에 확대, 회전 각도와 거리를 정의하기 위해 '교정 방향 매개 변수'탭을 선택합니다.
주 : 첫 번째 행은 X, Y, Z 방향으로 카메라 센서로 원점 보정 대상까지의 거리이다. 둘째 행은 X, Y, Z 축 주위 라디안 각도를 나타냅니다. 다음으로, 3 × 3의 데이터는 회전 행렬을 나타낸다. 그리고, 다음 두 행 (X)의 거리가 핀홀이며 - '추측'포인트가 검출 된 목표 지점과 일치되어 있는지 조사하기 위해 '검색'과 '쇼 초기 추측'을 클릭합니다.
- '추측'포인트까지 네 개의 뷰에 대한 단계를 반복 4.4 보정 이미지의 집합으로 정렬됩니다.
- 심문 볼륨의 방향을 재구성하는 '방향'을 클릭합니다.
주 : 보정 렌즈 왜곡과 아핀 변환을 조절함으로써 개선 될 수있다. 이제 조사 부피 보정 데이터를 처리 할 준비가된다. 자세한 설명은 저자의 논문 36 교정 과정을 참조하십시오.
5. 흐름 설정 / 데이터 수집
- 프레임 레이트에서 각 카메라 뷰 캡처 입자의 최대 양을 추정D 최대 유속. 이 예제에서, 기준 속도는 U ≈ 0.4 m / 초이며, 프레임 레이트가 550 Hz에서 ~ 4 × 4 × 4cm 3 질의 볼륨이다. 이 ~ 프레임 당 1,000 입자 발생했습니다.
- 3 단계에서 얻어진 최적의 설정으로 카메라를 켭니다.
- 시드 입자를 첨가하고, 정상 상태에 도달하는 흐름을 허용하는 여러 가지의 평균 체류 시간을 기다린다. 필요한 경우 더 많은 입자를 추가하지만 모호성이 발생할 수 있습니다 단계 5.1에서 추정 높은 파종 밀도를 피할 수 있습니다.
참고 :이 예에서, ~ 1.1 g / ㎤의 밀도의 100m 실버 도금 중공 세라믹 분야의 1.6 g의 유체 매체 (2 × 0.4 × 0.4 m 3)에 대한 시드로 사용된다. - 유동 이미지의 수를 기록한다.
주 :이 실험에서, 550 Hz에서 9000 사진 기록 소프트웨어를 사용하여 촬영 하였다. 카메라 및 / 또는보기 스플리터가 (심지어 약간의 모션이 많이 결과에 영향을 미칠 수 있습니다) 이동하면 반복 5.3 2.4 단계를 반복합니다.
6. 데이터 처리 (비아 OpenPTV 소프트웨어)
- 네 개의 독립적 인 서브 이미지로 단계 5.4에서 얻어진 원료 이미지를 나눈다.
- 네 개의보기에서 초기 이미지를로드하기 위해 '시작'탭에서 '초기화 / 다시 시작'을 클릭합니다.
- '실행'디렉토리를 마우스 오른쪽 버튼으로 클릭하고 카메라, 굴절률, 입자 인식, 대응을위한 시퀀스 이미지, 관찰 볼륨과 기준의 수의 '기본 매개 변수'제어를 참조하십시오.
- '일반'탭의 실험에 사용되는 카메라 (보기)의 수를 정의합니다. 이 실험에서는,도 4과 같이 카메라의 수를 설정한다.
- '굴절 지수'탭에서 카메라 뷰를 따라 굴절률을 정의합니다.
- '입자 인식'탭 아래 네 개의 뷰에서 입자 검출의 수를 최적화하기 위해 최소 및 화소 검출뿐만 아니라 회색 값 임계 값의 최대 수를 정의합니다. 픽셀 드의 최소 및 최대 수tection 회색 임계 입자 검출 화소 사이즈와 밝기 레벨을 결정한다. 그것은 초점이 소음과 입자를 제거합니다.
- '시퀀스 처리 파라미터'아래 과정 이미지 수를 정의한다.
- '관측 볼륨'탭에서 관찰 볼륨을 정의합니다.
- 스테레오 매칭을위한 총 밴드 매개 변수 (mm)를 포함하여 '대응 기준'에 따라 대응의 상관 관계를 정의합니다.
- '미리 처리'탭에서 '하이 패스 필터'를 클릭하십시오. 이 네 개의 뷰에서 입자 빛의 산란을 강화.
- 네 개의 뷰의 서브 픽셀 레벨에서 검출 된 입자의 중심을 결정하는 '입자 검색'을 클릭합니다. 반복 단계 5.1에서 계산 된 입자의 예상 번호와 유사한 검출 입자의 수까지 6.2 및 6.3 단계를 반복합니다.
- 입체 corresponden을 설정하는 '대응'을 클릭각 뷰에서 CES.
주 : 검출 입자의 3 차원 위치를 재구성하도록, 대응 세 전망으로부터 적어도 결정되어야한다. - 교정에 기초하여 검출 입자의 3 차원 위치를 얻기 위해 '3D 포지션'을 클릭합니다.
- 단계의 모든 이미지 시퀀스에 대한 6.4-6.7의 과정을 반복 '을 표시하지 않고 순서'를 클릭하십시오.
참고 :이 탭으로 구분 된 텍스트 파일 형식으로 프레임에서 검출 된 입자의 요약을 포함하는 각 이미지 세트에 대한 'rt_is'파일을 만듭니다. - '실행'디렉토리를 마우스 오른쪽 버튼으로 클릭하고 추적 후보 입자를 검색하려면 (mm / 프레임 예를 들어 dvxmin 및 dvymin) 반경 영역의 매개 변수를 정의하는 '추적 매개 변수'를 클릭합니다.
- 단계 6.7에서 얻어진 복원 된 입자의 입자 식별 (ID)를 정의하기 위해 '표시하지 않고 추적'을 클릭합니다.
주 : 4 컷 predic를 사용 추적 인접 프레임의 시퀀스를 상관토르 예측 - 보정 방식 24. 이 프로세스는 프레임에서 검출 된 입자의 추적 정보를 포함하는 각 이미지 세트에 대한 ptv_is 파일을 생성; 첫 번째 두 열은 각각 이전 프레임과 다음 프레임에서 입자 ID를 나타낸다. - 각 카메라 뷰에서 궤도를 시각화 '쇼 궤적'을 클릭합니다.
7. 후 처리 (희구 법)
참고 : 범위와 후 처리의 유형은 개인의 필요에 따라 그것을 따라서, 사용자 정의 할 수있다. 여기서,베이스 포인트 계산 간단히 예를 들어 설명한다.
- (matlab에를 통해) 라그랑 프레임의 데이터를 얻습니다.
- ptv_is 파일에서 각 입자 및 관련 ID의 3D 위치에 압축을 풉니 다. 그것은 궤적을 재구성 이미지 시퀀스 중 발견 입자를 연결 할 수 있습니다.
- 각 궤적에 대해, 지정된 프레임 속도에서 속도와 입자의 가속을 계산합니다.이 예제에서, 입자의 속도 및 가속도는 이동 큐빅 스플라인 (34), (37)로 위치 신호를 저역 통과 필터링함으로써 계산된다.
- 3D 위치, 3D 속도, 3D 가속 및 시간 스탬프뿐만 아니라 각 궤도의 궤도 ID를 포함하는 필드 구조체 배열 형식을 확인합니다. 이 데이터 형식, 구조체 배열의 길이는 궤적의 수를 나타낸다.
- (matlab에를 통해) 오일러 프레임의 데이터를 얻습니다.
- 각 입자의 타임 스탬프를 이용하여 시간적 하나 구조체 배열한다 (단계 7.1.3) 변환. 이 단계 7.1.3에서 얻은 유사한 구조체 배열 구조를 생성하지만, 구조체 배열의 길이는 현재 실험에서 9,000 프레임 번호를 나타낸다.
- 오일러 좌표 순간 속도 필드를 얻기 위해 각각의 시간 프레임에 대한 삼차원 격자로 배열 구조체 시간적 보간. 이 데모에서, griddata 재미매트랩 ction 특 보간을 수행하는데 사용된다.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
사진 및 설치의 개략도는도 1 및도 2에 나타내었다. 보정 대상은 기점 마크가도 3에 도시 된 뷰 스플리터 및 3D 재구성 보정에 반영. 인식 교정 대상의 RMS는 유선 방향 X, 스팬 방향 Y, 및 Z 방향의 깊이 7.3 ㎛의 5.7 ㎛ 내지 1,417 ㎛의 것이다. Z의 -coordinate에 대하여 높은 RMS가 다른 방향 X 및 Y 좌표를 비교하여 Z 축 네 개의 뷰 비교적 작은 각도에서 이들에 대한 감소 대상 점에 기인한다. 임의의 주어진 순간에 네 뷰 각각 검출 된 입자는 103의 순서였다. 검출 된 입자 중에서 성공적인 3D 재구성 수가 사실로 인해 절반으로 감소된다 교차 영역 AR에만 입자전자 캡쳐. VIDEO 1은 네 시청 스플리터로부터 제트류의 고속 비디오 샘플을 나타낸다.
중간 필드 영역에서 네 대표 입자 궤적 샘플 주위 반경 거리에서 X / D H = (16) 평면과 교차는 R / D 시간 = V0을, 제트 코어 1.5, (3)은도 4에 도시되어있다. 예상대로, 주어진 시간 간격 (의 Δt ≈ 1 초)의 긴 궤적은 제트 코어에 관찰된다. 제트 (R / D h를 ≥2)의 가장자리에, 추적 입자. 짧고 더 복잡한 궤적을 나타내는도 5를 X / d의 시간 = 16 평면을 넘어 모든 성공적으로 복원 된 입자의 궤적을 그림. 선택된 영역에서 입자의 속도는 약 0-0.6 U J와 U J에 이르기까지 넓은 분포를 나타낸다 시간 = 16 평면을 교차하는 입자의 경우 6 그림. 6 B, 6 C도 및도 6은 정규 시간의 함수로서 입자 궤도, 속도, 가속도의 3 요소를 표시 거라고. 로컬 입자 가속 수회 표준 중력 수 있다는 강조 가치가있다. 입자 궤적 소위 천을 통해 입자 궤도의 특정 기능을 얻기 위해 허용ENET-Serret 프레임입니다. 그것은의 (t)를 따라 직교 벡터 (접선, 정상, 법선)의 변경 사항에 대해 설명합니다. 특히 관련성 ρ, 곡률 반경의 역수이며, 다음과 같이 정의되는, κ, 곡률 :
어디에 
= DR / DS는 궤적의 접선 단위 벡터이며, R은의 함수로서 기록 될 수있는 시간의 함수, 즉 R (S) = R (t 같은 입자 (유클리드 공간)의 위치 벡터이다 (에스)). 곡률은, κ, x/d의 시간 = 16 x/d의 시간이 = 17 평면을 넘어 모든 입자에 대해 계산된다. 평균 곡률 
제트 코어 (R)으로부터의 거리의 함수로서 다음과 같이 계산된다 : 
ΔR =는 0.2㎜의 시간이 여기에 사용된다. 7 설명도 일 시간에 의해 정규화 = F (R). 비교적 낮은 거의 일정한를 도시 관의 단면이 원형, R / D의 시간으로 정의 ≤0.5 지역. X / 일 시간 = 16 평면에서의 제트 코어에서 더 큰 거리에서, 단조 증가한다. 유사한 경향 X / D = 17의 H면에서 획득하지만 함께 감소된다 제트 코어 (R / D의 시간 ≥0.5) 외부. 그것은 강조 가치가이 흐름 FEA진짜야은 3D-PTV 기술로 추론 할 수있다. G 및 B의 설정의 다양한 레벨에 기초하여 데이터의 품질은 표 1에 나타낸 3 차원 재구성 입자의 나머지 부분에 결합 입자의 비율의 관점에서 평가된다. 높은 링크 비율 300 500 G 및 B의 설정에서 관찰된다.
오일러 유동 특성 차원 입자 화상 속도계 (3D-PIV)를 본뜬 격자 보간에 의해 달성 될 수있다. 이 프레임의 상당히 높은 수는 오일러 설명을 진정으로 모방 PIV 품질에 필요 때문에 각각의 시간에 추적 비교적 낮은 입자 있다는 점 유의하는 것이 중요합니다. 이는 고차 통계 (예를 들어, 난류 강도 레이놀즈 응력)의 추정에서 더욱 중요하다. 다양한 G 및 B 레벨 제트 핵심 유선 속도는도 8에 도시되어있다. 측정은 theoret 비교iCal의 행동 :
U 0 (X)이 제트 핵심 유선 속도이고, B ≈ 6 정수이고, X는 0 가상 원점 (38)이다. 도는 G 및 B 레벨 설정의 관련성을 나타낸다.도 9는 X / D = 16의 H면에 분사의 평균 속도 분포를 도시한다.
마지막 위치에서 유선 속도의 대형 동작의 분광 분포 φ (F)를, h는 = X / D H = (10) 평면에서 0, 0.6, 1이도 10에 도시된다 / D R. 버터 워스 저역 통과 필터는 속도 시계열에 도포차단 주파수, f를 C = 200 Hz에서와.
그림 1 :. 실험 장치의 도식 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 2 :. (왼쪽 위) 상위 뷰 (왼쪽 아래) 카메라와 뷰 분할 시스템의 다시보기 (: 실험 셋업이 다양한 카메라 뷰와 네 개의 이미지보기 스플리터, 수로 및 심문 볼륨을 보여줍니다 상단 중간, 하단 가운데) 전체 실험 장치의 측면도, (오른쪽) 줌 - 제트 흐름의 시드 입자의보기. <A HREF = "https://www.jove.com/files/ftp_upload/53745/53745fig2large.jpg"대상 = "_ 빈">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 3 : 교정 :. () 교정 대상,보기 스플리터에서 보정 대상의 (b)의 이미지 세트, (c) 상기 보정 대상의 기준 마크의 3D 인식 이의 더 큰 버전을 보려면 여기를 클릭하십시오 그림.
그림 4 : r/d의 시간 = 0에서 선정 된 입자의 궤적, 1. 5, 3. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 5 :. 속도는 컬러 레벨로 표시되는 X / 일 시간 = 16 비행기를 건너 그림 심문 볼륨을 입자 궤적이 (x)를 /d의 시간 사이에 포함 
(14.5,18.5), y/d의 시간 (-2,2), 및 h를 z/d.JPG "/> (- 2,2), 여기서 (X, Y, z) = (0, 0, 0) 제트 원점의 중심에 위치한 각각의 궤도를 따라 속도, 부피 유속에 의해 정규화. U 0, 컬러 레벨로 설명된다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 6 :. () 입자의 궤적, (b)의 변위, (c)는 속도, 임의의 입자의 (d)에 가속 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 8 : (x)를 /d의 시간 내에서 제트 핵심 유선 속도 다양한 G & B 레벨 (15, 18). 세 G & B 수준 (300 500 (최적), 300 및 250, 100 및 250)에 포함되어 있습니다.745fig8large.jpg "대상 ="_ 빈 ">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 9 :. X / 일 시간 = 16에서 유선 속도 성분의 비 차원 분포 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
도 10 : R / D에 H에 위치한 지점에서 유선 속도 성분의 파워 스펙트럼 φ (F) = 0 X / d의 시간에 (제트 코어), 0. 6, 1 = 16면. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
비디오 1 : 네보기 스플리터에서 제트 흐름의 비디오 샘플, 550 fps로 얻은 실제 속도보다 10 배 느린 ( 오른쪽 다운로드 클릭 ).
티1 수 :. 다양한 G & B 수준에서 3 D-재구성 입자의 나머지 부분에 링크 된 입자의 비율이 세 가지 G & B 수준 (100 250, 300 250, 300 및 500)에 포함되어 있습니다.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
3D-PTV는 많은 다른 사람의 사이에서 대동맥 토폴로지 (28) 등의 낮은 대기 (25), 실내 공기 분배 (26), 또는 타악기 흐름에 큰 규모의 난류 운동 난류 흐름의 다양한 복잡한 물리학을 해명하기 위해 큰 잠재력을 가지고있다. 그러나, 장점 및 제한뿐만 아니라 경험 이해 가능성을 최대화하기 위해 중요하다. 시행 착오 예비 시험 및 프레임 레이트, 조명 원, G 및 B 레벨 및 이미지 필터링 방법을 포함하여 최적의 설정을위한 철저한 반복은 직접 (예를 들어, 추적) 입자의 집합의 라그랑 경로를 재구성하는 능력과 상관된다. 여기에 설명 된대로 중요한 프로토콜 단계, G & B 레벨의 조정 및 FOV의 조명 (할로겐 스팟 조명, 돋보기 렌즈의 조합 수로의 바닥에서 반사 거울) 참고하는 것이 중요하다.
ENT "> 이러한 조정은 네 개의 뷰에 대한 조사에서의 광 산란을 최적화하는 것을 돕는다. 고성능 측정 실험 설정을 식별 한 후, 철저한 수정 및 문제는 프레임 레이트에 기초하여 정확한 궤도의 최대 값을 계산하기 위해 만들어 져야한다 포획 입자의 개수가 더 높은 프레임 속도로 증가 될 수 있지만, 카메라의 해상도 및 조사 볼륨의 크기. 그것은 3D-PTV에서 추적 된 입자의 개수가 PIV에 비해 훨씬 낮다는 것을 알아 차리고 가치가있다. 3D의 큰 잠재적 -PTV 여러 입자 라그랑 경로를 기술의 독특한 능력이다.이 예제에서보기 스플리터 설정 여러 넓은 카메라를 사용하지 않도록 구현하였으나,이 셋업 높은 카메라를 필요로한다는 것이 중요하다 해상도는 샘플 볼륨의 크기를 제한한다.본 연구에서는 원형 제트의 중간 필드 기능은 항문이다3D-PTV 기술로 yzed. 방법은 오일러 및 라그랑 프레임으로부터 유동 중요한 기능을 얻는시켰다. 특히, 반경 방향 거리의 함수로서 입자의 평균 곡률 입자 궤적 라그랑 기능을 사용하여 두개의 단면 평면에서 처음 특징으로한다. 인식 교정 대상의 RMS는 유선과 스팬 방향 141.7 ㎛, 7.3 ㎛의 사이의 범위이다. 인해 z 방향에서보기 작은 각도로 스팬 방향이 높은 상대 오차 극복 완전히되지 않을 수 있지만, 추가로 각종 위치에서 2D 캘리브레이션 타겟을 사용으로 z 방향에서보다 목표 지점을 추가로 감소시킬 수있다 (멀티 플레인 교정).
전반적으로, 3D-PTV는 시간에 따른 플로우 또는 활성 스칼라들의 역학을 포함하는 다른 많은 문제에 적용 할 수있는 유용한 기술이다. 예를 들어, STU 것이 매우 유용수생 환경에서의 난류과 종 사이의 상호 작용을 DY.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
저자는 더 경쟁 재정적 관심이 없음을 선언합니다.
Acknowledgments
이 작품은 레오나르도 P. 차모로어의 스타트 업 패키지의 일부로 기계 공학, 어 바나 - 샴페인 일리노이 대학,학과에 의해 지원되었다.
Materials
| Name | Company | Catalog Number | Comments |
| Mikrotron 4CXP 4 lanes CXP-6 CoaXPress | ImageOps | CAMMC4082 | High-speed camera |
| Active Silicon FireBird CoaX Frame Grabber | ImageOps | FBD-4XCXP6 | Frame Grabber |
| 100 μm silver-coated hollow ceramic spheres | Potters Industries LLC | AG-SL150-30-TRD | Seeding Paritcles |
| StreamPix6 | Upstate Technical Equipment CO.,INC | MISNOR-STP-6-S-CL | Camera appliation |
| Four-view splitter | Photrack AG | Customized part and necessary if performing 3D-PTV with one camera | |
| 250 Watts Spotlight Halogen | General Electrics | 23719 | Light source |
| OpenPTV (Software) | OpenPTV (http://www.openptv.net) | Open source particle tracking software (Note: available as a service for anyone who wants to use it without all the installation mess or computer power availability problems). |
References
- Wygnanski, I., Fiedler, H. Some measurements in the self preserving jet. , Cambridge university press. (1968).
- Rajaratnam, N. Turbulent jets. , Elsevier. (1976).
- Panchapakesan, N., Lumley, J. Turbulence measurements in axisymmetric jets of air and helium. Part 1. Air jet. J Fluid Mech. 246, 197-223 (1993).
- Hussein, H. J., Capp, S. P., George, W. K. Velocity measurements in a high-Reynolds-number, momentum-conserving, axisymmetric, turbulent jet. J Fluid Mech. 258, 31-75 (1994).
- Yule, A. Large-scale structure in the mixing layer of a round jet. J Fluid Mech. 89, 413-432 (1978).
- Yule, A., Chigier, N., Ralph, S., Boulderstone, R., Venturag, J. Combustion-transition interaction in a jet flame. AIAA Journal. 19, 752-760 (1981).
- Quinn, W. Upstream nozzle shaping effects on near field flow in round turbulent free jets. Eur J Mech B-Fluid. 25, 279-301 (2006).
- Mi, J., Nathan, G. J., Luxton, R. E. Centreline mixing characteristics of jets from nine differently shaped nozzles. Exp Fluids. 28, 93-94 (2000).
- Karlsson, R. I., Eriksson, J., Persson, J. LDV measurements in a plane wall jet in a large enclosure. DTIC [Internet]. , Available from: http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADP008905 (1992).
- Liepmann, D., Gharib, M. The role of streamwise vorticity in the near-field entrainment of round jets. J Fluid Mech. 245, 643-668 (1992).
- Oh, S. K., Shin, H. D. A visualization study on the effect of forcing amplitude on tone-excited isothermal jets and jet diffusion flames. Int J Energ Res. 22, 343-354 (1998).
- Cenedese, A., Doglia, G., Romano, G., De Michele, G., Tanzini, G. LDA and PIV velocity measurements in free jets. Exp Therm Fluid Sci. 9, 125-134 (1994).
- Wang, H., Peng, X., Lin, W., Pan, C., Wang, B. Bubble-top jet flow on microwires. Int J Heat Mass Tran. 47, 2891-2900 (2004).
- Shestakov, M. V., Tokarev, M. P., Markovich, D. M. 3D Flow Dynamics in a Turbulent Slot Jet: Time-resolved Tomographic PIV Measurements. 17th Int Symp on Applications of Laser Techniques to Fluid Mechanics. , (2014).
- Bridges, J., Wernet, M. P. Measurements of the aeroacoustic sound source in hot jets. AIAA [Internet]. , Available from: http://arc.aiaa.org/doi/abs/10.2514/6.2003-3130 (2003).
- Scarano, F., Bryon, K., Violato, D. Time-resolved analysis of circular and chevron jets transition by tomo-PIV. 15th Int Symp on Applications of Laser Techniques to Fluid Mechanics. , (2010).
- Holzner, M., Liberzon, A., Nikitin, N., Kinzelbach, W., Tsinober, A. Small-scale aspects of flows in proximity of the turbulent/nonturbulent interface. Phys Fluids. 19, 071702 (2007).
- Holzner, M., et al. A Lagrangian investigation of the small-scale features of turbulent entrainment through particle tracking and direct numerical simulation. J Fluid Mech. 598, 465-475 (2008).
- Chang, T. P., Wilcox, N. A., Tatterson, G. B. Application of image processing to the analysis of three-dimensional flow fields. Opt Eng. 23, 283-287 (1984).
- Racca, R., Dewey, J. A method for automatic particle tracking in a three-dimensional flow field. Exp Fluids. 6, 25-32 (1988).
- Maas, H. G., Gruen, D., Papantoniou, D. Particle tracking velocimetry in three-dimensional flows. Exp Fluids. 15, 133-146 (1993).
- Kasagi, N., Matsunaga, A. Three-dimensional particle tracking velocimetry measurement of turbulence statistics and energy budget in a backward-facing step flow. Int J Heat Fluid Fl. 16, 477-485 (1995).
- Virant, M., Dracos, T. 3D PTV and its application on Lagrangian motion. Meas Sci Technol. 8, 1539 (1997).
- Willneff, J. A spatio-temporal matching algorithm for 3 D particle tracking velocimetry. , Mitteilungen- Institut fur Geodasie und Photogrammetrie an der Eidgenossischen Technischen Hochschule Zurich. Zurich. (2003).
- Rosi, G. A., Sherry, M., Kinzel, M., Rival, D. E. Characterizing the lower log region of the atmospheric surface layer via large-scale particle tracking velocimetry. Exp Fluid. 55, 1-10 (2014).
- Fu, S., Biwole, P. H., Mathis, C. Particle Tracking Velocimetry for indoor airflow field: A review. Build Environ. 87, 34-44 (2015).
- Kolaas, J., Jensen, A., Mielnik, M. Visualization and measurements of flows in micro silicon Y-channels. Eur Phys J E. 36, 1-11 (2013).
- Gülan, U., et al. Experimental study of aortic flow in the ascending aortavia Particle Tracking Velocimetry. Exp Fluids. 53, 1469-1485 (2012).
- Mei, R. Velocity fidelity of flow tracer particles. Exp Fluids. 22, 1-13 (1996).
- Tropea, C., Yarin, A. L., Foss, J. F. Springer handbook of experimental fluid mechanics. 1, Springer Science & Business Media. (2007).
- Melling, A. Tracer particles and seeding for particle image velocimetry. Meas Sci Technol. 8, 1406 (1997).
- Hering, F., Leue, C., Wierzimok, D., Jähne, B. Particle tracking velocimetry beneath water waves. Part I: visualization and tracking algorithms. Exp Fluids. 23, 472-482 (1997).
- Biferale, L., et al. Lagrangian structure functions in turbulence: A quantitative comparison between experiment and direct numerical simulation. Phys Fluids. 20, 065103 (2008).
- Lüthi, B., Tsinober, A., Kinzelbach, W. Lagrangian measurement of vorticity dynamics in turbulent flow. J Fluid mech. 528, 87-118 (2005).
- Hoyer, K., et al. 3d scanning particle tracking velocimetry. Exp Fluids. 39, 923-934 (2005).
- Kim, J. -T. Three-dimensional particle tracking velocimetry for turbulence applications. , UIUC. http://chamorro.mechse.illinois.edu/3d.htm (2015).
- Lüthi, B. Some aspects of strain, vorticity and material element dynamics as measured with 3D particle tracking velocimetry in a turbulent flow. ETH Zürich. , Nr. 14893 (2002).
- Pope, S. B. Turbulent flows. , Cambridge university press. (2000).
