CFD中文网

@李东岳 懂了，谢谢李老师

各位老师好，迎风格式的示意图（红线所示）见的很多，但是中间缺的那部分我用绿线表示了（根据每个网格内变量线性变化），不知道我这样理解对不对，如果是这样，那么在界面左右两端取值发生突变，这个合理吗？

问题已解决，忘了这是1维算例了 宽度方向改成一个网格就行了

@李东岳 李老师好，表面张力是0.07，我看到的文献[1][2]都是数值结果，原始文献[3]下载不到，计算域高度是15m，宽度为1m，网格长度为0.025m。我已经尝试将宽度缩小为0.1m, 并考虑接触角，但是紊乱的问题仍然存在，压强也是衰减的。

[1]Sun H., Sun Z., Liang S., et al. Numerical study of air compressibility effects in breaking wave impacts using a CIP-based model. Ocean Engineering[J]. 2019,174: 159-168.
[2]Ma Z. H., Causon D. M., Qian L., et al. A compressible multiphase flow model for violent aerated wave impact problems. PROCEEDINGS OF THE ROYAL SOCIETY A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES[J]. 2014,470(2172).
[3]Dias F, Brosset L. 2010 Comparative numerical study: description of the calculation case. InProc. 20th Int. Offshore and Polar Engineering Conf., Beijing, PR China, 20–25 June.

各位老师好，我用compressibleInterFoam模拟了一维重力作用下液体活塞运动基准算例，文献中底部压强是震荡且保持稳定的幅值，但是我模拟出来的压强是逐渐衰减的


查看alpha场发现在底部气体压缩的过程中，气液界面紊乱了，推测是这个原因导致衰减。

算例二维，上下左右4个边界均是wall边界，无其他特殊设置。改变边界类型为滑移和无滑移没有大的变化，大佬们知道有什么方法可以解决界面变形这个问题吗？

感谢回答，尝试添加了fvc.H, 只报错以下内容

/home/dell/OpenFOAM/OpenFOAM-4.x/src/finiteVolume/lnInclude/cyclicAMIFvPatch.H:39:10: fatal error: cyclicAMILduInterface.H: 没有那个文件或目录
 #include "cyclicAMILduInterface.H"
          ^~~~~~~~~~~~~~~~~~~~~~~~~

不懂在哪里用到了cyclicAMILduInterface.H，先不管，它在src/meshTools路径下，在options里面添加路径，但是编译后又回到了一开始的报错

twoPhaseMixtureThermo.C:94:9: error: no matching function for call to ‘Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>::GeometricField(const char [4], Foam::tmp<Foam::Field<double> >)’
         )

求助大佬，为了在compressibleInterFoam求解器中添加有效粘度muEff

surfaceScalarField muEff
    (
        "muEff",
         mixture.muf()      
      + fvc::interpolate(rho*turbulence->nut())
    );

需要定义muf()函数，mixture是twoPhaseMixtureThermo类型，在对应的H文件和C文件进行更改

virtual tmp<surfaceScalarField> muf() const;     //muf()声明
Foam::tmp<Foam::surfaceScalarField> Foam::twoPhaseMixtureThermo::muf() const   //muf()定义
{
    const surfaceScalarField alpha1f     
    (
        min(max(fvc::interpolate(alpha1_), scalar(0)), scalar(1))
        
    );

    return tmp<surfaceScalarField>    
    (
        new surfaceScalarField   
        (
            "muf",
            alpha1f*(thermo1_->rho())*fvc::interpolate(nuModel1_->nu())          
          + (scalar(1) - alpha1f)*(thermo2_->rho())*fvc::interpolate(nuModel2_->nu())
        )
    );
}

编译时主要报这些错误 该如何解决？为什么类似的代码在incompressibleTwoPhaseMixture中可以编译？

/home/dell/OpenFOAM/OpenFOAM-4.x/src/OpenFOAM/lnInclude/GeometricField.H:99:40: error: invalid use of incomplete type ‘class Foam::surfaceMesh’
         typedef typename GeoMesh::Mesh Mesh;
/home/dell/OpenFOAM/OpenFOAM-4.x/src/OpenFOAM/lnInclude/GeometricField.H:103:48: error: invalid use of incomplete type ‘class Foam::surfaceMesh’
         typedef typename GeoMesh::BoundaryMesh BoundaryMesh;
twoPhaseMixtureThermo.C:83:17: error: ‘fvc’ has not been declared
         min(max(fvc::interpolate(alpha1_), scalar(0)), scalar(1))

举手手2718068983@qq.com

@ibelief 请问是视频教程么，可否分享一下，最近在自学这本书，有视频学起来应该效果更好