선택한 단어 수는 입니다.
최소 단어 이상 선택하여야 합니다.
최소 단어 이상 선택하여야 합니다.
선택한 단어 수는 30입니다.
최대 10 단어까지만 선택 가능합니다.
최대 10 단어까지만 선택 가능합니다.
문제 정의
- 본 연구에서는 이러한 항공기 노즐 열유동장, IR 신호 해석기법 및 IR 저감 설계기법을 검증할 목적으로 먼저 지상에서 실제 시험이 가능한 마이크로 터보제트 엔진에 적용하였다. 축소 모형인 마이크로 터보제트 엔진 노즐의 열유동장 예측과 시험은 시간과 예산을 아낄 수 있는 방법이기 때문이다.
- 본 연구에서는 항공기 IR 저감 설계 연구를 위하여 마이크로 터보제트 엔진 노즐의 열유동장 특성을 CFD 전산코드로 분석하였다. 마이크로 터보제트 엔진 성능분석을 통하여 노즐 입구에서의 열유동 조건을 도출하였다.
- 그리고 마이크로 터보제트 엔진 노즐에 대한 CFD의 해석값과 최신 무인전투기 중에 하나인 Dassault Aviation의 nEUROn 급 항공기의 full-scale 추진기관 노즐의 열유동장과 특정 파장대에서의 IR 신호 수준을 비교하였다. 이를 통하여 full-scale 전투기와 마이크로 터보제트 엔진 노즐 유동 및 IR 신호 특성의 상관 관계를 분석하고 궁극적으로 지상에서의 마이크로 터보제트 엔진을 통한 축소 모형 시험의 타당성을 점검하고자 하였다.
가설 설정
- 설계된 노즐은 지상조건의 외부 대기 조건하에 작동되며, 노즐 출구의 형상과 세장비 그리고 노즐의 곡률을 변화시키며 해석하였다. 노즐 입구로 유입되는 유동은 Soot가 존재하지 않는 완전연소를 가정하여 N2 74%, H2O 13%, CO2 13%로 몰분율 기준에 따라 설정하였다.
제안 방법
- 마이크로 터보제트 엔진 성능분석을 통하여 노즐 입구에서의 열유동 조건을 도출하였다. Full-scale에 관한 이전 연구와의 비교를 위하여 nEUROn의 노즐의 형상 변수(AR 및 곡률) 형태를 벤치마킹하여 여러 개의 노즐을 설계하였다.
- 노즐출구를 기준으로 5 m의 위치에서 고도각과 방위각에 대하여 계산하였다. IR 신호를 계산하고자 하는 지점에서 노즐 중심까지의 직선거리를 1000등분 하여 각 점에 대한 온도, 압력, 분압을 도출했다. IR 신호 계산 결과 같은 지점에서 계산 하였을 때, Fig.
- UCAV와 마이크로 터보제트 엔진의 상사성에 관한 이전연구[13]를 확장하여 운용조건 및 스케일 차이에 따른 무차원수 분석을 수행하였다. 본 연구에서는 노즐 열유동장 및 IR 신호의 스케일 효과를 분석하기 위해 에너지 보존 편미분 방정식의 대류(convection)항과 복사(radiation)항의 비인 Bo(Boltzmann number)[14], 전도(conduction)항과 복사항의 비인 NR(Stark number)를 도입하였다.
- nEUROn급의 엔진을 획득하여 시험하는 것은 많은 시간과 예산의 문제로 어렵기 때문에, 효과적이고 다양한 조건에서의 검증을 위해 본 연구에서는 노즐 직경이 13 센티미터 정도인 소형 마이크로 터보제트 엔진을 고려하였다. Fig.
- 먼저 x에 대한 y의 좌표를 식 (1)을 이용하여 먼저 구하고 x, y가 포함된 단면적의 타원 넓이와 타원식을 이용하여 z좌 표를 구한다. 그 다음 hyperbolic tangent를 노즐 단축의 가장 윗선과 1/4선에 적용하여 두 가지(소, 대)의 곡률을 생성하였다. 설계된 형상변형 노즐을 Fig.
- 그리고 M(Mach number)와 Re(Reynolds number), τ(optical thickness)를 도입하여 분석하였다.
- 이를 위해 본 연구에서는 상용 CFD 프로그램인 ANSYS의 FLUENT[7]를 이용하여 열유동장 해석을 수행하고 노즐 출구의 온도 및 압력, 추력 등 그 특성을 분석하였다. 그리고 마이크로 터보제트 엔진 노즐에 대한 CFD의 해석값과 최신 무인전투기 중에 하나인 Dassault Aviation의 nEUROn 급 항공기의 full-scale 추진기관 노즐의 열유동장과 특정 파장대에서의 IR 신호 수준을 비교하였다. 이를 통하여 full-scale 전투기와 마이크로 터보제트 엔진 노즐 유동 및 IR 신호 특성의 상관 관계를 분석하고 궁극적으로 지상에서의 마이크로 터보제트 엔진을 통한 축소 모형 시험의 타당성을 점검하고자 하였다.
- 엔진의 곡률을 변형하기 위해 nEUROn과 유사한 곡선을 가지는 hyperbolic tangent 함수를 적용하였다[8]. 노즐의 단면적이 바뀌면 그 특성이 바뀌기 때문에, 최대한 같은 면적을 가지는 노즐을 이용하여 계산하기 위해, Fig 2의 원형 노즐을 기준으 로 하여 변형 노즐을 설계 하였다. 먼저 엔진 성능 분석을 통해 얻어진 원형 축소 노즐을 x축 방향으로 8등분을 하였다.
- 본 연구에서는 항공기 IR 저감 설계 연구를 위하여 마이크로 터보제트 엔진 노즐의 열유동장 특성을 CFD 전산코드로 분석하였다. 마이크로 터보제트 엔진 성능분석을 통하여 노즐 입구에서의 열유동 조건을 도출하였다. Full-scale에 관한 이전 연구와의 비교를 위하여 nEUROn의 노즐의 형상 변수(AR 및 곡률) 형태를 벤치마킹하여 여러 개의 노즐을 설계하였다.
- 마이크로 터보제트 엔진의 변형 노즐을 해석하기 위해서 형상 변형을 위한 기준 형상을 설정하였다. 기준 형상으로 Dassault Aviation에서 개발한 nEUROn의 노즐을 참고하였다.
- 추진 시스템에 대한 연구는 엔진 성능 분석을 통한 노즐 입구의 초기조건 계산과 엔진 배기 노즐 형상 설계로 이루어진다. 본 연구에서는 네덜란드 AMT의 마이크로 터보제트 엔진의 성능분석을 통하여 얻은 노즐 입구의 조건들을 이용하였다. 설계된 노즐은 지상조건의 외부 대기 조건하에 작동되며, 노즐 출구의 형상과 세장비 그리고 노즐의 곡률을 변화시키며 해석하였다.
- UCAV와 마이크로 터보제트 엔진의 상사성에 관한 이전연구[13]를 확장하여 운용조건 및 스케일 차이에 따른 무차원수 분석을 수행하였다. 본 연구에서는 노즐 열유동장 및 IR 신호의 스케일 효과를 분석하기 위해 에너지 보존 편미분 방정식의 대류(convection)항과 복사(radiation)항의 비인 Bo(Boltzmann number)[14], 전도(conduction)항과 복사항의 비인 NR(Stark number)를 도입하였다.
- 본 연구에서는 네덜란드 AMT의 마이크로 터보제트 엔진의 성능분석을 통하여 얻은 노즐 입구의 조건들을 이용하였다. 설계된 노즐은 지상조건의 외부 대기 조건하에 작동되며, 노즐 출구의 형상과 세장비 그리고 노즐의 곡률을 변화시키며 해석하였다. 노즐 입구로 유입되는 유동은 Soot가 존재하지 않는 완전연소를 가정하여 N2 74%, H2O 13%, CO2 13%로 몰분율 기준에 따라 설정하였다.
- 최종적으로 축소 모델인 마이크로 터보제트 엔진과 full-scale nEUROn 급의 노즐 해석결과를 비교하여 스케일링 특성을 분석하였다. 에너지 보존에 기여를 하는 복사 등 여러 물리적 항들의 상대적 값 크기 비교를 통해 플룸의 광학적 특성에 관한 Boltzmann 수와 Stark 수를 도출하였다. 원칙적으로 이러한 무차원수들은 비슷한 크기로 유지하는 것이 이상적이나, 이를 달성하기에는 현실적인 제한이 많다.
- 4에 해석에 사용된 노즐 주변의 격자계를 나타내었다. 열유동장 해석을 위한 경계조건은 Table 2와 같으며, 노즐 입구의 입력 조건은 마이크로 터보제트 엔진 성능분석을 통해 도출된 결과를 적용하였다.
- 이를 위해 본 연구에서는 상용 CFD 프로그램인 ANSYS의 FLUENT[7]를 이용하여 열유동장 해석을 수행하고 노즐 출구의 온도 및 압력, 추력 등 그 특성을 분석하였다.
- 5에 나타내었는데, 노즐 표면과 중심에서의 압력분포가 실험값과 정성적으로 일치하는 것을 확인하였다. 이를 통해 본 연구에서 사용된 전산해석 기법의 신뢰성을 확보하였다.
- 전산해석 결과를 바탕으로 narrow-band 모델을 이용한 IR 신호 계산을 수행하였다. 관측 위치는 실제 실험 상황과 유사하게 설정을 하였다.
- 최종적으로 축소 모델인 마이크로 터보제트 엔진과 full-scale nEUROn 급의 노즐 해석결과를 비교하여 스케일링 특성을 분석하였다. 에너지 보존에 기여를 하는 복사 등 여러 물리적 항들의 상대적 값 크기 비교를 통해 플룸의 광학적 특성에 관한 Boltzmann 수와 Stark 수를 도출하였다.
- Narrow-band 모델은 산란은 하지 않고 흡수와 방출만을 하는 경우의 복사 열전달 방정식을 기반으로 계산하며, 흡수계수에 대한 통계표를 사용하여 계산된다. 플룸 IR 신호 해석수행 시, 각 지점의 온도, 압력, 혼합물의 물질에 따른 분압을 도출하여 파장별 강도를 도출하였다.
- 온도 변화에 관한 해석결과는 아래 Table 3에 나타내었다. 한편 형상 변형 노즐의 온도 감소뿐만 아니라, 각종 페널티에 대한 영향도 같이 고려해야 하므로, 노즐출구의 압력, 속도, 질량 유량을 이용하여 각 노즐 형상별 추력을 계산하였다. 그 결과 형상 변형 노즐에서의 추력이 원형 노즐 대비 최대 12%의 감소를 보였으며, 추력의 감소는 가로세로비가 크거나 곡률이 클수록 더 크게 나타났다.
대상 데이터
- 사용된 수치코드 FLUENT의 노즐 열유동장 해석에 관한 신뢰성을 확보하기 위해 해석코드 검증시 널리 사용되는 JPL 노즐 실험값[9]을 고려하였다. 공기의 전 압력, 전 온도를 각각 1.032 MPa, 550 K로 설정하였다. 그 결과를 Fig.
- 사용된 수치코드 FLUENT의 노즐 열유동장 해석에 관한 신뢰성을 확보하기 위해 해석코드 검증시 널리 사용되는 JPL 노즐 실험값[9]을 고려하였다. 공기의 전 압력, 전 온도를 각각 1.
이론/모형
- SST 모델은 두 모델의 장점을 살릴 수 있도록 벽 근처에서는 k-ω 모델이 사용되고 그 이외의 영역에서는 k-ε 모델이 사용된다.
- 마이크로 터보제트 엔진의 변형 노즐을 해석하기 위해서 형상 변형을 위한 기준 형상을 설정하였다. 기준 형상으로 Dassault Aviation에서 개발한 nEUROn의 노즐을 참고하였다. 엔진의 곡률을 변형하기 위해 nEUROn과 유사한 곡선을 가지는 hyperbolic tangent 함수를 적용하였다[8].
- 먼저 엔진 성능 분석을 통해 얻어진 원형 축소 노즐을 x축 방향으로 8등분을 하였다. 나누어진 총 9개의 단면적을 형상 변형 노즐에 똑같이 적용하기 위하여 hyperbolic tangent 함수와 타원 방정식을 이용하였다.
- 난류 해석을 위해 k-ε 모델과 k-ω 모델을 가중함수를 사용하여 결합한 SST 모델을 사용하였다.
- 노즐 내부 및 외부의 효과적인 해석을 위해서는 문제의 조건에 적당한 수치해석 기법을 적용하는 것이 중요하다. 본 연구에서는 열유동장 및 열전달을 해석하기 위해서 Navier-Stokes-Fourier 구성법칙에 기초한 질량, 운동량, 에너지 방정식, 화학종 방정식, 난류 방정식, 복사 열전달 방정식을 사용하였다. 상용 프로그램인 ANSYS FLUENT[7]를 사용하였으며, 주요 수치기법으로 밀도기반, 셀기반 유한 체적법, 내재 시간 전진법, 공간 이산화로 2차 풍상 차분법을 사용하였다.
- 본 연구에서는 열유동장 및 열전달을 해석하기 위해서 Navier-Stokes-Fourier 구성법칙에 기초한 질량, 운동량, 에너지 방정식, 화학종 방정식, 난류 방정식, 복사 열전달 방정식을 사용하였다. 상용 프로그램인 ANSYS FLUENT[7]를 사용하였으며, 주요 수치기법으로 밀도기반, 셀기반 유한 체적법, 내재 시간 전진법, 공간 이산화로 2차 풍상 차분법을 사용하였다. 플럭스 계산기법으로는 Riemann 근사해법에 기초한 Roe 기법을 사용하였다.
- 기준 형상으로 Dassault Aviation에서 개발한 nEUROn의 노즐을 참고하였다. 엔진의 곡률을 변형하기 위해 nEUROn과 유사한 곡선을 가지는 hyperbolic tangent 함수를 적용하였다[8]. 노즐의 단면적이 바뀌면 그 특성이 바뀌기 때문에, 최대한 같은 면적을 가지는 노즐을 이용하여 계산하기 위해, Fig 2의 원형 노즐을 기준으 로 하여 변형 노즐을 설계 하였다.
- 열유동장 해석결과와 narrow-band 모델을 사용하여 플룸의 IR 신호를 계산하였다. 각각의 AR과 곡률에 대하여 파장별 강도를 분석한 결과 AR이 커질수록 파장별 강도가 약하게 발생하는 것을 확인하였다.
- 플룸 및 후방동체에서 발생하는 파장별 복사신호는 항공기 추진기관 IR 신호분석에 있어서 중요하다. 이전 연구[8,10-12]에서 사용된 근사기법 narrow-band 모델을 이용하여 플룸 IR 신호를 계산하였다. 이 모델은 각 파장대에서의 복사 물성치를 평균을 내어 나타내는 밴드(band) 모델의 일종으로 많은 계산시간이 요구되는 line-by-line 밴드 모델보다 효율적 계산이 가능하다.
- 상용 프로그램인 ANSYS FLUENT[7]를 사용하였으며, 주요 수치기법으로 밀도기반, 셀기반 유한 체적법, 내재 시간 전진법, 공간 이산화로 2차 풍상 차분법을 사용하였다. 플럭스 계산기법으로는 Riemann 근사해법에 기초한 Roe 기법을 사용하였다. 난류 해석을 위해 k-ε 모델과 k-ω 모델을 가중함수를 사용하여 결합한 SST 모델을 사용하였다.
- 화학종의 질량 확산은 Schmidt 수(Sct = μt/ρDt)와 선형 Fick의 법칙을 사용하였다.
성능/효과
- IR 신호를 계산하고자 하는 지점에서 노즐 중심까지의 직선거리를 1000등분 하여 각 점에 대한 온도, 압력, 분압을 도출했다. IR 신호 계산 결과 같은 지점에서 계산 하였을 때, Fig. 9과 같이 AR이 커질수록 최대 파장별 강도 값이 감소하는 경향을 보였고, 곡률이 커질수록 파장별 강도 값이 감소하는 경향을 보였다. 또한 고도각, 방위각에 따라 파장별 강도의 값이 변화하는 것을 확인하였다.
- IR 신호의 결과 값들에서 공통적인 특징은 배기가스 조성으로 인해 CO2, H2O의 고유진동 수에 관련된 2.5-3.5 μm, 4-4.5 μm, 6-7 μm 에서의 파장별 강도 값이 증가한다는 것이다.
- 열유동장 해석결과와 narrow-band 모델을 사용하여 플룸의 IR 신호를 계산하였다. 각각의 AR과 곡률에 대하여 파장별 강도를 분석한 결과 AR이 커질수록 파장별 강도가 약하게 발생하는 것을 확인하였다. 그 이유로는 정 후방에서의 측정되는 고온 부분이 AR이 커질수록 최대 온도가 낮아지고 플룸 범위가 줄어들기 때문이다.
- 이를 통해 스케일은 다르지만 마이크로 터보제트 엔진과 nEUROn 열유동장의 특성이 세장비 및 AR 설계 변수에 따라 유사한 경향성을 보이는 것을 확인할 수 있다. 결과적으로 마이크로 터보제트 엔진 노즐의 열유동장 해석결과로 nEUROn의 열유동장 해석결과를 정성적으로 예측하는 것이 가능함을 의미한다고 할 수 있다.
- 한편 AR과 곡률이 증가할수록 온도 저감 효과가 있지만, 추력의 페널티도 함께 고려를 해야 한다. 계산 결과 원형 노즐 대비 최대 12%의 추력 감소를 확인 하였다. 이와 같은 결과를 통해 IR 신호 저감을 위한 노즐 설계 시 온도 저감 뿐만 아니라 추력의 페널티도 함께 고려를 해야 한다는 것을 알 수 있다.
- 한편 형상 변형 노즐의 온도 감소뿐만 아니라, 각종 페널티에 대한 영향도 같이 고려해야 하므로, 노즐출구의 압력, 속도, 질량 유량을 이용하여 각 노즐 형상별 추력을 계산하였다. 그 결과 형상 변형 노즐에서의 추력이 원형 노즐 대비 최대 12%의 감소를 보였으며, 추력의 감소는 가로세로비가 크거나 곡률이 클수록 더 크게 나타났다. 이와 같은 결과는 노즐 출구에서의 압력이 형상변형으로 인해 감소하고 노즐 출구 방향 변화로 인해 추력 방향 성분이 줄어들었기 때문으로 판단된다.
- 032 MPa, 550 K로 설정하였다. 그 결과를 Fig. 5에 나타내었는데, 노즐 표면과 중심에서의 압력분포가 실험값과 정성적으로 일치하는 것을 확인하였다. 이를 통해 본 연구에서 사용된 전산해석 기법의 신뢰성을 확보하였다.
- 9과 같이 AR이 커질수록 최대 파장별 강도 값이 감소하는 경향을 보였고, 곡률이 커질수록 파장별 강도 값이 감소하는 경향을 보였다. 또한 고도각, 방위각에 따라 파장별 강도의 값이 변화하는 것을 확인하였다. Fig.
- 먼저 같은 곡률을 적용한 경우 AR이 커질수록 플룸이 측면으로 넓게 형성 되고 플룸의 온도 또한 중앙이 아닌 측면에서 높게 나타남을 알 수 있다. 또한 곡률이 커질수록 고온부의 넓이가 축소되고 플룸의 유동이 y-방향으로 기우는 것을 확인할 수 있다. 이를 통해 스케일은 다르지만 마이크로 터보제트 엔진과 nEUROn 열유동장의 특성이 세장비 및 AR 설계 변수에 따라 유사한 경향성을 보이는 것을 확인할 수 있다.
- 전체적으로 무차원수의 값의 크기(order)는 큰 차이가 나지 않아 무차원수가 비슷한 크기로 유지되는 이상적인 상황에 근접하다고 할 수 있다. 또한 앞서 확인하였던 열유동장과 특정 밴드의 IR 신호에는 무차원수의 일부 차이가 그리 큰 영향을 끼치지 않음을 확인할 수 있다.
- 전산 해석결과 노즐출구에서 최대 33%, 노즐 벽면에서 최대 5%의 온도 저감 효과가 있는 것을 확인하였다. 또한 플룸의 온도 분포가 AR과 곡률이 증가할수록 고온부의 영역이 좁아지고 양 끝단으로 고온 부분이 집중되는 현상을 확인하였다. 이는 AR과 곡률이 증가하면 고온의 배기 플룸이 얇아지고 폭이 넓어지며 차가운 대기와의 혼합 효율이 훨씬 증가하여 생기는 현상이다.
- 모든 IR 신호의 결과를 살펴 보았을 때, 약 4-4.5 μm와 6-7 μm 밴드에서 파장별 강도가 다른 파장에 비해 두드러지게 나오는 것을 확인하였다.
- 주요한 다섯 개의 무차원수에 대한 결과를 full-scale 노즐과 마이크로 터보제트 엔진 열유동장에 대해 Table 4에 요약하였다. 온도 및 광학적 특성에 관한 무차원 수인 Bo은 온도에 민감하며, 플룸의 광학적 특성에 관한 NR은 상대적인 크기는 차이나지만 그 값이 매우 작음을 확인하였다. 전체적으로 무차원수의 값의 크기(order)는 큰 차이가 나지 않아 무차원수가 비슷한 크기로 유지되는 이상적인 상황에 근접하다고 할 수 있다.
- Scale의 크기와 IR 신호의 관측 거리 등이 서로 다르지만, CO2, H2O의 특성은 동일하게 작용을 하고, AR과 곡률이 증가 할수록 파장별 강도의 값이 낮게 나타나는 특성을 보여 준다. 이를 통해 마이크로 터보제트 엔진 결과로부터 full-scale의 결과를 정성적으로 유사하게 예측하는 것이 가능함을 보여준다.
- 또한 곡률이 커질수록 고온부의 넓이가 축소되고 플룸의 유동이 y-방향으로 기우는 것을 확인할 수 있다. 이를 통해 스케일은 다르지만 마이크로 터보제트 엔진과 nEUROn 열유동장의 특성이 세장비 및 AR 설계 변수에 따라 유사한 경향성을 보이는 것을 확인할 수 있다. 결과적으로 마이크로 터보제트 엔진 노즐의 열유동장 해석결과로 nEUROn의 열유동장 해석결과를 정성적으로 예측하는 것이 가능함을 의미한다고 할 수 있다.
- 전산 해석결과 노즐출구에서 최대 33%, 노즐 벽면에서 최대 5%의 온도 저감 효과가 있는 것을 확인하였다. 또한 플룸의 온도 분포가 AR과 곡률이 증가할수록 고온부의 영역이 좁아지고 양 끝단으로 고온 부분이 집중되는 현상을 확인하였다.
- 원칙적으로 이러한 무차원수들은 비슷한 크기로 유지하는 것이 이상적이나, 이를 달성하기에는 현실적인 제한이 많다. 축소 마이크로 터보제트 엔진 및 full-scale 노즐 모델에 관한 전산해석 결과를 비교해 보았을 때, 전체적으로 무차원수의 값의 크기에는 두드러진 차이가 나지 않아 두 모델 사이의 상사성이 존재함을 알 수 있었다. 특히 이 무차원 상수들의 일부 차이에도 불구하고, IR 신호 파장대 특성과 세장비 및 곡률에 따른 저감 또는 최적 경향성은 거의 동일하게 유지되어 지상 축소 모델을 통한 IR 신호 저감 설계 개념이 full-scale 모델에 적용이 가능함을 확인하였다.
- 축소 마이크로 터보제트 엔진 및 full-scale 노즐 모델에 관한 전산해석 결과를 비교해 보았을 때, 전체적으로 무차원수의 값의 크기에는 두드러진 차이가 나지 않아 두 모델 사이의 상사성이 존재함을 알 수 있었다. 특히 이 무차원 상수들의 일부 차이에도 불구하고, IR 신호 파장대 특성과 세장비 및 곡률에 따른 저감 또는 최적 경향성은 거의 동일하게 유지되어 지상 축소 모델을 통한 IR 신호 저감 설계 개념이 full-scale 모델에 적용이 가능함을 확인하였다.
- 외부 유동장의 크기는 노즐 입구 지름 D를 기준으로 반경 방향 30배, 축 방향 50배로 설정하였다. 플룸이 형성될 것이 예상되는 부분에 격자의 밀집도와 갯수를 증가시켜서 해석 결과의 정확도를 높였다. Fig.
- 해석 결과 원형 노즐에 비하여 형상 변형 노즐의 출구 평균온도가 최대 33% 감소하였고, 노즐 벽면 근처의 평균온도가 최대 5% 감소하였다. 온도 변화에 관한 해석결과는 아래 Table 3에 나타내었다.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.