[论文解读] Example Setups of Navier-Stokes Equations with Control and Observation: Spatial Discretization and Representation via Linear-quadratic Matrix Coefficients

Maximilian Behr, Peter Benner+1 more|arXiv (Cornell University)|Jul 27, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 3被引用 12

一句话总结

本文提出了一套系统化的不可压缩Navier-Stokes方程空间半离散化框架，结合控制与观测，采用矩阵表示以实现高效的数值模拟。该方法引入基于张量的非线性项，并将雷诺数缩放分离，使结果可直接用于线性代数软件包。主要贡献在于提供了一个可复现的、开源的数据集，包含预组装的系统矩阵，用于两个基准流动问题——顶盖驱动涡旋与圆柱绕流，支持快速模拟与控制系统设计，无需有限元工具箱。

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.

研究动机与目标

• 为不可压缩Navier-Stokes方程提供一种空间半离散化、基于矩阵的表示形式，并引入输入与输出，使其可立即用于数值线性代数软件包。
• 通过解耦雷诺数依赖性并采用非线性项的张量表示，实现流体动力学问题的高效模拟与控制系统设计。
• 提供一个可复现的、开放获取的数据集，包含预组装的系统矩阵、边界条件与控制算子，用于两个典型流动基准问题：顶盖驱动涡旋与圆柱绕流。
• 通过将非线性项表示为矩阵张量H，使得H(v,v) = H·v⊗v，以支持模型降阶与系统理论分析。
• 通过惩罚型Robin方法实现分布式与边界控制，确保动量方程中线性输入算子的一致性。

提出的方法

• 采用速度与压力有限元空间的混合有限元方法进行空间半离散化，使用节点基函数。
• 通过标准Galerkin有限元公式组装系统矩阵A（扩散）、M（质量）、J（散度）与F（强迫项）。
• 通过展开的三线性形式将非线性对流项(v·∇)v表示为张量H，使得(H·v)⊗v可实现纯矩阵运算。
• 将扩散算子A中的1/Re缩放分离，通过缩放矩阵A即可实现不同雷诺数下的模拟。
• 采用惩罚型Robin方法处理非齐次Dirichlet边界条件与自然边界条件，实现通过线性算子一致建模输入。
• 通过独立的输入算子实现分布式与边界控制，边界控制通过Robin惩罚参数α实现。

实验结果

研究问题

• RQ1 如何将非线性Navier-Stokes方程以完全基于矩阵的形式表示，以适用于数值线性代数软件包？
• RQ2 在空间半离散化的Navier-Stokes系统中，如何最有效地解耦雷诺数依赖性，以实现高效模拟与模型降阶？
• RQ3 如何一致地将分布式与边界控制输入整合到半离散系统中，同时保持结构完整性并支持系统理论分析？
• RQ4 能否为基准流体流动问题建立统一的开源数据格式，以同时支持模拟与控制系统设计？
• RQ5 该非线性项的张量表示在实现由Navier-Stokes方程导出的双线性-二次系统高效模型降阶方面，其有效性如何？

主要发现

• 非线性项(v·∇)v成功表示为矩阵张量H，使得H(v,v) = H·v⊗v，实现高效纯矩阵运算，无需显式存储张量。
• 通过缩放扩散矩阵A，成功将雷诺数依赖性清晰分离，仅通过缩放A即可实现不同Re值下的模拟。
• 惩罚型Robin方法在边界条件中实现了Dirichlet条件与控制输入的一致建模，边界控制通过线性输入算子实现。
• 数据集包含两个基准问题（顶盖驱动涡旋与圆柱绕流）在多个网格分辨率下的系统矩阵（N=1,2,3用于圆柱绕流；N=1,2,30用于顶盖驱动涡旋），自由度范围为NV=5812至19512。
• 数值示例展示了在Re=40, 90, 800与1200下的稳态与瞬态模拟成功，使用Paraview兼容的.vtu文件准确可视化了速度与压力场。
• 提供的示例脚本（Python与MATLAB/Octave）可完整复现所有模拟，包括分布式与边界控制情形，输出包含选定时间点的速度与压力结果。