CFD中文网

@一二 我的意思是这主要是多孔界面的问题 对于扩散方程来说 基本都是基于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里面变形调整 这个重新生成的速度是相当快的

连续方程 如果你是在模拟一个工程实际设备 这种不能达到-4的情况是比较常见 通常在-2左右 这代表在你的设备设计中间 真正的运行的时候 设备并不能保证一个进气和出气均衡的稳态（比如内部有一部分涡旋致使一部分进气永远在这里堆积） 而其他方面的原因上述几位已经讲解过了 对于工程模拟来说 连续性不达标的问题是比较常见的 你可以看看你的设备内部是不是存在涡旋流 或者是因为本身你的边界条件的设定是与实际不符合的 比如你在一个明显不均匀的出口面取了一个平均压力出口值 这样有的时候你需要扩充出口 直到真正符合工况的外界条件

TP_N 输出的是在一个颗粒簇中间（一般dpm都是用质量流量来算的）的真实颗粒数量 但是这些东西在17.0版本之前好用一点

物理速度和这不相关 物理速度描述多孔的统计了孔隙率的速度 但是实际上除非你是固定孔道的多孔板式 才能够很好的预计内部速度 与理论计算速度可能比较符合 因为中间真实速度也不好测量 我在某吊桥线缆上有过项目实测经验这种方式预测得还不错 而关于浓度连续性 你的设想是在一个连续的流道内他的浓度是连续变化的 那么这个连续流道实际上已经属于需要对多孔介质直接建立模型 这种方式你可以具备连续性 但是fluent处理是经验公式性的 不是真正求解通道的 所以实际上你的浓度当在一个cell zone求解时 当进入下一个多孔介质孔隙率跳跃变化的cell zone时 他直接会利用这个cell zone的porosity求解来得出这个cell zone cell内部的值 它不需要保证连续性 实际上fluent不太可能实现预测准确分布 但是如果你的结构很大并且分布比较均匀的时候 fluent能够大概预测一下趋势

CFD models comparative study on nanofluids subcooled flow boiling in a vertical pipe
你看看这篇文章 然后你看看他对于纳米部分的经验化就了解了

如果说从长远的技术来讲 我更看好实际颗粒推积的真实内部场的格子法模拟 虽然这个方面目前尺度有限 但是如果之后计算机算力飞速上升的话 那么这块可以做到部分可控 同时包括目前的扫描技术实现内部多孔建模的技术也日趋成熟 德国的Math2公司这块也已经推向应用化阶段 对于你的这个研究课题你首先应该考虑这种方式 如果你作为一个研究课题的话 fluent这块几乎多孔内部分布没有任何前途

其实真正严格研究这个问题 目前还不可解 你采用fluent可以不用管这些 因为fluent 如果你选的默认superficial选项 求解多孔介质的时候 不考虑孔隙率 只有在非稳态项才考虑 那么实际流场跟孔隙率就无关 然后切向的速度用适当的该方向的阻力系数来模化 一般可以近似成为无滑移的条件 因为你用fluent那么还是精度和稳定的折中了 这种办法是最佳收敛手段 表观速度这个词最早是针对颗粒层床的经验数据 就是流量除以过滤面积 也是不考虑孔隙率的
关于这方面真做研究目前存在介面连续和不连续两种模型 但是我试过几乎都难收敛 实际上真实情况介面浓度更接近不连续 你如果有兴趣真的做深这部分我给你提供参考 但是我个人建议不要深入搞这块 严格来说 fluent除开压降 其他的都是错的 但是他能算整体设备 其他方式复杂几何目前都不能收敛flowporous.doc

换热根据经验公式得出体积分数相关的纳米流体导入率影响 再通过求解particle in cell 计算体积分数 再编UDF吧

这种属于相变模型 实际上fluent内置已经做得比较主流了 但是如果你需要更高精度的模型 还是得编UDF

这个模型正在做 这个UDF实际上类似树叶系数沙尘的模拟 一般来说 具备相间的mass传质 然后mixture相的动量和湍流因子 能量等的模化
哪一项有变化就加哪一项 实际上你可以写在同一个UDF内部 根据不同条件来判据这种物性变化就行 采用phase指针来切换

fluent使用的是交错网格加重叠网格技术