看很多文献中都是将计算域的长、宽、高和二分之一的高度$\delta$做比值等于多少倍的$\pi$,这其中的原理是什么?
yt
帖子
-
-
@fubianhanshu
在particleIO.C下添加一行代码,注释原先的输出可以输出颗粒全局坐标。os << position() << token::SPACE << celli_;
先采用最笨的方法,把保存时间步设置小,遍历读取坐标,其他的再慢慢实现,看着lagrangian粒子程序太复杂了
-
@星星星星晴 谢谢哈,这个昨天我也看到了,我再看看他的这个实现思路。另外应该还有两种思路,一是:src/functioObjects/lagrangian/cloudInfo可以在计算过程总输出颗粒质量随时间的文本信息(这正是我想实现的效果),应该也可以将particle的坐标输出的程序代码到这里新建类,但是对于底层代码理解的不透,尝试了还没编译通。二是:至少v1706版本以后的functionObject/field/中都有particleDistribution库函数,也可以将它改写到of版本里实现,这个工作应该更大些,涉及太多的依赖头文件和关键字写法的不一样。
-
DPMFoam中将是不是可以实现颗粒位置信息按照时间顺序输出到postProcessing中的一个文本中?正常情况颗粒位置是保存在每个输出的时间步文件夹下有positions文本中的,我计算颗粒比较少,想要输出颗粒连续时间的位置,若是是每个时间步保存文件大小太大,且后期通过shell脚本遍历每个结果路径下的positions内容写入单独文本中比较麻烦。是不可以改写src/lagrangian/basic/IOPosition中的远程序实现,不知道具体怎么改,有前辈可以指点下吗?
-
@李东岳 这个实现了就能很好的解决关于颗粒的后处理问题了。一种是要想知道某个颗粒归属哪个网格就一直追踪这个颗粒的坐标,然后再确定网格;另一种是要知道某个网格内有哪几个粒子就像你说的遍历一遍网格确定粒子,会很慢。但可能还会遇到下面的问题,我在做CFD-DEM的时候在paraview中选取的粒子发现有两个编号,一个是:ID,一个是id,而且在多核计算的时候发现id会在跨越block时进行重新编号(后处理显示时),DPMFoam中应该不存在。
-
@sungda 好的,谢谢哈。我使用的5.X版本,跟liggghts耦合用,想要给定随距离变化的速度梯度作为初始场,我安装试下,有问题了向你多请教
-
@sungda 你好,请问你是在openfoam哪个版本实现的呀
-
@李梦瑶 粒子可以循环,你试下就知道了。
-
目前是用的MPPICFOAM
-
@星星星星晴 我都是合成后处理的,选定某个粒子,整个时间序列ID号应该是不变的,多核分割后怎么处理界面的也不太清楚,我注意到分割后并行计算时,每个块内的粒子是动态变化的
-
@upc_ngh 选完后显示全域,装上Python后,用find找也很方便,就像楼下道友说的
-
@bestucan 你是行家呀,厉害厉害,终于知道原因了
-
要是虚拟机的话,装一个版本临时看看结果,处理下,后期处理我感觉转成VTK格式在Windows版本上流畅很多,速度也快
-
@fubianhanshu 上面的图没输出好,贴个清晰点的
-
@fubianhanshu 自己先回复下吧,在paraview中可以处理,我先贴出来一个图,某个固定ID粒子的Z坐标随时间的变化曲线图,具体的还在摸索,不是很清楚,希望得到大家更好的建议
-
@zhangxc0223 我不清楚你怎么做的,不过我刚发现可以在paraview中可以显示粒子ID,随着时间步推移,好像每个粒子的ID号是不变的,这样就可以通过后处理输出某个固定ID的各个时间步的信息了,可以交流下
-
楼主你这个问题解决了吗?
-
欧拉拉格朗日方法的MPPICFoam求解器计算的颗粒输运,该怎么输出或处理粒子的位置信息,比如:给每个粒子编号,追踪某个固定编号(id=100)的粒子在每个时间步的空间位置、速度?纠结纠结······ 跪求懂得前辈指点
-
-
请问目前CFDEM最高支持到哪个版本?
-
@胖小帅 速度是不是要设置成slip?你试试
-
@液固两相流 我是在计算域内初始化一定量的颗粒,输出计算过程中每个颗粒的位置和速度信息,这个可以实现吗?
-
大佬,最近我也在做固液的模拟,输出的positons文件中的内容怎么看,您知道吗?
-
文献有点大,上传不了,麻烦东岳老师下载吧(1)http://dx.doi.org/10.1016/j.ijsrc.2017.09.006(2) 10.1002/2013JF002911文献里都是用DEM方法模拟的,绝大多数文献里MPPIC都是用于气固或气液,想着试试用它做液固,跟别人的对比看看如何,增加些模型什么的,目前刚开始接触欧拉拉格朗日法,没什么思路
-
@东岳 东岳老师,意思就是MPPICInterFoam现在没办法吗?我感觉计算水沙MPPICInterFoam比MPPIC好点,想转MFix了
-
在计算颗粒在连续相中运动时,使用MPPICFoam和DPMFoam对U和P文件中的进口和出口使用“cyclic”边界,颗粒能够在出口和进口周期性运动。但是,在使用MPPICInterFoam时,对U和P文件中的进口和出口使用“cyclic”边界,颗粒从出口逸出后不再从进口流入,无法进行周期性运动,有知道的大佬这可能是什么原因吗?
-
@东岳 我在速度文件力设置好内部速度为定值(0.6),压力文件里设置内部压强(1e5)通过一段时间计算后,流速越来越大。像我需要指定初始速速,只需要设置速度文件,压力文件中设置为0就可以,或者将指定的初始速度转化为相应的压力添加在压力文件中,速度文件设置为0呢!
-
@东岳 李老师,内部流速也不需要设置吗?若是这样,那流动的初始流速是怎么来确定大小值的呢?
-
@bestucan 好的,谢谢
-
想做两相混合的DNS模拟,应该怎么下手
-
@东岳 在 twoLiquidMixingFoam中物性参数中“Dab”、“alphatab”分别代表什么意思? 中说:
考虑过driftFluxFoam么?加什么模型?
加两相间的扩散模型,具体还在学习琢磨中。我刚才看了下driftFluxFoam感觉跟我的更贴切,可能都不用改方程了,我先试试看,现在算了ke湍流,下一步算les,弄完后可以话后面我把例子分享出来
-
@东岳 在 twoLiquidMixingFoam中物性参数中“Dab”、“alphatab”分别代表什么意思? 中说:
你在模拟什么?这个求解器太简单了,我周末就能弄完。也没啥需要太大的改动,就是符号需要改一下
打算找点算例,更新个非tutorials的算例
感觉这个还行:Stratified flows and associated shear instabilities modelling over an inclined plan
谢谢大佬,下去看下提供的资料,再更新的时候可以把施密特数的定义加上(湍流黏性与湍流扩散系数的比值,为无量纲参数),层流和湍流扩散系数其实是有量纲的,下意识会觉得系数应该是无量纲参数,今天在这里迷惑了一会,希望其他人不要走老路
-
异重流,水沙两相流动,目前算了几个工况,发现nu和alphatab对两相混合影响很大,这两个参数敏感度很高。后面要把动量方程改一下,加点模型,不知道能不能完成,个人感觉很慌
-
自己又看了看,结合东岳流体网上的twoLiquidMixingFoam解析,Dab应该表示层流扩散系数,与Dc对应;源程序里定义的alphatab是施密特数的倒数,alphatab*湍流粘度vt,为湍流扩散系数Dt,其量纲与湍流粘度相同,与东岳流体代码解析中相对应。不会上图片,感觉应该说清楚了,希望对其他人有帮助吧
-
根据后面量纲看,猜测Dab代表混合项的粘度,alphatab代表两相混合的程度是这样子吗?自己在改变alphatab的数值后(自己加了ke湍流后,把案例里设置的1过渡到0),两相的混合程度确实会降低,相当于削弱湍流强度似的,向层流过度。感谢有了解的前辈能够解答下
[图片]
湍流计算域物理尺寸一般怎么考虑?
输出particle坐标按时间排序到指定文本下
输出particle坐标按时间排序到指定文本下
输出particle坐标按时间排序到指定文本下
粒子与网格归属问题
初始速度场给定
初始速度场给定
关于MPPICInterFoam的周期性边界的疑问?
关于MPPICInterFoam的周期性边界的疑问?
MPPICFoam中如何输出某个颗粒在各个时间步的位置信息?
MPPICFoam中如何输出某个颗粒在各个时间步的位置信息?
多版本OF共用paraview
多版本OF共用paraview
MPPICFoam中如何输出某个颗粒在各个时间步的位置信息?
MPPICFoam中如何输出某个颗粒在各个时间步的位置信息?
颗粒的全局ID输出问题?
颗粒的全局ID输出问题?
MPPICFoam中如何输出某个颗粒在各个时间步的位置信息?
CFDEM安装
CFDEM安装
求教对于一个四周固壁的开敞水面边界条件设定问题,谢谢~(仅仅只有水,不涉及空气)
关于颗粒属性的输出
关于颗粒属性的输出
关于MPPICInterFoam的周期性边界的疑问?
关于MPPICInterFoam的周期性边界的疑问?
关于MPPICInterFoam的周期性边界的疑问?
请教一个周期性边界的设置问题
请教一个周期性边界的设置问题
twoLiquidMixingFoam计算时,设置成laminar,就是DNS模拟吗?
twoLiquidMixingFoam计算时,设置成laminar,就是DNS模拟吗?
twoLiquidMixingFoam中物性参数中“Dab”、“alphatab”分别代表什么意思?
twoLiquidMixingFoam中物性参数中“Dab”、“alphatab”分别代表什么意思?
twoLiquidMixingFoam中物性参数中“Dab”、“alphatab”分别代表什么意思?
twoLiquidMixingFoam中物性参数中“Dab”、“alphatab”分别代表什么意思?
twoLiquidMixingFoam中物性参数中“Dab”、“alphatab”分别代表什么意思?