CFD中文网

李东岳

@evensun 可以有两种方法在动量方程中添加项。一种是在速度方程中，一种在压力方程中。如果加在了速度方程中，那就不需要处理压力方程。反之亦然

李东岳

@evensun 排版好评。

我想在simpleFoam的动量方程中增加一个湍流应力项

已经包含湍流应力了啊。如果你只是想在栋梁方程上添加一个$\nabla\cdot\tau$，最简单的就是直接写在栋梁方程上就好了。不需要考虑压力方程。

李东岳

ugks, dugks和LBM都是解得动力学方程，但是算法还是有极大区别的

李东岳

@cresendo 那就是纯对流方程了

李东岳

Ubuntu 下的安装

手头没有 Ubuntu 的机子，所以这里在虚拟机中演示；考虑到用 Ubuntu 的用户中很多在终端环境（命令行）下工作，这里就尽可能在终端中完成各个步骤。

需要管理员权限的地方会明确的用 sudo 指出来。

挂载安装文件

打开终端。首先，创建挂载的位置：

sudo mkdir -p /tmp/texlive

切换到下载文件所在的位置，比如你将 texlive.iso 下载到了 ~/Downloads 下，那么就首先切换到这里：

cd ~/Downloads

然后将安装包挂载到前面所创建的位置上去：

sudo mount -o loop texlive.iso /tmp/texlive

然后切换到 /tmp/texlive 目录下，就可以看到安装包里的所有内容了：

cd /tmp/texlive
ls  # 可以看到安装目录下的内容，其中 install-tl 是安装脚本

如果你不想只用终端，则在图形界面（Ubuntu 下的 GNOME 桌面）中，打开文件管理器，对着安装包点右键，也会有挂载的选项；此时一般会把盘挂载到 /media/xxx/TeXLive2021 位置（xxx 是你的用户名），在文件管理器的左侧栏中可以看到，进入即可。

运行安装脚本

接下来就可以运行安装脚本 install-tl 了。这是一个 Perl 脚本，需要有 Perl 解释器才能运行，不过 Ubuntu 上默认装了有，不担心。假如你的 Ubuntu 上没有 Perl 环境，请先安装一下：

sudo apt install perl

然后就可以执行脚本了。需要考虑一下你的需求：

如果你想给当前用户安装，安装在你的家目录（~，或者说 /home/xxx）下，那么不需要 sudo；
如果你想给所有用户安装，安装在系统目录（比如 /usr/local、/opt/ 等等）下，那么就需要 sudo。

对新手来说，不用考虑太多，越安全越好。所以就不要 sudo 吧！直接安装在你的家目录里：

./install-tl

如果发现找不到 install-tl，请检查一下所在位置是否正确！

配置安装选项

和 Windows 不同，Linux 下执行安装脚本时的配置界面限制在终端内，没有图形界面，如下图所示。

注意一下上面用蓝字标出来的，由于前面启动安装脚本的时候没有 sudo，所以无法安装在它的默认安装位置 /usr/local/texlive/2021 下。因此，首先要做的事情是输入 <kbd>D</kbd> ，回车，修改安装目录，切换到下面这个界面：

根据提示，按 1 再回车修改 TEXDIR 变量的值，这里我们把 /usr/local 整体替换为家目录 ~ 就好了：

输入完成之后，回车，可以看到配置中的各项路径都跟着修改了：

再输入 R，回车，返回最开始的界面；此时所有配置都没问题了，所以直接输入 I，回车，开始安装！

跟 WIndows 上的情况差不多。可以注意到这里的模块数是 4164，比 Windows 上的还多几个。和 Windows 一样，耐心等待安装结束，过程中不要关掉终端（或在终端中杀掉当前进程）。等到出现「欢迎进入 TeX Live 的世界！」（或者其英语形式）的时候，说明安装成功，此时该进程应当就自动终止了。

如果安装过程因意外中断，请首先清理安装目录：

rm -rf ~/texlive

然后再重新安装。

配置环境

安装结束时，会出现欢迎语；与 Windows 系统不同，请注意最后一段：

这一段是说，请手工配置以下三个环境变量：

