@Rachel0096 可以看看interisoFoam作者Johan Roenby及合作者Henning Scheufler最新的论文及代码，里边有一个compressibleInterFlow求解器，应该就是你需要的。
论文：https://arxiv.org/abs/2103.00870
代码：https://github.com/DLR-RY/TwoPhaseFlow

大家好，我来更新这个问题的相关进展。最近两天我查了些资料，基本上解决了这个问题，感谢大家的帮助，尤其是bestucan大佬。
参考帖子关于相邻cell我写了个遍历进行暴力求解。
（1）代码思路：对所有网格做遍历，找到气液界面网格，针对每一个气液界面网格都对所有网格再做一次遍历，计算两网格之间的距离，对满足距离要求的网格在markRegionNearInterface这个场上+1。最后，可以利用markRegionNearInterface>=1来筛选出所有满足条件的网格。
（2）这段代码计算起来实在是费时间，因为我的需求(见贴名)不需要每一个时间步都运行这段代码，所以我在代码最外层套了一个时间步判断，可以实现多个时间步运行一次该代码。
（3）我从其他地方抄了点代码写了两个接口，分别从transportProperties读取 加密时间步间隔 和 加密范围，其中加密范围是依据气泡半径写的无量纲范围，所以需要配合前面代码(未在本段代码中给出)求出的气泡半径radius.value()定义标记区域。

代码如下

//read timestep interval required for the following code, as the following code is very time consuming so I needn't it run for every timestep.
const label calculateTimestepInterval = transportProperties.get<label>("calculateTimestepInterval");
//read information from transportProperties so that we can modify the size of marked region easily.
scalar markRegionNearInterfaceSize = transportProperties.get<scalar>("markRegionNearInterfaceSize");

  if (mesh.time().timeIndex() % calculateTimestepInterval == 0)//judge whether the timestep satisfy the timestep interval criterion.
     {
        markRegionNearInterface.field()=0;//reset the field to be zero first. Of course such field is claimed in createFields.H first.
        forAll(mesh.C(), i)//Cycle all meshes
           {
              if (alpha1[i]>0.0001&&alpha1[i]<0.9999)//judge whether the cell contains gas-liquid interface
                 {
                    vector centerOfSurfacePoint = mesh.C()[i];//store center information of such a interface cell
                    forAll(mesh.C(), cellI)//Cycle all meshes again
                    {
                           scalar offset = mag(centerOfSurfacePoint - mesh.C()[cellI]);//distance between any cells and the interface cell
                           if (offset <= markRegionNearInterfaceSize*radius.value())//judge the criterion, note the radius is calculated before with code which is not shown here
                              {
                                 markRegionNearInterface.field()[cellI] += 1;//any cells which meet the criterion will be marked with value>0
                              }
                    }
                 }
           }
     }

运行结果如下

BTW
（1） 因为我是个beginner，代码是东抄抄西抄抄弄出来的，上述代码可能有些冗余，如果大佬们有时间，想麻烦大佬们帮我看看我这些代码能从哪些角度优化优化。代码中的变量类型scalar vector label之类都是我一次次试错试出来的，也不知道这些变量类型用的合理不合理。
（2）我在代码中写了两个链接到transportProperties的接口读取 加密时间步间隔 和 加密范围，测试中发现程序只会在刚开始运行时读取这些常量。请问各位大佬，有没有什么办法能否让程序运行中也能读取这些常量呢？因为假设运行过程中我想改动这两个常量，目前只能把计算停下来改了后重新继续算，不能边算边改。

@香柏树 好的，多谢楼主分享

@香柏树 楼主您好，请问这个问题您找到解决方法了吗

@bestucan 太谢谢大佬了！我研究研究

@bestucan 对了大佬，我还碰到了一个小小的问题：

@ir77 在 标记气液界面附近一定距离的网格用于自适应网格加密 中说：

const labelList& neighbour = mesh.cellCells()[i];
scalar sumforinterface(0);
forAll(neighbour, cellI)
{
if (alpha1[neighbour[cellI]]>0.0001&&alpha1[neighbour[cellI]]<0.9999)
sumforinterface += 1;
else
sumforinterface += 0;
}
if (sumforinterface>0)
regionforinterface.field()[i]=1;
else
regionforinterface.field()[i]=0;

我把这段代码第一句中的mesh.cellCells()[i]分别改成了mesh.pointCells()[i]和mesh.edgeCells()[i]，运行结果却非常糟糕（如下图），似乎只有mesh.cellCells()[i]才可以正常使用。
大佬您能从经验上帮我判断一下可能是哪里出错了吗

@bestucan bestucan 大佬您好，我提到的

只要目标网格的任意一个相邻网格alpha位于标记范围内则该目标网格标记为1

我通过下面这段代码解决了

const labelList& neighbour = mesh.cellCells()[i];
          scalar sumforinterface(0);
          forAll(neighbour, cellI)
          {
           if (alpha1[neighbour[cellI]]>0.0001&&alpha1[neighbour[cellI]]<0.9999)
	      sumforinterface += 1;
           else
	      sumforinterface += 0;
          }
        if (sumforinterface>0)
	   regionforinterface.field()[i]=1;
        else
           regionforinterface.field()[i]=0;

