@李东岳 谢谢李老师夸奖,我现在课题改为离心泵等旋转机械了,关于这方面以后还要向李老师以及中文网的各位大能学习!
悬浮的猪
帖子
-
-
哈哈哈哈好了,总算善始善终
-
@李东岳 抱歉啊李老师,失踪了一段时间,主要是前段时间在谋未来,刚刚有了还算满意的着落,我现在可以把硕士写的三篇小论文都发上来,希望还来得及
-
有需要相关代码的朋友可以私聊我,我会酌情提供
-
关于本专题讨论已发表了一篇文章,在知网可以查询到,小论文是基于低雷诺数κ-ε模型冲击射流传热的数值模拟,大论文的第五章有更详细的讨论,大论文题目是三种典型贴壁射流的数值模拟,过段时间应该也可以检索到。有错误之处还请大家指出,欢迎讨论谢谢大家,希望能帮助有需要的伙伴。
-
@李东岳 东岳老师,关于这个主题帖子,4月份会有一篇相关用OpenFoam完成计算的论文刊出,有兴趣的小伙伴可以参考着看一下,里面加入了一个直接求解热通量的双方程标量方程,也许传热相关的小伙伴会感兴趣,如果有小伙伴需要,我可以选择性提供模型代码,以及调试好的计算算例。谢谢东岳老师,因为马上毕业,不想自己弄出来的一点东西就没有了,还希望广大网友完善,谢谢!
-
@King_RVM 没错的问题就在乱用了tmp指针,我已经解决了。关于tmp指针,我在网上搜索到过问题,祝好!
-
@东岳 李老师我后来也用LS模型,计算了一下,对于湍动能的设置可能李老师还需要改进一下,我计算的结果要比这个稍好。这个结果我用原来设置的参数计算得到过,但我觉得不太完美,就修改了湍动能和湍流耗散的设置得到了更高的结果
-
@东岳 李老师我取的是雷诺数10400下的射流工况,w取0.01423为特征长度
-
@东岳 李老师我这个求解器没什么特别,只是在simpleFoam中加了一个温度场,相关的文献我给老师发过来
-
-
@东岳 李老师,怎么样提供算例呢,我非常乐意提供的,只是前段时间有事没有来得及上网,对此很抱歉
-
最近写了一个基于simple算法的磁流体求解器(假设感应磁场影响不存在,用电势泊松方程求解),并且写了一个基于低雷诺数ke系列的磁流体湍流模型,在运行方管算例的时候,单核运算没有问题,但考虑到更大的网格,想尝试多核并行计算,结果报出长串的错误,下面是正常运行时的迭代图以及残差图。针对这样的问题,我不清楚是求解器的问题还是模型的问题,命令以及分核文件应该都没有问题,希望有路过的老师能指点一下,小子感激不尽。
-
@东岳 东岳老师,openFoam最新的版本里面simpleFoam代码是不是就是采用的simplec算法,多出的rauUt似乎就是相比较simple算法多的东西
-
@linhan-ge 哈哈,老师你好,这个只不过是我个人的一点误解,当不得真的
-
@东岳 文章正在给导师以及相关的老师看,目前我在研究相关磁流体领域类似的传热现象,期望能再研究生生涯的末期再写一篇文章,恐怕到时候又是一连串问题要请教东岳老师
-
@东岳 抱歉东岳老师最近才看到回复您,我整理了一下数据内容,疫情期间写了一篇关于冲击射流传热的文章,主要是探讨不同湍流热扩散系数求解方案对预测结果的影响,最近正在修改文章
-
@东岳 李老师,那换句话说,LamBremhorstKE没有这一项,是否代码有问题呢
-
谢谢李老师
-
最近由于疫情关在农村,想把一些由标准k-E模型衍生出来的低雷诺数模型写进openFoam里面往后可以运用,现在openFoam中有两个以标准k-E模型为基准的低雷诺数模型,分别是LaunderSharmaKE以及LamBremhorstKE模型。在文献上来看,这两个低雷诺数模型上的区别只是在于阻尼函数以及源项上的不同,但我发现在epsilon耗散和k湍流方程中LaunderSharmaKE多出了这两项,标准K-E也是(加粗倾斜项)
- fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_) - fvm::SuSp((2.0/3.0)*alpha()*rho()*divU, k_)
然而LamBremhorstKE模型则没有这两项,我翻阅了陶老师的数值传热学也并没有这两项的说明,希望大神解惑。
-
这也说明我往里面写的模型大体正确
-
李老师这是我修改后的结果
-
李老师给出的建议一针见血,刚刚看到李老师的建议我就去进行修改了,我把改正后的对比图放上来,细节决定成败,李老师今天的给的建议让我感觉很多细节上我都没有到位还需要和老师们多多学习
-
@东岳 李老师必须是严格的一层网格嘛?那我之前是有问题的,我做的是2d数值试验
-
@cccrrryyy 谢谢老师的建议,我准备用都有的Laundershmark做一下试验
-
@东岳 李老师这样做是不正确的吗?of进行2d计算我的师兄师姐们是这样教我的,将不用的那个面empty。小白问题还请李老师指点!见笑了老师
-
@东岳 好的谢谢李老师,我尽快将我的算例以及模型贴上来,网格的数量在十几万左右,比较小是一个二维的网格
-
@cccrrryyy 老师您好!我准备在原有的rans低雷诺数模型上写一个自己稍加改动的模型,加上一些自认为算是有所改善的尺度或者说是一些阻尼函数,那么我的工作意义也就是在于对这些新兴的尺度或者说阻尼函数在特定的模型上的适用性,那么第一步就是现在of上还原这个模型,并且能够与一些权威的软件上的湍流模型相验证,然后再在其上做出调整,做一些数值试验。我现在的这个算例算出来的k与epsilon的分布差异还是蛮大的,所以现在我认为这与我自己的某些细节设置存在关系,最后,谢谢老师抽出时间回答。
-
@dzw05 谢谢您的建议,我的湍流模型是拿来验证的,基本按照其出处文献的说明来写,我会再看看比较一下的,谢谢您的建议。
-
@cccrrryyy 你好,请问您验证的是什么样的算例呢?我现在反思了一下或许我这边的问题还要出在边界条件上,亦或是温度方程上,湍流模型基本上与fluent保持一致。
-
-
同样的网格,也近乎相同的边界条件,最后都达到了相应的参差需求,然而对比实验数据,of的结果惨不忍睹,连大致的趋势都没有正确,相反fluent却得到了更令人满意的结果。至于湍流模型是我自己添加的,尽量仿照文献公式编译,与fluent中自带的湍流模型保持一致,具体来说式低雷诺数k-e系列的,AB模型,研究的现象是冲击射流的流向努塞尔数变化,话不多说下面式对比图
-
仅仅运行一步直接报错,这是具体的错误提示:
Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0667125, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0804578, No Iterations 3 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0500194, No Iterations 4 time step continuity errors : sum local = 0.879188, global = -0.272208, cumulative = -0.272208 smoothSolver: Solving for T, Initial residual = 1, Final residual = 0.0597181, No Iterations 2 smoothSolver: Solving for epsilon, Initial residual = 0.418304, Final residual = 0.0304297, No Iterations 2 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.00011942, No Iterations 1 smoothSolver: Solving for epsilont, Initial residual = 1, Final residual = 0.0345744, No Iterations 2 bounding epsilont, min: -1800.64 max: 30630.5 average: 1357.93 --> FOAM FATAL ERROR: tmp<N4Foam14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE> deallocated From function const T& Foam::tmp<T>::operator()() const [with T = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>] in file /opt/openfoam6/src/OpenFOAM/lnInclude/tmpI.H at line 278.
-
根据各位老师的建议这次我将代码直接贴上来,补充说明问题,有老师说我的错误提示应该在tmp的使用不当上面,但是我现在对具体出错的地方还是不甚明白。我之前对于tmp的理解停留在一个智能的指针分配管理者上,知道它能减少内存的调配,提高效率,但实际所知并不多,对它的使用也停留在仿照源代码湍流模型上的类似公式。这次的问题我直接贴上来,承接上次帖子贴出来的问题:
tmp<fvScalarMatrix> epstEqn ( fvm::ddt(epsilont_) + fvm::div(phi_, epsilont_) - fvm::laplacian(alphaEff(), epsilont_) == Cp1_*Gt*epsilont_/kt_ + Cp2_*G*epsilont_/k_ - Cd1_*(sqr(epsilont_)/kt_) - Cd2_*((epsilont_*epsilon_)/k_) ); epstEqn.ref().relax(); epstEqn.ref().boundaryManipulate(epsilont_.boundaryFieldRef()); solve(epstEqn); bound(epsilont_, epsilontMin_); tmp<fvScalarMatrix> ktEqn ( fvm::ddt(kt_) + fvm::div(phi_, kt_) - fvm::laplacian(alphaEff(), kt_) == - Gt//速度产生项 - fvm::Sp(2*epsilont_/kt_, kt_) ); ktEqn.ref().relax(); solve(ktEqn); bound(kt_, ktMin_); correctalphat();
-
@bestucan 谢谢老师指点,这个问题我已经解决,老师有空看一下,我运行这个湍流模型时出现的错误,万分感谢老师
-
@白礼耕 老师您好,老师能具体说一下用info检查错误的方法嘛?谢谢老师抽出时间回答我的问题,谢谢老师
-
@wwzhao 恩恩谢谢老师抽出时间回答我,老师我想问一下这个具体应该怎么去解决呢,希望老师能提供一个关于tmp的解决思路,万分感谢老师
-
@星星星星晴 老师您好,因为我写的是新的两个标量方程,类似于求解脉动速度的k与epsilon,尝试的时候我先采用的就是类似于k和epsilon的边界条件,用的也是k与epsilon的壁面函数,我想先算一下,有必要的时候我再添加自己的壁面函数。我的计算是用的simpleFoam,所以时间步用的是一秒一步,slover我编译了一个我自认为相匹配的求解器,读取了一下由kt和epsilont求解更新的湍流热扩散系数。数学公式的话我是采用的文献上的写的代码,我觉得老师说的除以0以及负数开根的情况很有可能,我再仔细阅读一下文献。
-
@星星星星晴 恩恩,谢谢老哥
-
@bestucan 恩恩,老哥你的说法应该是对的,谢谢,但是现在又有了一个新的问题,你不妨看一下,万分感谢
-
-
谢谢老哥,但你可能误会了我的意思,我在湍流模型中就需要用到这个固有参数,难不成在里面写一个这样的函数吗
-
@东岳 只能是厚颜无耻地再问问李老师了,万分感谢
-
还有一个待解决的问题就是怎么在湍流模型中读取算例文件中constant文件夹下的transportProperties文件,如果有老师看见并指点一下的话,真是不胜感激
-
是否需要用到指针?新编译的湍流模型需要用到求解器中的层流普朗特数参数,而相应的求解器中的温度方程也需要用到湍流模型新解出的湍流热扩散系数
之前在simpleFoam中看到的参照是turbulence->divDevReff(U)指针解决的动量方程的求解,但是我的新温度方程中的湍流热扩散系数应当如何更新呢?是应该利用类似于turbulence这样的指针嘛?归根结地就是更新这个alphaEff项
-
谢谢李老师
-
@东岳 老师这一项我看来应该不是简单的温度梯度,之前我有写如果是温度梯度的话这一项就变成了矢量,但按照等式的运算,这一项必定式标量,恳请李老师能指点一下代码实现的思路
-
等式的两边按照道理来说都应该是标量,标量之间进行运算。但是我在写湍流模型的时候,等式右边的第一项感觉让我很为难,翻阅了几遍编程指南也没有找到合适的解决方案。总结一下问题,等式右边的第一项在of中应该如何用代码表示,谢谢各位老师,希望知道的老师可以指点一下 -
有同学碰到过这样的错误嘛?跪求指点,卡了好多天了
-
这张图是报的错误
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
调用新编译的湍流模型进行运算,只计算了一步,出现了这样错误提示,有老师碰到过吗?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
关于并行计算报错的问题
一维不可压缩流体如何计算压力变化?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
OF中标准K-E模型是否做过调整
OF中标准K-E模型是否做过调整
OF中标准K-E模型是否做过调整
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
fluent总能比of得到更满意的结果?
tmp类使用问题,承接上次的帖子问题
tmp类使用问题,承接上次的帖子问题
湍流模型中添加新的方程式编译过程中出现加减操作符调用提示错误
调用新编译的湍流模型进行运算,只计算了一步,出现了这样错误提示,有老师碰到过吗?
调用新编译的湍流模型进行运算,只计算了一步,出现了这样错误提示,有老师碰到过吗?
调用新编译的湍流模型进行运算,只计算了一步,出现了这样错误提示,有老师碰到过吗?
湍流模型中添加新的方程式编译过程中出现加减操作符调用提示错误
是否需要用到指针?
调用新编译的湍流模型进行运算,只计算了一步,出现了这样错误提示,有老师碰到过吗?
是否需要用到指针?
是否需要用到指针?
是否需要用到指针?
是否需要用到指针?
关于标量求偏导在openFoam中的表达代码,希望有大佬指点
关于标量求偏导在openFoam中的表达代码,希望有大佬指点
关于标量求偏导在openFoam中的表达代码,希望有大佬指点
湍流模型中添加新的方程式编译过程中出现加减操作符调用提示错误
湍流模型中添加新的方程式编译过程中出现加减操作符调用提示错误