CFD中文网

@学流体的小明 我想请问一下，这里的体积源项假定了一个平均流速 这个平均流速是怎么计算出来的，可以和重力挂钩嘛

@李东岳 谢谢李老师回复，可是全场加密下1000w网格跑的比局部加密67w网格快是不是不太合理

网格是需要接近1：1：1加密吗？整个区域尺寸的比例约为1：300：100

各位老师好，我在计算三维多相流流动时采用界面局部加密的动网格，发现时间步长特别小，界面局部加密后的网格67w时间步长1e-7 1e-8左右；
而在全场加密1000w网格情况下计算是正常的时间步长1e-5；
请问是局部加密非结构化网格的原因导致步长很小？ 四面是周期性边界条件，因此采用了cyclicAMI的边界。

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh   dynamicRefineFvMesh;

// How often to refine
refineInterval  1;

// Field to be refinement on
field           alpha.water;

// Refine field in between lower..upper
lowerRefineLevel 0.001;
upperRefineLevel 0.999;

// If value < unrefineLevel unrefine
unrefineLevel   10;

// Have slower than 2:1 refinement
nBufferLayers   1;

// Refine cells only up to maxRefinement levels
maxRefinement   2;

// Stop refinement if maxCells reached
maxCells        20000000;

// Flux field and corresponding velocity field. Fluxes on changed
// faces get recalculated by interpolating the velocity. Use 'none'
// on surfaceScalarFields that do not need to be reinterpolated.
correctFluxes
(
    (phi none)
    (nHatf none)
    (rhoPhi none)
    (alphaPhi0.water none)
    (ghf none)
);

// Write the refinement level as a volScalarField
dumpLevel       true;


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2212                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;


hN  1.01e-4;

L   #calc "100*$hN";
D   #calc "1.0*$hN";
H   #calc "3.0*$hN";
z   #calc "M_PI*$L";

nX  100; 
nY  50;
nZ  100;

vertices
(
    (0  0  0) //0 
    ($L 0  0) //1
    ($L $H 0) //2
    (0  $H 0) //3

    (0  0  $z) //4
    ($L 0  $z) //5
    ($L $H $z) //6
    (0  $H $z) //7
);

blocks
(
    hex (0 1 2 3 4 5 6 7) ($nX $nY $nZ) simpleGrading (1 1 1)
);

edges
(
);


boundary
(
    movingWall
    {
        type wall;
        faces
        (
            (0 1 5 4)
        );
    }
    inlet
    {
        type                  cyclicAMI;
        neighbourPatch        outlet;
        transform       translational;
        separationVector ($L 0 0);
        faces
        (
            (0 4 7 3)
        );
    }
    outlet
    {
        type                  cyclicAMI;
        neighbourPatch        inlet;
        transform       translational;
        separationVector (-$L 0 0);
        faces
        (
            (1 5 6 2)
        );
    }
      atmosphere
    {
        type patch;
        faces
        (
            (3 2 6 7)
        );
    }
    front
    {
        type cyclicAMI;
        neighbourPatch        back;
        transform       translational;
        separationVector (0 0 $z);
        faces
        (
            (0 1 2 3)
        );
    }
    back
    {
        type cyclicAMI;
        neighbourPatch        front;
        transform       translational;
        separationVector (0 0 -$z);
        faces
        (
            (4 5 6 7)
        );
    }
);

mergePatchPairs
(
);

// ************************************************************************* //

各位老师好，在求解NS方程的时候，理论上压力只是一个过程量，不需要设置边界条件，只需要给定一个参考点压强就可以。openfoam能否实现这个操作呢？

比如一个自由界面上的边界，只设置计算网格上一点压力为0

blockMesh里面的edge可以解决

各位老师好，请问openfoam需要用弯曲的壁面，比如正弦函数，用什么工具或者方法

@学流体的小明 是的，因为p_rgh不是周期的。你的方法我也想尝试一下

@李东岳 看到李东岳老师的推送了 没有这么邪乎哈哈哈

@学流体的小明 不是，我直接在interFoam求解器里修改，动量方程里把rgh单独拿出来计算了。

@User0221 你具体做的什么东西我不太清楚，不好评价，但是个人觉得和网格细化没有多大的关系。

@User0221 入口速度流量小 就应该变薄吧

问题已解决， 是interFoam和周期性边界条件的问题，参考：
https://www.cfd-china.com/topic/6577/interfoam-周期边界条件的压力问题/4

更新，已解决。 将重力在interFoam中当成体积力拿出单独计算，结果如图所示

各位老师好，最近在研究在有重力作用下的interFoam周期性边界条件，如气泡上升，液膜降落等，
看到几个论坛交流的帖子，有fixedJump设置压力梯度的，但是不适合所有的情况 https://www.cfd-china.com/topic/1992/cyclic周期性边界条件

在这篇论坛里提到
https://www.cfd-online.com/Forums/openfoam-solving/79300-interfoam-cyclic-boundary-condition-issue.html
interFoam求解器求解的时候把重力放到了压力计算的过程中，$p_{\rho g h}=p-\rho gh$。 请问有老师做过相似的问题可以交流一下吗

你把速度放大试试

@李东岳 太感谢李老师了！https://core.ac.uk/reader/36634570
interFoamHF在这篇论文里有提及。74页提及了一些参数设定。

@李东岳 我的程序就是最简单的interFoam，求解器文献中有对比，没有液膜厚度。但是我算的结果液膜没有流动，实在找不到原因了。文献中优化后的interFoamHF有厚度方程。
https://core.ac.uk/reader/36634570

@李东岳

各位老师好，我最近在做一个液膜受重力在平板上流动的模拟，示意图如下

初始薄膜厚度$h_0=200\mu m$的速度为半抛物线速度分布$u_{water}=\frac{g}{\nu}(h_0y-y^2/2)$,$u_{Air}=\frac{g}{\nu}\frac{h_0^2}{2}$,初始速度示意图如下
初始相场：

初始条件设置如下


运行结果如下
film.avi
结果好像液膜没有按照重力方向沿着x流动，左右两边向中间收了。重力设置确定没有问题，有老师可以帮忙看看哪里出了问题吗。非常感谢！
参考文献：Rohlfs, W., Scheid, B., 2015. Phase diagram for the onset of circulating waves and flow reversal in inclined falling films. J. Fluid Mech. 763, 322–351. https://doi.org/10.1017/jfm.2014.684
我的程序代码
ultracyclic4.zip 再次感谢！