최근 선박으로부터 배출되는 입자상 물질(Particulate Matters)이 지구의 온난화 효과는 물론 북극해의 해빙에도 영향을 미친다는 보고가 나오고 있다[1]. 특히, IPCC의 기후변화 보고서[2]에 의하면 PM 특히 블랙카본의 온난화지수는 0.2W/m² 으로 보고하고 있으며, 미국의 연구진은 0.9W/m²의 한 차수가 큰 값까지도 제시하고 있다[3]. 이로 인해 최근 국내에서는 선급을 비롯한 산학연 협동으로 선박으로부터의 매연이 온난화에 미치는 영향 및 선박에서 발생하는 매연의 기본적인 구조 등을 조사하기 위한 연구가 진행 중이다.

국제해사기구(International Maritime Organization, IMO)에서도 선박으로부터의 PM을 비롯한 매연(soot)의 배출규제를 강화하려고 하고 있다[4].

일반적으로 선박에서 배출되는 PM을 제어하기 위해서는 매연의 생성에 영향을 주는 연료, 온도 및 체류시간 3대 인자의 영향에 대하여 알아볼 필요가 있다. 이에 본 연구에서는 가스상 에틸렌 확산화염을 이용하여 매연의 생성 특성에 알아보고자 한다.

매연은 일반적으로 오염 물질로만 잘 알려져 있지만, 화염 내에서는 중요한 복사 매체로서의 역할을 한다. 이로 인해 매연 생성에 관한 연구는 오랜 기간 동안 많은 연구자에 의해서 연구되고 있다[5]-[7]. Ku et al. [5][6]은 지상중력 및 마이크로중력 화염 내에서 매연 생성에 대한 부력의 영향에 대하여 연구를 수행하였다. Atreya et al. [7]은 실험과 수치해석을 통해 마이크로중력 환경 내에서는 복사 열전달이 매연 생성에 중요하다는 것을 파악하였다. 특히 최근에는 부력에 따른 매연 생성의 연구가 활발히 이루어지고 있다[8]-[10]. Fujita et al. [8]은 지상중력과 마이크로중력환경에서의 실험을 통해 부력의 유무에 따른 열영동력(Thermophoretic force)이 화염 내에서 매연 생성에 미치는 영향에 대하여 조사하였다. Hermanson et al. [9]는 층류, 천이 및 난류 제트 확산 화염 내에서의 매연 생성에 대한 부력의 효과에 대하여 화염의 높이, 모양 및 산화제의 효과에 대한 다양한 연구 결과를 제시하였다. Jeon et al. [10]은 6m의 낙하탑(Drop tower)을 이용하여 다양한 부분 중력(Partial gravity) 조건에서 확산 화염 내의 매연 생성 특성에 대한 중력의 영향을 조사하였다.

본 연구는 이러한 연구들의 일환으로 매연 생성에 대한 부력의 영향에 대하여 알아보기 위하여 수치해석을 수행하였다. 특히 1g의 지상중력과 0g의 무중력에서의 에틸렌 확산 화염을 대상으로 하여 확산 화염 내 온도 및 속도의 유동 특성과 매연 분포를 조사하였다. 본 연구의 수치해석은 부력과 비부력 화염에서의 매연 특성의 차이점에 대하여 알아보기 위하여 1g와 0g 조건에서의 물리적 및 화학적 과정의 상호작용에 대하여도 고려하여 수행되었다. 또한 확산 화염 내에서 매연 생성 특성 및 탄소나노물질의 합성 연구 등에 있어 이들의 채취를 위한 적정 위치에 대해서 고찰하였다.

