CFD中文网

各位老师，近期在研究立管涡激振动问题，想采用切片法，主要思路是针对具有大长细比特征的立管，我想沿展长方向设置多个切片，通过求解每个切片上二维圆柱的升力和阻力来获得立管不同高度处的外力，然后再将这些外力作为输入项，自己编一段有限元程序来求解立管结构动力学问题。
现在存在的主要问题和疑惑是：我如果想设置100个二维切片，应该采用什么样的方式实现循环求解这100个二维圆柱涡激振动问题，以及如何将这100个切片上求解得到的外力与有限元程序耦合实现流固耦合求解呢？
想请问各位老师有什么思路吗？

@yuky 可以参考这个文章https://www.bilibili.com/read/cv6292980/

@yuky 是的，用画网格软件画的，这个网格最初是来源于NASA官网

@李东岳 对，我们也没想到影响会这么大，有时候这类原因都想不到

@李东岳 是的李老师，之后就是用您调的结果算的

@luomuliunian 谢谢大佬

不知道是不是可以在这个版块问

导师投稿后，我收到一个合作者关系确认的邮件，请问一下是否需要我也注册一下，然后确认才能继续审稿流程啊？搞不清楚了，各位是怎么处理的呢？
“Elsevier asks all authors to verify their co-authorship by confirming agreement to publish this article if it is accepted for publication.

Please read the following statement and confirm your agreement by clicking on this link”

@King_RVM 我解决了我的问题，是乱用了tmp这个内存管理工具，删掉就编译通过了

@Cr 非常非常感谢你的分析和回复，我的编译通过了，编译通过之后运行算例有新的问题，原因是tmp这个内存管理小工具用得不好，我最新的通过编译的部分修改如下，供大家参考

@悬浮的猪 我想请问你解决了你的问题么？我是修改了komegaSST湍流模型，编译后运行自己的算例，结果报错，跟你的错误很类似

@Cr 非常感谢这部分资料的提供，我没想过源代码会有问题

@Cr 我尝试过改成一样的类型，所有都改成vol##Type##Field，但是这样sij会报错

在这里 dev(symm(tgradU()()))是张量除以一个标量还是张量，还是报错。
dev(symm(tgradU()()))的灵感来自于李老师的东岳流体网站http://www.dyfluid.com/docs/tensor.html

@Cr 我发现问题来自于单位矩阵的相乘，当我用更简洁的形式去表达sij来避免单位对角矩阵I的表达，sij没有报错，但是V3关于I 的部分报错了。其实我想添加的表达式是

其中sij为除以omega_来做到无量纲化，而图中omegeij为
除以omega_来做到无因次化
我现在修改后的定义为

V3中的I不可避免，因为需要将那一项从标量转为张量，现在的报错如下

@Cr 是的，我定义I1 I2的时候调用了转置函数，但是我现在发现实际上我不需要，于是删去了转置计算。关于量纲的问题，你说的很对，虽然sij omegeij、标量I1 I2和张量V1 V2 V3没有量纲，但是定义过程中S 与 Omega有量纲，为了减少量纲的问题，我定义如下

sij 与 omegaij的无量纲化是通过除以omega_实现的。现在的报错是与定义的张量I 有关，在sij中，场离散的项symm(tgradU())不能直接与单位对角矩阵I相加减，该怎样把I转换成一样的类型呢？我想通过volTensorField定义I，但是报错。我定义单位对角阵是因为定义sij中1/3*partU/partx需要乘以狄拉克函数，实际计算中就是单位对角阵。
目前报错如下:

@Cr 谢谢你的回复！我把跟sij的量都修改成了：：Internal的形式，这几个量编译倒是没有报错，但是出现了新的错误，提示是openfoam里面源代码的错误，如下

这是怎么回事呢？DimensionedField是定义单位的，我定义的标量矢量I1 I2 V1 V2 V3都是无量纲的量，最后的量纲是乘以密度与湍动能决定的

@bestucan 谢谢跟我讨论！divU的定义是标量

S跟I 都是张量，我觉得应该可以相加减啊