将 ~/texlive/2021/texmf-dist/doc/man 加入到 MANPATH 中；
将 ~/texlive/2021/texmf-dist/doc/info 加入到 INFOPATH 中；
将 ~/texlive/2021/bin/x86_64-linux 加入到 PATH 中（如果你安装在不是 x86_64 架构的机子上，最后一部分可能不太一样）

前两个无所谓，大家现在都习惯于上网查资料；但若不配置最后一条，就无法从命令行正确启动 LaTeX 软件。所以，请通过以下命令将 TeX Live 的程序目录加入到环境变量中去：

echo "export PATH=$PATH:~/texlive/2021/bin/x86_64-linux" >> ~/.bashrc

如果你用的是 zsh 等其他的 shell，就加到它们自己的配置文件里去。关掉当前的终端，再重启；或者在当前的终端内 source ~/.bashrc。然后，检查一下 PATH 的值：

可以看到，TeX Live 的安装目录已经出现在了 PATH 中。尝试运行一下 xelatex 程序：

一切正常！现在就可以在命令行中使用各种 LaTeX 程序了。让我们用两句话生成一个 PDF 文档：

echo '\documentclass{article}\\begin{document}Hello, world!\end{document}' > test.tex  # 生成一个最简单的 LaTeX 文档
pdflatex test.tex  # 编译这个文档

看看结果：

这样就表明，我们的 TeX Live 已经在 Ubuntu 上正确安装了！安装过程到此结束。

需要说明的是，TeX Live 在 Linux 上安装时不会附带安装 Windows 上有的 TeXWorks IDE。如果你想要配置一个用起来顺手的 LaTeX 编辑器，就需要自己再花些心思。这个属于安装以外的范畴了，不多讲。

附录 1：安装过程常见问题

欢迎各位读者补充！

安装路径中可以有中文吗？

答：最好不要有，在 Windows 系统下尤其如此。能避免就避免！安装过程中可能一切正常，但装好之后就可能出现各种问题，不得不通过重装解决。
以前安装的版本可以直接升级吗？

答：通常不行，一般的建议是：要么保留旧版本，安装新版本，然后把环境变量 PATH 指向最新的版本；要么先把旧版本删干净，然后再重新安装。如开头所言，TeX 软件历史悠久，所以很多方面的发展没有跟上时代潮流；有人喜欢这样，也有人不喜欢——尤其是新手们。
怎么配置编辑器 / IDE？推荐用哪些？

答：这个可以参考 install-latex-guide-zh-cn 中的第 5 节：使用编辑器，那里的介绍非常全面。

附录 2：其他安装方式简介

除了以上所介绍的通过安装包安装 TeX Live 的方式之外，还有以下一些方式：

安装 MikTeX：直接到 MikTeX 官网 https://miktex.org/ 下载安装包，然后走一遍安装流程即可，这个过程要比 TeX Live 的安装容易很多、快很多；但是，正如前面所说的，MikTeX 的特点是：刚开始只安装很少的一些模块（总共几百 MB 左右），在以后的使用过程中自动安装其所发现的缺失模块——这个过程通常很麻烦，比如你的某一文档中有几十个模块没有安装，则 MikTeX 可能要反复运行好几次才能把所有的模块装齐，有时候可能还找不到相应的模块（比如你调用了某个字体模块，用的并不是 \usepackage 命令，MikTeX 就识别不了）。所以，除非你经验丰富，并且希望用运行时安装模块的时间来换一次性装齐所有模块的硬盘空间，否则不要选择这种方式。
安装 ProTeXt：ProTeXt 是一套整合起来的工具，包括 MikTeX（提供 LaTeX 程序）、TeXStudio（IDE）、Sumatra PDF（阅读器）。只能在 Windows 上使用。安装过程基本上是图形界面，操作起来比较容易；但是使用的人似乎不多。下载地址：http://mirrors.ctan.org/systems/win32/protext/protext.zip
通过 Linux 系统的包管理器：包括 Ubuntu、Fedora 在内的主流 Linux 发行版中都有 TeXLive 的软件包，比如 Ubuntu 上可以直接通过这个命令安装 TeX Live 完整版：
```
sudo apt install texlive-full
```
Fedora 上则是：
```
sudo dnf install texlive
```
这种方法的问题在于：发行版提供的 TeX Live 软件包一般比较旧，比如 Ubuntu 现在提供给 20.04 LTS 的版本是 TeX Live 2019，即使是以新著称的 Fedora 也只提供到 TeX Live 2020。如果你使用 Linux 系统，嫌以上的安装过程繁琐，或者想统一管理软件，可以采用这种方式。

