东岳,记录美好算法
walten
帖子
-
-
@JamesCheng
括号的是变量, tp是个指针,应该指向某一相, 1 应该是这一相的第二个组分,从0开始计数。 -
@东岳 我这就让他去试试 @IO-Turbulent
-
@东岳 哈哈哈 没事他已经找我代他发帖了。
-
大家好,我需要植入EMMS模型到multiEulerFoam中,其中Hd的系数需要针对液相体积分数的不同而有所改变,因此我需要写入一个分段函数到Scalar.H文件中,函数是这样写的
inline Scalar between(const Scalar s1, const Scalar s2) { return (s1>0 && s2<=0) ?1:0; }
但是我在编译的时候调用between函数会出现数据类型不对的错误
我想知道应该如何定义才能定义出符合需求的数据类型。
另外,我想详细了解openfoam中src/OpenFOAM/primitives中的数据的正确定义、调用的方式,请问有什么资料可以学习,谢谢 -
李老师辛苦了,手动点赞!
-
差不多情况。。。泪目
-
昨天因故错过 东岳脱口秀第二季,有没有哪位大侠录屏了的,球球分享下,谢谢了!
邮箱:walten0503@126.com -
@zousiyu
一段程序可以处理所有的数据对吧,微调就可以了,改明儿试试 -
绘图matlab代码好评,已收藏,感谢楼主的干货分享。这种类型的图,excel是不是也能做出来,万年excel选手路过。。。。。。
-
今天的我,悄悄地来,宛如昨天的我。我挥一挥衣袖,对明天的我说一声,你丫快点的!!!
-
当年,,人人网还叫校内网的时候。。。。我也特别喜欢这个图,我还用这个图做过封面。。做过头像。。
-
-
@东岳
不不不,您才是大佬,我是仰慕小菜鸡 -
@东岳
人在上海,鞭长莫及,不过fluent装起来也简单,无脑skip,最后LM安排一下就好了,实在不行我微信给你传个安装教程视频。ezpez -
@东岳
让我们一起,为服务器窒息 -
@东岳
老贾是贾跃亭么,终于回国了? -
@西湖冷月
nice兄弟,d=====( ̄▽ ̄*)b -
@bestucan 开的是双精度的,之前用pbm都是自己写UDF的,这次用的自带的,第一次遇到这个问题。。。回头有空搞搞,这感觉是软件自身的问题,我装个低版本的c试试,现在装的是 2012,改明儿装个2010,我还是感觉2010比较稳定。
-
@Samuel-Tu
学校的ip才能用么?出了学校能用么? -
@bestucan
试了试,并不是udf的问题,而是fluent自带pbm的问题,我试了14.0和19.2两个版本,合并模型用luo的,破碎模型选择luo和 lehr的都会出现该问题。合并和破碎都用常数或者合并用湍流模型,不会出现该提示,其他模型暂时没有试,所以我感觉就是fluent在调用luo和lehr的合并破碎模型的时候出的问题,究竟为啥,咱也不知道。。也没法知道。。。 -
@闻久STU
乐观点, 六年 这样才学得深入 -
@bestucan
我也这么想过,看来也只能如此逐行破解了,谢谢老师,我回头试试,不过只是在启动软件首次开始算的时候会出现,暂停计算,然后继续,就不会提示了。 -
支持,也是为了预防ansys被卡脖子
-
@bestucan
之前从来没遇到这个问题,这个怎么解决呢老师?我需要修改c还是fluent启动项?而且我UDF里并没有用到int型,,用的是real。。。 -
前排打卡,老板大气!
-
@walten
没有老师遇到过么。。这是c的问题还是udf的问题啊。。。 -
给位老师同学,我最近有个case一计算就提示我数据类型的问题,但是计算并没有发散。我用的版本是14.0,计算用了欧拉双流体,开了传质和反应,开了PBM模型,加载了UDF。百度了,也查了fluent手册,都没有说明这个提示是什么意思。是c语言的提示么?有没有同学老师遇到过类似的问题?这是为什么?
-
@关伯兰在笑 逐行排查吧,所有分母都添加一个很小的量。
-
横坐标小r是指画图点的位置,R是反应器半径,比值是指该点位于反应器的径向位置,比如,0 就是值反应器中心,-1,就是左壁面,1 就是右边壁面,能迅速定位到径向位置。
-
@杰克逊 只要你知道颗粒团聚和破碎的机理,可以写UDF啊。fluent里自带的合并和破碎机理一般是气液系统的。
-
@东岳 我现在接触的还很浅。。以后碰到这方面问题再来请教~~
-
@东岳 就是还是那个问题,书写习惯中log是不是一般是指以10为底的对数。如果这样的话,您文章里应该写成ln,避免引起误会。
-
@walten 东岳老师,我试出来了。还是用的您给的公式,
只不过这里的αd不是指床层的固含率,而是指 所有填充颗粒占整个反应器的体积分数,这样的话文献中的初始固含率0.6经过换算应该为0.033。然后代入N公式,最后计算出来的值与文献基本一致。
感谢东岳老师百忙之中的指点,我终于解开了这个谜团了。。。
另:我又发现了您分享文献里的几个笔误。
1:公式94下边第一句话,standard 写成了dtandard。
2:公式96分母右边括号多了一半。 -
@东岳 东岳老师,C语言中的语法是 log 是以e为底,log10 是以10为底。。但是平时书写log是以10为底,以e为底是不是习惯写成ln。所以看到您文章里的log我就想当然的用了10为底的log,excel中log也是默认以10为底。。所以会有以上的计算误差。现在剩下的问题就是N的计算了。。
-
@东岳 东岳老师,我把计算excel表中A3和A4 两个数值的公式中的log改成了ln,得到的结果比较一致了,但是N仍然不对。您看下。
-
@东岳 是的,0阶矩是颗粒总数量,1阶是颗粒总长度,m1/m0应该是平均粒径,所以我计算时是不是哪里出问题了?我也还在排查,而且N的计算值也不对。
-
@东岳
文献中的一些物性参数。
-
@东岳 颗粒平均直径 d10 我用的米的单位,即颗粒粒径是235微米。
-
@东岳
东岳老师您好,您分享的文章中N值按理说就是文献中矩阵标准化之前的M0值,即单位体积的粒子数量,我按照您文章中的算法计算了N 与文献中M0差距了5个数量级。感觉M0还是没有算对。然后对计算值和文献值的矩进行标准化后,也还是存在差别。标黄的部分是B列的计算公式,您看下那里出问题了。。。感激不尽 -
@东岳 东岳老师,您17年AIChE文章里的公式和您刚分享的文章里的公式不太一样,不知道哪个是笔误。
-
@东岳 感谢感谢!这个问题困扰我好久了,就是不知道初始矩是怎么算出来的。。
-
@东岳 谢谢东岳老师的耐心指教,我现在的问题就是不知道4个m的初值是怎么得到的。。。。卡在这里了。。。
-
@东岳 了解了一点点,但是东岳老师,我现在的问题就是每单位体积的粒子数量这个数是根据填料量自己算出来的吗? 还有就是,我还是算不出来文献里的那些数字,,您贴的图是您自己算的吗? 东岳老师能再指点一点么,,我感觉自己快从泥淖里出来了。。。。
-
太难了,,感觉自己智商不够搞矩方法。。。。
-
@东岳 东岳老师 ,我试了,,算不出来。。。难道是我方法不对么。
按照您的方法,m0是1,但是文章里 m0/m0才是1,那么m0=4.85e9是怎么来的呢?根据填料量算出来的吗?但是这个文献参考文章的填料量是10kg的颗粒。但是m0如果是颗粒数量,乘以颗粒平均质量的话,得到的是76kg。。。这就是我疑惑所在了。而且公式里 d10是知道的,平静粒径。d32并不知道。是要迭代几次吗?
-
没有老师解答么。。。
-
@lcj4ever 感谢提醒,我看到了。试试看吧,不知道跟PD算法有什么不同~
-
以前用Fluent做过颗粒增长方面的模拟,用到过QMOM法。有一个问题就一直困扰着我,不知哪位老师给指点指点。
颗粒增长模拟采用的是气固流化床,床内初始填充了颗粒,模拟床内颗粒的增长过程。可参见文献:LIU, Si-Si; XIAO, Wen-De. CFD–PBM coupled simulation of silicon CVD growth in a fluidized bed reactor: effect of silane pyrolysis kinetic models. Chemical Engineering Science, 2015, 127: 84-94.文献中图3给出了六个初始矩的值,然后根据我对手册和文献的阅读,六个矩的初值是根据PD算法求得。但是我的问题是,在知道颗粒粒径、密度和填料量后怎么计算六个初始矩呢? 按照Fluent手册介绍,0阶矩代表总数量,那我是否可以用填料总质量除以单个颗粒平均质量得到0阶矩的数值。但是我按照刚才文献中的数据反推,并得不到相等的结果。。还望各位老师解惑一下,QMOM方法中矩的初值是怎么算来的,需要用到哪些数据?谢谢大家!~~
五一了,送李老师一份礼物,请笑纳
fluent DPM模型UDF使用求助
OpenFORM植入EMMS模型时函数变量不匹配的问题
C币运行机制测试
OpenFORM植入EMMS模型时函数变量不匹配的问题
OpenFOAM国内下载地址,满速
昨天东岳老师办的talk show 有没有哪位大侠有录像,求分享下
昨天东岳老师办的talk show 有没有哪位大侠有录像,求分享下
Daniele Marchisio的CFD Talk登记
k-omega 模型的网格无关性验证
k-omega 模型的网格无关性验证
C币运行机制测试
一个照片
C币运行机制测试
有人周四9月17号有空不,来我们这帮我们装个fluent?
有人周四9月17号有空不,来我们这帮我们装个fluent?
有人周四9月17号有空不,来我们这帮我们装个fluent?
有人周四9月17号有空不,来我们这帮我们装个fluent?
分享Matlab读取Fluent数据代码
fluent非错误提示
ansys 2020学生版
fluent非错误提示
如何CFD自学成才
fluent非错误提示
如何CFD自学成才
fluent非错误提示
CFD巨著开始印刷!
fluent非错误提示
fluent非错误提示
求助UDF并行计算的问题
数据归一化
有关Fluent里面的PBM求解
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
多相流与矩方法的数学模型
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
关于QMOM矩初值的给定
多相流与矩方法的数学模型
关于QMOM矩初值的给定