Introduction and Scope | Browse Archives | Support Funding Sources |
Sorry.
You are not permitted to access the full text of articles.
If you have any questions about permissions,
please contact the Society.
죄송합니다.
회원님은 논문 이용 권한이 없습니다.
권한 관련 문의는 학회로 부탁 드립니다.
[ Original Paper ] | |
The KSFM Journal of Fluid Machinery - Vol. 19, No. 2, pp. 5-10 | |
Abbreviation: The KSFM Journal of Fluid Machinery | |
ISSN: 2287-9706 (Print) | |
Print publication date Apr 2016 | |
Received 11 Aug 2015 Revised 05 Jan 2016 Accepted 11 Jan 2016 | |
DOI: https://doi.org/10.5293/kfma.2016.19.2.005 | |
다점 피토관을 이용한 기체 유량 측정의 불확도 평가 | |
양인영*, † ; 이보화*
| |
*한국항공우주연구원 엔진시스템연구팀 | |
Uncertainty Assessment of Gas Flow Measurement Using Multi-Point Pitot Tubes | |
Inyoung Yang*, † ; Bo-Hwa Lee*
| |
*Engine system research team, Korea Aerospace Research Institute | |
Correspondence to : †E-mail : iyyang@kari.re.kr | |
Funding Information ▼ |
Gas flow measurement in a closed duct was performed using multi-point Pitot tubes. Measurement uncertainty was assessed for this measurement method. The method was applied for the measurement of air flow into a gas turbine engine in an altitude engine test facility. 46 Pitot tubes, 15 total temperature Kiel probes and 9 static pressure tabs were installed in the engine inlet duct of inner diameter of 264 mm. Five tests were done in an airflow range of 2~10 kg/s. The flow was compressible and the Reynolds numbers were between 450,000 and 2,220,000. The measurement uncertainty was the highest as 6.1% for the lowest flow rate, and lowest as 0.8% for the highest flow rate. This is because the difference between the total and static pressures, which is also related to the flow velocity, becomes almost zero for low flow rate cases. It was found that this measurement method can be used only when the flow velocity is relatively high, e.g., 50 m/s. Static pressure was the most influencing parameter on the flow rate measurement uncertainty. Temperature measurement uncertainty was not very important. Measurement of boundary layer was found to be important for this type of flow rate measurement method. But measurement of flow non-uniformity was not very important provided that the non-uniformity has random behavior in the duct.
Keywords: Gas Flow Measurement, Multi-point Pitot Tubes, Measurement Uncertainty, Boundary Layer 키워드: 기체 유량 측정, 다점 피토관, 측정 불확도, 경계층 |
한국항공우주연구원에서는 가스 터빈 엔진의 고공 환경 성능을 시험적으로 평가하기 위한 엔진고공환경시험설비(altitude engine test facility; AETF)를 운용하고 있다. 가스 터빈 엔진 중 터보젯 및 터보팬 엔진의 경우 엔진의 성능을 평가하기 위해서는 엔진으로 유입되는 공기의 유량을 평가하는 것이 반드시 필요하다. 이에 따라 AETF에는 벤투리 파이프 형식의 유량계를 설치하여 사용하고 있다. 그러나 벤투리 파이프가 대형, 대용량(공기 유량 2~10kg/s)임으로 인해 교정의 어려움 등이 있어 다른 공기 유량 측정 방법을 사용하여 상호 검증하는 것이 필요하였다.
벤투리 파이프 이외의 유량 측정 방법으로는 초음파 유량계, 열선식 유량계 등 다른 유량 측정 방법이 있다. 이 연구에서는 다점 피토관을 이용한 기체 유량 측정 방법을 적용하고자 하였다. 이것은 시험 대상물인 엔진 입구에서의 시험 조건(온도, 압력)을 정확히 정의하기 위해 기왕에 다점 피토관을 설치하여 사용하는 것이 필요하기 때문이었다. 이 연구에서는 이러한 기체 유량 측정 방법을 사용하는 데 있어서 측정 불확도를 평가함으로써 이 방법이 실효성을 가질 수 있는 측정 범위를 파악하고 이 방법의 측정 불확도를 향상시킬 수 있는 방법을 모색하고자 하였다.
다점 피토관을 사용한 기체 유량 측정 및 불확도 평가 방법에 관해서는 ISO 3966 표준(1)을 참조할 수 있으나, 기존 연구 사례 중에서는 이를 충실히 따른 것이 많지 않다. 피토관의 가장 대표적인 용례인 유동 속도 측정에 관해서는 최진철(2), 오세윤(3) 등의 연구가 있으나 관내에서의 유량 측정과 연결시킨 연구는 아니었다. 임재명(4), Nguyen(5) 등의 연구는 유량 측정에 관한 연구이기는 하나 압축성을 고려하지 않은 저속 유동을 대상으로 하는 연구였다. 이번 연구에서는 저자 일동이 기존에 수행한 연구(6)의 연장선상에서 이상과 같은 연구 사례에서 다루지 않았던 압축성 기체 유동의 관내 유량 측정에 관한 연구를 수행하고자 하였으며, 앞서 언급한 ISO 표준에 따라 실제 시험 데이터를 적용해 보고자 하였다.
엔진고공환경시험설비에서는 엔진 공급 배관 상류에 벤투리 파이프가 설치되어 있다. 이 연구에서는 이에 더하여 엔진 입구 배관에 다점 피토관을 설치하여 공급 공기 유량을 측정하고자 하였다. 엔진 입구 배관은 내경 264mm인 원형 배관이다. 공기 유량 범위는 2~10kg/s이며 이 때 관내 유속은 최대 90m/s로서 압축성이 고려되어야 하는 조건이다.
Fig. 1은 엔진 입구 배관의 형상이며 유동 방향은 그림의 오른쪽에서 왼쪽이다. 그림에서 “Stilling chamber”로 표시된 부분의 상류에 벤투리 파이프가 설치되어 있으며 “09”로 표시된 section의 하류에 엔진이 설치되어 있다. 다점 피토관은 “05”로 표시된 section에 장착된다. Fig. 2는 05 section에 설치된 다점 피토관의 배관 내 반경 방향 위치를 나타낸다. 피토관은 레이크의 형태로 설치되어 있는데, Fig. 2와 같이 반경 방향으로 16개의 피토관이 배치된 레이크가 배관 안에 원주 방향으로 120˚ 간격으로 3세트 설치되어 있다. 단, 배관 중심에 위치한 피토관은 그 중 1세트에만 있으며 나머지 2세트에는 없으므로 나머지 2세트에는 15개씩의 피토관이 있다(총 46개 피토관). 이와 별도로 공기 유동의 전온도를 측정하는 Kiel type 프로브 레이크가 있는데, 역시 3세트이며, 원주 방향으로 120˚ 간격으로 설치되어 있고, 각각 반경 방향으로 5개의 프로브에서 전온도를 측정한다(총 15개 프로브). 또한 05 section의 벽면에 40˚ 간격으로 9개의 측정용 hole을 설치하여 각각의 벽면 정압력을 측정한다.
기체 유량 측정 방법은 이보화 등(6)에 상세히 기술되어 있는데, 여기에서는 논문의 전개를 위하여 간략히 인용하고, 부록에 ISO 3966 표준(1)과의 차이를 기술하였다.
공기 유량 W05는 식 (1)로 계산한다. 식 (1)에서 공기 유속 V05는 배관 내 유동의 마하 수 MD,05와 유동의 음속 a05를 구하여 식 (2)와 같이 계산한다. 음속은 식 (3)과 같이 계산하며, 여기에서 정온도 Ts,05는 식 (4)로 계산한다. 식 (1)에서의 밀도 ρ05는 식 (5)로 계산한다.
(1) |
(2) |
(3) |
(4) |
(5) |
여기에서 ps,05는 원주 방향으로 여러 위치에서 측정한 벽면 정압력 MD,05,k의 산술평균값이다. 그런데 MD,05나 Tt,05는 반경 방향으로도 일정하지 않기 때문에 각 지점에서 측정한 압력과 온도를 사용한 평균값을 적용하여야 한다. 위 식 (2) 및 식 (4)에서 MD,05와 Tt,05는 각 프로브가 대표하는 배관 내 면적의 가중 인자를 고려한 면적 평균으로 구한다. 즉,
(6) |
(7) |
여기에서 i, j는 각각 전압력 레이크, 전온도 레이크 상의 프로브의 번호이며 rpt,i, rTt,j는 그 각각의 프로브의 반경 방향 위치이다. 전온도 Tt,05,j는 직접적인 측정값이며 MD,05,i는 식 (8)과 같이 구한다.
(8) |
정리하면, 공기 유량은 각 프로브에서 측정한 전온도, 전압력, 정압력과 나머지 인자(덕트 직경, 기체 상수, 기체 비열비, 프로브의 위치)의 함수이다. 이 나머지 인자는 아래와 같이 상수로 가정하였다.
측정 불확도 계산 방법은 ISO(7)의 방법을 따른다. 측정 불확도 계산에 있어서 R, γ, D05, rpt,i, rTt,j의 측정 불확도는 무시하였다. 결과적으로 기체 유량의 측정 불확도에 기여하는 변수는 pt,05,i, Tt,05,j, ps,05,k의 세 변수이다. 이 세 변수 각각의 감도 계수와 표준 측정 불확도를 계산하고 이를 합성하여 기체 유량의 측정 불확도를 계산한다. 식으로 나타내면 식 (9)와 같다.
(9) |
여기에서 감도 계수를 계산하기 위한 편미분은 각 측정 인자 1% 변화에 대한 W05의 변화량으로 수치적으로 근사화하여 계산한다. 측정 인자 pt,05,i, Tt,05,j, ps,05,k의 표준 측정 불확도는 실제 시험 데이터의 산포 및 각각의 측정 계통 상에서 발생하는 불확도 요인을 A형(7)과 B형(7) 방법으로 평가한다. 이 연구에서 pt,05,i 및 ps,05,k 측정 센서는 실제로는 대기압과의 차압을 측정하며, 그 측정값에 대기압 측정 센서의 측정값을 더하여 pt,05,i=pt,05,i,g+patm, ps,05,k=ps,05,k,g+patm과 같이 계산되므로 이를 합성하여야 한다. 또한 pt,05,i,g와 ps,05,k,g는 센서의 불확도 외에 추가 불확도 요인이 없으나 patm은 센서 불확도에 아날로그-디지털 변환(ADC) 보드의 불확도를 추가하여야 한다. Tt,05,j도 센서 불확도에 ADC 보드의 불확도를 추가하여야 한다. Table 1에 이러한 측정 불확도 요인을 정량화하여 정리하였다.
Measured Parameter | Uncertainty Component | Uncertainty | Remarks | |
---|---|---|---|---|
pt,05,i | A-type | 0.04-0.39%Rd | ||
pt,05,i,g | Sensor | 0.05%FS | FS = 103 kPa | |
patm | Sensor | 0.02%FS | FS = 80-110 kPa | |
ADC Board | 0.02%FS | FS = 32 VDC 0.07 kPa equiv. |
||
Tt,05,j | A-type | 0.1-0.2%Rd | ||
Sensor | 1.5 ℃ | - | ||
ADC Board | 0.02%FS | FS = 32 VDC 1.2 ℃ equiv. |
||
ps,05,k | A-type | 0.03-0.09%Rd | ||
ps,05,k,g | Sensor | 0.05%FS | FS = 103 kPa | |
patm | Sensor | 0.02%FS | FS = 80-110 kPa | |
ADC Board | 0.02%FS | FS = 32 VDC 0.07 kPa equiv. |
측정 불확도 평가를 위하여 실제 시험을 수행하였다. 시험은 정해진 배관 및 다점 피토관 형상에 대해 공기 유량 약 2~10kg/s 범위 내의 5개 공기 유량 값에 대해 수행하였다. 시험 데이터 및 그로부터 위 2.2절과 같은 방법으로 계산한 유량 데이터는 Table 2와 같다. 표에서 전압력, 전온도, 정압력 데이터는 편의상 1개 대표 센서의 값만 나타내었다. 파이프 레이놀즈 수는 ReD=450,000~2,220,000 범위이다.
case number | 1 | 2 | 3 | 4 | 5 |
---|---|---|---|---|---|
Test data | |||||
pt,05,1 (kPaA) | 104.7 | 115.1 | 132.3 | 154.6 | 185.3 |
Tt,05,1 (K) | 296.3 | 299.5 | 299.4 | 300.2 | 301.0 |
ps,05,1 (kPaA) | 104.2 | 112.8 | 128.0 | 148.2 | 176.6 |
V05 (m/s) | 30.4 | 54.5 | 71.4 | 80.5 | 86.0 |
ReD | 4.5E5 | 8.7E5 | 1.3E6 | 1.7E6 | 2.2E6 |
Air flow | |||||
W05 (kg/s) | 2.015 | 3.911 | 5.829 | 7.642 | 9.724 |
Wven (kg/s) | 1.998 | 3.872 | 5.855 | 7.776 | 9.866 |
W05 - Wven (%) | 0.85 | 1.01 | -0.44 | -1.73 | -1.44 |
Extended uncertainty | |||||
2 ×u(W05) (kg/s) | 0.123 | 0.077 | 0.059 | 0.061 | 0.074 |
2 ×u(W05) (%) | 6.1 | 2.0 | 1.0 | 0.8 | 0.8 |
2 ×u(Wven) (kg/s) | 0.038 | 0.034 | 0.045 | 0.057 | 0.072 |
2 ×u(Wven) (%) | 1.9 | 0.9 | 0.8 | 0.7 | 0.7 |
Uncertainty percentage contribution | |||||
19.6 | 15.9 | 29.5 | 16.3 | 18.0 | |
0.1 | 0.7 | 2.7 | 4.3 | 4.7 | |
80.3 | 83.4 | 67.9 | 79.4 | 77.3 |
W05는 벤투리 유량계로 측정한 Wven에 비하면 최대 1.73%까지 차이가 난다. 저유량에서는 Wven보다 W05가 많은 것으로 측정되고 고유량에서는 적은 것으로 측정되는데, 아래에 기술한 것처럼 저유량에서는 W05의 불확도가 다소 높기 때문에 고유량에서의 경향을 받아들이는 것이 옳다고 판단된다. Wven보다 W05가 적게 측정되는 이유에 대해서는 추가적인 분석이 필요하다.
2.3절의 방법으로 평가한 유량의 확장 불확도(extended uncertainty)도 Table 2에 나타내었다. 불확도의 절대값(2×uW05)을 보면 대체로 저유량에서 값이 크고 고유량에서 값이 작음을 알 수 있다. 따라서 상대 측정 불확도도 저유량에서 크다. Case 1(유량 2.015kg/s)에서는 상대 확장 불확도가 6.1%에 이르러 다점 피토관을 이용한 유량 측정을 적용하는 것이 적절하지 못한 영역이라 할 수 있다. 반면 case 2(유량 3.911kg/s) 이상에서는 불확도가 2% 이하이므로 다른 측정 방법과 비교해 나쁘지 않다고 할 수 있다.
case number | 1 | 2 | 3 | 4 | 5 |
---|---|---|---|---|---|
Air flow from table 2 | |||||
W05 (kg/s) | 2.015 | 3.911 | 5.829 | 7.642 | 9.724 |
Air flow neglecting boundary layer (sec. 3.3) | |||||
W05,nobl (kg/s) | 2.135 | 4.156 | 6.191 | 8.129 | 10.329 |
W05,nobl - W05 (%) | 5.3 | 6.3 | 6.2 | 6.4 | 6.2 |
Air flow neglecting flow non-uniformity (sec. 3.4) | |||||
max(σ(pt>,05,r)) (kPa) | 0.040 | 0.14 | 0.45 | 0.46 | 0.57 |
W05,1rake (kg/s) | 2.027 | 3.907 | 5.795 | 7.619 | 9.669 |
W05,1rake - W05 (%) | 0.6 | 0.1 | -0.6 | –0.3 | -0.6 |
불확도의 절대값이 저유량에서 큰 이유는 다음과 같다. 다점 피토관을 이용한 유량 측정 방법에서는 유동의 전압력과 정압력의 차(pt,05,i-ps,05)를 이용하여 유량을 계산하는데(식 (8)) 저유량에서는 이 값이 거의 0에 가까워지기 때문에 pt,05,i나 ps,05k의 감도 계수가 커져 측정 불확도가 커진다. 이번 연구에서의 데이터를 보면 case 1에서는 pt,05,i-ps,05=0.556kPa이고 case 2에서는 pt,05,i-ps,05=7.70kPa이다. 참고문헌 (1)에 따르면 이번 시험에서 pt,05,i-ps,05=3Pa 이상이면 피토관을 이용한 유량 측정 방법을 적용할 수 있다고 하였으나(해당 문헌 8.1절의 조건식) 실제로는 불확도 측면에서 만족할 만한 측정을 하기 위해서는 그보다 훨씬 높은 차압 조건이 필요하였다.
또한 Table 2에는 각 측정 인자가 공기 유량의 측정 불확도에 얼마나 기여하는지를 나타내는 불확도 퍼센트 기여도(uncertainty percentage contribution; UPC)를 함께 나타내었는데, 편의상 46개 전압력 측정 인자의 UPC의 합, 15개 전온도 측정 인자의 UPC의 합, 9개 정압력 측정 인자의 UPC의 합으로 나타내었다. 이 연구에서 사용한 데이터의 경우 정압력의 UPC가 나머지 인자에 비해 크게 높았으며, 다음으로 전압력, 전온도의 순이었다. 전온도의 UPC는 고유량에서 늘어나며, 전압력과 정압력의 UPC는 유량에 따라 뚜렷한 경향이 나타나지 않았다.
Fig. 3은 피토관 데이터를 이용하여 계산한, 각 시험 조건별 덕트 내 반경 방향 유속 분포를 나타낸다. 그림을 보면 case 1의 경우에는 경계층 두께가 약 12mm 정도이고 나머지 case의 경우에는 20mm 내외이다. 또한 유속이 높을수록 경계층 두께가 두꺼워지는 것도 관찰된다.
경계층을 무시하고 배관 중심에 가까운 5개의 피토관에서 측정한 전압력만을 사용하여 유량을 계산한 결과를 Table 3에 나타내었다. 경계층을 고려한 Table 2의 결과에 비해 유량이 5.3~6.4% 높게 계산되었다. 따라서 다점 피토관을 사용한 유량 측정의 경우 경계층을 잘 측정하는 것이 반드시 필요함을 알 수 있다.
이 연구의 경우 원주 방향으로는 3군데에서만 측정했으므로 측정 데이터를 사용하여 유동 불균일성을 논하기는 어렵지만 유동 불균일성의 영향을 일반화하여 논의하는 데에 참고값으로 삼기로 한다. 이 연구에서의 유동은 경계층이 있는 유동인 관계로 반경 방향의 불균일성은 필연적으로 존재하므로 원주 방향의 불균일성만 다루기로 하고 유동 불균일도를 “동일한 반경에 있는 3개의 전압력 측정값의 표준 편차(σ(pt,05,r))”로 정의하였다.
유동 불균일도는 유속에 따라서도 다르고 반경 방향 위치에 따라서도 다른데, 시험 데이터로부터 계산한 각 반경에서의 유동 불균일도의 최대값(max(σ(pt,05,r)))을 Table 3에 나타내었다. 표를 보면 유동 불균일도는 유량 계산의 참조값인 pt,05,i-ps,05과 비교하면 7.2~12%에 달하는 큰 값이며, 시험 조건에 따라서는 case 1에서 최소이고 case 5에서 최대이다.
유동 불균일도를 무시하고 측정을 수행하여 이를 바탕으로 유량을 계산하는 경우 유동 불균일도를 고려하는 것에 비해 유량 계산 결과가 얼마나 차이가 나는지를 알아보기 위해 시험 데이터 중 1개 레이크의 측정값만을 사용하여 유량을 계산해 보았다. 그 결과를 Table 3에 나타내었는데, Table 2와 비교하여 최대 0.6% 차이가 났다. 유동 불균일도에 비하여 유량 계산값의 차이가 적은 것은 이번 연구의 시험에서 유동의 불균일도가 어느 한 방향으로 편중되지 않았기 때문이다. 또한 유동 불균일도를 정의하는 데 사용된 pt,05,i의 공기 유량에 대한 UPC가 앞서 3.2절에서 기술한 것처럼 높지 않다는 점도 원인이다.
이상과 같은 연구 결과로부터 다음 결론을 얻었다.
본 논문은 2015년도 국가과학기술연구회 사업 재원으로 지원을 받아 작성된 논문임(CMPp.-13-5-KARI).
1. | ISO, (2008), Measurement of Fluid Flow in Closed Conduits – Velocity Area Method Using Pitot Static Tubes, ISO 3966:2008. |
2. | Choi, J. C., Lee, K. G, and Park, C. W, (2001), “Development of the Pitot-Tube to Measure Flight Altitude and Velocity,”, Journal of the Korean Society of Propulsion Engineers, 5(3), p19-2424. |
3. | Oh, S. Y, Lee, J. G, Kim, S. C, Kim, S. H, and Ahn, S. K, (2012), “An Experimental Study of Test Section Velocity Calibration for Low-Speed Wind Tunnel,”, Journal of The Korean Society for Aeronautical and Space Sciences, 40(3), p230-236.![]() |
4. | Lim, J. M, (2000), An Experimental Study on the Flow Measurement by Using Multi-point Pitot Tube, Master Thesis, Mokpo National University. |
5. | Nguyen, D. T, (2013), Evaluation of Uncertainty for Gas Flow Measurement with S Pitot Tube, Master Thesis, University of Science and Technology. |
6. | Lee, B. H, Lee, K. J, and Yang, S. S, (2010), “The Preduction of Air Flow and Pressure Loss at Inlet Duct,”, Journal of the Korean Society of Propulsion Engineers, 14(1), p48-55. |
7. | ISO/IEC, (2008), Uncertainty of Measurement – Part 3: Guide to the Expression of Uncertainty in Measurement, ISO/IEC Guide 98-3. |
이 논문에서 사용한 공기 유량 측정 방법은 ISO 3966 표준(이하 “표준”으로 칭함)과 거의 동일하며, 아래와 같은 점에서만 차이가 있다.