CFD中文网

理论上了解一些求解方法，软件没有用过 :chitang:

@东岳 我在controldict 试过一次 好像不管用

@搬运工 后来我发现我这个问题不是周期性问题，虽然结构具有周期性，谢谢你的解答:happy:

@random_ran :upset: 老师让做这方面的，没办法啊

shurufa diu shi
Could you please check the other terms? e.g., lift force

@李东岳 哦，明白，谢谢李老师

@邙山老王 看图

不同的语言在这块不太一样，估计是写惯了其它语言，转到c++就惯性思维了:duang:

@东岳 感谢东岳大哥

@LXJCFD 波浪问题很少用 LES，波流问题倒是可以用 RANS 或 LES，取决于是否有大分离流动。

液体射流喷射雾化后，在高温高压的环境中燃烧

不是很熟悉你这个物理过程。普遍性原则是你们是否需要某一内部坐标点的数量密度函数，如果需要就需要重构。比如蒸发过程就需要当内部坐标为0时候的权重，因此需要NDF重构。

0_1532484662934_6ade2d8e-6cbc-492d-9683-fcd95f369558-image.png

图片出自文献《AMESim and 3D CFD EOLE models for Diesel Injection》。

@昝磊 好像是做凝固，李喜老师那的:huahua:

@random_ran

不满足微分方程，但是满足积分方程的解就是弱解吧。本来俩方程应该等价的。

其实我觉得关键单不在于什么是弱解，而是这帮数学家为啥要引入弱解这个概念，我的理解是丫们发现物理上一些有意义的解不严格满足他们推导出来的微分方程，但是似乎又差得不多，然后整吧整吧搞出一个弱解的概念来。

有点感觉想和稀泥。。。。

今天青年节，下午放假 :chigua:

很普适性的问题。如果你说的是如何判定非稳态的收敛性的话，可以监控一些参数，如果稳定波动一定周期，就可以认为收敛了。

phi是用来调用判断速度方向的： surfaceScalarField phiBD(upwind<scalar>(psi.mesh(), phi).flux(psi));

@kasitefy
非常对。

另外一个问题是这种加和的数值积分忒慢了

@李东岳

常见的问题是求$\mathbf A \cdot \mathbf x=\mathbf b$中的$\mathbf x$，可以有直接解法（如LU），也可以用迭代法（如Gauss-Seidal），求逆矩阵肯定是下下之选，原因有几点： 操作太多，$~O(N^3)$; 稀疏矩阵的逆矩阵不是稀疏的； 你要求得$\mathbf x$还要再乘一次，$~O(N^2)$ 逆矩阵求了之后没有用，因为CFD方程组是非线性的，每次求的矩阵$\mathbf A$都不一样，所以你求得的逆矩阵只用一次，下次迭代又得重新求解。如果是针对线性问题，每次针对不同的$\mathbf b$求对应的$\mathbf x$，用逆矩阵还勉强说得过去。 现代的CFD程序似乎迭代法用得更多，JFNK, IDR(s)，最不济还有CG和BiCGSTAB.有的连矩阵都不需要组装，何苦去求逆...

Fluent我们不太清楚开发人员是如何做的数值稳定性工作。不过如果你看OpenFOAM源代码会发现OpenFOAM里面对湍流方程的守恒、有界、减少数值扩散、以及针对网格质量的自适应数值调整是非常非常多的。没有一个方程在openfoam中是直接求解的。So~