[논문 리뷰] Example Setups of Navier-Stokes Equations with Control and Observation: Spatial Discretization and Representation via Linear-quadratic Matrix Coefficients

Maximilian Behr, Peter Benner외 1명|arXiv (Cornell University)|2017. 07. 27.

Advanced Numerical Methods in Computational Mathematics참고 문헌 3인용 수 12

한 줄 요약

이 논문은 제어 및 관측이 가능한 비압축성 라우아-스토크스 방정식의 공간 이산화를 위한 체계적인 프레임워크를 제시한다. 매트릭스 기반 표현을 사용하여 효율적인 수치 시뮬레이션을 가능하게 하며, 텐서 기반 비선형 항과 레이놀즈 수 스케일링을 분리함으로써 선형 대수 패키지에서 직접 사용할 수 있도록 한다. 주요 기여는 라이드드 캐비티와 실린더 웨이크와 같은 기준 유동에 대해 사전 조립된 시스템 행렬을 포함하는 재현 가능한 오픈소스 데이터 세트를 제공하는 것이다. 이는 유한요소 도구상자 없이도 빠른 시뮬레이션과 제어 시스템 설계를 가능하게 한다.

ABSTRACT

We provide spatial discretizations of nonlinear incompressible Navier-Stokes equations with inputs and outputs in the form of matrices ready to use in any numerical linear algebra package. We discuss the assembling of the system operators and the realization of boundary conditions and inputs and outputs. We describe the two benchmark problems - the driven cavity and the cylinder wake - and provide the corresponding data. The use of the data is illustrated by numerous example setups. The test cases are provided as plain PYTHON or OCTAVE/MATLAB script files for immediate replication.

연구 동기 및 목표

• 수치 선형대수 패키지에서 즉각 사용 가능한 입력과 출력이 있는 비압축성 라우아-스토크스 방정식의 공간 이산화된 매트릭스 기반 표현을 제공하기 위해.
• 비선형성의 텐서 표현과 레이놀즈 수 의존성의 분리함으로써 유체역학 문제의 효율적 시뮬레이션과 제어 시스템 설계를 가능하게 하기 위해.
• 두 가지 표준 기준 유동인 라이드드 캐비티와 실린더 웨이크를 위한 사전 조립된 시스템 행렬, 경계 조건, 제어 연산자와 함께 재현 가능한 오픈액세스 데이터 세트를 제공하기 위해.
• 비선형 항을 매트릭스 텐서 H로 표현함으로써 모델 차수 축소와 시스템 이론적 분석을 지원하기 위해. 여기서 H(v,v) = H·v⊗v이다.
• 분포형 및 경계 제어의 구현을 촉진하기 위해 페널티가 부여된 루빈 접근법을 통해 운동량 방정식에서 일관된 선형 입력 연산자를 가능하게 하기 위해.

제안 방법

• 속도 및 압력 유한요소 공간에 노드 기반 기저 함수를 사용한 혼합 유한요소 방법을 통한 공간 이산화.
• 표준 갈레르킨 유한요소 공식을 사용하여 시스템 행렬 A(확산), M(질량), J(발산), F(외력)의 조립.
• 비선형 대류 항 (v·∇)v를 전개된 삼선형 형식으로 변환하여, (H·v)⊗v를 가능하게 하는 텐서 H로 표현함으로써 행렬 연산만으로 계산 가능하게 함.
• 확산 연산자 A에서 1/Re 스케일링을 분리함으로써, A의 스케일링만으로 다양한 레이놀즈 수 범위에서의 시뮬레이션 가능.
• 페널티가 부여된 루빈 접근법을 사용하여 비동차 딜레르흐 및 자연 경계 조건을 통합함으로써, 선형 연산자를 통한 일관된 입력 모델링 가능.
• 분포형 및 경계 제어를 별도의 입력 연산자로 구현하며, 경계 제어는 루빈 페널티 매개변수 α를 통해 실현함.

실험 결과

연구 질문

• RQ1 비선형 라우아-스토크스 방정식은 어떤 방식으로 수치 선형대수 패키지에 적합한 순수 매트릭스 기반 형태로 표현될 수 있는가?
• RQ2 공간 이산화된 라우아-스토크스 시스템에서 레이놀즈 수 의존성을 어떻게 효과적으로 분리하여 효율적 시뮬레이션과 모델 축소를 가능하게 할 수 있는가?
• RQ3 분포형 및 경계 제어 입력은 반산형 시스템에 어떻게 일관되게 통합될 수 있으며, 이를 통해 구조 유지와 시스템 이론적 분석이 가능하게 할 수 있는가?
• RQ4 기준 유체 유동 문제에 대해 시뮬레이션과 제어 시스템 설계를 모두 지원하는 통합된 오픈소스 데이터 포맷을 어떻게 확립할 수 있는가?
• RQ5 비선형 항의 텐서 표현 방식은 비선형-이차형 시스템에서 유도된 라우아-스토크스 방정식에 대해 효율적 모델 차수 축소를 얼마나 잘 가능하게 하는가?

주요 결과

• 비선형 항 (v·∇)v는 H(v,v) = H·v⊗v를 만족하는 매트릭스 텐서 H로 성공적으로 표현되어, 명시적 텐서 저장 없이도 효율적인 행렬 연산만으로 계산 가능.
• 확산 행렬 A의 스케일링을 통한 레이놀즈 수 의존성의 청결한 분리로, A의 스케일링만으로 다양한 Re 값 범위에서의 시뮬레이션 가능.
• 경계 조건에 대한 페널티가 부여된 루빈 접근법은 딜레르흐 조건과 제어 입력을 일관되게 모델링할 수 있게 하며, 경계 제어는 선형 입력 연산자로 실현됨.
• 데이터 세트에는 라이드드 캐비티와 실린더 웨이크에 대해 다수의 메쉬 해상도(N=1,2,3 for cylinder wake; N=1,2,30 for driven cavity)에서의 시스템 행렬이 포함되어 있으며, 자유도 수는 NV=5812에서 19512까지 다양함.
• 수치 예제는 Re=40, 90, 800, 1200에서의 안정 상태 및 일시적 시뮬레이션을 성공적으로 구현하였으며, Paraview 호환 .vtu 파일을 통해 속도 및 압력장의 정확한 시각화를 가능하게 함.
• 제공된 예제 스크립트(Python 및 MATLAB/Octave)를 통해 분포형 및 경계 제어 사례를 포함한 모든 시뮬레이션을 완전히 재현 가능하며, 특정 시간점의 속도 및 압력 출력 제공.