Distributed optimal control problems governed by poroelasticity equations

기술 요약: 다공성 탄성 방정식에 의해 지배되는 분산 최적 제어 문제

문제 정의
본 논문은 비오(Biot)의 압축성 다공성 매질 모델에 의해 지배되는 다공성 탄성(poroelasticity) 내의 분산 최적 제어 문제의 수치적 근사를 다룬다. 이 물리적 시스템은 간극 유체 흐름과 다공성 매질의 변형 사이의 상호작용을 기술하며, 이는 지구물리학, 생물학적 조직 역학 및 탄소 격리 분야의 응용에 있어 매우 중요한 모델이다. 특정 제어 문제는 고체 변위($u$)와 유체 압력($p$)이 목표값($u_C, p_C$)으로부터 벗어나는 정도를 측정하고, 제어 변수의 크기에 페널티를 부여하는 비용 범함수(cost functional)를 최소화하는 것을 목표로 한다. 제어 변수는 영역 내에 분포된 유체 소스/싱크 항($m$)으로 식별된다. 본 연구는 가속도 효과가 무시되는 준정적 압축(quasi-static consolidation)에 초점을 맞추며, 특히 비저장 계수($s_0$)가 0에 가까워져 유체와 탄성 매질이 거의 비압축성이 되는 '다공성 탄성 잠김(poroelastic locking)' 현상과 관련된 수학적 및 수치적 과제를 다룬다.

방법론
저자들은 주요 미지수가 고체 변위와 유체 압력인 **새로운 2-장 대칭 정식화(two-field symmetric formulation)**를 제안하고 분석한다. 방법론은 다음과 같은 단계로 진행된다:

1. 연속 정식화 및 적정성(Well-Posedness):

   • 저자들은 운동량 균형의 시간 미분을 취함으로써 지배 방정식을 재정식화하여 시스템을 미분 방정식 집합으로 변환한다.
   • 적절한 함수 공간($V$는 변위, $Q$는 압력)을 정의하고, 세들 포인트 이론(saddle-point theory)과 레이디젠샤나-바부슈카-브레지(LBB) 조건을 사용하여 약해(weak solution)의 적정성을 확립한다.
   • 비용 범함수의 볼록성과 상태 연산자의 선형성을 활용하여, 직접 변분법(direct method in the calculus of variations)을 통해 최적 제어의 존재성과 유일성을 증명한다.
   • 1차 필요 및 충분 최적 조건은 아드조인트 상태(adjoint state) 접근법을 통해 유도되며, 이는 제어를 위한 변분 부등식으로 이어진다.

2. 이산화 스킴(Discrete Scheme):

   • 공간 이산화: 영역은 inf-sup 조건을 만족하는 유한 요소($V_h \times Q_h$, 예: MINI-element 또는 Taylor-Hood element)를 사용하여 이산화된다.
   • 시간 이산화: 후방 오일러 스킴(구체적으로 dG(0) 방법)이 사용되어 완전 이산 시스템을 생성한다.
   • 제어 이산화: 본 논문은 제어 변수에 대해 변분 이산화 접근법(Hinze가 제안한 방식)을 활용한다. 즉, 제어 공간을 명시적으로 이산화하지 않고, 이산 제어를 허용 가능한 제어 집합 위로 투영된 이산 아드조인트 상태로부터 얻는다. 이는 제어 공간을 이산화할 때 흔히 발생하는 근사 오차를 방지한다.

3. 오차 분석:

   • 저자들은 완전 이산 스킴에 대한 **사전 오차 추정치(a priori error estimates)**를 도출한다.
   • 분석은 전체 오차를 (고정된 제어에 대한) 이산화 오차와 연속 및 이산 최적 제어 간의 차이로 인한 섭동 오차로 분리한다.
   • 전방 상태(forward state) 및 후방 아드조인트 시스템 모두에 대한 안정성 추정치를 확립하여 이러한 오차를 제한한다.

주요 기여
본 논문은 문헌에 세 가지 주요 기여를 한다:

• 새로운 대칭 정식화: 비오 모델 내의 최적 제어 문제를 위해 특별히 설계된 2-장 대칭 정식화(변위 및 압력)를 도입하며, 이는 기존의 4-장 또는 혼합 정식화와 구별된다.
• 존재성 및 유일성: 새로운 정식화의 적정성을 엄격하게 증명하고, 최적 제어의 존재성과 유일성 및 1차 최적 조건을 도출한다.
• 사전 오차 추정치: 완전 이산 스킴에 대한 완전한 사전 오차 분석을 제공한다. 특히, $L^\infty(0, T; H^1(\Omega))$ 및 $L^2(0, T; L^2(\Omega))$ 노름에서의 변위와 압력, 그리고 $L^2(0, T; L^2(\Omega))$ 노름에서의 제어 변수에 대한 오차 추정치를 도출하며, 표준적인 정칙성 가정 하에서 최적의 수렴율을 달성한다.

결과

• 이론적 결과: 본 논문은 수치 스킴이 $O(h^s + \Delta t)$의 차수로 수렴함을 증명하며, 여기서 $h$는 공간 격자 크기, $\Delta t$는 시간 단계, $s$는 해의 정칙성과 관련된다. 구체적으로, 해가 충분히 정칙적이면($s=1$), 공간과 시간 모두에서 선형 수렴($O(h + \Delta t)$)한다.
• 수치적 결과: 단위 정사각형 영역에서의 수치 실험을 통해 이론적 발견을 검증한다. 오픈 소스 라이브러리인 FEniCS를 사용하여 저자들은 전방-후방 알고리즘(forward-backward algorithm)을 통해 최적성을 시스템을 해결하였다. 결과는 제안된 스킴이 최저 차수의 MINI 및 Taylor-Hood 요소를 사용할 때 변위, 압력 및 제어 변수에 대해 예측된 최적 수렴율을 달성함을 보여준다.

의의
저자들은 이 연구를 다공성 탄성 시스템에 의해 지배되는 최적 제어 문제의 수치 해석을 위한 기초적인 단계로 위치시킨다. 저자들은 (비저장 계수가 0인 가정을 바탕으로 한 Bociu 등의 연구와 같이) 변형 가능한 다공성 매체 내의 유체 흐름에 대한 최적 제어는 연구되어 왔으나, 일반적인 비오 모델에 대한 이러한 문제의 수치적 근사는 이전에 조사되지 않았음을 언급한다. 새로운 대칭 정식화의 적정성을 확립하고 완전 이산 스킴에 대한 엄격한 오차 추정치를 제공함으로써, 본 논문은 문헌의 공백을 메우며 생물학적 조직 모델링 및 지구물리학적 시뮬레이션과 같은 응용 분야에서 다공성 탄성 과정을 제어하기 위한 견고한 이론적 틀을 제공한다.