附录 3：怎么查看 TeX Live 自带的文档

在环境变量已经正确配置的情况下（见前面的安装说明），在 Windows 系统下打开命令提示符 / 在 Linux 系统下打开终端，执行如下命令：

texdoc xxx

其中 xxx 是你要查找的模块名称，可以是：

一个程序的名字，比如 texdoc xelatex 就会打开 xelatex 程序的说明书；
一个文档类的名字，比如 texdoc thuthesis 就会打开清华大学学位论文 LaTeX 模板 thuthesis 的文档；
一个宏包的名字，比如 texdoc geometry 就会打开用来调页面边距的 geometry 宏包的参考文档（这样你就不需要再在网上搜它的参数怎么设置了！）；
一个字体的名字，比如 texdoc fandol 就会打开 Fandol 中文开源字体的说明文档。

当然可能还不限于上面几类。最常用应该是第三种，因为大家用 LaTeX 时经常为各种宏包的选项犯难，这时候通过 texdoc 命令查文档就比去网上漫无目的地搜索「经验」要快得多、好得多。

当然，绝大多数文档都是英文的，这对于工程技术人员来说应该不算一个大困难。如果你一定想读中文文档，可以到这个地方找现有的一些文档的中文翻译版本。

李东岳

@stan同学 reconstructParMesh

李东岳

有几个人联系我了。东西大体上就是网上的东西凑到一起。丝毫看不到用心所在。让我怀疑这种方式领取邀请码是不是浪费我的审核时间。看来为了利益做工作，就是应付。只有真心想写出教程方便大家，才能写出好东西。

如果给我发邮件，没回复。就表示你写的东西我没感到热度。

李东岳

@joann 那个太难debug了跑不动..

李东岳

方程44 http://dyfluid.com/interFoam.html

但实际上比这个要复杂。因为这个方程使用MULES进行求解。

另外，这个方程里面，没法关闭对流项，不可能看扩散情况。

李东岳

实在没空弄了我看看有没有人给你弄 https://www.cfd-china.com/topic/4861/领取任务免费赠送cfd中文网注册资格

李东岳

我帮你处理下

李东岳

@wying 要非常高的精度，比如普方法、有限差分那种7 8 9阶，有限体积法大体都不行

李东岳

@zhy2022 我琢磨下帮你处理下这个问题

李东岳

很多人都说不满足注册条件。现在发布一些任务，如果你可以完成，那么就送你注册资格。

Latex安装方法（win版本）

任务要求：在CFD中文网发布window下的Latex安装方法。并通过测试人员的安装。

Latex安装方法（ubuntu2004版本）

任务要求：在CFD中文网发布ubuntu2004下的Latex安装方法。并通过测试人员的安装。

ANSYS在ubuntu2004下的安装

任务要求：在CFD中文网发布ubuntu2004下的ANSYS安装方法。并通过测试人员的安装。

OpenFOAM-8-ICC编译方法

任务要求：在CFD中文网发布使用ICC对OpenFOAM-8进行编译安装的方法。并通过测试人员的安装。

测试人员

任务要求：依据他们写的安装流程，实测能否顺利安装。中间可能你们要彼此沟通，修修改改。

上面的任务，邮件联系，我会给你CFD中文网测试账户，使用测试账户撰写相关内容。经过修修改改之后，就可以获得邀请码。

OpenFOAM虚拟机

任务要求： 参考连接：http://dyfluid.com/docs/install.html 里面有很多虚拟机是20G的，你需要把20G的虚拟机，换成100G的。

注意，你需要做3个虚拟机，获取一个注册资格。

CFDEM虚拟机

任务要求： 制作CFDEM虚拟机。

foam-extend 3.2 虚拟机

任务要求： 制作foam-extend 3.2 虚拟机。

如何做：做好虚拟机后，邮箱发我即可。注意，虚拟机内不要有任何无关的软件。要一个纯净的系统。