본 연구의 대상인 동축류 에틸렌 확산화염에 대하여 수치 시뮬레이션을 하기 위해서 2차원의 비정상계산을 수행하였다. 지배 방정식은 질량, 운동량, 에너지, 화학종 성분 및 기체 상태 방정식이며, 반응속도정수는 Arrhenius형 반응식을 이용하여 계산을 수행하였다[11]. 여기서 1차 총괄 반응을 위한 계수로서 선지수 함수(pre-exponential factor)는 5.8E+12 g/cm³s 이며, 활성화 에너지(Activation energy)는 125.6 kJ/mol을 적용하였다. 무중력상태의 조건을 맞추기 위해서는 운동량 방정식에서 중력 부분을 제외하여 시뮬레이션을 수행하였다. 화학반응은 1차의 총괄반응이며 따라서 낮은 마하수 가정(Low mach number approximation)을 적용하였고, 계산의 단순화를 위해서 Soret-Dufour 효과, 압력구배에 의한 확산효과는 무시하였다. 벽면 조건은 버너 부분을 제외한 벽면에는 미끄럼이 존재하지 않는 No-slip　condition을, 그리고 출구 조건은 대류경계 조건(Convective boundary condition)을 적용하였다. 계산에 필요한 열역학적 물성치들과 열전달계수는 Chemkin Ⅲ와 Transport package를 이용하여 해석하였다. 수치 해석을 위한 지배방정식 및 상세한 수치 방법은 선행 연구인 참고 문헌[12]을 참조 하였다.

또한 본 연구에서는 매연의 수밀도(Number density) 및 체적 분율(Volume fraction) 등을 파악하기 위하여 Moss et al. [13]에 의해 제안된 매연 지배 방정식을 도입하여 이에 대해서도 계산이 이루어졌다. 이 지배 방정식은 매연의 핵 생성(Nucleation), 표면 성장(Surface growth), 응집(Coagulation) 및 산화(Oxidation)에 대한 소스항이 포함되어 있다. 아래에 매연 수밀도(1)와 체적분율(2)에 대한 방정식을 나타낸다.

이 식에서 n은 매연 수밀도, f_v는 매연 체적분율, v_t,xi는 매연의 열영동속도를 나타낸다. ω_n은 매연 수밀도의 소스항이며 ω_fv는 체적 분율의 소스항이다. 소스항 ω_n과 ω_fv에 대한 상세한 방정식 및 값 들에 대한 내용은 선행 연구[14][15]를 참조하였다.

또한 방정식 (1)과 (2)는 대류와 열영동에 대한 항을 포함하며 특히 열영동 속도에 대해서는 아래의 식으로 정의되는 값을 취하였다[16].

여기서 ρ은 밀도, μ는 점성계수, ᐁT는 온도구배, T는 온도를 나타낸다. 열영동 효과에 대해서 본 연구에서는 반경방향의 r(속도 u)과 축방향 z (속도 v) 모두에 대하여 적용하였다.

복사의 고려에 대해서는 에너지 방정식에 매연, CO₂와 H₂O에 의한 복사 효과를 포함시켰으며, 방정식 (4)의 식으로 나타내었다.

여기서, q_r은 열복사 플럭스, σ는 스테판 볼츠만 상수 및 a_overall은 전체 플랭크 평균 흡수 계수(Overall Planck mean absorption coefficient)를 나타낸다.

이러한 전체 플랭크 평균 흡수 계수는 매연에 대한 값에 CO₂와 H₂O 혼합기체의 계수의 값을 모두 합한 값을 사용하였다[15].

Figure 1은 수치 시뮬레이션을 위한 모델 개략도를 나타낸다. 그림에서는 축대칭좌표를 이용하여 계산하였으므로 그림에서 점선으로 표시된 부분(Calculation zone)에 대해서만 해석을 수행하였다. 연료는 에틸렌(C₂H₄)를 사용하였으며, 연료와 산화제의 유속은 동일하게 5.1 cm/s 이다. 버너의 내경은 7mm이나 실제의 계산은 반경인 3.5mm에 대해서만 계산을 수행하였다. 그림에서 반경방향(r) 속도가 u, 축방향(z)의 속도는 v이다. 전체 계산 영역은 3cm(r) × 10cm(z)이다. 격자수는 128(r) × 512(z)이며 균일 간격으로 비정상 계산을 수행하였다.

