关于wedge边界的使用

[xuhuashi]

@韬智tz 1. 浮点数溢出导致发散了。发散的问题可能有很多，比如网格质量不好，比如某个变量计算变化太大（例如负压负温）。你可以试试减小松弛因子，或者先用低阶的离散格式试试，或者给个合适一点的初始参数。
2.气泡初始压力大于周围流体压力才会膨胀，你看看是不是气泡初始条件給错了，或者是不是气泡还没开始变化就已经发散了。

[韬智tz]

@xuhuashi 谢谢指点，我试试大佬的解决办法。

[韬智tz]

@xuhuashi 按您的解决办法我加密了网格（直接在blockMeshDict调x，z方向的网格数），没有解决。
松弛因子的话，我这个算例是OF自带的，fvSolution里面并没有设置relaxationFactor，是要在fvSolution里面加上松弛因子吗，我看tut里compressibleInterFoam那一类里面都设置为

relaxationFactors 
{ 
     equations 
     { 
             ".*"        1; 
     } 
} 

这个我不确定，所以问问您。
由于是自带算例，离散格式是这样的，我也没改，您看怎么改能解决发散的问题呢。

ddtSchemes
{
    default         Euler;
}

gradSchemes 
{
    default         Gauss linear;
}

divSchemes 
{
    div(phi,alpha)  Gauss interfaceCompression vanLeer 1;

    div(rhoPhi,U)  Gauss upwind;
    div(rhoPhi,T)  Gauss upwind;
    div(rhoPhi,K)  Gauss upwind;
    div(phi,p)      Gauss upwind;
    div(phi,k)      Gauss upwind;

    div(((thermo:rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes 
{
   default         Gauss linear uncorrected;
}

interpolationSchemes 
{
    default         linear; 
}

snGradSchemes
{ 
    default         uncorrected;
}

[韬智tz]

@xuhuashi 我加了这个，结果能算，不指定α和α=1不是一样的效果吗？为什么能算了呀

relaxationFactors 
{ 
      equations 
     { 
             ".*"        1; 
     } 
} 

加了这个能算，但是结果还是气泡没有变化，我的初始压力是大于周围流体的，之后又加大了两倍，也没用。

[xuhuashi]

@韬智tz 要不咱私聊吧，qq:1144231085

[李东岳]

@韬智tz $\alpha=1$收敛性更好。不指定alpha，可能会非对角占优。

气泡初始压力大于周围流体压力才会膨胀，你看看是不是气泡初始条件給错了，或者是不是气泡还没开始变化就已经发散了。

就像 @xuhuashi 大佬说的，看一下你的压力场，压力场也要setFields

[Prometheus10]

@yu_tian 你也来了 老哥 咱俩又见面了 之前这玩意查了好久 一直不知道咋整

[李东岳]

@prometheus10 你的ID让我想起了普罗米修斯

[Prometheus10]

@李东岳 东岳老师好~
是的 我以前添加您微信的时候 ID也叫这个

[韬智tz]

@xuhuashi 谢谢大佬指点，还想问下，楔形网格在指定区域加密用topoSet里面的source是用cylinderToCell，还是其他的呀。我这里用的cylinderToCell，但是出现下面的错误，请大佬看下！

[Prometheus10]

@韬智tz 你好！韬智tz
我用boxToCell 能加密成功
不妨一试：

actions
(
    {
        name        c2;
        type        cellSet;
        action      new;
        source  boxToCell;
        box     (0 -0.01 0) (0.002 0.01 0.0075);
    }
);

[韬智tz]

@prometheus10 谢谢，把区域选好了就直接在终端里输入topoSet，是吧，不用其他的操作了吧

[韬智tz]

@prometheus10 你好，用box加密成功了，也能算了，但是结果中出现这种啮齿状的：

这样太丑了，我觉得是网格问题，所以在同一个区域有加密了一次，出现了这样的网格：

这里请教大佬两个问题：1.这种啮齿状是不是网格的问题。2.加密两次出现这种网格，怎么让他更好看。

[Prometheus10]

@韬智tz
你好！韬智tz
在判断问题点之前，先咨询一下：

你的命令顺序是怎么样的？
是不是每次直接就接着做了refineMesh?
还有加密的方向对不对？
如果方便的话 给看一下代码和命令行

[韬智tz]

@prometheus10 就是先blockMesh，topoSet，两次refineMesh，然后setfields，之后就并行算了用的compressibleInterFoam。

topoSetDict如下

actions
(
    {
        name    c0;
        type    cellSet;
        action  new;
        source  boxToCell;
        sourceInfo
        {
          box (0 -0.015 0.008) (0.005 0.015 0.015); 
        }
    }
);

refineMeshDict如下

set             c0;

coordinateSystem global;

globalCoeffs
{
   e1              (1 0 0);
   e2              (0 1 0);
   e3              (0 0 1);
}

directions      (e1 e3);

useHexTopology  yes;

geometricCut    no;

writeMesh       no;

[Prometheus10]

@韬智tz
我手里有个现成的一个topsetDict，具体的各行作用暂时回答不上，要不你先试一下：

set             c0;

coordinateSystem global;

globalCoeffs
{
    tan1            ( 1 0 0 );
    tan2            ( 0 1 0 );
}

patchLocalCoeffs
{
    patch outside;  // Normal direction is facenormal of zero'th face of patch
    tan1 (1 0 0);
    tan1 (0 0 1);
}

directions      ( tan1 normal );

useHexTopology  yes;

geometricCut    no;

writeMesh       no;

[韬智tz]

@prometheus10 好，谢谢