입자 열관성에 의한 입자 라플라스 베르누이 시스템 대류 발현 지연

초록

본 연구는 입자 열관성, 즉 입자와 유체 사이의 비특이 열용량 비율 ε가 입자 라플라스‑베르누이(pRB) 시스템의 선형 안정성에 미치는 영향을 분석한다. ε가 증가할수록 임계 레일리 수(Raₙ) 가 상승하여 대류가 억제되고, ε≈1 정도에서 효과가 포화한다. 이는 입자와 유체 간 열교환이 기반 온도 구배를 완만하게 만들어 부력 구동을 약화시키기 때문이다. 무거운 입자와 가벼운 입자 모두, 유량·속도·주입 온도 변화에도 동일한 안정화 경향이 관찰된다.

상세 분석

본 논문은 희석된 입자 현탁액을 Eulerian 2‑fluid 프레임워크로 기술하고, 연속 방정식, 운동량 방정식, 에너지 방정식을 입자와 유체 각각에 적용한다. 입자에 대한 질량 보존식(2)과 입자 속도 방정식(4)에서는 Stokes drag, added‑mass, buoyancy를 포함하고, 입자 온도 방정식(6)에서는 열관성 시간을 τ_T = (3‑β)/(2β) ε d_p²/(12 κ) 로 정의한다. 여기서 β는 수정된 밀도비, ε = c_p^p / c_p^f 는 입자와 유체의 비특이 열용량 비율이다. 차원화 과정에서 Rayleigh 수 Ra, Prandtl 수 Pr, Galileo 수 Ga, 입자 크기 비 Φ 등 9개의 무차원 파라미터가 도출된다.

기저 상태는 무류동( U=0 )이며, 입자는 일정한 체적 유량 J 로 한쪽 벽에서 주입되고 반대쪽에서 즉시 배출된다. 열교환 항 α 12 Φ²(Θ‑Θ_p) 가 존재함에 따라 입자와 유체 온도 프로파일이 서로 영향을 주고받는다. ε가 작을수록 열교환이 빠르게 일어나 Θ≈Θ_p 가 되지만, ε가 커질수록 입자 온도가 유체 온도에 대한 응답이 지연되어 Θ와 Θ_p 사이에 유의한 차이가 발생한다. 이 차이는 특히 주입벽 근처에서 Θ의 기울기를 감소시켜, 부력에 의한 상승 흐름을 약화시킨다.

선형 안정성 분석에서는 정상 상태에 작은 3‑차원 파동형 섭동을 가정하고, 고유값 문제를 수치적으로 풀이한다. 결과는 ε가 증가함에 따라 가장 불안정한 모드의 성장률이 감소하고, 임계 Rayleigh 수 Ra_c 가 크게 상승함을 보여준다. 특히 ε≈1 에서 Ra_c 가 포화하는 현상이 관찰되는데, 이는 입자와 유체의 체적 열용량이 동등해져 추가적인 열관성 효과가 사라지기 때문이다. 무거운 입자(β<1)와 가벼운 입자(β>1) 모두에서 동일한 경향이 나타나며, 입자 부피 유량 J, 주입 속도 W*·, 주입 온도 Θ*_p 를 변화시켜도 ε 의 영향은 지배적이다.

기계적 결합(드래그, 부력)만을 고려한 기존 연구와 달리, 본 연구는 열적 결합을 명시적으로 포함함으로써 입자 열관성이 대류 억제에 핵심적인 역할을 함을 밝힌다. 특히 가벼운 입자는 기계적으로는 부력을 강화시켜 불안정성을 촉진할 수 있으나, ε 가 충분히 크면 열교환에 의해 온도 구배가 억제되어 전체 시스템은 오히려 안정화된다. 이러한 복합 효과는 입자‑유체 시스템 설계 시 열관성 파라미터를 조절함으로써 대류 억제 또는 촉진을 정밀하게 제어할 수 있는 새로운 설계 자유도를 제공한다.