udf离壁面距离函数

Question

问题遇到的现象和发生背景

遇到的现象和发生背景，请写出第一个错误信息

```c
BOUNDARY_FACE_GEOMETRY(f,t,A,ds,es,A_by_es,dr0);
tauwall1=C_MU_EFF(c,tcell)(du/ds);
wallfricv1=sqrt(tauwall1/C_R(c,tcell));
if (yplus<11.25)
{
tauwall2=C_MU_EFF(c,tcell)C_STRAIN_RATE_MAG(c,tcell);
wallfricv2=sqrt(tauwall2/C_R(c,tcell));
}
else
{
wallfricv2 = (1./0.41)log(yplus9.);
}
ds=C_WALL_DIST(c,tcell);
wallfricv4=(VISCyplus)/(dsC_R(c,tcell));

```用代码块功能插入代码，请勿粘贴截图。 不用代码块回答率下降 50%

这里面的ds好像是以前版本的不太正确，现在想用网格中心到壁面距离乘2来表示，但我在帮助文档里没找到求距离的函数，请问有人知道吗？求指导

Answer 1

在 OpenFOAM 中，有一个函数叫做 boundaryDistance，可以用来求出网格中心到壁面的距离。你可以使用以下语句来获取网格中心到壁面的距离：

scalarField ds = boundaryDistance(mesh, c);

其中 mesh 是你的网格对象，c 是你的网格单元对象。这个函数返回的是网格中心到壁面的距离的标量域。如果你想要在你的计算中使用距离的平方，你可以使用以下语句：

scalarField ds2 = sqr(ds);

如果你想要使用网格中心到壁面距离的两倍，你可以使用以下语句：

scalarField ds2 = 2*ds;