CFD中文网

@浪迹天大 老师您好，按照这个方法，我需要更改controlDict中的maxDeltaT。但是会发现文件更新不起作用。

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version         2.0;
    format          ascii;
    class           dictionary;
    location        "system";
    object          controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     interFlowVb1;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         0.3;

deltaT          0.00005;

writeControl    adjustableRunTime;

writeInterval   0.05;

purgeWrite      0;

writeFormat     binary;

writePrecision  8;

writeCompression off;

timeFormat      general;

timePrecision   8;

runTimeModifiable yes;

adjustTimeStep  yes;

maxCo           0.5;
maxAlphaCo      0.5;
maxCapillaryNum 10; // stability criteria is 1

 maxDeltaT       2e-5;//controlDict.a;
// maxDeltaT       1e-5;//controlDict.b;
// maxDeltaT       5e-6;//controlDict.c;


functions
{
    fileUpdate1
    {
        type              timeActivatedFileUpdate;
        libs              ("libutilityFunctionObjects.so");
        writeControl      timeStep;
        writeInterval     1;
        fileToUpdate      "$FOAM_CASE/system/controlDict";
        timeVsFile
        (
            (-1 "$FOAM_CASE/system/controlDict.a")
            (1e-4 "$FOAM_CASE/system/controlDict.b")
            (2e-4 "$FOAM_CASE/system/controlDict.c")
        );
    }
}

三个文件只有maxDeltaT有差别，请老师捞捞@李东岳 @浪迹天大

运行后结果提示：

timeActivatedFileUpdate: copying file
"~/orifice3Da/system/controlDict.c"
to:
"~/orifice3Da/system/controlDict"

Courant Number mean: 2.0171174e-08 max: 0.3369988
Interface Courant Number mean: 2.978204e-10 max: 0.22052861
Capillary Number: 2.4526254
deltaT = 2e-05
Time = 0.00022

文件有更新，并且maxDeltaT也更新为了1e-5; 但是它没有发挥作用，可以看到文件更新后的deltaT 依然是 2e-05，如果是我自己运行过程中手动改小maxDeltaT，它会马上从deltaT = 2e-05，便为deltaT = 1e-05。

@好运来辣 我以前测试的时候，记得paraview不支持中文路径，注意检查这个影响因素。
包括paraview的安装路径（盘符有中文名称不影响），以及仿真文件的所在路径，并且文件名中的特殊符号最好也避免一下。

@好运来辣
直接去paraview的官方网页下载
paraview官网
最好下载绿色版的直接解压就能用，win服务器没试过，需要你自己试一试。可以的话回来跟帖和大家分享一下
具体的见图

没看懂你具体问的是啥问题，是不是问windows下的wsl2的具体地址？

你在wsl2命令行输这个命令试试：

explorer.exe .

@李东岳 谢谢东岳老师，我用的就是v2106，然后我看了下我按照v2106的211215的补丁，用wmake -build-info看了下补丁版本是0，自己尝试手动打补丁（替换文件），但是最终没起效果。

最后干脆用了简单粗暴的办法，把openfoam-OpenFOAM-v2106.211215这个版本下载下来之后直接编译。
依然感谢东岳老师。

各位老师，最近在GitHub发现一位作者开源了它的求解器，其中明确提到需要用OpenFOAM v2106
我尝试用东岳老师提供的OpenFOAM全系列虚拟机下载中的v2106编译，一直报错，通过不了。

后来查询到作者有这么一句提示：
If using a pre-compiled version of OpenFOAM v2106, use pre-compiled patch 211215.

我对应去看了这个v2106的211215的补丁

请问一下各位老师，怎么打这个补丁，甚至是如何查看自己的版本是否已经打过这个补丁了？

谢谢各位老师 这个打补丁 完全不知道怎么下手
感谢东岳老师的全系列全集 非常好用@李东岳

@chengan-wang 投影不清楚怎么弄，缩放行不行？用transform把x轴的缩放因子改成0.00001，你看看效果？

