CFD中文网

@xuqiming 噢噢，这个问题解决了，需要厚度刚好覆盖这一层网格。

@yzwang 直接toposet就行，然后进行积分

constrainPressure主要是更新压力第二类边界条件，公式如下：

$$\begin{array}{}\text{(1)}& {\left(\mathrm{\nabla }{p}_{rgh}\right)}_f\cdot {\mathbf{n}}_f=\frac{({\mathbf{H}\mathbf{b}\mathbf{y}\mathbf{A}}_f^{\ast }-\frac{1}{A_{\mathrm{P},f}^n}(\mathbf{g}\cdot \mathbf{h}\mathrm{\nabla }\rho {)}_f-{\mathbf{U}}_f)\cdot {\mathbf{S}}_f}{|{\mathbf{S}}_f|\frac{1}{A_{\mathrm{P},f}^n}}\end{array}$$
上述公式与代码并不一致。在OpenFOAM中，constrainPressure为

forAll(pBf, patchi) { if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi])) { refCast<fixedFluxPressureFvPatchScalarField> ( pBf[patchi] ).updateCoeffs ( ( phiHbyABf[patchi] - rho.boundaryField()[patchi] *MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi]) ); } }

多乘了一个密度。应该改为：

forAll(pBf, patchi) { if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi])) { refCast<fixedFluxPressureFvPatchScalarField> ( pBf[patchi] ).updateCoeffs ( ( phiHbyABf[patchi] - MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi]/rho.boundaryField()[patchi]) ); } }

在非常老的OpenFOAM版本中，看起来是正确的，与公式一致

setSnGrad<fixedFluxPressureFvPatchScalarField> ( p_rgh.boundaryField(), ( phiHbyA.boundaryField() - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField()) )/(mesh.magSf().boundaryField()*rAUf.boundaryField()) );

可以试试看symmetry边界

诱导区长度很大，反应物的活性很差，很难起爆，起爆能量高了又会受初始起爆区的影响很久胞格数目和初始的高压区数目一样，没法过渡到自身特征的胞格尺寸

这个在算法上可以理解为这样么：

诱导区长度很大（计算域太长），反应物的活性很差（化学反应设置的不好），很难起爆（开始计算的非常小时间内，反应放热不够，能量上不来）

@李东岳 老师，我试过了，还是发散的，但是换更高质量的网格也可以收敛，（不收敛的网格icem中的quality为0.6，收敛的为0.75）这种情况是不是应该去尝试找一些新的收敛性更好的算法来实现呢