Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
phi，phiHbyA是通量，通量的边界是固定的，不需要插值
phi，phiHbyA
@李东岳 李老师您好，我最近也对这个问题感到好奇，看了您的回答，也就是说phiHbyA在边界处直接赋值为边界处的通量phi即可，然后压力泊松方程里还有Ap这一项也要插值到面上，还想请教那么Ap在边界处是如何处理的呢？
Ap是体场，在边界不存在值，对边界不影响
@李东岳 李老师，那意思是说在如图公式构建压力泊松方程的时候，边界上的Ap并不做出贡献么，公式右侧的外部面求和不用考虑边界面么？我对该过程的理解有限，还希望李老师解答。
Ap在边界处是如何处理的呢？
http://dyfluid.com/piso.html 方程15，Ap本身是定义在vol上面的但是要插值到面上，就是简单的插值
边界上的Ap并不做出贡献么
有贡献，但不是phi那种直接定义出来的贡献，是插值过去的贡献
我冲说一下。Ap是定义在vol上的，是对角线系数，在边界处不存在值。但下一步需要定义Apf，存在边界值（包含通量的贡献）。
@李东岳 李老师，Ap定义我大概了解了，意思是计算域内部面的Apf是直接插值得到，而实际上Ap本身是对当前vol所有面如图公式求和得到的， ，那么在边界处的Apf也就考虑在边界面处的 和 就行了么。
那么在边界处的Apf也就考虑在边界面处的
Ap是UEqn.A()，Ap可以理解为一个scalarField，仅在内部点有值，但是为了一致性处理为volScalarField。Ap的boundary不影响计算。
UEqn.A()
Apf是fvc::interpolate(UEqn.A())，是Ap的插值，就这么个区别。
fvc::interpolate(UEqn.A())
@李东岳 @李东岳 李老师，我可能还是有点绕糊涂了，实际上我是在自己写代码尝试复现有限体积法的simple算法，那么在构建如图压力泊松方程的时候， 我的迷惑点就在于在边界处单元对该项 进行所有面求和时，边界面的Apf该如何处理，是直接不考虑其贡献不加它，还是把其取值为边界处的 + ，还是说是其它的处理方法？
Apf还有$\rd t$的贡献，就是http://dyfluid.com/piso.html 方程15
哦，如果你稳态，那就没有，你写的对
@李东岳 不好意思，李老师，工作原因回复慢了一些，我目前只在尝试稳态simple方法的实现，那就先只考虑边界面 + 的贡献么，我去修改代码再调试一下，非常感谢！
是的，你试试
，那么在边界处的Apf也就考虑在边界面处的
和
就行了么。
我的迷惑点就在于在边界处单元对该项
进行所有面求和时，边界面的Apf该如何处理，是直接不考虑其贡献不加它，还是把其取值为边界处的
+
，还是说是其它的处理方法？
+
的贡献么，我去修改代码再调试一下，非常感谢！