@lyc 你好，我很少遇到。即便遇到都是重启就能解决的那种。

@lwjetmann 老师您好，再请教一下您。

1. 用wallshearstress后处理算壁面粘性剪切力，然后在paraview中计算出在壁面各点法向和切向的投影；
2. 基于wallshearstress functionObject进行修改，输出投影后的结果

按照您提到的上述方案1，对于圆球绕流（球体界面是“noSlip”边界条件），我顺利的输出了粘性剪切力，结果也验证的上。

想请您再帮忙看看另外一种情况：绕流气泡，即将边界条件改为Symmetry：
因为wallshearstress只对wall类型的边界计算粘性剪应力，因此不能再通过【方案1】输出界面上的粘性剪应力结果。
关于设置Symmetry边界条件的原因是源于这里：Stokes流中的圆球绕流的阻力系数，怎么设置才能算准。

对此，针对绕流气泡的结果（边界条件为Symmetry）：

我尝试手动在结果中将Symmetry边界条件改为Wall，然后执行一次wallshearstress；
通过对比文献中气泡的压差阻力系数、法向粘性应力和切向粘性应力结果，输出的结果并不对

除此以外，是否只有像您提到的【方案二】基于wallshearstress functionObject进行修改，才能输出“绕流气泡”的粘性剪切力结果。
如果您方便的话，能否请您指点一些修改思路，谢谢老师。

@小泽同学 我此前还在想是否需要重新考虑从右下角的CD/DVD中加载一个适用于Ubuntu14.04版本的linux.iso，里面有vm-tools的安装包，相当于也是用里面写好的脚本重装VMware-tools。

emm... 总之，能用上就太好了，祝收敛！

@小泽同学 你好，不好意思今天出差回来，所以现在才回复。
目前可以确定关于剪切板的问题基本上就是出在这个open-vm-tools上，
因为open-vm-tools的各式各样的问题总是出现，也不好描述具体的处理方案。
因此，我建议按照以下代码重装一下open-vm-tools：
然后关闭虚拟机，接着重启再试试！

sudo apt-get autoremove open-vm-tools
sudo apt-get autoremove open-vm-tools-desktop
sudo apt-get update
sudo apt-get install open-vm-tools
sudo apt-get install open-vm-tools-desktop

@小泽同学 能否大致描述下你的执行过程，以及目前的情况？

@lwjetmann 感谢老师，组里同学大家都不太会OpenFOAM，我马上也看一下。再次感谢。

@李东岳 好的，谢谢东岳老师。
我也自己看了下[forces.C、forces.H，forceCoeffs.C、forceCoeffs.H]四个文件，

正如您得教材里所说的，在forceCoeffs.C中，总阻力的计算是：

$Cd=\frac{F_{drag}}{(\sum|\mathbf{S}_f|)\frac{1}{2}\rho|\mathbf{U}|^2}$

顺着这个我在forces.C中找了下各个力的计算公式：

其关系到的几个变量fN、fT、fP以及Md，在forces.H最先出现，

但是整个文档都看了下，好像都没有给出这几个量的定义，类似的还有pf、pm等：

东岳老师，各位大佬，能否指点一下这几个变量通常要去如何理解它的物理意义？
谢谢！

各位老师好！

最近计算圆球绕流，通过在controlDict中添加functions，能在每步输出周围流体对界面施加的总阻力系数。

由于总阻力系数是“压差阻力系数”和“粘性阻力系数”两个分量的总和，进一步开启log选项后，能分别显示出这两项，如下：

也就是说：
总阻力系数$C_D=C_{D,p}+C_{D,\tau}$

其中压差阻力系数$C_{D,p}=-4\int_0^\pi p\cos\theta\sin\theta d\theta $

粘性阻力系数$\begin{aligned}
C_{D,\tau}& =C_{D,\tau1}+C_{D_{\tau2}} \
&=4\int_0^\pi\tau_{r\theta}\sin\theta\sin\theta d\theta+4\int_0^\pi\tau_{rr}\cos\theta\sin\theta d\theta
\end{aligned}$

