CFD中文网

李博继续多出好文！

重叠网格这一块投资咱没有涉及。知识盲区目前是

@李东岳

常见的问题是求$\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.有的连矩阵都不需要组装，何苦去求逆...

@CFDngu 这个运算符重载真是。。。放右边结果是左乘。。。
无论是左乘还是右乘，都没有交换律。所以，应该都不是一个角度分两次旋转、旋转顺序的问题。

然后

6DoF定义的orientation叫做“rotation tensor”

这个，是因为，应该是这样的：所有船体与外界的交互，都是基于世界坐标系的。所以要把世界坐标系的力应用到船体上，就要把力之类的各种矢量张量旋转到船体坐标系去和船体交互。

还有另外一个猜测。两个坐标系之间的转换，可以通过三个旋转矩阵达到。也可以通过一个旋转轴和一个旋转矩阵达到。

我猜测，可能这个

0.5Rx-0.5Ry-Rz-0.5Ry-0.5Rx

是在组那一个旋转矩阵，或者找那一个旋转轴。如果有了这个特殊的旋转轴，那么船体的"orientation"就确实是“rotation tensor”了

以上都是我猜的_(:3」∠)_

哦，对了，这句"所以要把世界坐标系的力应用到船体上，就要把力之类的各种矢量张量旋转到船体坐标系去和船体交互。"不是猜的。是代码注释：

Angular momentum stored in body fixed reference frame. Reference orientation of the body (where Q = I) must align with the cartesian axes such that the Inertia tensor is in principle component form. Can add restraints (e.g. a spring) and constraints (e.g. motion may only be on a plane).

The time-integrator for the motion is run-time selectable with options for symplectic (explicit), Crank-Nicolson and Newmark schemes.

我应该理解对注释了吧:134:

Cp不是常数的情况下，焓和温度不是线性的关系吧。还是说你要找某个特定流体的焓随温度变化的关系？

最近用pointwise中的非结构网格生成了一个风力机旋转域中的block，检查质量发现最小角度0.14，网格质量特别差，但目前只会调整domain的网格质量，block的没看到相关教程，有大神可以指导一下吗？谢谢！

@random_ran

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

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

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

做二维后向台阶模拟，怎么产生不了流动分离呢？以为是边界条件有问题，在三维里面尝试相同的边界条件是没有问题的，是数值方法的问题吗？之前用的KEqn模型，是模型不支持二维的吗？有没有其他模型的推介

在Latex里面公式可以交叉引用，一般重要的公式我都可以起个名字例如\label{mom}，但是经常有一些不重要的公式，一篇文章公式很多，这些公式很难起名字也很难记住，你们怎么处理的额？

半新手，看其他论坛上可以实现这样的操作，即把一个复杂的几何体分成多部分分别划分网格之后再进行拟合，但是我自己试了下建块merge后导入到fluent中出错，请问该如何操作呢？不胜感激

有没有试过改project type？

我只做过把domain投影到database上面，在这个选项里面可以修改投影方式（版本：Pointwise V18.2 R1）

project.png

@jpzhao 有monitor和report的功能，具体看看相应版本的手册哈
捕获.JPG

0_1510344430862_草图.png

这本书我2015年买的，花了118欧，还让外导签了名。今天双11，琢磨亚马逊有没有什么书可以买，忽然发现这本书65欧了，便宜了一半啊！

这几年买国外书估计都花了1000欧了 :angry: 因为喜欢写写画画，打印的话A4纸太重太厚。

是OpenFOAM吗？
如果是的话参考
https://www.cfd-online.com/Forums/openfoam-meshing/61596-cyclic-bcs-pointwiseopenfoam-export.html

@李东岳 +1

时间离散格式应该是计算时间导数的格式，模板没有限制；而推进格式应该要求模板是单侧偏心的。http://www.scholarpedia.org/article/General_linear_methods是一般的时间推进格式，偏心模板最多也就是隐式格式。

计算域是否对称，可否切出1/3生成网格后旋转

谢谢，我刚把这个导入到程序中，但结果仍然不理想，解出得P'偏大，在修正速度时会导致速度变成负的，而且是进口速度得负几十倍，所以总感觉压力方程和速度方程不匹配，肿么破？？:crying:

@李东岳 网格是用fluentMeshing画的，严格按照周期性边界设置的，导入OpenFOAM之后，当matchTolerance很小的时候，使用createPatch设置周期性边界的时候会报错，

以下是 createPatch 代码

patches ( { name cyclic1; patchInfo { type cyclic; neighbourPatch cyclic2; transform rotational; rotationAxis (0 0 1); rotationCentre (0 0 0); matchTolerance 1e10; } constructFrom patches; patches (cyclic_1); set f0; } { name cyclic2; patchInfo { type cyclic; neighbourPatch cyclic1; transform rotational; rotationAxis (0 0 1); rotationCentre (0 0 0); matchTolerance 1e10; } constructFrom patches; patches (cyclic_2); set f0; } );

以下是报错信息

--> FOAM FATAL ERROR: More than one patch accessing the same transform but not of the same sign. patch:cyclic2 transform:0 sign:1 current transforms:1(1) From function Foam::label Foam::globalIndexAndTransform::addToTransformIndex(Foam::label, Foam::label, bool, Foam::scalar) const in file lnInclude/globalIndexAndTransformI.H at line 210. FOAM exiting

问题在于，点在网格单元中的位置不确定，考虑相邻网格间可能有较大速度差异，如果点不在单元正中则单元对应的速度不能代表某个具体点的速度。

这个不太好处理，有限体积法里面REA的average就是对每个网格的速度做了平均出一个体心速度。如果任意空间位置都要做插值，这是高精度有限体积法。