OpenFOAM验证算例

任务要求：参考 dyfluid.com 的相关内容，比如这个段塞流算例 http://dyfluid.com/interfoam.html#testcase ，提供相关的内容素材。注意，算例不能是自带的tutorials

如何做：发送给我相关资料，我看看是否有意思，有意思的话，我就放上去。会标记你的贡献。

ICEM结构网格案例

任务要求：提供有意思的、网格划分技巧高超、质量高的ICEM结构网格。需要包括tin文件、blk文件。

如何做：把网格截几张图发给我，我看看。可以的话，提供测试账号。你撰写一下关键技术。修修改改。就可以获得一个注册资格。

CFD界公众号撰写

任务要求：参考CFD界以往公众号的内容风格，自选素材，撰写相关文章。

一篇文章可以获取一个注册资格。

如何做： 发送给我相关资料，我大体评估下能不能发，能发的话就可以。

回答技术问题

回答下面的技术问题。提供专业的讨论：

https://www.cfd-china.com/topic/4559 关于multiphaseEulerFoam中固相摩擦黏度的计算问题
~~https://www.cfd-china.com/topic/4829 tecplot处理wall边界上的速度不为0~~
https://www.cfd-china.com/topic/4795 关于高分辨率高精度格式的算例验证

邮箱： li.dy@dyfluid.com
请务必评估好自己能不能做。节约沟通时间。邮箱不接受答疑。

李东岳

细网格算例网格数多少？

李东岳

@micro

第一个是中心格式吧？跟格式有关。

扩散项两个对角的元素是因为边界条件的作用。

李东岳

@wying 对，非常老练 :-)

李东岳

find src applications -name "*.L" -type f | xargs sed -i -e 's=(YY_FLEX_SUBMINOR_VERSION)=YY_FLEX_MINOR_VERSION < 6 && \1='

李东岳

@cresendo 太复杂了 1 2个小时写不完

李东岳

@hongjiewang 代码什么样的

李东岳

@wying laminar+高精度格式就是DNS，laminar+低精度格式就是quasi-DNS

是等同的

李东岳

细网格看起来像是发散了... 原因不详

李东岳

@zhy2022

phi是一个surfaceScalarField，和volScalarField不矛盾吗？

参考 http://dyfluid.cn/无痛苦ns方程笔记.pdf 6.17节

后面那个用phi

另外，你那个0.5找到理论依据了么

李东岳

@wying 你确实可以用potentialFoam初始化一下流场，但是实际操作中，我感觉还不如费事的，各种设置...我都是用simpleFoam上来直接就算，先用simple，再用simpleC。

李东岳

@chonp 如果是周期性边界，就不存在进口了啊

李东岳

@yu_tian libCoxVoinovAngle看起来是一个自己写的库，貌似你没有这个库程序

李东岳

非常感谢！已更新

最后是都忽略了，这样就没问题了吧

李东岳

感谢分享代码！

但是去debug别人的代码很费事...尤其这种算两步就发散的问题

李东岳

原来这样。

最近在算好多LES算例，直接把服务器当主机用了。日常用没问题。但是发现最大的问题就是word文件的编辑。因为我这面很多都是别人发来的word文件我在上面继续编辑。但是linux下面用各种word编辑软件都会出现格式问题。有些文档就两页，但是导出pdf的时候格式就不对。有些是商务文件，格式对不齐，感觉不是很专业。

我先试试wps，看看会不会好一些。

李东岳

我有空试一下。

老铁知道linux这面有什么方便处理word的东西么。默认的libreoffice，稍微复杂一点的文档就格式有问题。office online有点卡，格式有时候也出现问题。linux下开虚拟机开word更卡..

李东岳

meanDiameter
    {
        type            coded;
        libs            ("libutilityFunctionObjects.so");
        name            error;
    
        codeExecute
        #{
            const volScalarField& d =
                mesh().lookupObject<volScalarField>("d.alpha.oil");

            scalar d32 = d.weightedAverage(mesh().V()).value();
    
            if (Pstream::master())
            {
                std::ofstream file;
                file.open ("d32", std::ofstream::out | std::ofstream::app);
                file << mesh().time().timeName() << " " << d32 << "\n";
                file.close();
            }
        #};
    }