我想运行的时候输出切向粘性阻力系数$C_{D,\tau1}$ 和法向粘性阻力系数$C_{D,\tau2}$
或是从运行完之后的流场结果中计算出来。

forces.C和forces.H，forceCoeffs.C以及forceCoeffs.H文件还不太看得明白，希望能有老师指点，谢谢！

@Qing 感谢老师！您说的特别清楚，这对我特别有帮助。再次谢谢您！

@Qing 您好，能否请教一下您，请问在曲面上（比如一个半球面），滑移边界条件和对称边界条件有什么区别？

@tidedrinker 老师您好，还是不行。
相当于我的困境是：

用Symmetry，T文件夹内，无法给目标边界指定固定的温度值fixedValue 373。

然后想到官方文档也提到过Slip和Symmetry基本上就是数值等效的
但因为我的目标边界是个曲面，不是平面，
验证对比发现Slip边界条件 和 Symmetry边界条件应用到曲面边界得到的速度结果（以及基于此得出的阻力系数）并不一致

因此，另一个困境是：
用Slip边界，T文件夹内，虽然能给目标边界指定固定的温度值fixedValue 373，但是Slip边界条件 和 Symmetry边界条件算出来的流场信息不一样，现在拿不准哪个对。

@cresendo 老师您好，您说得确实是正确的，官方也提到过Slip和Symmetry基本上就是数值等效的

Notes on Computational Fluid Dynamics: General Principle

但是，能否请您帮忙看看这个帖子（您可以直接看第#54楼的评论回复）：
Slip边界和Symmetry边界应用于弯曲界面的结果对比

Slip边界条件 和Symmetry边界条件应用到曲面边界得到的速度结果（以及基于此得出的阻力系数）并不一致，这也是我目前较为疑惑的地方。

@pengdi @李东岳
东岳老师， PengDi老师，您们好

我目前仍然有一点没有想明白，想再请两位老师指点：

首先是Slip边界条件的定义：
【在界面处，法向速度为0，同时切向速度的梯度为0.】

沿用知乎博主的解析神秘色彩OpenFOAM 边界条件系列解析—Slip边界

COMSOL的帮助文档也是一致的内容
Theory for the Wall Boundary Condition

再有是东岳老师在#44楼所提到的：

因为slip是衍生来自symmetry，因此切向方向二者处理是一样的，都是零梯度。

那么，slip和symmetry按道理都应该同时符合这几个公式
$\mathbf{U}_{b\perp}=(\mathbf{U}_b\cdot\mathbf{n})\mathbf{n}=0$
$\mathbf{U}_b{|}=\mathbf{U}_c{|}$
$\mathbf{U}_c{|}=\mathbf{U}_c-\left(\mathbf{U}_c\cdot\mathbf{n}\right)\mathbf{n}$

但是，当Slip和Symmetry应用到弯曲边界的时候：为什么速度有这么大的差异，并由此得出的阻力系数也大不相同。
雷诺数Re=1时，阻力系数：
CD(slip)=13.49,
CD(symmetry)=17.6

从以上两个图结果上来看：
Symmetry的速度明显是和主流速度之间渐变过渡的（左图蓝线，大概从U/2逐渐发展到U），也就意味着流体在界面处是受到切应力的；
而Slip的速度是和主流速度一致（左图蓝线，Slip界面处就是和主流速度一样的U），因此我认为：

symmetry条件应用到弯曲边界的时候，并不会对弯曲边界的切应力做限制，因此弯曲边界的symmetry条件只是以下表达式：
$\mathbf{U}_{b\perp}=(\mathbf{U}_b\cdot\mathbf{n})\mathbf{n}=0$

此时Symmetry已经不再等同于Slip条件，其真实意义已经为：

Symmetry：【在界面处，法向速度为0，同时切向速度的梯度可以不为0.】
而Slip依然是：【在界面处，法向速度为0，同时切向速度的梯度为0.】

不知道这样理解是否正确？