CFD中文网

@李东岳
是的，我确实还需要多加学习！谢谢李老师

@李东岳
内部我就想模拟为一个装满水的容器，弧面我视作固体壁面，然后通过映射把力映射到结构求解器中，在结构求解器中给予材料属性来计算结构变形，然后把变形映射回OpenFOAM。底部边界我也看作壁面，就如同水槽中底部边界那样设置。我现在把边界条件设置成这样alpha.water p_rgh pointDisplacement U, 这样的边界条件无法计算，报错是这样的log.interFoam 。

@李东岳
李老师，如果我采用分区域耦合的话，整个模拟分成Fluid-outer, Solid, Fluid-inner三部分。那我的Fluid-inner的边界条件该如何给定？Fluid-inner有两个边界弧(interface)和底部边界。弧(interface)是需要用来耦合交换数据的。

@李东岳
我用interfoam计算，算例会很快发散。查看结果发现结构物附近的p_rgh值有异常的变化, 而且结构物上的液面有流速，位置正是液面锯齿处。但是我用wave2foam来计算的话，计算不会发散，可以计算完20s.我把两个算例放在了github(https://github.com/NUGUY2022/OpenFOAM_revise.git). 请李老师给点意见

@李东岳
之前计算的单一流体域的是计算的弹性板

您的意思是让我只计算流体域的上半部分看看吗？把圆弧结构内部填实吗？

@李东岳
流体域只有一个的时候是可以很好的计算的，我之前一直计算的单一流体域。

@李东岳
算变形
我想知道我的流体求解器设置是不是有问题。

@李东岳
李老师，我这个网格有大大的问题啊。结构物上方的初始液面呈现锯齿状, 然后我单独计算Openfoam可以计算，但是流固耦合的时候会发散。这应该是网格导致的问题把, 李老师可以给我一些建议来解决这个问题吗？
我的算例放在了github(https://github.com/NUGUY2022/OpenFOAM_Calculix_membrane.git)

@李东岳
感谢李老师的回复。
两个区域是分开的，我用其它的结构求解器来计算结构的变形，然后传递到这两个区域，但流体的两个区域不会互相影响。
我想做二维的模拟，于是我把前后面设置为empty，这时checkMesh出现下面的问题。
log.checkMesh
该如何解决这个问题？

Fluent3DMeshToFoam将fluent网格转换成openfoam网格时出现这样一条内容“ <<Writing 158278 conflicting points to set nonManifoldPoints”，想请问一下大家nonManifoldPoints会对计算产生影响吗？整个网格checkMesh显示的是Mesh OK。网格如图所示

另外还想问一下，我中间结构将流体域分成了内外两部分，网格这样划分是否合适？
polyMesh和fluent.msh也都上传了，希望大家可以提供一些建议。fluent.rar polyMesh.rar

@李东岳 看到这些质量不好的网格点，好想像hypermesh一样去手移

请教一下各位老师，利用snappyHexMesh贴合STL模型时，有什么办法可以让我STL模型的边界网格点（在不改变细化等级和边界网格总数的前提下）等距分布吗？如图所示。

@李东岳 正在学习，谢谢李老师！

@李东岳 明白了，谢谢李老师！

@李东岳 李老师，下面那一行代码alpha()*rho()*betaStar_*omegaInf_*kInf_和alpha()*rho()*betaStar_*omega_()*k有什么区别吗？
我看上面代码中有对omegaInf_和kInf_定义，

    kInf_
    (
        dimensioned<scalar>::lookupOrAddToDict
        (
            "kInf",
            this->coeffDict_,
            k_.dimensions(),
            0
        )
    ),
    omegaInf_
    (
        dimensioned<scalar>::lookupOrAddToDict
        (
            "omegaInf",
            this->coeffDict_,
            omega_.dimensions(),
            0
        )
    )

那是否可以理解成omegaInf_和kInf_就是0？不明白alpha()*rho()*betaStar_*omegaInf_*kInf_这串代码的意义，请教一下老师，谢谢！

@李东岳 公式没问题，只是左边红框内代码fvm::Sp(alpha()*rho()*epsilonByk(F1, tgradU()), k_)不知道是什么意思，是公式里面的哪一项？刚开始学of编程，求李老师指点一下，十分感谢！

@李东岳 李老师您好，我目前用的版本是v1812

编程小白请教各位老师两个不成熟的问题。。。

图一中湍流模型方程左右都乘了一个相分数，那如果在两相流模型中，都是空气时，相分数为0，并且相分数的时间导数和空间导数也是0，这方程怎么解？

图二中左边的代码和右边的公式感觉不匹配，求老师们指点一下，十分感谢！

最近在做风浪相互作用的模拟时遇到一个问题：在气液剪切过程中，湍流的研究是非常重要的，我在加入湍流模型之后，计算结果耗散特别严重，但是没有气体作用（nowind）时，SSTk-w湍流模型能够较准确的模拟，从下面的水位图中也可以看出来。并且湍动能k和湍流粘度nut在出口outlet逐渐向计算域内部扩散导致波高衰减（见下图），但是不知道具体原因。单纯的平底地形结果也没问题，但是遇到到这种斜坡地形就有衰减。
目前我已经尝试过of自带的几乎所有RANS类的湍流模型，计算结果均耗散严重（以SSTk-w为例），也尝试过对边界条件和壁面边界层的各种处理方式，还是效果不理想，想请教一下各位大佬，想解决这个问题，我还应该从哪些地方入手？
二维模型如图：

水位图：

不同时刻湍动能k分布：



不同时刻湍流粘度nut分布：

@李东岳 是的李老师，我理解成了这个边界用在入口时相当于固定值，由inletValue指定，用在出口相当于零梯度，应该是我理解错了。。。