CFD中文网

fluent可以做FFT，将监测某点的速度时间曲线导入其中即可

为什么我这个推导的反了？应该是大于0 。

0_1529644803725_微信图片_20180622131807.jpg

2422ba2b-f9b8-402e-936e-8075325aae10-image.png 例如我生成一个part名字叫多“cavity1”，但是导入到fluent里面会多出“cavity1：007，cavity1：008，cavity1：008shadow，”这几个边界fb6ad2d6-0c7d-4266-a7b4-fa752482f0a6-image.png 。相应的网格也被分裂出好多细块。微信图片_20190126110946.jpg 微信图片_20190126110950.jpg

可以在join/intersect里面查找一下看看。

@Noob
也不可行，因为标准 kEpsilon 某些是所谓的 High Re 湍流模型，不能用于求解边界层内的流动。假设按照你说的，y+很小，k 和 nut 都用 lowReWallFunction，并且epsilonWallFunction 开启 lowRe correction，也是错的。你可以画一个平板边界层的算例，通过调整网格的首层厚度，得到 y+不同的网格，然后你用标准 kEpsilon 模型来计算不同网格下的壁面剪切力，就会发现 y+ 小于 30 以后，结果就不对了。
但是kOmegaSST就可以在不同的 y+ 下都得到较好的壁面剪切力结果。

@ljpgbl
这个是全貌
0_1539505430619_Clip_20181014_162251.png

http://advances.sciencemag.org/content/2/6/e1501869

@赵一铭

你好👋

@两月三年 good ，好主意。我试试哈。

short communication类的，找了几个都是做的生物、实验什么的，看的不明不白，打算研究研究文体和思路。谢谢

666eed00-e17c-4be5-8a8c-0ec0e66e1275-image.png 是用哪个UDF函数能输出转速，能给出个例子吗？

大都忽略重力源项，哪位大神给解释下为什么忽略重力源项

重力对流体的作用就是引起水动力学压力，如果是密闭的流体，其不会引起流动变化。因此可以忽略重力项，这样压力边界条件也可以好给出。

一般考虑重力的主要是重力为一种重要的体积力，如密度不均引起的浮力，其中密度不均进一步由温度导致。

还有一种情况要考虑重力是模拟流体的自由落体。

其他情况可能跟具体的应用有关，如旋转系、海洋热风那种普遍都添加重力。

如果不忽略，怎样恰当的处理？

如果考虑重力，通常求解的压力不是先前的压力p，而是压力p减去水力学压力$\rho g h$。从算法层面理解，重力引致下落，压力对这种效果抵消，采用这种压力计算的时候，可以省去这种来回下落-抵消的求解过程。求解出$p-\rho gh$之后，真实的压力即为$p-\rho gh+\rho gh$

参考coalChemistryFoam（OpenFOAM4）即可。
注意在修改cloud的名称，引用，以及在各个方程中添加对应的源项。

需要注意的是引用顺序，一定要在求解flow之前，不然两个cloud见到的Euler场是变化的。
尤其在2-way coupling中非常重要。不然第二个cloud见到的flow是被改变的。

同时，这个方法也有一个问题，也就是两个cloud相互之间是不可见的。
如果考虑4-way的话，或者涉及到考虑整个cell中的theta()的话，这个方法是不正确的。
因为两个cloud是相互独立的，求解的受力也都是独立的。

@东岳 好的，麻烦您啦！我在一点一点仔细研究下。如果能研究明白，万一以后也有人问，我再做个帖子出来♪(^∇^*)

type codedFixedValue; value uniform (0 0 0); // 初始速度 name UM2; // 自定义边界条件名称 codeInclude #{ #include "fvMesh.H" #include "polyPatch.H" #include "fvcMeshPhi.H" #include "addToRunTimeSelectionTable.H" #}; code #{ // 避免重复计算 if (this->updated()) { return; } // 获取网格和边界数据 const fvPatch& p = patch(); const polyPatch& pp = p.patch(); const fvMesh& mesh = internalField().mesh(); // 仅在网格运动时计算（静态网格跳过） if (mesh.moving()) { // 1. 计算上一时刻的面中心坐标 const pointField& oldPoints = mesh.oldPoints(); vectorField oldFc(pp.size()); forAll(oldFc, i) { oldFc[i] = pp[i].centre(oldPoints); // 基于上一时刻节点坐标 } // 2. 计算几何速度（位置变化/时间步长） scalar deltaT = mesh.time().deltaTValue(); vectorField Up = (pp.faceCentres() - oldFc) / deltaT; // 3. 通量修正确保连续性 const volVectorField& U = static_cast<const volVectorField&>(internalField()); scalarField phip = fvc::meshPhi(U, p.index()); // 网格运动产生的通量 vectorField n = p.nf(); // 边界法向量 scalarField magSf = p.magSf(); // 面面积 tmp<scalarField> Un = phip / (magSf + vSmall); // 通量对应的法向速度 // 4. 合成最终速度：几何速度 + 法向修正 // 公式：修正速度 = 原始速度 + 法向*(通量法向速度 - 原始法向分量) vectorField::operator=(Up + n*(Un() - (n & Up))); } // 调用父类方法完成更新 fixedValueFvPatchVectorField::updateCoeffs(); #}; codeOptions #{ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I$(LIB_SRC)/meshTools/lnInclude #}; codeLibs #{ -lfiniteVolume \ -lmeshTools #};

