[논문 리뷰] Refined radial basis function-generated finite difference analysis of non-Newtonian natural convection

Miha Rot, Gregor Kosec|arXiv (Cornell University)|2022. 02. 16.

Numerical methods in engineering인용 수 7

한 줄 요약

이 논문은 2차원 및 3차원 복잡한 기하구조에서 비뉴턴 대류를 시뮬레이션하기 위해 개선된 라디얼 기저 함수 기반 유한차분(RBF-FD) 메쉬리스 방법을 제안한다. 산점된 노드를 사용하고 국소 기울기에 기반한 조각별 선형 노드 밀도 함수를 활용한 적응형 균열을 통해, 정규 메쉬화 대비 동일한 정확도를 확보하면서도 최대 5배 적은 수의 노드를 사용할 수 있으며, 인위적 확산을 줄이기 위해 수치적 안정화 기법을 사용하지 않는다.

ABSTRACT

Code and data for https://doi.org/10.1063/5.0257896 published in Physics of Fluids 37, 033130 (2025) Refined radial basis function-generated finite difference analysis of non-Newtonian natural convection In this paper, we present a refined radial basis function-generated finite difference solution for a non-Newtonian fluid in a closed differentially heated cavity. The non-Newtonian behavior is modeled with the Ostwald–de Waele power law and the buoyancy with the Boussinesq approximation. The problem domain is discretized with scattered nodes without any requirement for a topological relation between them. This allows a trivial generalization of the solution procedure to complex irregular three dimensional (3D) domains, which is also demonstrated by solving the problem in a two dimensional (2D) and 3D geometry mimicking a porous filter. The results in 2D are compared with two reference solutions that use the finite volume method in a conjunction with two different stabilization techniques, where we achieved good agreement with the reference data. The refinement is implemented on top of a dedicated meshless node positioning algorithm using piecewise linear node density function that ensures sufficient node density in the center of the domain while maximizing the node density in a boundary layer where the most intense dynamic is expected. The results show that with a refined approach, more than five times fewer nodes are required to obtain the results with the same accuracy compared to the regular discretization. The paper also discusses the convergence with refined discretization for different scenarios for up to 2e5 nodes, the impact of method parameters, the behavior of the flow in the boundary layer, the behavior of the viscosity, and the geometric flexibility of the proposed solution procedure.

연구 동기 및 목표

• 복잡한 기하구조에서 비뉴턴 대류를 위한 차원에 관계없고 메쉬리스 수치적 접근법을 개발하기 위해.
• 계산 영역에서 메쉬 생성 및 위상적 연결성의 필요성을 제거하기 위해.
• 상향 또는 QUICK과 같은 안정화 기법을 피하여 수치적 확산을 줄이기 위해.
• 국소 스텐실 및 기저 함수 제어를 통해 효율적인 hp-적응 가능성을 제공하기 위해.
• 기본 2차원 및 복잡한 3차원 다공성 필터 유사 영역에서의 정확도와 기하 유연성을 입증하기 위해.

제안 방법

• 이 방법은 위상적 제약 없이 산점된 비정규 메쉬 노드를 사용하는 개선된 RBF-FD 이산화를 적용한다.
• 전용 메쉬리스 노드 위치 결정 알고리즘이 조각별 선형 노드 밀도 함수를 사용하여 경계층 및 고기울기 영역에 노드를 집중시킨다.
• 오스탈트-데와일 전력법 모델과 부우신스끄 근사법을 적용한 라우지-니비에르 방정식을 창킨의 투영 방법을 통해 압력-속도 결합을 해결한다.
• 명시적 시간 스텝을 사용하며, 물리적 정확성을 유지하기 위해 인위적 안정화 기법(예: 상향 또는 QUICK)을 적용하지 않는다.
• 각 노드에서 스텐실 크기, 다항식 보정 순서, RBF 유형을 변경 가능하게 하여 hp-적응 가능성을 지원한다.
• 수식은 차원에 민감하지 않아, 최소한의 코드 수정으로 2차원에서 3차원으로의 원활한 확장이 가능하다.

실험 결과

연구 질문

• RQ1 개선된 RBF-FD 메쉬리스 방법은 수치적 안정화 없이도 비뉴턴 대류에서 높은 정확도를 달성할 수 있는가?
• RQ2 국소 기울기에 기반한 적응형 노드 균열이 균일한 노드 분포에 비해 효율성과 정확도를 어떻게 향상시키는가?
• RQ3 이 방법은 메쉬 없이도 다공성 필터와 같은 복잡한 3차원 기하구조를 얼마나 잘 처리할 수 있는가?
• RQ4 점성도가 비선형 점성 특성에 의해 변화함에 따라, 제한된 영역 내에서 유동 구조와 속도 프로파일은 어떻게 영향을 받는가?
• RQ5 확장성과 유지보수성을 고려하여 차원에 독립적이고 객체 지향적인 방식으로 이 방법을 구현할 수 있는가?

주요 결과

• 2차원 기준 사례에서 개선된 RBF-FD 방법은 정규 이산화 대비 동일한 정확도를 유지하면서 최대 5배 적은 수의 노드를 사용하였다.
• 2차원 사례 결과는 상향 및 QUICK 안정화를 사용한 기준 FVM 해와 뛰어난 일치를 보이며, 방법의 정확도를 검증하였다.
• 수치적 안정화의 부재로 인해 낮은 수치적 확산이 발생하여 더 선명한 유동 및 온도 기울기를 유지하였다.
• 3차원 다공성 필터 사례에서는 고속도 채널 및 벽과 좁은 통로에서 강한 점성 감소 효과를 성공적으로 포착하였다.
• 비선형 점성도 변화는 명확하게 시각화되었으며, 가장 낮은 점성도는 벽과 좁은 통로 내 고기울기 영역에서 관찰되었다.
• 이 방법은 강력한 기하 기하학적 탄력성을 보이며, 최소한의 코드 수정으로도 복잡한 3차원 영역의 정확한 시뮬레이션을 가능하게 하였다.