CFD中文网

http://dyfluid.com/pinn3.html 李老师您好，我想请问两个问题：
对于这段代码：

auto cont1 = dudx + dvdy;
        auto mom1 = 2.0*u_pred*dudx + v_pred*dudy + u_pred*dvdy 
          - nu*(dudxx + dudyy) - nu*(dudxx + dvdxy) + dpdx;
        auto mom2 = u_pred*dvdx + v_pred*dudx + 2.0*v_pred*dvdy 
          - nu*(dvdxx + dvdyy) - nu*(dudxy + dvdyy) + dpdy;
        auto loss_cont1 = crit(cont1, 0*cont1);
        auto loss_mom1 = crit(mom1, 0*mom1);
        auto loss_mom2 = crit(mom2, 0*mom2);

1. mom1中u_pred*dudx为何要*2
2. 粘度项为何写成 - nu*(dudxx + dudyy) - nu*(dudxx + dvdxy)，这好像并不太符合NS方程中粘度项沿X方向的分解？

请李老师指点一下

假如边界是 wall，是否有边界条件能设置液滴和边界的接触角？来体现wall的亲水性和疏水性

我看了mutiPhaseInterFoam,设置的是不同流体之间的接触角，没找到液体和边界的接触角设置。

@李东岳

李老师能不能把ML部分讲的最后一个例子的完整代码分享一下？
当时着急赶车最后一个没听完~
多谢李老师

邮箱：
wang.qi@hit.edu.cn

@李东岳 李老师，在写fvm::laplacian项之前已经通过体积分数alpha计算过nu了。因为我觉得直接修改turbulence::类里面的nu很麻烦。所以我就自己定义了提及标量场nu然后正常使用alpha计算之后带入laplacian，所以获得了与文献一致的结果。

最近在尝试改不同的两相流求解器，遇到了一个问题。
当前OF的两相粘度项是通过湍流模型计算的。
turbulence->divDevRhoReff(rho, U)
当我的case不需要湍流使用laminar时，结果和文献结果不匹配：

当我改成直接用拉普拉斯求解：

• fvm::laplacian(nu*rho, U)
  结果：
  

此时与文献结果一致，请问李老师这是什么原因啊？

RT，想算个液滴结冰。 谢过

@李东岳 收到 多谢李老师

@李东岳 OF勉强会用，简单的求解器都能写或者能改。太复杂的目前也没需求，基本满足当下的需求啦已经~

看了李老师的libtorch, 但是我个人对python熟悉一点，C++真的不太会，请问如何将pytorch挂载到0F呢，期望李老师能出个简单的例子或者说明~~

您能解释一下为什么会有rhoCpPhi=rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2吗？这个是什么物理意义？

有相变的温度方程形式不应该是很简单的这种吗：

这个rhoCpPhi是怎么回事？方程的形式又是怎么回事？
请知道的朋友帮忙解读一下谢谢

@李东岳 多谢李老师，可惜我现在人在国外向日葵貌似用不了。但是我跟您说下最终测试结果。我们可以使用./Allwmake -j -k把那一步跳过去。这样编译出来的还是能用的。从日志上来看，貌似是这个FoamtoFireMesh这里
3b01921f-39c2-4d31-a36d-b8310fcd497a-image.png
而且，如果我不下载代码包，直接用官网的命令在线安装是没有任何报错的。
我把报错的日志发给您看下，您看看是否有头绪。

链接：https://pan.baidu.com/s/1RcbLj7qFiBwl3yK4eEhDjA
提取码：wo84

@李东岳 是的，我测试了一下，相同的环境，我安装v2012(我之前最常用的版本)是和原来一样没有问题的。但是如果我安装比较新的，比如v2206和v2306都会出现相同的问题。并且，出现这个问题之后，虽然输入blockMesh也能出现预想中的信息，但是直接跑case的话会显示找不到求解器。所以还是未能编译完成。

PS，我编译都是下载压缩包解压之后 ./Allwmake的方法。可能是更新之后哪里和原来不一样了

@李东岳 好难受啊 我也第一次遇到，以前都没遇到过。不知道怎么解决

之前用的都没问题，新买了个电脑，官网下载了新的linux。按照东岳流体中得步骤，先
sudo apt update &&
sudo apt-get install -y build-essential flex bison git-core cmake zlib1g-dev libboost-system-dev libboost-thread-dev libopenmpi-dev openmpi-bin gnuplot libreadline-dev libncurses-dev libxt-dev

安装完之后对OpenFoam进行了编译，结果报错如下：

虽然报错，但是输入blockMesh 仍然能够获得正确的输出，请问有人知道这个错误是什么原因吗？看着别扭

@李东岳 我想做沸腾过程，下面是空气，上面是液体，其实是需要求解整个场的电势方程。奈何气液电导率差距过大，所以一直解不准

@李东岳 多谢李老师，请问openFoam中，您提到的这个多计算域如何实现？有什么关键词吗

@李东岳 对，就是想把界面的以及界面以上的液相电势全部设置为0

@李东岳 应该不是李老师。具体问题是这样，扩散项的数值扩散问题太大，我实在是无法解决（我看到了Cfd界发的那个文章，但是对于我并不适用）。 但是通过我和外导的数值推演之后，我们想到了另一种解决办法。 就是希望气液界面及界面以上的数值都是0，只解开界面下方，气相的扩散方程，这样就可以避免气液两侧扩散系数差距太大导致的数值扩散。
但是在of中我目前还没找到合适的办法实现它