Figure 1:  
Schematic diagram of computational domain

Figure 2는 지상중력인 1g에서의 동축류 에틸렌 확산화염의 온도 분포를 나타낸다. 비정상 상태의 화염온도분포의 변화를 나타낸다. 임의의 시간에서 왼쪽에서부터 오른쪽으로 0.1초 간격으로 변화하고 있다. 전체적인 최고 화염온도는 2152 K 정도를 나타낸다. 전체적인 화염의 형상은 화염의 온도와 주위공기의 밀도차에 의한 부력의 영향으로 축방향으로 길게, 반경방향으로는 폭이 좁은 것을 알 수 있다. 또한 흥미로운 점은 그림에서 알 수 있듯이 본 연구의 대상인 화염은 지상 중력인 1g에서 어떤 시간 간격을 두고 불안정한 형태로 진동이 발생하고 있는데 순간적으로 화염의 크기가 줄어들면서 다시 커지는 것을 알 수 있다.

이러한 조건에서는 확산 화염 내에서 생성되는 매연의 흐름의 관찰을 비롯해 매연의 생성 및 채취를 하기에는 매우 불리한 여건이다.

또한 최근 많이 이루어지는 연구[17] 중 화염을 이용한 확산 화염 내 탄소나노튜브의 채취를 하기에는 매우 불리한 조건으로 판단된다.

Figure 2: 
Temperature distribution of coflow ethylene diffusion flames in 1g

Figure 3은 지상중력인 1g에서 마이크로 중력으로 변화하는 동축류 에틸렌 확산화염의 온도분포를 나타낸다. 비정상 상태의 화염온도분포의 변화를 나타낸다. 임의의 시간에서 왼쪽에서부터 오른쪽으로 0.1초 간격으로 변화하고 있다. 그림에서 알 수 있듯이 1g에서 0g로 중력이 변화함에 따라서 화염온도는 2150 K에서 조금씩 계속적으로 낮아져 2123K를 나타내고 있다. 실제적으로 무중력으로 변화함으로써 화염의 온도가 많이 낮아질 것으로 예상되었으나 본 연구에서는 동축류를 대상으로 하고 있으므로 주위에서 충분한 산화제를 공급하고 있기 때문에 예상만큼 화염온도의 급격한 저하는 이루어지지 않았다. 반면에, Figure 4에 무중력 환경 하에서 정지 분위기 상태에서에 온도분포를 나타내고 있는데 화염 온도가 많이 낮아진 것을 알 수 있다.

Figure 2에서 나타내듯이 1g에서의 불안정한 화염의 상태는 시간이 지날수록 안정된 화염이 되어 가고 있는 것을 알 수 있다. 1g에서 0g로 순간적인 중력의 변화가 일어날 때는 화염의 축방향 길이가 작아졌다가 무중력 상태에서 순차적으로 화염이 커지고 있으며, 화염의 반경은 1cm 정도까지 커지고 있는 것을 알 수 있다. 이는 무중력 상태에 들어서면서 부력이 없어짐으로 인한 화염 내로의 산화제 공급량이 적어져 화염이 주위의 산화제를 공급 받기 위해 커지는 것이다.

