哦，

diag() = saveDiag;//恢复diag

是在括号之前的，所以每次的修改不会叠加。

之前理解错误的一点是所有的边界有关的系数都是放在internalCoeffs和boundaryCoeffs中的，不仅仅是coupled BC。从初始化时候的size可以看出来。

template<class Type> Foam::fvMatrix<Type>::fvMatrix ( const GeometricField<Type, fvPatchField, volMesh>& psi, const dimensionSet& ds ) : lduMatrix(psi.mesh()), psi_(psi), dimensions_(ds), source_(psi.size(), Zero), internalCoeffs_(psi.mesh().boundary().size()), //和边界的face数大小一样。 boundaryCoeffs_(psi.mesh().boundary().size()), faceFluxCorrectionPtr_(nullptr) { //...

不过非耦合的应该是internalCoeffs放比例系数，boundaryCoeffs放源项，而耦合的BC是internalCoeffs放owner的系数，boundaryCoeffs放的另一侧单元的系数。从addBoundarySource()的源代码可以看出。

template<class Type> void Foam::fvMatrix<Type>::addBoundarySource ( Field<Type>& source, const bool couples ) const { forAll(psi_.boundaryField(), patchi) { const fvPatchField<Type>& ptf = psi_.boundaryField()[patchi]; const Field<Type>& pbc = boundaryCoeffs_[patchi]; if (!ptf.coupled())//非耦合边界，只有owner，另一侧没有单元，pbc存的是边界源项。 { addToInternalField(lduAddr().patchAddr(patchi), pbc, source); } else if (couples)//耦合边界，如果couples==true, 另一侧有单元，pbc存的是另一侧单元的系数。 { const tmp<Field<Type>> tpnf = ptf.patchNeighbourField(); const Field<Type>& pnf = tpnf(); const labelUList& addr = lduAddr().patchAddr(patchi); forAll(addr, facei) { source[addr[facei]] += cmptMultiply(pbc[facei], pnf[facei]); } } } }

好吧，既然没有人回答，我就把我的笨办法贴出来
foamToVTK-->paraview-->打开算例文件和粒子Cloud文件-->File/Save Data-->得到每个时间步每个粒子信息.csv-->用python提取有用信息

18年我会搞这个...

在输运方程里加源相

在传输方程里面加源相？不就是加个字符么？还是怎么了

fvm::ddt(T) == S

并没有尝试过这个安装包。我推荐你采用bash自己装。http://dyfluid.com/install.html 在windows10中安装OpenFOAM

sampleDict 如何提取一个圆柱固壁上的所有速度啊

sigma是表面张力

需要给真实的压力，可以从单位判别出来，如果是压力/密度的，就不需要，如果是压力的单位，建议给真实的压力。

请问下这两者是不是同一种东西两个写法呢?

fvc::flux(U)可以选定各种格式，如果你指定他的格式是linear，那么fvc::flux(U)=linearInterpolate(U) & mesh.Sf()

@100yearsalone 我只试了一个算例，就这个算例来看，pimpleDyMFoam两个版本没有明显的速度差别

paraview对计算区域内部的一个面进行积分求平均流速
选了一个面slice
然后用plot Overline
最后用了integrateVariable积分，发现求得的平均流速差很多, 所以怀疑这个工具是否靠谱？有用过这个的吗？

//- Update moments manually //{ //const scalar& deltaT = mesh_.time().deltaTValue(); //scalarField& MIf = M_[kth]; //const scalarField& M0 = M_[kth].oldTime(); //MIf = 0.0; //fvc::surfaceIntegrate(MIf, mFlux_[kth]); //MIf = M0 - deltaT*MIf; //M_[kth].correctBoundaryConditions(); //}

但是这样没法驱动流动

如果开始就有速度的话，并且壁面没有摩擦力，可以一直流下去

可以通过在求解器中加一个源项在进口驱动流动

可以添加一个固定的压力梯度项

但是如果加入了紊流模型，紊动能在进口的的这个边界怎么给呢？仅仅是周期性边界好像不行吧？

给不了。周期性边界条件给不了随即的湍流波动。

@dzw05 好的，谢谢你的建议！我再试一试。

um. I dont understand. 我这是个2为网格，每个网格单元应该有4个内部面面，但是两个empty边界面上却具有通量？

你用的什么进口条件？你的设置符合常规设置么？see Ferziger and Peric page 315

@mohui 谢谢哈，

感觉有人手把手带着学OPenfoam是每个openfoam初学者的梦想。。所以最近国内有没有相关机构或个人组织培训？:big_mouth:

好的，谢谢李老师！

请问有谁测试过OpenFOAM的并行效率？在使用酷睿或至强的工作站或集群上并行最多可以达到单核性能的多少倍？我用I7-4790和I7-3920XM测试的结果，用阿姆达尔定律推测最高是2.5倍不到，而且线程数超过物理核心之后就是负优化了，感觉不太可能，无法解释那些上亿网格是怎么处理的。现在正在鼓动老板至少搞一台高性能工作站，同等价位下是主频优先还是线程优先？优先考虑酷睿还是至强？

0_1512907622643_2.jpg
fluent中设置是这样的：
Boundary A, B and C, D
are set as periodic boundary conditions. Boundary E is set as no slip
wall. Boundary F is set as opening pressure boundary.
不知道在OF怎样设置好

FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 100000; boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedFluxPressure; } atmosphere { type totalPressure; p0 uniform 100000; U U; phi phi; rho rho; psi none; gamma 1; value uniform 100000; } } dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedValue; value uniform (0 0 0); } atmosphere { type pressureInletOutletVelocity; value uniform (0 0 0); } }

体积项用setField，但这样不对，求大神帮忙:big_mouth: