CFD中文网

@桑葚熟了 目前没有用这些集群了，天河没用过，北京超算云那边感觉服务和价格都还不错，EasyCAE这边很久之前用的了之前价格有点小贵，现在不清楚了。

之前试了好多次自己搞都没搞出来，现在能直接用了！

@Zhy2022

有没有什么进展或者文章见刊啊？

试试reconstructParMesh之后再reconstructPar？
如果直接处理的话可以运行 touch case.foam && paraview case.foam或者paraFoam -builtin，但是paraview没法直接处理拉格朗日场，可以试试用tecplot

@leilei operation pressure是0的话，只要有压力下降就会有负数

我又回来了，这次发现虽然写了用坐标去计算不同点的阻力系数，但是实际计算中发现还是均一的，以为是定义的变量被重新赋值了，重写了所有的中间变量发现也没办法解决，目前认为是坐标计算的不对或者openfoam porosityModel不支持非均一的阻力系数设置
如果想用场自定义阻力需要如何操作呢？直接改求解器加入源项，然后在0文件夹写初始场？

        forAll(cells, i)
        {
	    gamma_[zoneI][i] = Zero;
	    scalar x = mesh_.C()[i].x();
	    scalar y = mesh_.C()[i].y();
	    scalar z = mesh_.C()[i].z();
	    scalar dir = y;

            if (xxx)
            {
                gamma_[zoneI][i].xx() = xxx;
                gamma_[zoneI][i].yy() = xxx;
                gamma_[zoneI][i].zz() = xxx;

            }

            D_[zoneI][i] = Zero;
            D_[zoneI][i].xx() += pow(gamma_[zoneI][i].xx(), 1/3) * dXYZ_.value().x();
            D_[zoneI][i].yy() += pow(gamma_[zoneI][i].yy(), 1/3) * dXYZ_.value().y();
            D_[zoneI][i].zz() += pow(gamma_[zoneI][i].zz(), 1/3) * dXYZ_.value().z();

            D_[zoneI][i] = coordSys_.R().transformTensor(D_[zoneI][i]);

            // leading 0.5 is from 1/2*rho
            F_[zoneI][i] = Zero;
            F_[zoneI][i].xx() += 0.5 * fXYZ_.value().x();
            F_[zoneI][i].yy() += 0.5 * fXYZ_.value().y();
            F_[zoneI][i].zz() += 0.5 * fXYZ_.value().z();

            F_[zoneI][i] = coordSys_.R().transformTensor(F_[zoneI][i]);
        }

@TINGHAIK 下载了还没用过，感觉算气动应该挺不错的了，看里边的案例做的还挺有意思的，不过界面上应该还可以下点功夫，不涉及工程计算用现有的工具可以了。

@Samuel-Tu 最近几天看了看salome，感觉照着官方的tut操作一遍，使用上应该没啥问题了，除了操作有点繁琐，其他感觉用起来还挺顺手的，我先学学，学会了可能在b站录个系列教程。
salome的第一个tut就可以创建只有一层的网格，楼主有需要可以自己去看看。

@xpqiu u 好的，感谢！

@Samuel-Tu 好滴，gmsh和salome哪个好一些呢？

preCICE很amazing，用好了应该比自己改求解器简单高效 https://www.bilibili.com/video/BV1iT4y1P7e7/
（现在发帖竟然需要间隔300s ）

cfmesh好像能定义不同方向不同的加密

想在计算完成之后输出ddt(U)和div(phi,U)
直接运行postProcess -func 'ddt(U)', 提示ddt is not supported with the postProcess utility FOAM Warning : functionObjects::ddt ddt(U) failed to execute.
运行pisoFoam -postProcess -func 'ddt(U)' 提示in command line pisoFoam -postProcess -func div(phi,U) Placeholder value is <fieldName>
想用coded输出也一直没成功，因为已经有计算结果了，不想动求解器：

libs    ("libutilityFunctionObjects.so");
type    coded;
name    ddt;
codeWrite
#{

const volVectorField& U = mesh().lookupObject<volVectorField>("U");
//const volScalarField& phi = mesh().lookupObject<volScalarField>("phi");

// Write
volVectorField ddt("ddt", fvc::ddt(U));
ddt.write();

#};

@Stan同学 这个看起来比IBM啥的简单一些，我看有人也用它算心脏的流动，坐等出教程，哈哈哈

@Stan同学 可以看看这个流固耦合的接口precice

preCICE is the leading open-source coupling library for partitioned multi-physics simulations, including, but not restricted to fluid-structure interaction and conjugate heat transfer simulations.

Partitioned means that preCICE couples existing programs/solvers capable of simulating a subpart of the complete physics involved in a simulation. This allows for the high flexibility that is needed to keep a decent time-to-solution for complex multi-physics scenarios.

The software offers convenient methods for transient equation coupling, communication, and data mapping.

visualisation of how preCICE couples different solvers

今天解决了这个问题
先在$WM_THIRD_PARTY_DIR编译paraFoam并执行wmRefresh
到$WM_PROJECT_DIR/modules/visualization路径下执行 ./Allwclean && ./Allwmake && wmRefresh
转到$WM_PROJECT_DIR路径下执行 ./Allwmake && wmRefresh

命令：

cd $WM_THIRD_PARTY_DIR
./makeParaView
wmRefresh

cd $WM_PROJECT_DIR/modules/visualization
./Allwclean && ./Allwmake && wmRefresh

cd $WM_PROJECT_DIR
./Allwmake && wmRefresh

你的paraFoam没有成功编译，打开的其实是paraview，在paraiew中可以使用VTK格式查看拉格朗日场，运行一下foamToVTK

@李东岳 李老师，可能是您python版本的问题，我用的是python3.8，您应该也有python3但是系统调用的是python2.7，可以试试python3 plot.py运行。
坐等大佬改我的代码