OpenFOAM中partialSlip边界条件中valueFraction的设置问题
-
最近在研究一个问题，涉及到壁面滑移条件partialSlip设置：
partialSlip中滑移速度的计算如下：
\begin{equation} u_w=\lambda {(\frac{\partial u}{\partial y})}_w \end{equation}
其中$\lambda$是滑移长度。
partialSlip中滑移速度实现的效果如下：
在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相关的，需要从流场中获取，我不清楚该怎么做？
请各位大佬帮着指点一下思路。
- 我了解到
-
在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：
我了解到valueFraction 的type是field 类型，我不知道是否可以定义个类型也为field 类型的
𝛼
，然后赋值给它？这个思路可以的。
wallShearStress，OpenFOAM有个后处理函数可以生成wallShearStress，它就是存在于patch上的。你可以遍历wall边界的面，先从wall边界上读取wallShearStress，然后有函数找这个面对应的内部网格，也就是你壁面上第一层的网格，这样给$\alpha$幅值就可以了。或者直接自己计算这样可以从wall边界上读取wallShearStress，算一个偏导数再乘以动力粘度就可以了。 -
@学流体的小明 在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：
在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：
我了解到valueFraction 的type是field 类型，我不知道是否可以定义个类型也为field 类型的
𝛼
，然后赋值给它？这个思路可以的。
wallShearStress，OpenFOAM有个后处理函数可以生成wallShearStress，它就是存在于patch上的。你可以遍历wall边界的面，先从wall边界上读取wallShearStress，然后有函数找这个面对应的内部网格，也就是你壁面上第一层的网格，这样给$\alpha$幅值就可以了。或者直接自己计算这样可以从wall边界上读取wallShearStress，算一个偏导数再乘以动力粘度就可以了。！感谢大佬，关注你好久了
OK，那我先试试。
