谢谢！！已好

欧拉欧拉？欧拉拉格朗日？欧拉拉格朗日比较好处理。欧拉欧拉的话，目前能想到的就是多孔介质。因为欧拉欧拉压力方程是公用的。速度可以单独设置无穿透。但是压力不能单独设置。 如果有大佬看到相关的sci，感谢分享

你把nu放在括号外试试

这个估计很难了。你这个涉及流固耦合的问题。CFD这边做不鸟。但是你可以这么做，流场获得你需要的压强时程曲线，然后你在固体场中把压强加上去。这样是可以的。

波浪得求解器很多得。一般三个，waves2Foam，olaFlow, ihFoam。我现在可能开始用olaFlow

去源代码里面找了一下caseConstant()对于串行和并行，都可以返回case里面的constant的地址，我试了一下，达到我想要的效果了

@东岳 感謝东岳老師的回覆.首先說個題外話,最近在英國這邊一個公衆號上看到您發的一片勸解博後選擇的文章.因爲是轉發的,起初不知道是您發的,後來越看越像您的口吻,最後在下面看到"轉自李东岳"的時候,心裏很是激動.謝謝东岳老師對後輩們的勸解! 我又做了一個用MPPICInterFoam在沒有粒子的情況下,繞柱(光滑壁面)的模擬.這次用了LES湍流模型.效果很好.並且並行運算時,LES比K-omega SST運行更快.大概是因爲這種情況更符合實際情況吧.跟實驗數據比較,甚至比之前interFoam做的還要好.圖中我做了MPPICInterFoam和interFoam的對比.上面的是MPPICInterFoam的.在柱後水流形成的湍流旋渦非常明顯.在不同高度的對比,效果也很好.至少說明,在水動力這塊兒,這個求解器配合LES,可以得到不錯的效果.當然,邊界條件也很重要. 有問題的大家一起討論啊.

删除算例目录下除了m开头的所有文件，注意备份！！注意会把constant文件夹内的也会删除 find . ! -name 'm*' -exec rm -f {} \;

in file db/functionObjects/functionObjectList/functionObjectList.C at line 278 看起来你改了functionObjectList.C的代码？

@Samuel-Tu 非常感谢您的回复。讲的非常明白，我之前没有想到这些。其实，我想实现的就是一个热源点慢慢加热最终提高出口温度，或许我应该试试改成热源强度再试一试看。非常感谢您。

感谢分享！！！收录在此 https://www.cfd-china.com/topic/3499/openfoam小代码

gmsh可以，这个是开源的网格软件。

@xpqiu 感谢，我去查一下!不过目前我用break也能顺利终止程序，并行和串行都可以

没有实际做过，但是这个链接提供的方法看起来很靠谱。https://www.physixfan.com/如何计算湍流能谱/

@东岳 非常感谢东岳老师！按照您的指导，我设置了新的边界条件immersedboundaryEnthalpyFvPatchField，其他部分的设置目前都正常，但是想要使用浸没温度边界上的函数来更新焓值时会报错，我的程序部分如下所示： const label patchi = patch().index(); wordList type = thermo.T().boundaryField().types(); if(isA<immersedBoundaryFvPatchScalarField>(type[patchi])) { immersedBoundaryFvPatchScalarField& TIb = thermo.T().boundaryField()[patchi]; scalarField & Tref = TIb.refValue(); } 之后程序就会报错， error: invalid initialization of reference of type ‘Foam::immersedBoundaryFvPatchScalarField& {aka Foam::immersedBoundaryFvPatchField<double>&}’ from expression of type ‘const Foam::fvPatchField<double>’ thermo.T().boundaryField()[patchi]; 应该是无法将基类转换为派生类，想和老师请教一下遇到这个问题要怎么解决呢，麻烦老师了！！！

@东岳 微信回复您了，不然还是在微信说吧

fvm::ddt(rho,T) = $\frac{\p \rho T}{\p t}$, rho*fvm::ddt(T) = $\rho\frac{\p T}{\p t}$

问题背景是在simpleFoam计算中添加被动标量运输，因此在controlDict里添加了function的scalarTransport。但是在设置参数的时候发现resetOnStartUp。 在源代码中的discription是Reset scalar to zero on start-up。 现在比较疑惑的是，如果打开resetOnStartUp: １．是在每个时间增量计算开始时把标量reset为０ ２．还是一整个计算结束后，再次重新计算时把标量设为０ 大家有经验吗

经过一上午的排查，找到原因了。在我使用并行的时候，往scalarField里面加数据的时候超过了定义的大小

有现成的，可以看看6DOF模块的代码，里面又控制方程，也有写好的阻力

@Sloan 感谢回复，overset是另一个路子，我这里没有涉及

@东岳 这位是我后来一步一步检查，用checkMesh得到得具体的问题 如果感觉没有那么复杂的不规则边界，那么基本就是几何拓扑关系出了问题吧

pEqn里面用到了rAUc

感谢分享！

