CFD中文网

请教各位老师，在计算域顶部设置一层空气，将温度设为368K，水域温度也设置为368K，计算域底部针管注入373K的蒸汽，计算后发现计算域上方气液交界T面存在一根线，明明是相同的温度，为什么会出现这样的现象呢？

@bestucan 万分感谢您抽时间回答了这个简单的问题，我之前被报错提示误导了，后面找到了这个小问题！最后再次感谢您！

@abc 您好，我之前使用过这种方法，但是由于makeAxialMesh这个工具和openfoam的版本有很强的联系，所以我在后面转化为轴对称的时候失败了，显示无法将以前的网格去掉，请问您遇到过相似的情况吗？

@yu_tian 请问您解决这个加密不均匀的问题了吗？

感觉应该是格式问题，但是又看不出。

Generating mesh for region fluid
Creating block mesh from
    "/home/dyfluid/OpenFOAM/dyfluid-v1812/run/scrivenWedge/D-4mm-70-nophase-end/system/fluid/blockMeshDict"
--> FOAM Warning :
    Reading "/home/dyfluid/OpenFOAM/dyfluid-v1812/run/scrivenWedge/D-4mm-70-nophase-end/system/fluid/blockMeshDict" at line 63
    Too many closing ')' ... was a ';' forgotten?
--> FOAM Warning :
    Reading "/home/dyfluid/OpenFOAM/dyfluid-v1812/run/scrivenWedge/D-4mm-70-nophase-end/system/fluid/blockMeshDict" at line 120
    Too many closing ')' ... was a ';' forgotten?
--> FOAM Warning :
    Reading "/home/dyfluid/OpenFOAM/dyfluid-v1812/run/scrivenWedge/D-4mm-70-nophase-end/system/fluid/blockMeshDict" at line 123
    Imbalanced brackets


--> FOAM FATAL IO ERROR: 
"ill defined primitiveEntry starting at keyword 'blocks' on line 48 and ending at line 123"

file: /home/dyfluid/OpenFOAM/dyfluid-v1812/run/scrivenWedge/D-4mm-70-nophase-end/system/fluid/blockMeshDict at line 123.

    From function void Foam::primitiveEntry::readEntry(const Foam::dictionary&, Foam::Istream&)
    in file db/dictionary/primitiveEntry/primitiveEntryIO.C at line 238.

FOAM exiting
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.001;
//D0=3mm,W=1.5D0,H=10D0//
vertices
(
    (0 0 0)
    (2 0.1 0)
    (0 0 20)
    (2 0.1 20)
    (6 0.3 20)
    (20 1 20)
    (0 0 40)
    (2 0.1 40)
    (6 0.3 40)
    (20 1 40)
    (0 0 60)
    (2 0.1 60)
    (6 0.3 60)
    (20 1 60)
    (2 -0.1 0)
    (2 -0.1 20)
    (6 -0.3 20)
    (20 -1 20)
    (2 -0.1 40)
    (6 -0.3 40)
    (20 -1 40)
    (2 -0.1 60)
    (6 -0.3 60)
    (20 -1 60)    
);

blocks
(
    hex  (0 14 1 0 2 15 3 2) (10 1 50)
    simpleGrading  (1 1 1)
    hex  (2 15 3 2 6 18 7 6) (10 1 50)
    simpleGrading  (1 1 1)
    hex  (6 18 7 6 10 21 11 10) (10 1 50)
    simpleGrading  (1 1 1)
    hex  (15 16 4 3 18 19 8 7) (20 1 50)
    simpleGrading  (1 1 1)
    hex  (18 19 8 7 21 22 12 11)) (20 1 50)
    simpleGrading  (1 1 1)
    hex  (16 17 5 4 19 20 9 8) (70 1 50)
    simpleGrading  (1 1 1)
    hex  (19 20 9 8 22 23 13 12) (70 1 50)
    simpleGrading  (1 1 1)
);

patches
(
    wall faceWall
    (
      (5 17 20 9)
      (9 20 23 13)
    )

    wedge front
    (
      (0 1 3 2)
      (2 3 7 6)
      (6 7 11 10)
      (3 4 8 7)
      (7 8 12 11)
      (4 5 9 8)
      (8 9 13 12)
    )

    wedge back
    (
      (0 14 15 2)
      (2 15 18 6)
      (6 18 21 10)
      (15 16 19 18)
      (18 19 22 21)
      (16 17 20 19)
      (19 20 23 22)
    )

    slip sideWalls
    (
      (1 14 15 3)
      (15 16 4 3)
      (16 17 5 4)
    )

   patch outlet
    (
      (10 21 11 10)
      (21 22 12 11)
      (22 23 13 12)
    )

   patch inlet
    (
      (0 14 1 0)
    )

    empty axis
    (
      (0 2 2 0)
      (2 6 6 2)
      (6 10 10 6)
    )
);

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

你好，我也遇到相同的问题了，请问您解决这个问题了吗？

@Exthan 你好，请问你运行makeAxialMesh出现warinning的问题是怎么解决的呢？

@浪迹天大 感谢大佬回答，我是这样设置的，但是不行，具体原因我还在找。

请问各位老师我使用makeAxialMesh后便提示我这个Warning,然后我谷歌后将writePrecision提高到12也不行，请问是怎么回事呢？我是在openfoam2.4.0中使用的。

@李东岳 好的，谢谢东岳老师！

请问interthermalphasechangefoam求解器是只能在2.4.0版本中编译吗？

@东岳 瞬态问题的话，是不是就是内残差作为收敛标准，初始残差不看的？

@东岳 使用pyFoamPlotRunner.py 这个命令，出来的只有初始残差。

@东岳 不好意思，老师！请问在openfoam中如何判断收敛呢？pyFoam可以查看其他的残差吗？还是在log文件中查看呢？确实是新手不太熟悉，网上也是说的不是很明白。

@freeman 是瞬态问题，那外残差是不是应该为0啊？

使用pyFoam查看残差，这个显示的是外残差吗？这个数值多少是收敛的呢？在fvSolution中设置的残差是10-6.请教各位老师。

@东岳
老师，请问您说的加密中线是指的气液交界这个面，也就是孔口的网络密度吗？入口很长是考虑针管的长度的影响。

@东岳 老师！您好，请问您知道问题出在哪儿了吗？

@东岳 这个模型是二维轴对称，在openfoam中转化为楔形的，只能整个文件夹发给您了，超过了这儿的上传大小了。

@东岳 已经通过邮箱发给您了。