shu-osher-wwc.tar.xz

shuosher.JPG

@xpqiu 哈哈哈必应第二条就是，谢谢邱老师！faf5e653-d711-4b05-aa92-a1af89def6f5-image.png

@capillaryFix 是的，我求解了Eqn.30 -> Eqn.34 ->Eqn.33，但是Eqn.33的结果用在哪里了？

@李东岳 在 如何根据stl模型定义cellZone？ 中说：

你需要看cellZone，不是看stl。paraFoam可以打开cellZone

老师，之前尝试过了，在一个snappyhexmeshdict中好像没有办法同时生成壁面边界模型和cellzone（即A.stl生成壁面模型，B.stl为cellzone区域），这个我在polymesh文件夹中检查了cellzone文件，的确为空。不过，后续我通过两次运行snappyhexmesh解决了，第一次只配置壁面边界模型，第二次只生成cellzone。

@xuluxi18 才看见，rhoCentralFoam能跑通，rhoPimpleFoam跑不通么

@chon01 这也只是我的猜测，因为我按了一下计算器，lg和ln函数算出来的结果差距确实比较明显，希望能帮到你吧。

@欧阳 我也尝试了一下，目前编译没啥问题

压力是没下去。原因不知道了，需要详细debug

搞出来了，一个小bug卡2小时

捕获.JPG

一说到这个库朗数，我看了一下，这面是用最大波速来计算的CFL。我用compressibleCourant对比了下。二者一个0.17，一个0.14。还有点差距。波速是求解之前预估的，phi是求解之后。我估计是差距的来源。不过没有细看代码，自己想的。

我今天在准备这个算例的时候发现，http://dyfluid.com/class.html#hyperbolic 特征速度和传输速度有2倍的差异。针对这种偏hyperbolic的方程可以用俩种方法计算库朗数：

特征速度 传输速度

所以，rhoCentralFoam应该就是用的特征速度计算的库朗数，而不是传输速度。如果想看传输度的库朗数，需要用compressibleCourant.H来计算。

CFD界的竞赛：https://mp.weixin.qq.com/s/0PjOnv09LerabVSZyp-g9Q

现在最新的机理中大部分都带有plog反应，形如44f5b470-8e44-46fd-a399-0d52188a215a-image.png ，但是OpenFOAM无法识别plog反应，因此没办法直接把CHEMKIN格式的机理转换为OpenFOAM格式的机理，想问问各位大佬怎么才能解决这个问题，我尝试过改机理把plog反应拟合成其他反应，但这样又会导致机理准确度较差，与原机理误差较大。

今天Ubuntu突然打不开了，错误信息如下所示：
89da0367-9d40-4dd4-ab7d-0948def2a515-image.png
在终端显示文件损坏
b95a8447-bf48-4559-ad9f-99568138d5da-image.png

在安装完Windows版ParaView之后, 也可以去WSL的.bashrc文件里面添加以下内容

alias paraview='"/mnt/ParaView/安装/路径/bin/paraview.exe"'

然后重新Source一下就可以使用paraview命令启动Windows版ParaView了, 可以直接打开WSL里面的文件夹, 无需在Windows目录下运行OpenFOAM算例

@李东岳 是的 我的意思是这个东西都能跑 所以应该不是case的问题 是那台电脑硬件的问题

@李东岳 谢谢李老师回复

让AI帮忙写了个脚本，供大家参考。原理是提前在当前的算例文件夹下创建好链接，链接到另一个硬盘已经创建好的时间文件夹。实测在/home/disk1/caseA目录下运行OpenFOAM求解器，可以将数据文件存到/home/disk2/caseB中。注意这个脚本得在/home/disk1/caseA目录下运行，不然ln命令总会有一些bug，有时候是没办法正确创建链接，有时候是会出现链接的嵌套。

#!/bin/bash # 配置路径 src="/home/disk1/caseA" dst="/home/disk2/caseB" # 创建目标目录 mkdir -p "$dst" # 处理单个 processor 的函数 process_proc() { local proc="$1" local dst="$2" local procname=$(basename "$proc") # 创建 processor 目录 mkdir -p "$dst/$procname" # 生成时间序列并创建时间文件夹（自动剪除尾随零） seq 8.0001 0.0001 8.5 | while read t; do formatted_t=$(printf "%g" "$t") mkdir -p "$dst/$procname/$formatted_t" # 创建符号链接，指向目标目录中的时间文件夹 ln -sfn "$dst/$procname/$formatted_t" "$procname/$formatted_t" done } # 导出函数以便并行使用 export -f process_proc # 收集所有 processor 目录 proc_dirs=("$src"/processor*) # 并行处理（使用所有CPU核心） parallel --jobs 0 --progress --bar process_proc {} "$dst" ::: "${proc_dirs[@]}" echo "processor 文件夹和时间文件夹创建完成。"

@李东岳 在 数据驱动CFD的意义是什么？ 中说：

博士后按道理来说应该已经看过很多sci了啊。现在很多sci都是啥用没有，纯粹发文章用...

太真实了，绝大部分的sci都是啥用没有

好像只要有网格尺寸分割切换的地方，中间过渡不够，都会出现这种问题，是不是还得尝试调整网格。

@xiaoyangfoamer
你试试这种方法：初始化的时候就在喷嘴后面设定一小段的气体柱，很短就可以。
之前用Basilisk的时候，看到过这样的处理，说是让求解器好收敛一些。不知道OpenFOAM是不是这样。