我最近也凑巧要用这个，我用的上面这种方式。

李东岳

你调用volScalarField p，然后计算一下fvc::grad(p)不就可以了么

李东岳

给你个思路，用codedFixedValue可以么

李东岳

@bestucan master pdf这个linux使用方便么这个okular不能跟adobe那样进行标注 okular的标注很奇葩搞出来是这种的：

2021-05-07 08-39-40屏幕截图.png

李东岳

这是模拟泡泡流么？好像跟我这个文章里面玩的东西是一个东西

https://www.jianguoyun.com/p/DfaaRYkQ9s3ZBhjG5fED

李东岳

你这个圆圈外面的那个圈放大一点就好了啊？

李东岳

pisoFoam不可以考虑重力

interFoam可以考虑中立，你这么操作就是对的，没有别的操作了。

李东岳

就是发散了，deltaT太小了，比如1e-15, 1e-16就没有步进

李东岳

@merlinizlf 在自动调节时间步长循环计算某一时间步中说：

一直在循环计算某一个时间

这不是闹鬼了么？时间旅行？一直循环在1-2秒？算完了2秒又回到1秒了？

李东岳

各位简直就是大佬中的大佬厉害厉害

李东岳

ipmitool -I lan -U ADMIN -H 192.168.1.18 sensor thresh FANB lower 100 150 200

李东岳


Courant Number mean: 0.143461 max: 0.780362
deltaT = 0.000740741
Time = 0.502963

Weighted Re = 9.08136, weighted breakage frequency = 0
MULES limit is done
MULES: Solving for alpha.oil
MULES: Solving for alpha.water
oil fraction, min, max = 0.003 0.000631032 0.00359543
water fraction, min, max = 0.997 0.996405 0.999369
Phase-sum volume fraction, min, max = 1 1 1
DILUPBiCGStab:  Solving for n, Initial residual = 0.000902233, Final residual = 1.49928e-09, No Iterations 1
    Inconsistency value, average, min, max = 2.52405e-11 -0.000617714 0.00282076
    Droplet diameter max, min, mean = 0.0002, 0.0002, 0.0002
MULES limit is done
MULES: Solving for alpha.oil
MULES: Solving for alpha.water
oil fraction, min, max = 0.003 0.000631032 0.00359535
water fraction, min, max = 0.997 0.996405 0.999369
Phase-sum volume fraction, min, max = 1 1 1
DILUPBiCGStab:  Solving for n, Initial residual = 0.000901165, Final residual = 1.49419e-09, No Iterations 1
    Inconsistency value, average, min, max = 2.51988e-11 -0.000617713 0.00282694
    Droplet diameter max, min, mean = 0.0002, 0.0002, 0.0002
Constructing momentum equations
oil min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00559932, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00565399, No Iterations 6
GAMGPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 3.20997e-08, No Iterations 15
DILUPBiCGStab:  Solving for epsilon.water, Initial residual = 0.000218679, Final residual = 6.3486e-09, No Iterations 2
DILUPBiCGStab:  Solving for k.water, Initial residual = 0.000310738, Final residual = 1.00525e-09, No Iterations 2
ExecutionTime = 91.01 s


Courant Number mean: 0.14345 max: 0.780081
deltaT = 0.000740741
Time = 0.503704

Weighted Re = 9.07954, weighted breakage frequency = 0
MULES limit is done
MULES: Solving for alpha.oil
MULES: Solving for alpha.water
oil fraction, min, max = 0.003 0.000631032 0.00359535
water fraction, min, max = 0.997 0.996405 0.999369
Phase-sum volume fraction, min, max = 1 1 1
DILUPBiCGStab:  Solving for n, Initial residual = 0.000900074, Final residual = 1.47821e-09, No Iterations 1
    Inconsistency value, average, min, max = 2.51619e-11 -0.000617711 0.00283307
    Droplet diameter max, min, mean = 0.0002, 0.0002, 0.0002