Figure 4는 무중력 환경하의 주위 산화제 속도가 0cm/s의 정지 분위기에서의 동축류 에틸렌 확산화염의 온도분포를 나타낸다. 앞서 설명한 바와 같이 화염의 온도는 2079K으로 1g에서의 동축류 에틸렌 확산화염보다 약 70K 정도 낮아진 것을 알 수 있다. 또한 화염의 형상은 축방향으로는 길이가 짧아졌으며, 반경방향으로의 화염 반경은 1g 에서의 화염 반경이 약 0.5cm 내외인 반면, 주위 산화제의 속도가 0cm/s의 정지 분위기 상태에서는 약 1.5cm 까지 커지고 있다는 것을 알 수 있다. 이것은 주위 유속이 0cm/s 이므로 부력의 영향은 물론 주위 유속에 의한 영향도 없으므로 순수한 연료의 유량에 의하여 화염의 길이가 결정되기 때문에 축방향으로는 짧아지고 반경으로는 커진다고 생각한다. 이러한 조건에서는 축방향 속도도 낮아져 화염 내체류 시간도 더 길어질 것으로 예상되며, 이로 인해 매연의 생성도 많아 질 것으로 생각한다.

Figure 3: 
Temperature distribution of coflow ethylene diffusion flames in 0g

Figure 4: 
Temperature distribution in 0g (oxidizer velocity is 0cm/s)

Figure 5는 지상중력(1g)과 무중력(0g)에서의 동축류 에틸렌 확산화염의 속도 분포와 정지 분위기에서의 확산화염의 속도 분포를 나타낸다. 그림에서 나타난 바와 같이 지상중력인 1g에서의 화염 내 최대 속도는 208cm/s의 매우 빠른 속도를 나타내고 있는 반면, 마이크로중력 환경에서 주위산화 제의 속도가 5.1cm인 경우 확산화염 내 최대 속도는 18.9cm/s로 11배나 느린 속도분포를 나타낸다. 더욱이 동축류가 없는 상태의 정지된 산화제 분위기에서는 최대 속도가 14.2cm/s를 나타내고 있다.

이러한 결과는 일반적으로 실제 화염 내에서 매연의 흐름의 관찰이나 매연의 채취를 위해서는 가능한 한 화염 내 속도가 느린 것이 유리하다고 할 때 1g에서의 관찰이나 매연의 채취는 매우 어려운 환경이 될 수 있다는 것을 의미하며, 반대로 무중력 또는 마이크로 중력에서는 유리한 조건이라는 것을 의미한다. 또한 화염 내에서의 유속이 느리다는 것은 화염 내에서 매연의 체류시간이 길어질 것이라는 것을 의미하며, 이는 매연의 생성량이 증가 할 수 있다는 것을 나타낸다. 이 결과는 무중력 또는 마이크로 중력에서의 화염은 화염 내 많은 양의 매연이 긴 체류시간을 가져 복사열 방출이 많아져 화염의 온도는 낮아질 것이라는 나타낸다.

더불어 탄소나노물질의 생성 및 채취에 있어서는 체류 시간이 그만큼 길어지고 그에 따라 탄소나노물질의 생성되는 량도 많아질 가능성이 있다는 것을 의미한다. 요약하면 지상중력 환경의 화염 내에서는 매연과 탄소나노물질의 채취에는 많은 제약이 있다는 것을 나타내고 있으며 마이크로 중력 환경에서는 그만큼 매연 및 탄소나노물질의 생성 메커니즘을 비롯한 여러 가지 물리적 현상을 파악하는데 매우 유리하다는 것을 보여주고 있다.

Figure 5:  
Velocity distribution in 1g and 0g

Figure 6: 
Number density and volume fraction of soot particles in coflow diffusion flames in 1g and 0g

Figure 7: 
Distributions of soot volume fraction, velocity and temperature in coflow diffusion flame under 0g

