Topology optimization of two-fluid turbulent heat exchangers: A Darcy flow-based multifidelity approach

본 논문은 두 유체 난류 열교환기 설계를 위해 계산 효율이 높은 다공성 매질 흐름 기반 저정밀도 모델을 고정밀도 RANS 모델과 보정하는 다정밀도 위상 최적화 프레임워크를 제시하며, 향상된 열전달과 관리 가능한 압력 강하를 균형 있게 조화시켜 기존 설계 대비 최대 22% 의 성능 향상을 달성합니다.

핵심

두 유체(예: 뜨거운 물과 차가운 물) 가 파이프를 통해 흐를 때 열적 악수를 하는 장치인 최고의 열교환기를 설계한다고 상상해 보세요. 목표는 유체가 통과하는 데 너무 많은 어려움을 겪지 않으면서 (이는 에너지 낭비로 이어질 수 있음) 가능한 한 빠르게 열을 교환하도록 하는 것입니다.

수십 년 동안 엔지니어들은 파이프 내부에 금속 리본을 꼬거나 핀을 추가하여 이러한 장치를 개선하려고 시도해 왔습니다. 하지만 이러한 방법들은 망치로 걸작을 조각하려는 것과 같습니다. 전통적인 제조 공정이 구부리고 꼬을 수 있는 것에 제한을 받기 때문입니다.

이 논문은 위상 최적화 (Topology Optimization) 라는 컴퓨터 "두뇌"를 사용하여 이러한 장치를 설계하는 새로운 방법을 소개합니다. 이는 파이프 내부에 들어맞는 한 상상할 수 있는 모든 형태를 조각해 낼 수 있는 디지털 조각가와 같습니다. 그러나 높은 속도에서 유체가 소용돌이치고 혼합되는 방식 (난류) 을 시뮬레이션하는 것은 허리케인 속에서 날씨를 예측하려는 것과 같습니다. 매우 정확하지만 슈퍼컴퓨터가 수년 동안 실행해야 할 정도로 시간이 많이 걸립니다.

문제: "완벽함" 대 "속도"

연구자들은 다음과 같은 딜레마에 직면했습니다:

1. 고정밀 (HF) 모델: 이는 "날씨 예보관"입니다. 복잡한 물리 (RANS 방정식) 를 사용하여 난류 유체의 거동을 정확히 예측합니다. 정확하지만 너무 느려서 최적의 설계를 찾기 위해 수천 번 실행하는 것은 불가능합니다.
2. 저정밀 (LF) 모델: 이는 "빠른 스케치"입니다. 유체를 스펀지를 통과하는 것처럼 취급하는 단순화된 수학 모델 (다르시 흐름) 을 사용합니다. 매우 빠르지만, 특히 유체가 잃는 압력에 관한 세부 사항에서 종종 오류를 범합니다.

스케치만 사용하면 실제 압력 하에서 붕괴될 수 있는 아름다운 파이프를 설계할 수 있습니다. 날씨 예보관만 사용하면 설계를 결코 끝내지 못할 것입니다.

해결책: "다중 정밀도 (Multifidelity)" 접근법

저자들은 다중 정밀도 접근법이라고 부르는 교묘한 2 단계 전략을 개발했습니다. 마라톤 훈련과 비슷하다고 생각하세요:

1. 훈련 주파 (최적화): "빠른 스케치"(LF 모델) 를 사용하여 수천 번의 연습 경주를 합니다. 설계를 수정하고 속도를 변경하며 다양한 형태를 시도하여 유망한 후보를 찾습니다. 스케치가 빠르기 때문에 수백 가지의 다양한 "만약에" 시나리오를 빠르게 탐색할 수 있습니다.
2. 보정: 훈련 주파를 시작하기 전에 그들은 스케치를 "보정"했습니다. 표준 파이프에 대한 날씨 예보관의 결과와 스케치의 결과가 일치하도록 수학에서 스펀지의 밀도를 조정했습니다. 이로 인해 스케치가 훨씬 더 똑똑해졌습니다.
3. 경기 당일 (평가): 컴퓨터가 빠른 스케치를 사용하여 흥미로운 설계들을 찾은 후, 상위 후보들을 "날씨 예보관"(HF 모델) 에 각각 한 번씩 통과시켰습니다. 이것이 실제로 어떤 설계가 승리하는지 확인하는 최종적이고 정확한 테스트입니다.

그들이 발견한 것

그들은 유체가 매우 빠르게 이동하는 (난류) "이중 파이프" 열교환기에 이 방법을 적용했습니다.

