GCC只要内核没问题 一般Linux下Fluent可以UDF直接编译 不需要额外的环境变量配置 否则你需要查看 对应的Fluent makefile文件
刚开始可以使用Ubuntu这样的桌面Linux系统先熟悉一下 之后再使用Journal file 因为一般性的UDF还是需要不少调试环节的 有一个图形显示可能会好点 个人不成熟的意见 仅当参考

@一二 我的意思是这主要是多孔界面的问题 对于扩散方程来说 基本都是基于AD的 这不是产生这个的原因 你要模仿的是一个多孔对于切向流的组分吸收的工况 我之前带过学生也搞过这个 实际上最后他还是放弃了侧向了

这已经是比较成熟的方法了 应该aixprocess在15年前就已经成型了 做木材洗涤工艺的除尘优化 如果说无脑的办法 还是直接写碰撞列表做加减法 也可以直接设置某一方质量半径都为0 然后另外一方做简单的质量累加 但是实际上真实的情况比这复杂得多 液固的碰撞的结局以及最后的补集情况是多样性的 这有很多相关文献 但是对于UDF编写来说 因为是不同injection 还是有一定难度的

fluent采用的是质量流量来做颗粒的 所以他的外形不是真实外形 对于这个 我建议还是用EDEM这类的试试吧 至于接触点 在fluent内部只是空间解析几何的问题 对于球形可以 你可以直接获取P_POS 然后根据直径 算出碰撞边界就行 接触算法我手写过 然后最后某985的在pt上还发表了两篇 还算不错 虽然当然我认为需要做particle_in_cell循环 他觉得太复杂就放弃了

显然这又是射流的一个严肃的问题 射流的边界条件实际上是非常困难的 有条件的话 最好走的路线是检测喷口处的压力变化曲线 但是实际上一般来说 国内能够检测的还真的不多 退一步就是根据出口压力的型线 简化成为升压-稳压-降压区间 这样通过理论可以计算出 大致的每一个时刻（简化成线性）的出口质量流量 但是因为可能压力比较小的时候你还没有choke 那么不同的区间你还得应用不同的方法 至于流量是恒定的值 你直接给质量流量入口就行了 但是基本上这跟实际情况还是差得比较远 因为你的储气罐 你的阀门 你的气力管道的变化 都会影响你的流量

这是目前fluent自带多孔的一个比较严肃的问题 实际上fluent自带的egrun式多孔还和真实的多孔相差非常远 因为实际中 这两层之间有一层称为过渡层的地带 在过渡层的里面 多孔介质的表面实际上会对组分有一定的吸收能力 但是fluent自带的多孔模型只是一个压力-动量源项而已 并不具备这样的功能 因此在这个设计里面 流体总是会走阻力最低的路线 自然你的多孔介质也就没有什么作用 如果你的多孔理论背景比较浅 我非常不建议你在目前做这块的工作

fluent应该是没有曲面这个说法的 因为你是有限体积 所以实际上你只有cell face这样的网格数据 那么获取一个网格面上的数据 基本上你先获取这个面的法向量 采用 F_AREA 然后通过获取直角坐标系的数据 C_U C_V C_W 就可以通过解析几何的办法 算出面法向上的速度分量 这种对于一些多孔介质的表面对比是有用的 虽然目前这块还是一个理论不完善的地带

射流 假设你是超音速流体 那么这是一个非常复杂的问题 这个问题是直接和初始条件相关的 对于瞬态 就是每一步的上部计算结尾的状态 因为它是对于初始条件非常敏感的 所以一个微小的积累可以导致很大不同的后续结果 实际上就是一句话 超音速这块实际上还是很难控制的 最好还是辅以实验 一般优化效果需要达到15%以上 才能说是设计真的有优势 如果要做严肃精确研究 fluent目前对于这种还不是太适合 但是fluent的优势是可以综合多场一起做出一个工程上可以大致量化的结论

vof模型只有一个混合相的温度 所以你这个方式其实难以实现 现在新出来了一个欧拉-欧拉-vof 可以尝试一下 实际上更好的办法是直接采用欧拉-欧拉同样算一遍 这样你可以方便的获取出口温度 实际上这样的方法也不会有太大的误差

除开物理上的不守恒性的考虑 还有计算时的不收敛考虑 你可以通过fixed values的方式来操作 这种方式对于某些燃烧反应器还是可以用用的

网格移动问题是非常复杂的问题 最简单的还是直接重新做网格 生成网格的时间和主频比较相关 可以买台高主频的计算机 也有通过网格节点移动宏来针对网格移动的UDF 比如做仿生机械鱼游动的仿真 但是这些东西你还是需要利用算法来通过局部的改变来进一步扩大优化临近的网格质量 如果变动比较小 或者是你的局部网格真的很粗糙 可能还是简单的 否则直接通过几个时间步变形到定位区域 光是这个理论还有每步的计算 可能比起重新做更没有效率 其实实际上 如果可能的话 是学术研究 你可以用ICEM做全结构化网格 然后再block里面变形调整 这个重新生成的速度是相当快的