MULES limit is done
MULES: Solving for alpha.oil
MULES: Solving for alpha.water
oil fraction, min, max = 0.003 0.000631032 0.00359527
water fraction, min, max = 0.997 0.996405 0.999369
Phase-sum volume fraction, min, max = 1 1 1
DILUPBiCGStab:  Solving for n, Initial residual = 0.000899041, Final residual = 1.47284e-09, No Iterations 1
    Inconsistency value, average, min, max = 2.51262e-11 -0.00061771 0.00283924
    Droplet diameter max, min, mean = 0.0002, 0.0002, 0.0002
Constructing momentum equations
oil min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00559881, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00565348, No Iterations 6
GAMGPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 3.20563e-08, No Iterations 15
DILUPBiCGStab:  Solving for epsilon.water, Initial residual = 0.000217858, Final residual = 6.31719e-09, No Iterations 2
DILUPBiCGStab:  Solving for k.water, Initial residual = 0.00030924, Final residual = 9.97882e-10, No Iterations 2
ExecutionTime = 91.12 s


Courant Number mean: 0.143439 max: 0.779868
deltaT = 0.000740741
Time = 0.504444

Weighted Re = 9.07807, weighted breakage frequency = 0
MULES limit is done
MULES: Solving for alpha.oil
MULES: Solving for alpha.water
oil fraction, min, max = 0.003 0.000631032 0.00359526
water fraction, min, max = 0.997 0.996405 0.999369
Phase-sum volume fraction, min, max = 1 1 1
DILUPBiCGStab:  Solving for n, Initial residual = 0.000897968, Final residual = 1.45572e-09, No Iterations 1
    Inconsistency value, average, min, max = 2.50956e-11 -0.000617709 0.00284538
    Droplet diameter max, min, mean = 0.0002, 0.0002, 0.0002
MULES limit is done
MULES: Solving for alpha.oil
MULES: Solving for alpha.water
oil fraction, min, max = 0.003 0.000631032 0.00359518
water fraction, min, max = 0.997 0.996405 0.999369
Phase-sum volume fraction, min, max = 1 1 1
DILUPBiCGStab:  Solving for n, Initial residual = 0.000896997, Final residual = 1.44647e-09, No Iterations 1
    Inconsistency value, average, min, max = 2.50645e-11 -0.000617707 0.00285154
    Droplet diameter max, min, mean = 0.0002, 0.0002, 0.0002
Constructing momentum equations
oil min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00559871, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00565306, No Iterations 6
GAMGPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 3.20133e-08, No Iterations 15
DILUPBiCGStab:  Solving for epsilon.water, Initial residual = 0.000217072, Final residual = 6.28692e-09, No Iterations 2
DILUPBiCGStab:  Solving for k.water, Initial residual = 0.000307752, Final residual = 9.91092e-10, No Iterations 2
ExecutionTime = 91.23 s

各位大佬，这是上面的log文件，我想从中提取几个关键的值：

Time = 0.000576911里面的0.000576911
Inconsistency value, average, min, max = -1.34441e-17 -0.000583605 0.000718129里面的0.000718129

目前我使用这个语句：

awk '{if($0~"Inconsist") print}' log

可以提取出来一部分，但是达不到我的需求。我最终想要的是两列数据，类似：

0.000576911  0.000718129

上面是一行的，第一列Time，第二列是后面那个

有没有大佬会搞这个？

李东岳

@小狗狗在关于多块网格并行剖分的负载平衡问题中说：

没有拓扑连接关系的分块网格呢？

这是什么意思没太懂

李东岳

@OItoCFD 在新写的程序单核计算无问题，但无法并行计算？中说：

mpi并行那些需要考虑吗？

一般不需要考虑

除非你的surface场进行一些特殊的插值

李东岳

@Micro 不，不是那个。我说的主要是alpha方程的子循环。你知道多相流这面alpha方程一般要显性离散，因此要加子循环的。

李东岳

感谢分享恭喜你debug成功

李东岳

@Samuel-Tu 不涉密可以把算例发上来玩玩

李东岳

最后网格做完了之后数量查的多么

你两个终端同时跑就应该能感觉出来那个步骤卡的慢了我看你blockMesh的edge大约是gmsh的2倍

李东岳

无厚度的面不应该是interface啊，interface是两个重叠的面，无厚度的面就是一个面。

李东岳

@hongjiewang 你这个interface是有两个面么，类似interface1和interface2，还是一个无厚度的面

Posts made by 李东岳