CFD中文网

@Linzd 多孔介质区域在toposet的时候本身只是选中了cell集合，只有在壁面上才能通过原始代码计算法向力和切向力。

http://www.mafengwo.cn/gonglve/ziyouxing/99502.html

毕超的《计算流体力学有限元方法及其编程详解》这本书里面详细介绍了方程各项的离散，不过使用MATLAB编写的，我现在主要就是参考这本书的

@liujm 对的，通常我们计算整个风机都是采用komegaSST湍流模型，毕竟这个湍流模型基本上是做风机的都会用的，很少用其他湍流模型。

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

大家好，估计你们都意识到广告条改了。嗯，因为几个月来，只收到几个3位数的广告，不足以支撑。索性自己全用了。

同行想做广告的可以撤了，没戏。我用心血做网站做这么几年，做到月几十万电机量，我不缺打广告那一顿饭钱，有的一顿饭钱都不想给

同行想做项目的可以试试，单打独斗年7位数没问题，前提是你混到CFD行业的塔尖。我个人每年只能做数量有限的项目，多了我也干不过来太占用时间

同行想做培训的也可以试试，请有点良心别瞎忽悠，好多机构只知道赚钱，就跟互联网公司跟风卖买肉夹馍

同行想做咨询的，不好说，这玩意不费事躺着就来钱。但做咨询的不是特别多。小客户这面不见得你竞争得过我。大客户那面我也拿不到。如果你还没开始做，那估计你也够呛

:papa:

感谢东岳老师，我尝试着改小了，问题得到了部分解决。

http://www.sciencedirect.com/science/article/pii/S0021999117306289

0_1505210410009_捕获.JPG

嘿嘿了解了一些，出口设置的是你估计的就好，它这个好像只在初始化的给定一个初始值而已，小白一个不知道其他大佬有什么看法:xiezuoye: :xiezuoye:

reFineMesh可以加密一个方向，这个需要指定区域，你试试topoSet+refineMesh呢》

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

把ubuntu装在了机械硬盘。

用了大约4年的固态硬盘，忽然发现难以忍受机械硬盘。用了1星期之后，又装在了固态硬盘。只有40G大小。

现在固态硬盘便宜了，250G的640元。

等径管里如果出现激波那就直接堵塞了吧

首先感谢李老师的热情帮助~~
这个PPT中的第9页给出公式（1.5.9）。如果这个公式是正确的，则我的疑问就没有了，因为我原图中（4）式的右端项，第二项与第三项抵消了。但可疑的，这个公式（1.5.9）似乎并不正确。贴上来再麻烦李老师帮我看看：
0_1502119393399_7fd74579-5b34-4c15-ac3c-1cc3aae56236-image.png

@coolhhh 是空的计算域，几何部分已经挖掉了

@chon01 NUMERICAL SIMULATIONS ON FLOW AROUND AN INCLINED CIRCULAR CYLINDER AT HIGH REYNOLDS NUMBER

替代文字

@xpqiu 感谢感谢

如标题，感谢大佬们