@李东岳 网格是用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就是对每个网格的速度做了平均出一个体心速度。如果任意空间位置都要做插值，这是高精度有限体积法。

@李东岳 现在社会语境下的“道友”几乎约等于“同志”:chigua2:

是的，Boussinesq是没有什么问题的，计算波形并没有衰减，而且计算结构物的波浪力也是基本准确的。 那我还想请教一下您，在波流耦合这方面有所了解吗，我知道的olaFlow造波流，但是我测试过了，有严重问题，比如逆向流作用下波浪高度会增加，这是基本的规律，但是olaFlow造出来的却不呈现这样的结果。然后我想在openfoam2112版本中开发波流耦合模型，但是目前感觉可以修改边界条件，那求解器和其它东西是否需要修改你，有思路的大佬们可否指点一下

QQ截图20240712161602.png

QQ截图20240716123826.png 参考了一篇文献，插值边界单元的主要作用是从其他网格接受流程信息，但我的算例中interpolated为0，请问是哪里出了问题呢？

@秋山结月 您好，请问您解决这个问题了吗，我目前也遇到了QQ截图20240712161602.png 应该怎么改呢

跟接触角、网格、时间步长都有关系。你可以尝试把Co缩放到最大是0.1试试

直接自己编程写。把每一个网格单元的面网格的值复制到另一套网格里面。

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

:146: :146: :146:

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

@李东岳 已经实现！感谢东岳老师回复

:146: :146: :146: :146:

@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:

目前只能判断粒子在哪个网格，然后输出某个网格的流场数据。稚嫩这样搞。

@江山不如凉竹 您好，你解决了吗，可以指点一下吗

@lwjetmann 好的，我试试

谢谢大佬，我再去修改lisa模型参数试试看，别的模型都没有出现问题

@郑学习 感谢老师！我有相关进展一定与您分享讨论

@派大星huangxiaojie 你好，windows下打开paraview解决了么，我打开contordict文件，但是有速度场压力场无法识别