@如果随风 pointswise可以画
一二
帖子
-
-
@micro 网格数量越多,网格尺度越小,计算肯定更加不稳定,您可以查一下CFL数
-
@micro 我用个G-S举例啊,如果对角不占优,就算不出正确的解,数字会越来越大,你可以试试。这是我5年前做的,当时也没有预处理等措施,BICG跟STABBICG经常算着算着就会发散。
-
@bestucan 我看了目前用到这个
-I.
的地方,包括reactingParcelFoam跟sprayFoam,都有子文件夹,比如sprayFoam下的egineFoam,应该跟compile子求解器有关 -
@bestucan reactingParcelFoam这个求解器的Make文件下有这一行,删除掉并没有什么变化
-
@bestucan 我不理解这个
.
是啥意思,如果是当前文件夹,那么所有程序应该都用到了当前文件夹的一些头文件,为什么只有少部分的Make文件里面有这一行呢? -
各位好,请问sprayFoam下面的make的options文件是什么意思?
EXE_INC = \ -I. \
-
@bestucan 这个混乱是可以解决的,在引用的库文件前面加入路径就行了
-L$(FOAM_USER_LIBBIN) \ -llagrangianIntermediate \
而且我在其他论坛上也看到,lagrangian的代码修改时必须把名字改成一致的,不然会出现很多warning,
-
@bestucan
我做的是droplet,属于颗粒的。组分说的是这个droplet是一种溶液,只是我不清楚这里面的溶液到底算不算组分,溶液里面的水是可以蒸发的,溶质又无法蒸发,溶质跟溶液到底算不算俩种phase,还是按照同一种phase的不同组分来处理 -
@星星星星晴 我对化学这方面不了解,假如我设置的droplet是个溶液,这个是算相同Phase的多组份吗?
-
等人回答我的问题
-
主要做droplet的喷射,不考虑化学反应,但是液滴是化学溶液,测试了一下好像只能做一种组分的
-
我在lagrangian里面添加新的颗粒受力,重新编译了intermediate这个文件夹,发现重新编译的库必须跟原来的名字一样。
比如:LIB = $(FOAM_USER_LIBBIN)/libmyLagrangianIntermediate
这样就不行,必须是
LIB = $(FOAM_USER_LIBBIN)/liblagrangianIntermediate
不然会出现一大堆warning,这些warning并不影响计算,最大问题是这个新添加的受力模型无法读入参数。
-
@l-j刘侃 我调整了Fluent的自带模块,把组分扩散改成dillute扩散,发现也会出现我这种问题。
-
@bestucan 您好!
这个是fluent源代码吗?可以给我也发一份吗? -
您好!
我要算的是燃料电池问题,目前加入了O2,H2的specie transport equation,从下图的进口进入流道,按说这些specie会从流道扩散到多孔介质中,但是扩散不进去,如果初始化把多孔介质中加入了这些specie,计算就会发散
-
UDF燃料电池时,气体几乎不进多孔介质层,有没有前辈帮忙解答下。 -
@东岳 李博,我准备试试,之前写过fluent多面体的网格转换程序,不过这个发不了论文
-
@东岳 假如一个很差的网格,跟一个收敛的结果,导入到Openfoam里面,openfoam会把这个结果算发散吗?
-
用商业软件作为初始值来增强OpenFOAM的收敛性,也能解决UDF自由度不够的问题
-
@Samuel-Tu 可行蓝色可以设置为internal面,顺便统计一下风量
-
@benqing 这是用什么画的呢
-
-
@东岳 您的邮箱多少,我给您发一下
-
@东岳 有Of的网格文件,可以吗
-
@东岳 不是大佬,工作上我还在请教您呢,这个是做肺部药物传输的
-
@东岳
-
@东岳 我有气管的几何
-
@多相流FOAM 热计算直接给定热流密度或者用对流换热算比较合适
-
@东岳 哈哈哈哈,直播时候还推荐了你的网站
-
【6月4日 20:00 B站直播预告】
1、直播主题:超算中的颗粒世界之CFDEM的运用与二次开发
2、直播讲师:王泽坤,北京大学力学系博士,颗粒流动与金属3D打印方向。长期从事OpenFOAM与CFDEM的二次开发。
3、直播大纲:
1)广泛存在的颗粒流动现场
2)颗粒流动基本理论
3)CFDEM介绍与算例实现
4)CFDEM代码解析
群里看到的 -
@Izumi 这个没添加,我回去试试
-
@东岳 我在pimple的option里面添加了这行代码
-L$(FOAM_USER_LIBBIN) \ -lmyPostProcessFunctionObject
编译通过了,但还是运行时侯找不到的我的functionObjcets,这么报错:
in file db/functionObjects/functionObjectList/functionObjectList.C at line 278 Cannot find functionObject file myPostProcess
是不是functionObjects这也要重新编译
-
@东岳 李博,这个具体咋操作啊,我加了这一行代码也没有效果:
libs ("libMyPostProcessFunctionObject.so");
-
我写了一个后处理用的functionObjects,在controlDict添加算的时候会找到这个新的functionObjects,但是用比如pimpleFoam -postProcess -func myPostProcess的时候就会找不到这个myPostProcess
-
@东岳 东岳老师,用LES的话模拟不了这种转捩吗?
-
@纵浪大化 按照kolmogorov的相似分析,在粘性底层应该是各向同性的,所以一般LES是基于这个假设做亚格子模型的
-
@yishanbb2018 具体怎么做,你可以搜一搜CHANNEL FLOW的文献,湍流边界层上对正弦扰动比较敏感,所以需要添加正弦扰动来生成湍流,还有一种方法就是用DNS的数据作为初始场,这样就自动发展成湍流了。反正不对湍流进行任何操作的话,肯定不对。
-
@yishanbb2018 recycle需要加入一个正弦扰动,就是rescale的步骤,不然没法生成湍流
-
你算的是槽道流吧,湍流需要人工进行生成的,一种是进口处就设置人造湍流,一种是recycle&rescale的循环方法,进口边界是LES的难点,不知道你怎么设置进口条件
-
@东岳 这个边界条件已经整合到exteranlHeatFluxTemperature里面了,算的也是温度梯度,也就是$q=-\kappa \frac{\partial T}{\partial n}$,这种边界按照FDS技术手册里面应该是DNS的计算方法。
LES计算时给定的方法是$q=h \Delta T$,$h$是对流换热系数。openfoam里面的边界条件要么是固定值要么是梯度,这种差值不知道在of里面能不能实现 -
我查看了openfoam的代码,里面有externalHeatFluxTemperature能直接设定heat flux,但是我看了代码,对于流体里面直接计算的是{$\lambda T=\frac{q}{\kappa_{eff}}$}:
case fixedHeatFlux: { refGrad() = (q_ + qr)/kappa(Tp); refValue() = 0; valueFraction() = 0; break; }
里面的{$\kappa_{eff}$}是有效导热系数
我参考FDS的技术手册,这种计算方式应该是在DNS里面采用的,LES用的是对流换热的那个公式{$q=\Delta Th$},不知道这种边界条件有没有朋友做过的 -
官方版本我搜了heat flux并没有这个边界条件,在cfd online上有使用fixedGradient来代替heat flux的讨论,但我想要的应该是固体表面对流换热的热流密度,应该是{$q=h\delta T$},有谁做过类似的问题吗
-
在fieldAverage的字典里面加入window这个命令,window后面的参数是5s就是最近5s的平均值。比如你算一个流场,然后写入"window 5",一共算了20s,输出的平均值就是15s~20s的平均值。
-
@东岳 是我搞错了,向量混合积跟微分算子不一样
-
@东岳 是的,这个负号应该是等于号吧
-
@东岳
公式(23)的混合积公式,负号应该是等于号吧
\begin{equation}
\vec{a}\cdot (\vec{b}\times \vec{c})=-\vec{b} \cdot(\vec{a} \times \vec{c})=-\vec{c} \cdot(\vec{a}\times \vec{b})
\end{equation} -
在模拟co2、水蒸气的对流扩散的时候,准备用[Fluent theory guide]{http://www.afs.enea.it/project/neptunius/docs/fluent/html/th/node128.htm}的species transport方程,其中对流项和扩散项都没有什么问题,主要是重力作用的加入,因为co2水蒸气都比空气重,在实际中重力作用很明显的
对于重力项之前别人推荐参考buoyantFoam,我也根据[东岳老师]{http://www.dyfluid.com/buoyantPimpleFoam.html}的这篇讲解,在OpenFOAM里面加入speices transport引起的重力变化。主要还是加入了密度场,通过不同密度对速度方程中的浮力项的作用,来实现重力的影响,并没有直接对co2对流扩散有作用。所以主要还是靠的速度变化,引起co2向下流动的现象。这种变化感觉非常不明显,比如人呼出的co2才占整体重量的4%,实际模拟中呼出co2基本不太向下聚集,Fluent在species transport中是不是也加入了源项?0_1519956843102_scalarTransportBuoyant.zip -
@李东岳 最后还是用了后处理的方法进行计算,又要做实验又要做模拟有点赶,用的是老板之前研究的方法几何平均粒径、几何标准差与浓度的关系,算了浓度再来得到颗粒数关系。不过布朗运动要加进去我感觉得用群体平衡方法。
请问除了ICEM外,还有什么工具可以划分混合网格吗?
关于OpenFOAM矩阵求解器
关于OpenFOAM矩阵求解器
openfoam的wmake的语法
openfoam的wmake的语法
openfoam的wmake的语法
openfoam的wmake的语法
在lagrangian库里面添加新模型时出现的问题
SprayFoam的颗粒可以作为多组分的吗?
SprayFoam的颗粒可以作为多组分的吗?
各位大佬们在干什么
SprayFoam的颗粒可以作为多组分的吗?
在lagrangian库里面添加新模型时出现的问题
Fluent UDF燃料电池时出现气体几乎不进多孔介质层的问题
UDF编程求助,NV_DS是什么函数
Fluent UDF燃料电池时出现气体几乎不进多孔介质层的问题
Fluent UDF燃料电池时出现气体几乎不进多孔介质层的问题
用fluent的结果作为openfoam的初始值能不能增加OpenFOAM的收敛性?
用fluent的结果作为openfoam的初始值能不能增加OpenFOAM的收敛性?
用fluent的结果作为openfoam的初始值能不能增加OpenFOAM的收敛性?
房屋开洞网格
这种图怎么做
类似这种血管几何有人有么
类似这种血管几何有人有么
类似这种血管几何有人有么
类似这种血管几何有人有么
类似这种血管几何有人有么
类似这种血管几何有人有么
buoyantPimpleFoam 求解器,热扩散为0
连CFD都开始搞地摊经济了?
连CFD都开始搞地摊经济了?
自己写的后处理的functionObjects怎么直接调用?
自己写的后处理的functionObjects怎么直接调用?
自己写的后处理的functionObjects怎么直接调用?
自己写的后处理的functionObjects怎么直接调用?
湍流模型计算层流转变为湍流
大涡模拟LES有各项同性假设吗?
请教:基于FLUENT的LES模拟算不出湍流原因
请教:基于FLUENT的LES模拟算不出湍流原因
请教:基于FLUENT的LES模拟算不出湍流原因
OpenFOAM里面有热通量heat flux的边界条件吗?
OpenFOAM里面有热通量heat flux的边界条件吗?
OpenFOAM里面有热通量heat flux的边界条件吗?
fieldAverage怎样对指定时间范围求时均?
东岳流体“张量基础”中的疑问
东岳流体“张量基础”中的疑问
东岳流体“张量基础”中的疑问
species transport方程的问题
pm2.5的计算应该采用什么模型?