我翻了一下tutorial中的案例，只找到了如下的一段controlDic中的function，是否是对其进行一定修改来达成目的呢？求各位大佬指教
1721204958853.png

@疏影横斜水清浅 您好，最近我也在学习这个帖子的问题，我想在已经o切的圆柱侧面通过patchinject注射颗粒，所以我要定义这个注射的边界面，您现在知道如何去用topoSet和createPatch去实现在圆柱侧面上画边界了吗

@没有时间出去玩 在 求助，更改求解器时的问题 中说：

但是发现现在的求解器不支持这个边界条件

是报错给你说没有这个边界条件吗？？那可能是版本原因，你这个版本没有这边界条件（我猜的），比如我实验室版本很老，没有noSlip条件，我们就用fixedValue（0 0 0）

信息量太大不好debug，我怀疑是不是正负号写反了。对单一网格离散和对整个计算域离散每个面的面矢量方向不一样。

我用系统自带的网格（ 1 75 1）可以算出来，但是我改成（30 2000 1）就算崩溃了，打算算高2m宽15.4mm的圆管截面

@李东岳 网格是用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:

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