然后我按照您建议的

因为你的 const labelList& neighbour2声明在 for 循环里，所以出了 for 循环无法使用
挪外面，和 neighbour 放一起声明，然后在 for 里赋值

将const labelList & neighbour2;放到了外边，但是报错error: ‘neighbour2’ declared as reference but not initialized
好像是说需要初始化，我不知道该怎么弄

后来我在cfdonline上找到了一条和我问题类似的帖子。参照18楼的代码完成了代码修改，实现了对气液界面附近的更多网格的标记。

代码放在下面，供大家参考，我是初学者，代码写的不规范，大家见谅

      const labelList& neighbour = mesh.cellCells()[i];
      labelHashSet setNBCells(1);
       for(int j=0; j<neighbour.size(); j++)
       {
      const labelList & neighbour2= mesh.cellCells()[neighbour[j]];
            setNBCells.insert(neighbour2);
       }
      labelList neighbour2=setNBCells.toc();
     scalar sumforinterface(0);
      forAll(neighbour2, cellI)
      {
       if (alpha1[neighbour2[cellI]]>0.0001&&alpha1[neighbour2[cellI]]<0.9999)
      sumforinterface += 1;
       else
      sumforinterface += 0;
      }

    if (sumforinterface>0)
   regionforinterface.field()[i]=1;
    else
       regionforinterface.field()[i]=0;

@bestucan 多谢大佬费心，我昨天也看到了这个帖子，昨天我参照这个帖子做出了一个初级版本，能够实现对气液界面附近三层网格的标记。
我采用的代码是这样的：

const labelList& neighbour = mesh.cellCells()[i];
          scalar sum(0);
          forAll(neighbour, cellI)
          {
            sum += alpha1[neighbour[cellI]];
          }
          sum /= neighbour.size();
        if (sum>0.0001&&sum<0.9999)
	   regionforinterface.field()[i]=1;
        else
           regionforinterface.field()[i]=0;

后来，我想给它添加一个嵌套进去，那样就可以加密气液界面附近五层网格了，采用如下代码：

const labelList& neighbour = mesh.cellCells()[i];
           for(int j=0; j<neighbour.size(); j++)
           {
          const labelList& neighbour2= mesh.cellCells()[neighbour[j]];
           }
          scalar sum(0);
          forAll(neighbour2, cellI)
          {
            sum += alpha1[neighbour2[cellI]];
          }
          sum /= neighbour2.size();
        if (sum>0.0001&&sum<0.9999)
	   regionforinterface.field()[i]=1;
        else
           regionforinterface.field()[i]=0;

但是编译的时候提示neighbour2未申明：
error: ‘neighbour2’ was not declared in this scope forAll(neighbour2, cellI)

bestucan大佬，能麻烦您帮我看下第二段代码吗，我不知道该怎么改
另外，bestucan大佬，我的第一段代码是参照东岳老师之前的一个帖子改的，是通过对目标网格所有相邻网格的alpha场做平均，如果平均值位于标记范围内则该目标网格标记为1。
我想，能否改一下代码，使得只要目标网格的任意一个相邻网格alpha位于标记范围内则该目标网格标记为1。感觉应该用if判断语句，但这里的判断对象又是多个网格中的alpha值，我不知道该怎么写这段代码，麻烦大佬指点一下迷津。
bestucan大佬费心了，谢谢大佬

bestucan大佬好，我能明白您的意思，脑子里也清楚整个流程，但是这一步*（给原始单层网格每个网格的左右+5的相邻网格都标记上）*我确实不知道该怎么用OpenFOAM语法表示出来

我能想到的只有我在第一楼最后提到的遍历网格，将0.001<α<0.999的网格赋值为1，但怎么将这个网格上下左右前后五个网格范围内的网格值也赋值为1，我不清楚该怎么弄

希望bestucan大佬能再花一点点时间帮我考虑一下这个问题

@李东岳
李老师好！
对于我的问题，我有两个思路，但限于水平太菜，两个思路都弄不出来：

1.标记气液界面附近的一定层数的网格，也就是前面说的“内部5层外部5层的网格”。我找到的比较相关的代码位于AMR源文件的dynamicRefineFvMesh.C里：

                // Extend with a buffer layer to prevent neighbouring points
                // being unrefined.
                for (label i = 0; i < nBufferLayers; ++i)
                {
                    extendMarkedCells(refineCell);
                }

这里的nBufferLayers是执行AMR时在dynamicMeshDict里填写的数字，可以通过这段代码控制缓冲层的层数。我想的是将气液界面标记成一个标量场（存在气液界面的网格赋值为1，其他网格赋值为0），并借鉴这段代码将标记区向内和向外“延拓”几个网格。但水平太菜，尝试了多次未能成功。

2.标记气液界面附近的的一定距离的网格。我想的是定义一个标量场，标量场的初始值都为0，然后对网格做遍历，只要网格到气液界面（0.001<α<0.999）的距离小于一定距离则该网格标量赋值为1。但是这种遍历很复杂我不知道该怎么写。

我的水平很菜，如果理解有偏差，表述有不清楚的地方望李老师指出 谢谢李老师。