一个小bug调试了一个月

[李东岳]

原来的源代码：

volVectorField Su 
(
    IOobject
    (
        "Su",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
	mesh,	
	dimensionedVector
	(
		"Su",
		g.dimensions(),
		vector::zero	
	)
	//zeroGradientFvPatchField<scalar>::typeName
);

开始的时候把SU受力场边界的零法向梯度注释掉了，弄到NS方程里面总是有问题。把零法向梯度加回去就好了

这个bug搞了一个月啊

https://coding.net/u/dyfluid/p/GPBEFoam_PUB/git/commit/acf25fd806169e30f4294016b7826fe2a372b3d3

这个看起来更直观，就这么2行用了我一个月时间。

[浪迹天大]

一般人早就抓狂了吧。
NO_READ 的量，如果不指定，默认边界类型就是 calculated 吗？

[李东岳]

对，然后不计算边界的话就是0