[논문 리뷰] An Introduction to Adjoints and Output Error Estimation in Computational Fluid Dynamics

Steven M. Kast|arXiv (Cornell University)|2017. 12. 03.

Computational Fluid Dynamics and Aerodynamics참고 문헌 9인용 수 6

한 줄 요약

이 논문은 계산 유체역학(CFD)에서 이웃 방법과 출력 오차 추정에 대한 종합적인 소개를 제공하며, 출력 기반 메쉬 적응에의 응용에 초점을 맞춘다. 이는 이산 및 연속 이웃 형식을 유도하고, 출력 오차 추정 및 메쉬 정밀화 유도에의 응용을 보여주며, 수치 예제를 통해 출력이 초수렴할 경우 고차 오차 추정이 최적 수렴률을 달성할 수 있음을 보여준다.

ABSTRACT

In recent years, the use of adjoint vectors in Computational Fluid Dynamics (CFD) has seen a dramatic rise. Their utility in numerous applications, including design optimization, data assimilation, and mesh adaptation has sparked the interest of both researchers and practitioners alike. In many of these fields, the concept of an adjoint is explained differently, with various notations and motivations employed. Further complicating matters is the existence of two seemingly different types of adjoints -- "continuous" and "discrete" -- as well as the more formal definition of adjoint operators employed in linear algebra and functional analysis. These issues can make the fundamental concept of an adjoint difficult to pin down. In these notes, we hope to clarify some of the ideas surrounding adjoint vectors and to provide a useful reference for both continuous and discrete adjoints alike. In particular, we focus on the use of adjoints within the context of output-based mesh adaptation, where the goal is to achieve accuracy in a particular quantity (or "output") of interest by performing targeted adaptation of the computational mesh. While this is our application of interest, the ideas discussed here apply directly to design optimization, data assimilation, and many other fields where adjoints are employed.

연구 동기 및 목표

• CFD에서 이산 및 연속 이웃의 기본 개념을 명확히 하여, 서로 다른 표기법과 형식으로 인한 혼동을 해결한다.
• 이웃 해를 이용한 출력 기반 오차 추정을 위한 엄밀한 프레임워크를 수립하여 관심 있는 양의 정확한 예측을 가능하게 한다.
• 이웃 기반 오차 추정의 유용성을 바탕으로 CFD 시뮬레이션의 정확도 향상을 위한 메쉬 적응을 유도하는 데에 그 기능을 보여준다.
• 이산 및 연속 이웃 접근법을 비교하여 일致성 및 수렴 성질을 강조한다.
• 출력이 초수렴할 경우 고차 오차 추정이 최적 수렴률을 달성할 수 있음을 보여주며, 특히 고차 방법에서 그러한 성능이 두드러진다.

제안 방법

• 행렬 역행산을 통한 이산 이웃 유도를 통해 특정 출력을 계산하기 위해 역행렬의 한 행만 필요하다는 것을 보여준다.
• 변분법과 약한 형태를 통해 연속 이웃를 도입하여 연속 함수 공간 내의 PDE에 대한 이웃 해를 가능하게 한다.
• 이중 가중 잔여법(DWR)을 적용하여 출력 기능에 대한 가중치를 유도한 이웃 방정식을 풀어 출력 오차 추정을 도출한다.
• 갈레르킨 방법을 사용하여 원본 및 쌍대 형태의 오차 추정 간의 등가성을 보여주며 일관성을 확보한다.
• 이차 및 삼차 오차 추정을 사용하여 출력을 보정하고, 다양한 초수렴 상황에서의 수렴률을 분석한다.
• 3D 날개 및 비정상 상태의 날개류 시뮬레이션을 통해 접근법을 검증하며, 이웃 기반 오차 지표에 기반한 메쉬 적응을 보여준다.

실험 결과

연구 질문

• RQ1 이산 및 연속 이웃 형식은 유도 방식과 일치성 측면에서 어떻게 다를까? 각각의 형식은 CFD에서 어떤 장점을 지니는가?
• RQ2 어떤 조건에서 이웃를 통한 출력 오차 추정이 최적 수렴률을 달성할 수 있으며, 고차 추정이 필요한 경우는 언제인가?
• RQ3 이웃 기반 오차 추정은 관심 있는 양의 정확도 향상을 위해 메쉬 적응을 효과적으로 유도할 수 있는가?
• RQ4 출력 기능의 초수렴 현상이 이차 및 삼차 오차 추정의 성능에 어떤 영향을 미치는가?
• RQ5 고차 유한요소법에서 신뢰할 수 있는 오차 추정 및 메쉬 적응을 보장하기 위해 이웃 일치성은 어떤 역할을 하는가?

주요 결과

• 출력이 2p+1의 속도로 초수렴할 경우, 삼차 오차 추정은 최적 수렴률을 달성하지만, 이차 오차 추정은 부분 최적 수렴을 보인다.
• 표준 수렴 속도(p+1)의 경우, 이차 및 삼차 오차 추정 모두 미세한 공간의 출력보다 더 빠르게 수렴하므로, 보정된 출력은 추정 순서와 관계없이 동일한 수렴 속도를 갖는다.
• 이웃 기반 메쉬 적응은 3D 날개 및 비정상 상태의 날개류에서 목표 출력(예: 양력 및 항력)의 오차를 성공적으로 감소시키며, 더 적은 요소로도 오차를 크게 감소시킨다.
• 이산 이웃 연산자는 신뢰할 수 있는 오차 추정을 보장하기 위해 원본 이산 시스템과 일致해야 한다. 일치하지 않으면 부분 최적 수렴이 발생한다.
• 이중 가중 잔여법은 출력 오차 추정을 위한 강력한 프레임워크를 제공하며, 갈레르킨 방법에서는 원본 및 쌍대 형태 간의 등가성이 유지된다.
• 비정상 문제에서는 이웃 해가 시간에 따라 변하는 변화에 대한 민감도를 포착하여, 진동하는 날개류와 같은 일시적 유동에 대한 동적 메쉬 적응을 가능하게 한다.