CFD中文网

@Samuel-Tu 会的

@东岳 感谢李老师！

看到官方的求解器，根本摸不清楚边界条件和有限体积是如何交互的。有人能解释一下这和我上述思路的差别和联系吗？

粗略看了下，这个求解器完全没用openfoam的思路来，采用的他自己的思路，就像你说的，他在开发这个求解器的时候，可能是有其他的代码（比如同组的Fortran）转过来的。并且完全是面向过程的思想。

如果你要用OpenFOAM求解的思路，还是看OpenFOAM的求解器比较好他这个就是把求解的过程自己实现了，然而没有使用OpenFOAM自带的方法。

他的main函数已经写的很清楚了：

int main(int argc, char *argv[]) { #include "setRootCase.H" #include "createTime.H" #include "createMesh.H" #include "setInputValues.H" #include "createFields.H" #include "readFluxScheme.H" /// Time step loop /// Posts the non-blocking send/recv of fields long int iter = 0; while( runTime.loop() ) { /// 构造通量 #include "constructFaceFlux.H" /// 有限体积离散 #include "sumFlux.H" /// 边界通量修正 #include "boundaryFlux.H" Info << "Iteration = " << ++iter << " "; /// 矩阵计算 #include "stateUpdateLTS.H" Info << " Max residue = " << rhoResidMax << endl; /// 输出结果 runTime.write(); } return 0; }

这几天在看有限体积法的书，大佬的论文、《数值传热学》和 《The FVM inCFD OpenFOAM..Matlab》都看了，逛论坛发现都在推荐《Computational Methods for Fluid Dynamics》这本；去SpringerLink一搜正好今年出了新版，网上资源好像还不多，在这里分享一下： 点击下载

@cccrrryyy 受教了，谢谢您的指导，您说得对极了，假如要真能收敛，应该最大默认的最大迭代步数是够用的了，问题还是出在格式上，我把格式调成PISO，压力调成standard，然后步长变成0.0001的一半之后，这个问题就解决了，cfl数在fluent里输出的问题我会自己琢磨的，谢谢您！

换个湍流模型
换个壁面函数

@supersoldier OpenFOAM大量使用了C++面向对象的设计思想，不能用一般面向过程的编程思维去阅读代码。了解面向过程的设计思想有助于去理解OpenFOAM的框架和具体实现。

好的，谢谢，我已经找到输出的地方了!

方程实际要求的是：
\begin{equation}
\frac{\partial \rho\mathbf{U}}{\partial t}+\nabla\cdot(\rho\mathbf{U}\mathbf{U})=0
\end{equation}
对应：

fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U)

连续性方程：
\begin{equation}
\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{U})=0
\end{equation}
对应：

- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

因此实际求解的就是：

fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) + turbulence->divDevRhoReff(rho, U) );

把- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)包含进来的原因为这有利于使求解变量有界以更好地收敛。

@cfd-china 我看你网站上有这个模型文档的翻译。你是专门研究这个的？

:kiss:

https://codepen.io/pkra/pen/EPeKjo

@HITSC30 icem不是很熟练，我去试一下

是不是可以设置一下relTol

@东岳 如此说来，除了加密网格，没有好的解决方案啦:zoule: :zoule:

你用什么做的网格？如果是icem的话改一下scale factor就可以。

我i比较喜欢用进行建模和mesh处理然后导出STL，你这个报错有点像组合的STL文件，直接导出STL的时候连接位置有缝，没有闭合，类似于一个圆柱，如果你整体mesh处理或者直接导出stl是可以的，但是你把圆柱炸开，形成三个面STL，之后不管是组合还是分别mesh导出STL或者直接导出STL，连接处都会有缝，在liggghts会报错

weblate ，或者 github。
debian中文手册翻译
但学习成本都挺高。

但是特性好，版本控制啥的，weblate 就是个轻量的版本控制系统。po 文件可以统一术语。后台就是这个软件 gettex

曾经想参与 debian 手册的翻译，没时间入门这些东西。:xinlei:

也有会用软件一个人带着其他人翻，用邮件列表报翻好的片段，然后统一上传。

对这种不是那么重的翻译，用专业的方法成本太大。那种翻译框架搭好了是照着往多种语言翻译的。用土法对个人压力过大，这种压力不只是时间紧迫性和工作量上的。没有好的工具，hold 不住相当碎片化的翻译数据，邮件勉强可以，版本控制器最适合（全历史记录）

可以扣搜扣搜有没有日文版和日文翻译者，取取经:chitang: 一般好多技术文档有国际翻译，八成就有日语翻译。
日本 OpenFOAM 大本营：googlegroups
组织上决定派你去踩点，加油:chigua2:

@noodles 开个玩笑，我觉得你这个问题过于简洁，不知道是周期性边界条件还是什么意思
udf可以实现

@东岳 不好意思，表述的有点问题，fluent 里面没有质量流量出口，只要质量流量入口，所以我把出口设置为负的质量流量入口边界。