CFD中文网

@木木枭 厉害呀，连续伴随、离散伴随都挺难的，我想试试PINN求反问题 想当论文的最后一部分，还不清楚可行性

@李东岳 清华最近有工作“2024-工热学报-基于物理信息神经网络的气膜冷却湍流模型反演学习”，用的Si-LU激活函数:）

俺才刚开始用就要关么，，

@李东岳 李老师，就是初始流场，比如均匀初始化或者用一个RANS收敛的流场初始化

看完大佬们讨论，对OpenFOAM+神经网络模型有初步了解了，还有些疑问，想请教下各位老师：
用神经网络来确定非线性模型中的一些参数，并在每一步计算都调用神经网络，对初场的敏感性高吗？如果初场不同，神经网络更新的参数有很大区别，感觉可能导致收敛到不同结果？

感谢感谢！

@木木枭 还做伴随吗，感觉很不容易，对PINN反演，确定方程特定参数比较感兴趣，想先跟着李老师课程摸索下PINN

应该写成fvc::div(turbulence->devRhoReff() & U)吧，最近准备研究对流换热问题，不太清楚OF rhoSimpleFOAM、rhoPimpleFoam中不考虑粘性热项影响大吗，打算加上这一项试一试

@jinlinna 您好，想问下为什么alphat初场会对计算结果有明显影响呢，在0时刻后，alphat不还是用nutSGS/prt计算么

我也想知道L怎么定义的，OpenFOAM2312，这么解释

“The order of elements of the integral scales for the vector-based condition is “(Lxu Lxv Lxw Lyu Lyv Lyw Lzu Lzv Lzw)” and for the scalar-based condition is “(Lxi Lyi Lzi)”, where \c x here indicates the streamwise components and \c y/z lateral components. The streamwise components (i.e. \c x) are Lagrangian time scales in units of seconds and the remaining components are length scales in units of meters.”

Lxi对应拉格朗日时间尺度，Lyi/Lzi对应长度尺度

采用LUST 介于Guass linear和limitedLinear 0之间

@李东岳 div(phi,U) Gauss limitedLinear 0;

div(phi,U) Gauss linear;

Q=1e6

老师好，我用混合方法计算后台阶，采用div(phi,U) Guass linear 和 div(phi,U) Guass limitedLinear 0，瞬时速度场差距不小，div(phi,U) Guass limitedLinear 0依然存在明显耗散，小尺度不如div(phi,U) Guass linear丰富，理论上Guass limitedLinear 0和Guass linear不应该一致吗

我之前算过高亚声速叶栅，用rhoSimpleFoam求解器收敛不了，换成LUSGS就很好收敛

受教了，exp()，log()都不能计算量纲场，求ln(A)运算，用log(A/Aref)*Aref

受益匪浅呀

也可能需要模版参数特化的某个函数，指定的模版参数有问题，希望有用哈，我遇到的问题是，在自定义基类中把：
BasicTurbulenceModel::correctNut();
换成了TurbulenceModel::correctNut();
解决方法是传入
BasicTurbulenceModel和TurbulenceModel两个模版参数，
或者BasicTurbulenceModel::correctNut();只在子类中出现

@吴琼 问题解决了吗，我最近也遇到类似的问题，貌似是和构造函数继承类有关，编译不会错，初始化就出问题

@李东岳 李老师，我意思是SA和RSM算得的压力是真实值，因为两者计算时没有2/3k参与，kOmega模型算得的压力跟真实值差了2/3ρk

感谢感谢，这应该对SA或RSM模型没影响，对kOmega有些许影响

OpenFOAM RANS方程中默认将Boussinesq假设的各向同性分量▽·(2/3ρkI)带入压力梯度▽p中，那算得的压力p是否应该减去2/3ρkI才是真实压力？OpenFOAM输出压力场时是否做了类似处理，未找到相关代码，望得到指正