Figure 6은 지상중력 1g와 무중력 0g에서 동축류 에틸렌 확산화염 내의 매연 수밀도와 매연 체적 분율을 나타낸다. Fig. 6에서 나타낸 결과는 앞서 설명한 매연 생성과 관련한 방정식(1)~(4)를 이용하였다. 결과에서 알 수 있듯이 지상중력에서의 수밀도는 무중력의 수밀도 보다 큰 값을 나타내고 있다. 그에 반해 매연의 체적 분율은 지상중력 환경에서는 최고 체적분율 값이 6.788×10^-7을 나타내고, 무중력 환경에서는 최고 체적분율 값이 1.005×10^-6이다. 무중력 환경에서의 매연 체적 분율이 약 1.48배 증가하였다는 것을 알 수 있다. 이러한 현상은 앞서 서술한 바와 같이 무중력 환경에서 축방향 속도의 감소로 인하여 화염 내에서 매연 생성을 인한 충분한 체류시간을 가지기 때문이다. 또한 무중력 환경에서 보다 큰 매연 체적분율 값을 나타내는 것은 지상중력환경에서는 빠른 유속으로 인해 1차 입자의 작은 매연들이 많이 분포한 반면 무중력에서는 느린 유속으로 인해 화염 내 생성되는 매연들이 응결을 위한 상대적으로 충분한 시간적 여유가 있기 때문에 직경이 큰 매연들이 많이 분포하게 된다.

이러한 결과는 확산 화염 내 생성되는 매연의 채취와 탄소나노물질의 채취에 있어 각각의 적정 위치를 파악하는데 있어 매우 중요한 수치 데이터로 활용이 가능할 것으로 사료된다.

Figure 7은 무중력 환경에서 동축류 에틸렌 확산화염 내의 매연체적분율, 속도분포 및 온도분포를 나타낸다. 앞서 나온 결과들을 요약한 것으로 동축류 확산 화염 내에서 탄소나노물질 등의 채취를 위한 최적의 위치에 대해서 나타내었다. 일반적으로 무중력 환경하에서 확산화염 내 탄소나노튜브의 채취를 위한 최적의 위치로는 그림에서 나타난 바와 같이 확산 화염 내 매연의 분포가 많고 유속이 느리며 탄소나노물질의 합성이 가능한 1400K 전후의 온도 분포를 가진 지역이어야 하는데 그 위치를 나타내면 온도분포 그림의 점선 사각형 범위로 생각할 수 있다.

본 연구에서는 1g의 지상중력과 0g의 무중력에서의 부력의 영향에 의한 에틸렌 확산 화염 내 매연의 생성 특성을 파악하기 위하여 수치 시뮬레이션을 통하여 확산 화염 내 온도 및 속도의 유동 특성과 매연 분포를 조사하였다.

그 결과 동축류 에틸렌 확산화염의 경우 지상 중력인 1g에서는 부력으로 인해 불안정한 화염 형상을 나타내고 있는 반면, 무중력인 0g에서는 안정된 화염 형상을 나타내었다. 화염의 온도는 무중력 환경에서 동축류 확산화염의 경우가 30K 정도 낮아졌으며, 확산화염 내 최대 속도는 11배나 느린 속도분포를 나타낸다. 이러한 결과는 확산 화염내의 느린 유속으로 인해 체류 시간의 증가로 인해 매연 생성이 유리하다는 것을 의미한다.

또한 매연의 농도를 조사한 결과 지상중력 환경에서는 최고 매연 체적분율 값이 6.788×10^-7인 반면, 무중력 환경에서는 최고 매연 체적분율 값이 1.005×10^-6로 무중력 환경에서의 매연 체적 분율이 약 1.48배 증가하였다. 무중력 환경에서 보다 큰 매연 체적분율 값을 나타낸다는 것은 지상중력환경에서는 1차 입자의 작은 매연들이 많이 분포된 반면 무중력에서는 축방향 속도의 감소로 인한 화염 내 매연의 충분한 체류 시간으로 인해 화염 내 매연들의 활발한 응결 작용으로 인해 직경이 큰 매연들이 많이 분포하고 있다는 것을 의미한다. 또한 확산 화염 내에서의 매연 및 탄소나노물질의 채취 및 생성 특성의 파악을 위해서는 무중력 및 마이크로중력 환경을 이용한 확산 화염 내에서 이루어지는 것이 유리하다고 판단된다.