限制函数 — CAE术语解说

主要是限制数值计算中物理量的行为。例如，在CFD中防止密度或压力出现负值，并防止冲击波附近的解产生剧烈振荡。具体来说，流量限制函数（Flux Limiter）很出名，它确保高阶精度格式的单调性。

很好的问题。1阶迎风差分保持单调性但数值扩散很大，解会变得迟缓。另一方面，中心差分等高阶精度格式会产生振荡（超调）。限制函数通过混合这两种流量来发挥"混合系数"的作用，在光滑区域保持高阶精度，在变化剧烈的区域接近1阶迎风。混合的程度由限制函数 $\phi(r)$ 决定，其中 $r$ 是上游和下游的梯度比。

有几种，常用的是MINMOD、Superbee和Van Leer。MINMOD最为保守，稳定性很高但解有点迟缓。Superbee能保留锐利梯度但有过度平坦化的风险。Van Leer介于两者之间，被许多商用软件采用（如Ansys Fluent的2阶迎风格式默认使用），是个实用的选择。

在对流项的空间离散化阶段融入。MUSCL（Monotone Upstream-centered Schemes for Conservation Laws）格式是典型例子。在对单元界面处变量值 $u_{i+1/2}$ 进行重构时，通过从单元中心值 $u_i$ 及其梯度出发，乘以限制函数来修正梯度。写成公式就是：$u_{i+1/2}^L = u_i + \frac{1}{2}\phi(r_i)(u_i - u_{i-1})$ 这样的形式。

原理相同，但实现会更复杂。在非结构网格中，跨越单元界面的"上游"和"下游"单元不一定对齐。因此需要评估界面法向方向上的梯度，根据该方向沿着上游和下游单元中心值的差来定义 $r$。OpenFOAM的`limitedLinear`格式等采用这个方法。

影响很大。限制函数是非线性的，在隐式求解器中线性化会导致雅可比矩阵变得复杂。因此，完全隐式求解会降低收敛速度。在实务中，有时会"冻结"限制函数部分（使用前一次迭代的值）来线性化，或者采用Ansys Fluent的做法——"部分松弛格式"，在迭代中途更新高阶项重构。

应该从软件的默认设置开始。Ansys Fluent的"Second Order Upwind"默认使用Van Leer限制函数。如果出现振荡或负密度，试试更保守的MINMOD系（在Fluent中接近"First Order"行为的"Modified HRIC"等）。反过来，在航空航天领域需要锐利捕捉冲击波时，可以临时尝试Superbee，但发散风险高，通常不用于最终结果。

关系非常密切。粗网格中梯度评估会不准确，$r$ 值往往变成极端值。此时采用激进的限制函数（Superbee）容易在错误位置起作用，扭曲解。实务的铁则是"首先充分细化网格"，然后在此基础上使用默认限制函数。在网格无关性研究中，也应检查改变限制函数时的敏感性。

有两个重要检查点。第一，在物理量的云图或等高线中，看极端梯度部分（冲击波、火焰面、自由表面）是否出现不自然的"阶梯状"台阶。这是过度限制（过度扩散）的信号。第二，残差历史是否振荡。特别是密度或能量残差周期性波动时，说明限制函数的非线性作用很强。此时应切换到更光滑的限制函数或改进网格。

在"Solution Methods"窗格的"Spatial Discretization"部分。例如在"Momentum"的格式选择"Second Order Upwind"后，其下的"Limiter"下拉菜单可选择"Default"、"Minmod"、"Superbee"、"Van Leer"、"QUICK"等。但默认值是Van Leer系的专有实现。对于多相流的VOF法，在"Interface Modeling"下有"Sharpening/Diffusion"设置，这里选择界面专用的特殊限制格式如HRIC或CICSAM。

Star-CCM+在"Continua"→"Models"→"Implicit Unsteady"或"Segregated Flow"的设置内有"High-Order Scheme"选项可选限制。默认使用"Venkatakrishnan Limiter"，这是个高级函数，根据梯度大小平滑地降到1阶精度。OpenFOAM在`fvSchemes`文件的`gradSchemes`或`divSchemes`中直接指定。例如 `div(phi,U) Gauss limitedLinearV 1.0;` 最后的`1.0`是限制函数系数（1.0是正常限制，0是不限制）。

直接不太用这个术语，但类似概念存在。例如，LS-DYNA的显式格式中，元素应变或体积变化变得极端时的"体积粘性（bulk viscosity）"或"沙漏模式抑制"与限制函数起相同"防止发散"的作用。Abaqus/Explicit的"体积粘性缩放系数"也一样。这些从防止冲击波和大变形区域数值不稳定的角度看，在哲学上与CFD限制函数是一样的。

作为第一步对策很有效。特别是湍流模型（k-ε等）的变量应为正值，限制太弱时容易出现负值。出现这个错误时，在"Solution Methods"中将"Turbulence"的离散化格式暂时改为"First Order Upwind"（相当于强限制）。收敛后再改回"Second Order Upwind"，并将"Limiter"设为"Minmod"等更强的选项。根本上是细化产生错误的区域网格最稳妥。

很可能。细化网格后冲击波位置仍偏移数个单元以上时，使用的限制函数过于扩散（如MINMOD），把冲击波面弄钝了。反之，Superbee这样锐利的限制函数在粗网格中会过高估梯度，使冲击波出现太早。对策是充分细化网格后使用标准限制函数（Van Leer或Venkatakrishnan），在冲击波前后至少保持10～20个单元。可用超音速基准问题（如NASA楔形问题）来验证设置。

有可能。特别在分离流和尾流预测中，限制函数导致的过度数值扩散会减衰湍流结构，有效增加粘性，结果分离区变小，阻力被低估。验证方法是用同样网格采用不同限制函数（如默认 vs MINMOD）计算，比较后流涡结构和分离线位置。差异大时，根本原因通常是网格分辨率不足。在DES或LES等湍流模型中，限制函数影响会更明显，需要特别注意。