CFD中文网

@tiankai 看不出来什么问题？问问李老师@李东岳

@tiankai

1. 查看/修改OpenFOAM执行wmake时采用的Gcc版本，OpenFOAM/wmake/rules/General/Gcc/c++

   CC          = g++$(COMPILER_VERSION) -std=c++14
   

2. pytorch官网下载libtorch
   
   2.3.0版本依赖C++17
   
   
   2.0.0版本依赖C++14
   
   
   1.10.0版本依赖C++14
3. 修改/home/***/.bashrc，设置环境变量

   export LD_LIBRARY_PATH=/home/***/OpenFOAM/libtorch/lib:$LD_LIBRARY_PATH
   

4. wmake你的库（比如李老师贴的代码）PINNFoam
5. 最后调用PINNFoam

@李东岳 在 OpenFOAM外挂libtorch 中说：

scalarField P0internal = P[0].primitiveField();
scalarField P1internal = P[1].primitiveField();

李老师，这个P[0]和P[1]分别是什么物理量？怎么得到呢？

@李东岳 可运行，👍👍👍谢谢李老师！

@李东岳 在 OpenFOAM外挂libtorch 中说：

#include "net.H"

李老师，不知道是否方便贴一下net.H？
感谢！

就是你说的，计算得到的结果是一样的。

https://www.cfd-china.com/topic/4627/后处理时候singlegraph可以提取两条线的值吗

你想在openfoam里面直接计算？有困难，虽然网格是结构化的形式，但是非结构网格求解器仍然不方便找到你想要的信息。
可能还是praview+宏，批量处理n条线，把数据提取出来会简单些。

我感觉比较简单的应该是，paraview用plotOverLine，将数据另存为文件，之后在对这条线上的量计算。

@李东岳 在 fluent用LES算水翼边界层选择中心差分引起的震荡如何消除 中说：

库朗数减小一些，比如0.3或者0.2左右来尝试纯粹的linear。

这个倒是找到一个很有意思的格式（"Blends two schemes based on the local face-based Courant number"），

    div(phi,U)      Gauss CoBlended
                    0.01                  // Co below which scheme1 is used
                    linear                // scheme1
                    0.05                  // Co above which scheme2 is used
                    LUST grad(U);         // scheme2

这个格式可以根据Co数自动调整格式，可以实现按区域混合各种格式（虽然不是按权重混合）。

@李东岳 在 fluent用LES算水翼边界层选择中心差分引起的震荡如何消除 中说：

权重是0.25与0.75，不过可以相对简单的改一下权重

这个权重0.25和0.75是可以调整吗？

@心里的小漩涡 在 fluent用LES算水翼边界层选择中心差分引起的震荡如何消除 中说：

有试过二阶中心差分格式中耦合小部分一阶迎风格式吗，比如给个0.02%的一阶迎风格式，之前用fluent算3900的圆柱绕流，用纯二阶中心出现过数值振荡容易发散的情况，改成混合之后收敛性好很多，精度也还行，没有很明显的下降

@李东岳 李老师，这个在OpenFOAM里面能直接设置吗？

空的是指真空还是空气？
空气就两相流（入口是water，计算域本来是air），真空就知识盲区了。

假设问题中流体不可压缩，左侧有流体流进+1，右侧当然有流体流出-1，瞬时完成。

我的意思是
1.我只试过SST k-omega湍流模型，这个没问题。 （）
2.大涡模拟里面没有k方程，也就没这回事了。

SST k-omega湍流模型应该可以这么实现

按我的理解，如果你滤波器的尺寸跟网格尺寸一样，可以不用显式地使用滤波器，因为有限体积法用的体心值取的就是单元平均（？），相当于做了滤波操作。
如果要加滤波器，滤波器的尺寸应该是根据问题要捕捉的涡的大小来确定的吧。

@chapaofan

这里面有矛盾的设置。
每条边上的单元数改为10试试（下图），再慢慢调整。

@chapaofan 在 blockmesh生成网格总报错没有相邻面 中说：

hex (14 11 12 13 29 30 31 32) (10 10 1) simpleGrading (1 1 1) //3

上面的顺序错了。
可能是下面这个：
hex (14 11 12 13 32 29 30 31) (10 10 1) simpleGrading (1 1 1) //3

@Prometheus10

请教，symmetry边界，直观上理解是没有界面法向的速度变化吧，没有速度变化，为什么有切应力？