OpenFOAM中partialSlip边界条件中valueFraction的设置问题

OYYO

最近在研究一个问题，涉及到壁面滑移条件partialSlip设置：
partialSlip中滑移速度的计算如下：
\begin{equation} u_w=\lambda {(\frac{\partial u}{\partial y})}_w \end{equation}
其中$\lambda$是滑移长度。
partialSlip中滑移速度实现的效果如下：
partial slip img
在OpenFOAM中实现partialSlip的效果用到了partialSlip的边界条件，具体设置如下：

<patchName>
{
    // Mandatory entries
    type            partialSlip;
    valueFraction   <Field>;

    // Optional entries
    refValue        <Field>;
    writeValue      <bool>;

    // Inherited entries
    ...
}

其中valueFraction 是设置partialSlip上的速度$u_w$相对于靠近壁面第一层网格中心$u_1$的大小，具体公式如下：
\begin{equation} u_w=(1-\alpha)u_1 \end{equation}
其中$\alpha$就是valueFraction 。
当$\alpha=1$，$u_w=0$就是无滑移边界。
当$\alpha=0$，$u_w=1$就是完全滑移边界。

实际问题中要是表面每个点的滑移长度$\lambda$不同，导致在模拟中一个表面网格上每个网格片上的$\alpha$不同，这该如何进行设置？

我了解到valueFraction 的type是field 类型，我不知道是否可以定义个类型也为field 类型的$\alpha$，然后赋值给它？
另外我的$\alpha$是与壁面网格的每个网格片剪切应力WallshearStress相关的，需要从流场中获取，我不清楚该怎么做？
请各位大佬帮着指点一下思路。

OYYO

笔误了几处：

其中valueFraction 是设置partialSlip上的速度$u_w$相对于靠近壁面第一层网格中心速度$u_1$的大小。
当$\alpha=0$, $u_w=u_1$就是完全滑移边界。

学流体的小明

在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题中说：
我了解到valueFraction 的type是field 类型，我不知道是否可以定义个类型也为field 类型的
𝛼
，然后赋值给它？

这个思路可以的。
wallShearStress，OpenFOAM有个后处理函数可以生成wallShearStress，它就是存在于patch上的。你可以遍历wall边界的面，先从wall边界上读取wallShearStress，然后有函数找这个面对应的内部网格，也就是你壁面上第一层的网格，这样给$\alpha$幅值就可以了。或者直接自己计算这样可以从wall边界上读取wallShearStress，算一个偏导数再乘以动力粘度就可以了。

OYYO

@学流体的小明在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题中说：

在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题中说：
我了解到valueFraction 的type是field 类型，我不知道是否可以定义个类型也为field 类型的
𝛼
，然后赋值给它？

这个思路可以的。
wallShearStress，OpenFOAM有个后处理函数可以生成wallShearStress，它就是存在于patch上的。你可以遍历wall边界的面，先从wall边界上读取wallShearStress，然后有函数找这个面对应的内部网格，也就是你壁面上第一层的网格，这样给$\alpha$幅值就可以了。或者直接自己计算这样可以从wall边界上读取wallShearStress，算一个偏导数再乘以动力粘度就可以了。

！感谢大佬，关注你好久了
OK，那我先试试。

CFD中文网

OpenFOAM中partialSlip边界条件中valueFraction的设置问题