libs ("libmyfvOptions.so"); 这句不能像一些官方提供的案例那样放在fvoptions文件里,得放在controldict文件下
麦迪文
帖子
-
-
用fvoption框架想外加一个阻力,扩展了src/fvoption/文件夹fv::buoyancyForce这个类,类里增加了两个volScalarField类型的私有成员(67行68行)
下面是构造函数初始化
成员函数里把这俩volScalarField类型相乘
结果串行没问题,但是并行却却报错,把调用弹出来后是并行相乘两个volScalarField的变量调用multiply函数,查这个函数又是个宏,所以就很神奇了
所以为什么串行可以并行却出错呢,按理说我以前经历的并行出错大多集中一些边界条件上数据传输上,求各位大神提点建议 -
Ansys Maxwell 3D如何输出洛伦兹力,并且将洛伦兹力导入Fluent中实现磁场和流场的耦合?如图中所示,Maxwell计算出四分之一圆角矩形块的洛伦兹力密度,如何输出四分之一圆角矩形块的mag文件?
-
@wying 邀请码在哪里发出呀
-
@aaaachenjian 模拟是在微重力环境下,上下壁面有个温差,在自由面上加个热毛细力,已经放弃用VOF方法了,像这种变形较小的直接挪网格,这是算法天生问题不太好解决
-
@东岳 稀相流夹杂物分布的想听一听fjgwd001@163.com
-
@dzw05 谢谢
-
@dzw05 您好,想问一下凝固后的固相 密度 和 粘度是怎么考虑的?因为看文章的方程计算域全部采用的vof,对于常见的气液都是流动状态很容易理解,但是,凝固相是不再运动的,还用vof的话感觉很奇怪,所以猜测是不是操作了密度或粘度使固相区域难以流动
-
不很精通也不常用fluent,个人尝试使用感觉fluent看似让cfd入门门槛降低,减少学习时间,但是对于技术壁垒高的学科缺少对基础的扎实掌握,遇到问题无法对问题有清晰和精确的定位...............谢谢各位老师回复
-
问题已解决 体积分数换QUICK格式,体积分数残差标准由2e-4降1e-5,单步迭代次数20提40
-
大佬真好
-
fluent的dpm倒是并行的,楼主可以44,颗粒多了也没慢下来太多 -
正在试,体积分数从一介迎风格式换成了QUICK格式
-
跟这个情况类似,我觉得说的有道理 -
没想到引来大佬关注,而自己并未深入了解过pbm,就提出这样初级的问题,感觉特别惭愧
-
几何模型:尺寸50cm×20cm,入口0.4cm,物理模型:欧拉两相流,子模型默认,湍流模型标准ke除多相湍流模型用dispersed其他默认,pbm使用Discrete,聚并破碎全用luo模型,最小直径0.1cm,最大直径2cm。边界条件入口全气速度1m/s湍流强度5%直径4mm,bin3为1其他为0,压力出口,气体回流100%,bin1其他0,初始值液相高度44cm,所有bin0,上方气相bin1其他0,初始时刻液相略微上涨,随后剧烈下降,未做pbm方向不了解这种原因的细节,但是了解一些欧拉两相流模型,未查出原因,感觉像是数值上的原因,猜测用qmom可能会改善?,附qmom边界条件设置可参考D. Li﹡, others, Comparison of Eulerian QBMM and classical Eulerian-Eulerian method for the simulation of poly-disperse bubbly flows, AIChE Journal
-
是The Discrete Method,其他的方法不太明白速度入口该怎么设置
-
就是简单的体验体验新版本fluent的Population Balance Model.。。。。。。 。。。之前检查了一圈不知道出的啥问题,关了pbm就不会下降
-
@麦迪文 纯属娱乐用丷丷
-
开始时候是多了一点点,可能是因为气相冲出来了,后来液相却是剧烈下降。。。。。。后来靠udf每隔一段时间步统计液相体积,小于初始值就在底面喷液相大于了就不喷
-
我担心算到最后成了单相流
-
#include "udf.h" #include "dpm.h" DEFINE_DPM_SCALAR_UPDATE(ceshi, c, t, initialize, p) { Injection *I; Injection *Ilist = Get_dpm_injections(); loop(I, Ilist) { Particle *p_new; p_new = new_particle(I, FALSE); p_new->part_id = get_next_part_id(); P_INJECTION(p_new) = I; p_new->I_id = 0; ((p_new)->stream_index) = 1; P_POS(p_new)[0] = 0; P_POS(p_new)[1] = 0; P_POS(p_new)[2] = 5; P_VEL(p_new)[0] = 0; P_VEL(p_new)[1] = 0; P_VEL(p_new)[2] = 0; P_DIAM(p_new) = 0.001; p_new->n_steps = 0; p_new->next_time_step = solver_par.flow_time_step; P_TIME(p_new) = solver_par.flow_time; memcpy((char *)&(p_new->init_state), (char *)&(p_new->state), sizeof(particle_state_t)); P_FLOW_RATE(p_new) = 0.001; p_new->time_of_birth = 0.0; alloc_pvars(p_new); append_particle_to_list(p_new, &(I->p), I);/* 问题函数 */ I->n_particles++; } }
append_particle_to_list(p_new, &(I->p), I); 这个函数总是出问题,有人知道怎么用么这个函数
-
文献里没说考没考虑,不过fluent里可以加入,有针对颗粒的相关湍流分散力模型
-
主要参考这篇文献
Cho S M , Thomas B G , Kim S H . Bubble Behavior and Size Distributions in Stopper-Rod Nozzle and Mold during Continuous Casting of Steel Slabs[J]. Isij International, 2018, 58(8):1443-1452.
现在想在DEFINE_DPM_INJECTION_INIT(init_bubbles,I) 这个宏里创建和删除颗粒实现颗粒聚并破碎效果,但是在injection里新增一个颗粒 的宏没在标准的帮助里,在dpm_*.h的文件中试了些个都不好用,想咨询下专家的建议
-
我要刷火箭
-
这个边界条件是在一个单独的文件目录下
-
@yuan_neu 是李宝宽老师的学生吗?
-
最近看了篇晶体融化凝固方面的文章
Numerical Simulation of Thermo-Solutal Marangoni Convection in the Floating-Zone under Microgravity Fields但是在OF里没找到相关的求解器,heatTransfer目录下的求解器都不合适,multiphase目录下也没找到,请各位大神指点一二
-
这个在三维就还复杂些,在二维,实际简单得多,热毛细力就是在帮助里的温度梯度*(方向扭转90度的相分数梯度)再乘一个系数,这些在帮助里都能查到
-
@东岳 下面这个实验的图效果好点
下面是fluent设置接触角 下壁面设置的110度,上壁面设置70度时数据
工作站有点问题,暂时用不了openfoam,过两天贴下of加接触角的结果
-
我这个没设接触角,不过在openfoam的tutorials
OpenFOAM-5.x/tutorials/multiphase/interFoam/laminar/capillaryRise/
这个case里,已经设置了接触角,然后将液相入口设为壁面条件,入口压力也设为fixedFluxPressure,结果如下
自由面附近出现了较大的速度
-
(热毛细对流反向加在自由面后,自由面上速度被扭转回来)
-
@Albert-Lee 我想到的一种解决办法是界面变形不大的情况下或许可以用动网格的方法解决
-
是是是,我也因为这个没啥解决办法,好在我导师不怎么管我,自己定方向,研究了一阵子没突破出来就整稀相流去了
-
我之前用的VOF法,只能加在计算域里,如果用普通方法,就加在边界上
-
首先,拿到温度梯度和界面梯度,两个变量做一点数学计算就得到沿界面的热毛细力了,一旦我发文章了会把代码贴文章后面的,数学公式也很简单,UDF不是很麻烦哈
-
之前的数据没找到,又重新画网格算了一下,也有这个问题。。。。但是,fluent把热毛细力的UDF插入进去后,虚假速度影响会减弱,自由面上速度方向会变化,但openfoam用fvoption方法加入热毛细力,没法抑制虚假速度(以上两图是没加热毛细力,纯自然状态的形状)
晚上会把fluent加热毛细力使自由面的速度反向图附上
-
是按照实验来的,周围是空气,上下是俩圆柱
-
回复上上个话题,同一个网格点的不同传输速度不会处理所以没有处理。。。。。。在网格内的颗粒速度认为一致。。。。。请问泰勒涡的这个不同算法对比的case的文章发表了吗?能不能发个这个文章的连接学习一下
-
附:
在OF5中 interfoam计算曲率在:
OpenFOAM-5.x/applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C中的414行-----530行 -
不好意思,回复晚了,没有什么进展,已经用了fluent解决此问题, 如李老师所说,应该是openfoam曲率计算不准导致,当在openfoam把重力设置为0时,表面不再变形,虚假速度也相继大幅度降低基本消失,虽然猜测是这个问题导致,但是由于数学基础不好,所以没有什么改进的措施
上传的case是在openfoam5_x版本,直接运行interFOAM即可(由于研究的是热毛细对流,其特性是使自由面上出现沿自由面的切向速度,刚好与虚假速度方向相同或相反,从而使计算结果误差巨大)
-
-
还有就是我觉得颗粒作为杂质也可能会有一些聚并与破碎现象,cfd-PBM改造后可能可以应用在这一领域上,但是需要一些实验数据支撑,实验台我们组已经搭好,就是实验着实难做,而且对openfoam-7中的pbm模型也是一头雾水
-
@麦迪文
1.就是在连续相的动量方程中去掉曳力项 ---------这个是指的是twophaseeularfoam的中进行的操作 -
这个稳态simpleFoam不可能出来,必须是瞬态的,瞬态流场耦合dpm,单向耦合,不考虑颗粒对流场的作用,就能出来三角形的PAS
水平太低,不能理解 “同一个网格点不同的传输速度” 指的是啥,我现在的简化思路有两个
1.就是在连续相的动量方程中去掉曳力项
2.就是用pimlefoam ,外挂两个方程
固相动量方程
固相分数也就是想用欧拉数值模化拉格朗日模型
感觉上可行,还没具体开始操作。。。。这两个想法在fluent里实现特别难(得用巨量的udf),也然而,在openfoam里实现却没有那么的麻烦,这也是openfoam比fluent的魅力所在啊 -
是的是这个,我自己写了个单相流无形变的液桥的求解器,只不过对壁面边界进行了一特殊点处理,结果跟实验也差距不是特别的大,然后用欧拉-拉格朗日模型也能做出这个效果,的确现在正思考用twophaseeularfoam简化成稀疏流求解器,做出欧拉-拉格朗日模型的效果
做自由界面追踪是想整点创新看能不能发个文章充数好毕业哈哈哈哈哈哈,结果是个坑。。。。
-
但是fluent计算VOF模型的计算速度远远赶不上相同配置的OPENFOAM,所以我一有时间还是在思考openfoam的问题出在哪,感觉上是压力的问题但又不知具体该怎么入手
-
惭愧,没有,用fluent做了上述case的二维对比试验,也有一些奇怪的速度,但是非常非常的小,几本不影响流场内温度的分布,openfoam在边界处的奇怪速度分布影响了它的正常温度传递。。。
附上一张维也纳工业大学该领域作者的研究的惊艳成果(在上述模型中加入小颗粒)
引自:Muldoon F H, Kuhlmann H C. Origin of particle accumulation structures in liquid bridges: Particle–boundary-interactions versus inertia[J]. Physics of Fluids, 2016, 28(7): 073305.
-
@Dingcy 哈哈哈哈,我就瞎看看,我主要做CFDEM方面
-
@Dingcy interfoam告一段落了,接下来跟你们学习点先进的,整点夹杂物搁进去了
在自编类里两个volScalarField变量相乘,串行可运行,并行报错
在自编类里两个volScalarField变量相乘,串行可运行,并行报错
Maxwell与fluent耦合求解
注册打卡
interFoam求解器的液滴的奇怪现象(虚假速度?)
Daniele Marchisio的CFD Talk登记
请问OF里有没有融化和凝固相关的求解器
请问OF里有没有融化和凝固相关的求解器
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
首个OpenFOAM的VS Code插件公开发布
DPM并行计算
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
fluent pbe模型 气液两相流里 液相越来越少什么情况...........
用ddpm模型 模拟气泡聚并与破碎
用ddpm模型 模拟气泡聚并与破碎
用ddpm模型 模拟气泡聚并与破碎
ThermoMPPICFoam
chtMultiRegionFoam共轭传热边界的换热原理
一个CFDer的日常
请问OF里有没有融化和凝固相关的求解器
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
求助:液池热毛细对流算例,已编译好求解器和算例,但就是和文献自编程序不一样
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
interFoam求解器的液滴的奇怪现象(虚假速度?)
twoPhaseEulerFoam 的升力项怎么不在 动量方程里
twoPhaseEulerFoam 的升力项怎么不在 动量方程里