请教一下各位大佬，我想要通过codedfixedvalue实现movingwallvelocity这个边界条件的效果，这个代码是我仿照movingwallvelocity边界条件弄得，我用自己写的这个边界条件算了涡激振动，算出来的结果和movingwallvelocity结果几乎没有差别。计算完成后，我将输出结果导入tecplot，可以正常导入，但是切片时会出现以下报错，目前没有找到问题出在哪儿，应该是我写的这个边界条件tecplot无法识别导致的，但我不知道问题出在哪儿。
641bf902-dad9-4f14-91ab-7f5bc7eda851-QQ20251011-140535.png

@lsprxd

$h=\frac{q}{T_w - T_b}$

其中：$q$为热流密度，$T_w$为壁面温度，$T_b$为参考温度

这样可以吗？

对于这个问题，作为一个重度OpenFOAM v4.1, 轻度OpenFOAM v1606+,和目前轻度但在不远未来的重度OpenFOAM extend-4.1的使用者，我很想去了解这其中的一些渊源。我想从一个最普通用户的角度，来谈谈这个问题。

OpenFOAM v4.x 这个系列的最早的开创者是 Henry Weller. 来源 他在 1989年将FOAM,也就是OpenFOAM的前身创造成OpenFOAM，他被认为是OpenFOAM的创始者之一。目前这个系列的代码由 OpenFOAM Foundation 有限责任公司所推动。网站的页角列出了三名主要的董事会成员:

Henry Weller Chris Greenshields Cristel de Rouvray

OpenFOAM+ 系列的基地则对应这里。 它是由ESI公司的子公司OpenCFD有限责任公司所推动。

OpenFOAM extend 对应的是一个“延展”计划，它的基地在这里 目前的管理员主要有 来源：

Bernhard Gschaider Henrik Rusche Hrvoje Jasak Håkan Nilsson Martin Beaudoin Robert Keser

以上是仅列举了一些主要三大分支的一些基本信息。

============================================================================================

下面我想谈谈基于这些基本信息的个人感受：

个人来讲，我使用OpenFOAM的最主要原因是Licence的问题。Fluent的费用是无法承担的。 可以这么说，如果研究组里能用钱能解决Licence的问题，我恐怕很难会像今天这么依赖OpenFOAM。

Henry Weller对于我来说，一直是一个神秘一般的存在。作为OpenFOAM这样一款优秀的C++库的作者之一，我在互联网上很难找到有关他的个人信息。如果有一天，当面见到他，我恐怕都无法认出他。对Hrvoje Jasak 则是我在这个帖子中才听说的一个人，可能是由于我过分依赖OpenFOAM vx.x 系列。但是，我简单搜索就能找到关于Hrvoje Jasak的个人信息。我不想也不愿意做过多的猜测Henry Weller为什么如此“低调”。

另外，我想谈谈自己对开源精神的理解。

我觉得开源最核心的思想是这样：当一个东西被开源出来，它就不属于与任何一个人，而存在于人们共同的某种“想象”。人们觉得在自己力所能及的范围内，贡献出自己的东西，能够让这个“想象”更好，是推动这个“想象”不断前进的动力。至于推动的动机何在，不同的人有不同的理解。商业运作并没有被排斥在外。维基百科的理念，我个人觉得就是对开源精神的最好实现之一。但是人是复杂的，社会是复杂的。一个人的想法都会随着年龄段的不同而呈现不同想法，更何况千千万万人组成的社会呢？Linux各种各样的发行版，何尝不是这种对开源精神的不同解读而产生的结果。

============================================================================================

adept

好的，我看看，文章不长，图很多，简单点。。。