@程迪

对于fvm算子生成的fvMatrix和fvc生成的volScalarField相加生成fvMatrix的时候，volScalarField会被自动乘以网格体积？

没错，参考 https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C#L1691

template<class Type> Foam::tmp<Foam::fvMatrix<Type>> Foam::operator+ ( const fvMatrix<Type>& A, const tmp<GeometricField<Type, fvPatchField, volMesh>>& tsu ) { checkMethod(A, tsu(), "+"); tmp<fvMatrix<Type>> tC(new fvMatrix<Type>(A)); tC.ref().source() -= tsu().mesh().V()*tsu().primitiveField(); tsu.clear(); return tC; }

感觉A是除以了体积的对角系数

没错，参考 https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C#L724

template<class Type> Foam::tmp<Foam::volScalarField> Foam::fvMatrix<Type>::A() const { tmp<volScalarField> tAphi ( new volScalarField ( IOobject ( "A("+psi_.name()+')', psi_.instance(), psi_.mesh(), IOobject::NO_READ, IOobject::NO_WRITE ), psi_.mesh(), dimensions_/psi_.dimensions()/dimVol, extrapolatedCalculatedFvPatchScalarField::typeName ) ); tAphi.ref().primitiveFieldRef() = D()/psi_.mesh().V(); tAphi.ref().correctBoundaryConditions(); return tAphi; }

好滴，谢谢啦~

@HydroJames 我暂时还查不到原因，测试过程中我还引入了其他差异，比如我进行过换硬盘的操作

forAll(walls.boundaryField(), patch) { walls.boundaryField()[patch] = L*hFlux[patch]/kappa; }

可以写如上代码计算，非常简单。上述代码没测试可能需要改动。

@李东岳

没，还是不知道那个公式咋来的，不过基本解决了flux limiter 和slope/gradient limiter的等价关系问题。

但是这玩意儿也是个大坑，对于多维网格带skew和non-uniform的情况，限制器设计和计算量都挺大。

OpenFOAM具体实现还没看太懂。有一个先有鸡还是先有蛋的问题，到底是interpolation先调用gradient呢？还是gradient调用interpolation？

我看的文章里是先leastSquare算gradient，不用Gauss的情况下是不用插值的。然后interpolation可以用一些limited gradient。

暂时只看懂了cellLimited, faceLimited没看懂。

@aufalpen 谢谢，公众号我早已经关注过了😄我去看看那个网页，多理解理解理论，还请多指教

边界条件不直接修改patchInternalField，而是通过改变离散后的系数矩阵起作用。

@李东岳 之前发过pdf，不过后来好像删了，你问问他。

单水滴破碎的倒是见过，多水滴碰撞的还没见过，楼主找到没有？

没用过gmsh...

网格体积：mesh.V()
网格面矢量 mesh.Sf()
网格面大小 mesh.magSf()

#问题1 自我解决，感觉像是bug.
0_1488141087886_streamline.png

@cfd-china 谢谢回复，我研究了源码已经搞懂了，顺便附上我总结的链接laplacian项离散解码

我这样写不对啊？ 不知道哪位可以给看看毛病在哪@administrators @李东岳

相对于计算发散，通常存在一个致命的因素可以解决。这种关于计算不准的东西，很难断定为什么，我经常计算不准。只能看sci，调节参数，调节网格，调节各种各种，毕竟人家算收敛了，收敛表明这个设置的结果，就是这样的。

Henry最近添加了部分描述：点击查看

确实，默认的rhoCentralDyMFoam没有使用correctPhi, 而且在其他可压流求解器中有使用correctPhi，不过需要注意的是，由于我在这里涉及到的是激波模拟，因此有压力间断，但是correctPhi似乎是通过压力进行通量修正，我不敢确定这是否会对压力间断产生影响，因此目前暂时还没有使用。我尝试过网格不运动，即直接使用rhoCentralFoam情况似乎是正确的。

只要原来的网格点都在就没问题

@100yearsalone
手头没有openfoam，但是在你使用dynamicMothion...代替solidBody....之后要定义solver，在你的字典里添加solver就可以了

solver xxx;

也有可能是：

solver { }

gedit打开文件很慢，可以换个编辑器，hoho