[论文解读] Shape sensitivity analysis of time-dependent flows of incompressible non-Newtonian fluids

Jan Sokołowski, Jan Stebel|arXiv (Cornell University)|Sep 14, 2011

Rheology and Fluid Dynamics Studies参考文献 23被引用 4

一句话总结

该论文严格建立了在有界平面区域中，由幂律模型控制的时间依赖、不可压缩非牛顿流体流动的形状微分性质。通过速度法与物质导数框架，推导出阻力型代价泛函的形状梯度，由于在亚临界区域解具有强正则性和唯一性，因此可在无小数据限制的条件下实现基于梯度的形状优化。

ABSTRACT

We study the shape differentiability of a cost function for the flow of an incompressible viscous fluid of power-law type. The fluid is confined to a bounded planar domain surrounding an obstacle. For smooth perturbations of the shape of the obstacle we express the shape gradient of the cost function which can be subsequently used to improve the initial design.

研究动机与目标

• 严格分析二维时间依赖、不可压缩非牛顿流体流动的代价泛函的形状可微性。
• 在剪切速率依赖的幂律粘度模型背景下，推导阻力泛函的形状梯度。
• 建立所考虑的非线性、非定常的纳维-斯托克斯型系统中物质导数的存在性与正则性。
• 在无小数据假设的前提下，实现流固耦合问题的基于梯度的形状优化。
• 为工程与生物医学应用中涉及非牛顿流体的数值形状优化提供理论基础。

提出的方法

• 采用光滑向量场 T 定义域扰动的速度假设法，进行形式化的形状敏感性分析。
• 状态问题由速度 v、压力 p 及剪切速率相关粘度 ν(|Dv|²) 的非线性时间依赖方程组建模。
• 通过变量变换与解在域变形下的渐近展开，推导出物质导数。
• 通过对状态问题的弱形式关于域扰动求导，推导出形状导数的线性化系统。
• 通过将代价泛函的变化分解为两部分（一部分涉及解的导数，另一部分涉及几何项），计算形状梯度。
• 利用弱收敛与紧致性论证，使扰动方程中的极限成立，从而确保物质导数与形状梯度的存在性。

实验结果

研究问题

• RQ1 能否为幂律粘度模型下的时间依赖、不可压缩非牛顿流体流动，严格推导出形状梯度？
• RQ2 在不依赖小数据假设的前提下，形状敏感性分析是否依然有效，而依赖于解的强正则性？
• RQ3 在非线性、非定常设置下，科里奥利项的存在对形状导数计算有何影响？
• RQ4 物质导数在具有变粘度的非线性PDE形状敏感性框架中起什么作用？
• RQ5 形状梯度能否以适用于形状优化算法数值实现的形式表达？

主要发现

• 阻力泛函的形状梯度以闭式表达式推导得出，涉及速度的物质导数与应力张量。
• 在速度场梯度有界的假设下，确立了物质导数的存在性，该条件在亚临界区域成立。
• 形状导数关于向量场 T 连续，确保了形状优化过程中的稳定性。
• 该分析无需小数据或短时间区间假设，而是依赖于弱解的唯一性与正则性。
• 该方法可扩展至其他形状泛函与边界条件，尽管当前设置为便于理解而采用齐次狄利克雷边界条件。
• 该理论框架为未来在非牛顿流体动力学形状优化中的数值实现提供了支持。