• 결과: 컴퓨터가 설계한 형태는 야생적이고 복잡하여 표준 파이프와 전혀 닮지 않았습니다. 그들은 유체가 강하게 소용돌이치고 혼합되도록 강제하는 정교한 내부 벽을 만들었습니다. 이는 소스를 더 빨리 식히기 위해 요리사가 격렬하게 저어주는 것과 매우 비슷합니다.
• 비교: 그들은 새로운 설계들을 열전달을 개선하는 일반적인 업계 기법인 "꼬인 테이프"가 있는 표준 파이프와 비교했습니다.
  • 꼬인 테이프는 열전달을 개선했지만 거대한 "교통 체증"(높은 압력 강하) 을 유발하여 전체적으로 비효율적이었습니다.
  • 새로운 컴퓨터 설계된 형태는 일반 파이프에 비해 열전달을 최대 **66%**까지 개선했습니다.
  • 결정적으로 그들은 "교통 체증"을 훨씬 더 잘 관리했습니다. 전체 점수 (열 획득 대 에너지 비용의 균형) 를 볼 때, 그들의 설계는 꼬인 테이프보다 최대 22% 더 우수했습니다.

결론

이 논문은 훌륭한 설계를 찾기 위해 허리케인의 모든 소용돌이를 시뮬레이션할 필요가 없음을 증명합니다. 가능성을 탐색하기 위한 빠르고 보정된 "스케치"와 승자를 검증하기 위한 느리고 정확한 "예보관"을 사용함으로써 엔지니어들은 전통적인 방법으로 현재 구축할 수 있는 것보다 훨씬 우수한 고성능 열교환기를 설계할 수 있습니다.

이 연구는 특히 이러한 설계들이 다양한 속도 범위에서 잘 작동한다고 지적하며, 이는 실용적이고 제조 가능 (저자는 이러한 복잡한 형태를 가능하게 하는 핵심 요소로 3D 프린팅을 언급함) 한 경우 실제 사용에 준비되어 있음을 시사합니다.

기술 요약

기술 요약: 두 유체 난류 열교환기의 위상 최적화

문제 제기
이중 관 열교환기 (DPHX) 를 포함한 두 유체 열교환기 (HX) 의 설계는 전통적으로 꼬임 테이프 삽입, 핀, 또는 주름관과 같은 수동적 향상 기술에 의존해 왔습니다. 이러한 방법들은 효과적이지만, 기존 제조 공정의 한계에 제약받으며 열전달 향상과 압력 강하 패널티 사이의 비최적 균형을 초래하는 경우가 많습니다. 위상 최적화 (TO) 는 혁신적인 기하학적 구조를 생성할 수 있는 경로를 제공하지만, 난류 두 유체 시스템에 TO 를 적용하는 것은 상당한 계산적 도전을 제기합니다. 난류의 비선형성, 극도로 미세한 메쉬의 필요성, 그리고 수치적 불안정성으로 인해 고충실도 (HF) 난류 모델 (예: Reynolds-Averaged Navier-Stokes, RANS) 을 최적화 루프에 직접 통합하는 것은 계산적으로 불가능합니다. 또한, 기존 두 유체 HX 에 대한 TO 연구는 주로 층류 또는 낮은 레이놀즈 수 영역에 집중되어 있어, 난류 유동의 최적화는 largely 다루어지지 않았습니다.

방법론
저자들은 계산 효율성과 정확성 사이의 균형을 맞추기 위해 설계 탐색과 고충실도 평가를 분리하는 다중 충실도 위상 최적화 프레임워크를 제안합니다.

1. 저충실도 (LF) 최적화 모델:

   • 최적화 과정은 완전한 난류 모델 대신 다르시 (Darcy) 유동 기반 모델을 활용합니다. 이 접근법은 "가난한 자의 접근법 (poor man's approach)"으로 불리며, Navier-Stokes 방정식을 선형적이고 계산 비용이 저렴하며 수치적으로 안정적인 다르시 방정식으로 대체합니다.
   • 단일 설계 변수 필드 ($\gamma$) 가 유체 1, 유체 2, 그리고 고체 벽의 세 가지 상의 분포를 나타냅니다.
   • 보간 기법: 각 유체의 투과도는 RAMP 함수를 사용하여 보간되어, 고체 벽 영역 (중간 $\gamma$ 값) 에서 유체 혼합을 방지하기 위해 투과도가 충분히 낮아지도록 합니다. 열전도도는 유체와 고체 특성 사이의 매끄러운 전환을 위해 가우스 함수를 사용하여 보간됩니다.
   • 보정: 단순화된 다르시 모델과 현실 간의 불일치를 완화하기 위해, LF 모델 매개변수 (최대 투과도 및 열전도도 경도) 는 직관형 파이프 구성에 대한 참조 HF 모델을 기준으로 보정됩니다. 이 보정을 통해 LF 모델의 압력 강하 및 온도 차이 예측이 HF 모델과 일치하도록 조정됩니다.