用并行计算的时候，想通过某个face的local index得到其global index，查了一下似乎是这个函数faceProcAddressing，但是写入solver的时候调用不了这个函数。说没有定义此函数，如何调用这个函数呢。。ｗiki上说globalIndex=mag（faceProcAddressing[localIndex]）-1

@东岳 东岳老师，decomposed后的patch里面的face是否还是按分解前的顺序排列呢。例如分解前的patch中有6个face编号为0，1，2，3，4，5。分解的时候刚好从中间分开了，分解成两个processor里面，则processor1和processor2里的face编号都变成了 0，1，2。假如processor2里面的face是未分解时候的4，5，6，那么现在face 0，1，2一定对应的是4，5，6吗，也就是说这个顺序会不会改变，会不会0，1，2对应成原来的5，6，4了

@东岳 感谢!！!！!！

@东岳 踩坑踩多了而已

朋友，均方根脉动速度和雷诺应力提取解决的怎么样了，我也存在同样的问题，希望予以解答，感谢了哈！

看起来很复杂的样子....

@东岳 好的，谢谢老师～～

@veen

@东岳 那截面云图总是有小的变动，这个影响结果分析吗？谢谢老师～

@xulei-zju 直接看trap不行吗？

@lin_lia 可以用吧，在DPM模型里面设置

@史浩 嗯嗯，好的，很详细了，谢谢您的回复！我看看mpi里的参数设置，琢磨下看看

@东岳 这个是出现问题前的最后一步，网格变形太厉害了，我在考虑是不是可以用圆尾的好一些

#include "udf.h" #include "dpm.h" DEFINE_DPM_SCALAR_UPDATE(ceshi, c, t, initialize, p) { Injection *I; Injection *Ilist = Get_dpm_injections(); loop(I, Ilist) { Particle *p_new; p_new = new_particle(I, FALSE); p_new->part_id = get_next_part_id(); P_INJECTION(p_new) = I; p_new->I_id = 0; ((p_new)->stream_index) = 1; P_POS(p_new)[0] = 0; P_POS(p_new)[1] = 0; P_POS(p_new)[2] = 5; P_VEL(p_new)[0] = 0; P_VEL(p_new)[1] = 0; P_VEL(p_new)[2] = 0; P_DIAM(p_new) = 0.001; p_new->n_steps = 0; p_new->next_time_step = solver_par.flow_time_step; P_TIME(p_new) = solver_par.flow_time; memcpy((char *)&(p_new->init_state), (char *)&(p_new->state), sizeof(particle_state_t)); P_FLOW_RATE(p_new) = 0.001; p_new->time_of_birth = 0.0; alloc_pvars(p_new); append_particle_to_list(p_new, &(I->p), I);/* 问题函数 */ I->n_particles++; } } append_particle_to_list(p_new, &(I->p), I); 这个函数总是出问题，有人知道怎么用么这个函数

@东岳 感谢李老师回复！圆柱周围的网格细化我明白，但是在整个水槽的自由液面处网格细化也是用snappyHexMesh做的吗，而且这里的自由液面区域和圆柱周围区域的衔接也很好，看上去像是一个整体加密区域。但是我总感觉自由液面处的加密像是通过blockMesh实现的，我在这里还是有点困惑，李老师您能详细的指点我一下吗。

@东岳 而且我也不要把其他场映射过去，只需要颗粒速度场

弱耦合是可以实现的，如果你做的是刚体运动，就用DEFINE_CG_MOTION，如果你做的是非刚体的变形（比如弯曲变形之类的），你就要用DEFINE_GRID_MOTION，你自己定义的结构求解模块放入这两个宏内，就会随着时间推进每一步求解一次。从FLUENT中读取节点、面气动力的方法手册里都可以查到。

@东岳 @Hunahum ，我在cfd-online上看到这样一段代码，想尝试着写出来，不过总是报错，请问你知道问题在哪儿吗

@东岳 谢谢东岳老师的关注，第一个报错应该是关于内存释放的，有时候还会报错free() invalid next size (normal).尝试了把EEqn.H能量方程中的relaxfactor项删除，就可以避免这个报错，但是会对计算结果有影响…… 还尝试过用gdb运行程序，就不会碰到这个报错，能量方程可以顺利求解（至少第一个时间步可以），但还是会提示下面那个网格数量不匹配的错误

@东岳 @xiezhuoyu 谢谢！

@东岳 谢谢！

@东岳 谢谢老师，我已经找到网格生成规律了，从左至右，由下至上地规律依次生词。满足误之前需求工作，reNumberMesh有需要我再去琢磨吧，再谢！

@东岳 ，好的，万分感谢

@西极若木 邮箱：mechanics_xzy@163.com 微信：xzy11235813 多交流多交流 共同进步

setFields里alpha.water=1的区域里就是water，性质是给你说的1项的。OpenFOAM基于FVM，一个cell里相分数alpha=1表示都是水，alpha=0表示空气。供参考