CFD中文网

@李东岳 感谢回复。我是想先写个简单的特定算例学习一下算法，格式边界条件啥的都直接写死，然后再慢慢做的通用一点，能搞出来肯定开源。

glgoo，我一直在用的

@vonboett

if I may implement two kinematic clouds in DPMFoam, one that accounts for particle collisions (larger stones but smaller than gridsize) and one that only accounts particle-wall collisions (many gravel particles). This way, I could model erosion..."

Dear Albrecht,

Thanks for your comments. The E-L algorithmn is expensive since it needs to track the lagrangian particles. So it depends on your particles' number. From what you've written, it looks you are on your right way doing this.

D

@匡之舒云 我在做射流燃烧，可以加个联系方式聊一下:xiezuoye:

As an application of this geometric formulation we prove that the solution map of the incompressible Euler equation, associating intial data in Hs(Rn) to the corresponding solution at time t > 0, is nowhere locally uniformly continuous and nowhere differentiable. by Hasan Inci (2013)

你是不是想把曲面展开为一个平面？https://blog.csdn.net/mrbaolong/article/details/106949340

teoria.png
如图所示，这是非常经典的Le Méhauté (1976)波浪理论适用图，但我有点迷糊，主要有两方面疑问：

图中红色部分看起来是Stokes三阶波/四阶波与五阶波的交接处，那么我应该选择用三阶波/四阶波还是五阶波呢？同样地，图中浅蓝色和五阶波也有交集，那么在既是浅蓝色又是五阶波的区域应该用cnoidal还是五阶波呢？ 我还注意到图中五阶波下面还有一个流函数理论，但我看olaFlow的用户手册，用流函数的话需要输入一大堆参数（如下图所示），这些参数应该怎么选取呢？本菜鸡在此请教一下各位:135:
2024-07-08 21-55-00 的屏幕截图.png
2024-07-08 22-01-59 的屏幕截图.png

不清楚商软那面怎么回事，本身我也没太关注这个，我查了一下，OpenFOAM本身这面的streamline看起来是你说的patiline：

This function object generates streamline data by sampling a set of user-specified fields along a particle track, transported by a user-specified velocity field.

后处理软件paraview的streamline是如何定义的还没找到。

https://www.esi-group.com/sites/default/files/resource/other/6395/abstract-roe-university-michigan-upwinding-past-present-future.pdf

We are pleased to announce that Prof. Philip Roe, best known for his work on high-resolution compressible flows based on approximate Riemann solvers and Godunov methods will deliver the keynote at this year’s conference.

请问大家为什么我用meshing划分流体域网格时,明明两个面没有重叠，模型树上还会出现图片上这种接触区域呢？![6392L5}X9$J}C{GTK_ZZ7F.jpg

居然还发公众号了

我模拟的对象是反应堆内的稳态流动传热，由于堆芯内组件数量庞大，精细建模复杂，因此采用多孔介质模型对堆芯区域进行简化。

为了验证这种等效的合理性和误差，我先对单个组件栅元进行多孔介质简化，并将多孔介质模拟结果与未简化的模型结果对比验证多孔介质模型的合理性。

单个组件栅元的截面图如下：
图片1.png

其二维轴对称模型如下图：
f510f5de-0d5e-4389-b577-ad2d3defe29c-image.png

模拟结果
(1)未采用多孔介质模型
我用了二维轴对称模型来模拟，设置固体流体内热源，入口质量流量和出口压力边界条件，流体和固体分别加载体积热源P1和P2（W/m3），流体和固体的导热系数均为常数，得到的稳态温度场和压损分布如下：
bca065f1-7cc4-4385-8113-d7e0acd6d23e-image.png
流体出口温度700℃，固体由于带内热源温度高于流体，固体最高温度为799℃，固体平均温度为756℃，组件压损为2400Pa

(2)采用多孔介质模型模拟
由于固体带内热源，流固存在较大温差，采用非热平衡模型
根据栅元几何结构指定：孔隙率γ=0.1 ，界面换热系数hfs=5000W/M2K（与未采用多孔介质模型算得的换热系数一致），界面面积密度Afs=6.71m-1

通过流速-压力实验关系确定粘性损失项系数和惯性损失项系数，只考虑x方向的流动损失，y方向系数放大100倍

别给多孔介质流体和多孔介质固体设置了内热源，分别为P1×孔隙率和P2×（1-孔隙率），初始化后检查内热源加载总量无误；流体固体的热导率输入值和上个模型一致

得到的模拟结果：
af1efa10-a95f-4e49-a69c-63b4b4357ca0-image.png
采用多孔介质非热平衡模型模拟得到的流体温度是正确的，压损也是正确的，但固体温度分布明显不对，固体域温度热点仅为725℃，平均温度为711℃，均远小于实际模型算得的固体温度值。

我自己的猜测：
多孔介质模型将单组件变为流固均匀介质后，不能模拟固体区域沿导热路径带来的固体温升？刚好我这个栅元模型中固体由于导热产生的温升非常高，如果是这样，说明多孔介质非热平衡模型不适用于我所研究的问题，因为简化后会显著低估固体温度热点。

请问我的理解是正确的么

@李东岳 在 关于LES模拟中采用大时间步长的问题 中说：

你可以放大库朗数，或者你可以用URANS更快。这些都没问题。真正分析的时候库朗数放小就行。并且库朗数一般是局部库朗数的最大值。如果1000个网格里面有1个网格特别小。会出现999个网格的库朗数都在1以下，但这1个网格在10左右。都不是问题。

好的，谢谢李老师！我明白了

我今天看一篇论文，一个简单的方形腔，通过添加金属颗粒来改变腔内相变材料的热性质。文中一个研究点是改变金属颗粒的体积分数，按理说相同的纳米颗粒体积分数下混合流体的热性质都是一样的，根据Gr数计算公式，Gr数不应该变化呀，但是相同的颗粒体积分数下文中又设置了10000与100000两个量级的Gr数，不知道怎么弄的，求大佬解惑。
对了，我找遍那篇文章都没发现方形腔的物理尺寸，不知道这个与Gr的控制有没有关系

@hhh 解决了，需要现在ICEM中设置周期性边界，之后就可以了
我参照的是https://blog.csdn.net/muyangzixue/article/details/108631943

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

https://www.cfd-online.com/Forums/openfoam-meshing/133850-snappyhexmesh-problems.html

格式不对噻