Skip to content
  • 国内唯一专注于CFD的讨论社区

    3 主题
    325 帖子
    H

    我写 固我在,嘿嘿

  • OpenFOAM交流区

    5k 主题
    32k 帖子
    R

    各位老师们好,目前我在OpenFOAM-v1912中尝试把一个压力基的拉格朗日两相反应求解器(基于coalChemstryFoam)耦合进detonationFoam(密度基)中。目前编译没有问题,可以顺利通过,但是在测试算例时运行会报错(使用的是NS_Sutherland模型):

    === Start Sensible Enthalpy Transport ==== --> FOAM FATAL ERROR: incompatible fields for operation [rhoE] - [h] From function void Foam::checkMethod(const Foam::fvMatrix<Type>&, const Foam::fvMatrix<Type>&, const char*) [with Type = double] in file /share/home/zhoufan/OpenFOAM/OpenFOAM-v1912/src/finiteVolume/lnInclude/fvMatrix.C at line 1337. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 void Foam::checkMethod<double>(Foam::fvMatrix<double> const&, Foam::fvMatrix<double> const&, char const*) at ??:? #3 Foam::tmp<Foam::fvMatrix<double> > Foam::operator-<double>(Foam::tmp<Foam::fvMatrix<double> > const&, Foam::tmp<Foam::fvMatrix<double> > const&) at ??:? #4 ? at ??:? #5 __libc_start_main in /lib64/libc.so.6 #6 ? at ??:? Aborted (core dumped)

    这里显示动量方程和组分方程求解都没有问题,在求解能量方程时遇到了赋值不统一的问题。
    这是原本detonationFoam中rhoEEqn的代码:

    surfaceScalarField sigmaDotU ( "sigmaDotU", ( fvc::interpolate(muEff)*mesh.magSf()*fvc::snGrad(U) + fvc::dotInterpolate(mesh.Sf(), tauMC) ) & fvc::interpolate(U) ); solve ( fvm::ddt(rhoE) + fvc::div(rhoEPhi) - fvc::div(sigmaDotU) == reaction->Qdot() ); e = rhoE/rho - 0.5*magSqr(U); e.correctBoundaryConditions(); thermo.correct(); rhoE.boundaryFieldRef() == rho.boundaryField()*(e.boundaryField()+0.5*magSqr(U.boundaryField())); solve ( fvm::ddt(rho, e) - fvc::ddt(rho, e) + thermophysicalTransport->divq(e) ); thermo.correct(); rhoE = rho*(e + 0.5*magSqr(U)); p.ref() = rho()/psi(); p.correctBoundaryConditions(); rho.boundaryFieldRef() == psi.boundaryField()*p.boundaryField(); Info<< "min/max(p) = "<< min(p).value() << ", " << max(p).value() << endl; Info<< "min/max(T) = "<< min(T).value() << ", " << max(T).value() << endl;

    这是我修改添加颗粒源项之后的代码:

    clock_t t_begin_h = std::clock(); Info << "=== Start Sensible Enthalpy Transport ====" << endl; auto dt = runTime.deltaT(); surfaceScalarField sigmaDotU ( "sigmaDotU", ( fvc::interpolate(muEff)*mesh.magSf()*fvc::snGrad(U) + fvc::dotInterpolate(mesh.Sf(), tauMC) ) & fvc::interpolate(U) ); solve ( fvm::ddt(rhoE) + fvc::div(rhoEPhi) - fvc::div(sigmaDotU) - fvm::laplacian(turbulence->mut()/Prt + turbulence->alpha(), e) == combustion->Qdot() + AlParcels.Sh(e) + radiation->Sh(thermo, e) ); e = rhoE/rho - 0.5*magSqr(U); e.correctBoundaryConditions(); thermo.correct(); rhoE.boundaryFieldRef() == rho.boundaryField()*(e.boundaryField()+0.5*magSqr(U.boundaryField())); solve ( fvm::ddt(rho, e) - fvc::ddt(rho, e) ); thermo.correct(); radiation->correct(); //add by vv rhoE = rho*(e + 0.5*magSqr(U)); p.ref() = rho()/psi(); p.correctBoundaryConditions(); rho.boundaryFieldRef() == psi.boundaryField()*p.boundaryField(); Info<< "min/max(p) = "<< min(p).value() << ", " << max(p).value() << endl; Info<< "min/max(T) = "<< min(T).value() << ", " << max(T).value() << endl; clock_t t_end_h = std::clock(); double elapsed_secs_h = double(t_end_h - t_begin_h) / CLOCKS_PER_SEC; Info << "=== Sensible Enthalpy Transport [Done]: " << elapsed_secs_h << " s ====\n" << endl;

    想请教一下大家有没有什么更改的建议。

  • 基本理论区

    420 主题
    2k 帖子
    Q

    连续体中选择拉个朗日多相流,选的是瞬态,可是没见到如何设置颗粒的时间不长和喷射时间,,请教各位大佬

  • Fluent交流区

    1k 主题
    5k 帖子
    L

    @李东岳 感谢指点

  • 网格交流区

    240 主题
    1k 帖子
    李东岳

    这种没法搞,注定就是深宽比很大

  • 有界的CFD!无界的讨论!

    494 主题
    2k 帖子
    李东岳
    差评获奖感言

    5d92bd29-ea41-4f41-abb9-d497c5fb5c53-image.png

    之前我有一个朋友,聘用了一个月嫂。他们对月嫂的评价,并不是很好。但是他们还发现,之前还有人给月嫂送锦旗。朋友很难理解,为什么有人送锦旗,有人还嫌弃?我当时跟朋友说,每个家的情况不同,因此每个月嫂对每家的适用情况不一样,因此褒贬不一,这很正常

    产品褒贬不一,这在世界上所有的产品都存在的。顺丰有理赔,小米汽车有索赔,华为手机有粉有黑,就算是货币,在汇率低的时候都会被卖。我从来没见过任何一个产品,具有100%的好评。

    我自己非常能够理解产品的好评差评一定是共生的。然而自己的产品设计,遇到不好的评价,要用一生去治愈。

    捕获.JPG

    这个老铁23年9月来的。10月份回去后私下还微信问我几个问题。后来好像把我删了,公众号也没关注,上面是他唯一一条留言。不过他的评论,应该是真实的。

    除了这个课之外,我写的CFD公众号、《无痛苦ns方程笔记》、包括 dyfluid.com, 甚至 cfd-china.com ,也有人褒贬不一。尤其是CFD界公众号,隔三差五就有人骂人。骂人的我直接删了,感觉素质特别低。

    顺着这个话题,我们从技术的问题上,讨论一下潜在的差评来源。

    我的原因 课程没有PPT 课程不承诺给代码 课程结构设计 教室环境

    我的解决方案:

    没有PPT没有解决方法。我只能手推方程。想看PPT,我的课不合适。

    今年会给出一个方案。实在写不出来,你可以迂回的拿走源代码。但课程不给代码的原因,是我要让你们去写代码,而不是我给你们代码!!这个我强调很多遍了。要别人的代码不是学CFD的态度。方程你要自己推,代码你要自己写!

    我每年都会投入非常多的时间去课程设计。但众口难调。我这面风格非常严重。比如方程就是要手推,代码就是要你们自己写。如果你们想直接看方程,直接要能用的代码。我的课不适合你。

    针对教室环境这个。每年会议室都要花个几万块钱,这已经很贵了。我尽力找性价比高的

    学生的原因 本身主动性比较低

    主动性比较低。遇到卡壳你不问我,你就会跟不上。跟不上效果就会很差。

    我一直都在回答各种代码的问题,协助解决。但是有一些人真的是也不听我讲课,自己做跟课程完全无关的事情。因为我一直在教室下面来回走,有人举手我就解决,我能看到大家都在干什么。这我怎么说,我说你为什么不听课?我不可能说。

    主动性低可能有很多原因,一些可能是导师让来的,自己也不太愿意学(我见过一些)。一些可能是感觉学着学着感觉太难了,后面放弃了(较少)。

    最后说一个奇怪的现象

    另外,做这么10多年来,我见过2个课堂上睡觉的。

    只有2个。印象太深了。一个是北京某高校做可压缩的老铁。一个是青岛某研究院的老铁。青岛的老铁拿了个摄像机来录制。我阻止后,后续就趴桌子睡觉不听了。

    我也想不太明白。你花了这么多钱,为什么要来课堂上睡觉?这个我也不好处理,让他们走?人家花钱了。那咋整,那就睡吧。

    之前还有一次课程占了1天中秋假期。有学生留下匿名差评,说课程不考虑假期竟然安排在中秋节,谁中秋节还学习? 这种我属实控制不了。有些时候课程的时间不太好安排,只能这样了。

    正常的课堂效果

    正常的参加我的课程之后,你的效果应该是:

    之前学习OpenFOAM一年了感觉也不会编程,算法也看不太明白。课程5天时间里,主流的CFD算法也能听懂个80%。对算法有全局到细节的认识,一些东西可能还不是很懂,这需要回去消化。代码方面,你能写出来课堂上90%的求解器。虽然一些求解器没跟上。但自己写的时候遇到了上百个bug,问老师后都解决了。5天时间之后,自信心爆棚!回去跃跃欲试想自己植入自己的算法。

    另外,每个求解器,只要写出来的,我都会发礼品。因此我大体知道有多少人能写出来。那种热身的,基本都能写出来。

    复杂的求解器,估计有30%-50%的人能写出来。那剩下的人没写出来,一些人去做自己的不相关的科研去了,一些人临时有事,一些人纯主动性不足不想写这个特定的求解器。

    只要你主动的去吸收,主动问我,100%都能给你解决各种bug,所有求解器都能写出来

    在高中,流弊的老师,带100个学生,考上清华也只有2个。一个班级同样有排名靠后的。要做成任何一件事,都需要双方的努力。我好好备课,你好好去学。才能获得最好的效果。缺一不可。