2. 다중 충실도 전략:

   • 단일 최적화를 실행하는 대신, 프레임워크는 시딩 (seeding) 매개변수 전략을 사용합니다. 보정된 LF 모델을 사용하여 입구 속도와 트레이드오프 매개변수 (열전달 대 압력 강하 균형) 를 변화시키며 여러 독립적인 최적화를 수행하여 다양한 후보 설계 세트를 생성합니다.
   • 고충실도 (HF) 평가: LF 과정에서 추출된 최적화된 설계들은 RANS 기반 HF 모델 (특히 벽 함수를 사용한 $k-\omega$ 난류 모델) 을 사용하여 평가됩니다. 이 단계는 최적화 루프 외부에서 수행됩니다.
   • 선정: HF 모델은 전체 열전달 계수를 압력 강하 패널티로 정규화하는 성능 평가 기준 (PEC) 을 기반으로 모든 후보의 성능을 평가합니다. 가장 우수한 성능을 보이는 설계는 이러한 정확한 평가를 바탕으로 선정됩니다.

3. 수치적 구현:

   • LF 모델은 밀도 필터와 헤비사이드 프로젝션을 사용하여 명확한 유체 - 고체 경계를 보장하는 COMSOL Multiphysics 를 사용하여 해결됩니다.
   • HF 평가는 최적화된 설계 변수에서 생성된 바디 피팅 (body-fitted) 메쉬를 활용하며, 난류에 대해서는 $k-\omega$ 모델을 사용하고 결과를 경험적 상관관계 (꼬임 테이프에 대한 Dittus-Boelter, Blasius, Manglik & Bergles) 와 비교하여 검증합니다.

주요 기여

• 프레임워크 개발: 본 연구는 주로 층류 유동이나 단일 유체 시스템에 초점을 맞춘 기존 문헌의 공백을 해소하는 두 유체 난류 열교환기를 위한 다중 충실도 TO 프레임워크를 확립합니다.
• 보정 전략: 두 유체 난류 시나리오에서 다르시 기반 LF 모델과 RANS 기반 HF 응답을 정렬하기 위한 새로운 보정 전략이 도입되어, 압력 강하 및 온도 예측의 일관성을 향상시켰습니다.
• 정량적 검증: 논문은 광범위한 레이놀즈 수 (최대 13,000) 에 걸쳐 최적화된 설계와 기존 참조 설계 (꼬임 테이프 삽입이 있는 매끄러운 파이프) 간의 엄격한 정량적 비교를 제공합니다.

결과

• 설계 생성: LF 모델은 유체 혼합을 촉진하고 열교환 표면적을 증가시키면서도 유선형 흐름 경로를 유지하는 다양하고 복잡한 벽 구조를 성공적으로 생성했습니다.
• 성능 향상: 매끄러운 파이프 참조 설계와 비교하여, 최적화된 설계는 전체 열전달 계수를 최대 **66.7%**까지 향상시켰습니다.
• PEC 향상: 꼬임 테이프 참조 설계 (표준 향상 기술) 와 비교할 때, 최적화된 설계는 성능 평가 기준 (PEC) 값을 최대 **22%**까지 높였습니다. 이는 최적화된 기하학적 구조가 압력 강하 패널티에 비해 우수한 열전달 향상을 제공함을 나타냅니다.
• 모델 불일치: LF 모델은 온도 분포를 합리적으로 예측하는 것으로 확인되었으나, 점성 전단 및 난류 모델링의 부재로 인해 압력 강하를 과소평가하는 경향이 있었습니다 (유체 1 의 경우 평균 오차 약 72%). 그러나 다중 충실도 접근법은 HF 모델을 통해 후보를 필터링함으로써 이를 효과적으로 보완했습니다.
• 강건성: 더 높은 레이놀즈 수 (예: Re = 13,000) 에서 최적화된 설계는 광범위한 작동 조건 전반에 걸쳐 높은 PEC 값을 유지하며 강건성을 입증했습니다.

의의 및 주장
본 논문은 제안된 다중 충실도 접근법이 난류 조건 하에서 두 유체 열교환기를 최적화하는 것과 관련된 계산적 장벽을 효과적으로 극복한다고 주장합니다. 탐색 단계에는 계산 비용이 저렴한 다르시 모델을, 평가 단계에는 고충실도 RANS 모델을 활용함으로써, 난류 모델을 사용한 직접 최적화로는 발견이 불가능했을 고성능 구성을 찾을 수 있게 합니다.

저자들은 최적화된 구성이 열전달과 압력 손실 간의 트레이드오프 측면에서 특히 기존 꼬임 테이프 삽입보다 우수하다고 강조합니다. 연구는 다르시 기반 LF 모델이 절대적인 압력 강하 예측에는 한계가 있지만, 다중 충실도 프레임워크가 고성능 이중 관 열교환기 설계를 위한 실현 가능하고 효과적인 전략임을 결론지으며, 복잡한 난류 열 - 유체 문제에 대한 이 접근법의 잠재력을 보여줍니다. 이 작업은 시딩 매개변수의 수동 선택과 현재 연구에서의 제조 제약 조건 부재와 같은 한계를 인정하며, 데이터 기반 자동화와 제조 가능성 제약 조건의 포함을 향한 향후 방향을 제안합니다.