相对边界参数的灵敏度

基于边界条件修改的目标优化需要该目标相对于边界参数的灵敏度,即 d L d Q b

考虑以下形式的成本函数 L

1. EQUATION_DISPLAY
L = L ( Q ( Q b ) , Q b )
(5100)

其中, Q 为以下非线性离散方程组的求解:

2. EQUATION_DISPLAY
R ( Q ( Q b ) , Q b ) = 0
(5101)

R 为流体方程的残差, Q b 为施加的边界条件。

由于此方程组的求解 Q 取决于施加的边界条件 Q b ,因此 L Q b 不仅具有显式相关性,还存在隐式相关性。 目标是量化此相关性。

可通过与网格相关性类似的方式实现该相关性的量化,即通过对残差方程求微分并将结果代入目标微分中:

3. EQUATION_DISPLAY
d L d Q b = L Q b ( L Q ) R Q 1 R Q b
(5102)

可通过使用以下公式简化此方程:

4. EQUATION_DISPLAY
R Q T ( L Q ) = Λ Q
(5103)

其中,可通过 Eqn. (5085) 获得流体伴随 ΛQ 的求解。

此简化会产生:

5. EQUATION_DISPLAY
d L d Q b = L Q b Λ Q T R Q b
(5104)

Eqn. (5104) 中可以看到,成本函数相对于边界条件的灵敏度计算包含两个步骤:

  • 成本函数对所施加边界条件的显式相关性的微分。
  • 考虑由 Eqn. (5104) 中第二项给出的隐式相关性,其需要获得全局伴随求解。