Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. 秋山结月
    秋
    • Profile
    • Following 0
    • Followers 3
    • Topics 9
    • Posts 27
    • Groups 0

    秋山结月

    @秋山结月

    214
    Profile views
    27
    Posts
    3
    Followers
    0
    Following
    Joined Last Online

    秋山结月 Unfollow Follow

    Latest posts made by 秋山结月

    • RE: DNS求解湍流混合层出口处出现异常压力问题

      @李东岳

      李老师您好。

      目前出口是采用如下方程
      3.png

      这里Uc尝试过:

      1. 使用上下层速度的平均值1.5。
      2. 在X方向上采集最靠近边界的三个点的速度值进行平均,作为Uc。

      如果采用2的话,Uc是一个随着Y坐标以及时间变换的函数Uc(y,t)

      Phi是压力泊松方程构造的中间量,这个查了好多文章都说的很模糊,我也就直接用了零梯度。

      posted in Algorithm
      秋
      秋山结月
    • DNS求解湍流混合层出口处出现异常压力问题

      2.png

      各位老师好。目前我所模拟的是两层不同速度流体的混合问题。
      N-S方程的求解采用的是Fractional Step Method,时间步进为三步龙格库塔法。
      上下速度边界采用的是零梯度,左右采用周期性边界条件,进口指定一个速度Profile,出口采用的是对流流出的边界条件。

      在计算过程中我遇到了一个问题,当压力结构到达出口边界时,会发生严重的扭曲。

      P.0000.png
      P.0010.png P.0020.png

      这里主要有几个疑问:

      1. 差分格式,库朗数,时间积分方式会是造成这个现象的原因吗?
      2. 出口的速度边界采用了对流流出边界条件,出口的压力(以及构造压力泊松方程时出现的标量:速度散度)应该采用什么样的边界条件?
      3. 我试着统计了一下压力泊松方程的残差,发现即使迭代100步,残差值 r = Ax-b 依然在0.1这个量级,这正常吗?
      posted in Algorithm
      秋
      秋山结月
    • 未来的路该怎么走呢?

      自己在读本科的时候,大概是对科研一点想法也没有的吧。从小到大被教育好好学习,上了大学就轻松了,于是我就真的大学轻松了四年,之后在找个工作,拿一份能够养活自己的薪水,大概这就是自己的人生了。然而事实却告诉我,并不能这样,2014年恰好遭遇行业严冬,因为学校牌子多少还算是个985,所以虽然自身条件不怎么样,我也顺利进入了一所国企,但是实际的薪资却是差强人意,面对高额的房价再估算一下自己工资未来的涨幅,无奈之下只能工作之余复习考研。

      当时的自己特别想出国读书,看看外面的世界,然而遭到了家里的强烈反对,家长纷纷表示按我本科的表现想要出国读书实在是好高骛远,后来我又表示最少换个专业换个更好的学校,至少避免重蹈本科工作的覆辙,又被家里拒绝了。我个人性格又软(可能这就是我一直做不好的事情的原因),于是就听家里的话老实报回了本校本专业。当时或许也是觉得已经没有退路,反而激发了我的斗志,最后考出了400多分的成绩,让家人后悔不已,没有听我当初的话让我拼一把,不过当时我是觉得本校也有本校的好处,毕竟朋友都在这里生活上也比较开心。于是开开心心地区复试了。

      复试前便是选导师,由于本科没有实验室经历,也就凭本科上课的感觉选了现在的导师。导师年纪不大,但是却非常有为,每年论文稳定产出,学术称号也是一堆,让我觉得我的选择是正确的。当时刚进实验室时,心里是想着虽然硕士没能出去,至少博士能够出去读。然而很可惜,事实却很骨感,导师初始对我很热情,然而得知我出国读博而不是继续跟他的想法之后,我立即变成了“别人的学生”,进入不管不问偶尔打杂的时间。无奈下我开始自己找课题自己弄,也就是这个时期了解到了OpenFOAM。当然最初只是单纯觉得这个做出来图“很好看”,外加许多编程的地方比较多,自己也多少有点编程的天赋,于是就入了这个坑,顺便还能再结合本专业的内容进行研究。因为周围也没有搞这个的,也就自己摸索,外加的确Openfoam在外网开源的资料比较多,也就对照的看一下。由于之前只写过一些面向过程的程序,记得刚看Openfoam这种大量运用面向对象特性的程序感到无从下手,不过毕竟只要愿意下功夫,还是会有收获的。自己也算磕磕绊绊地上手了,虽然自己没有涉及过CFD和FVM底层的东西,但最后也能自己根据自己需要写一些功能,修改一下求解器,也算是能够上手了。这个过程或许是非常痛苦的,然而当有所进步时,也同时是最开心的,自己也是在这个过程中第一次有了做科研很有趣的想法,后来自费参加了东岳老师的暑期培训,也逐渐思考着未来自己能不能向纯CFD算法靠拢。然而个人的能力毕竟是有限的,尤其是对于我们这种工科非理科,除非是算法上很大创新,否则光有数值模拟没有试验很难发出什么高质量的文章。结果满打满算,自己即将毕业时不过手握两篇EI水会,虽然凭借自己努力考出了比较高的语言成绩,然而这是博士并非硕士,语言成绩已经不是什么重要因素,自己的论文又实在不够看,在申请季收到的回复寥寥无几,要么是“实验室已满”,要么是“我没有钱如果你自费我们再考虑你”,于是自己现在还是前路未卜。

      现在离毕业只有3个多月,我却手里还没有收到什么像样的offer,加上春招由于疫情缘故也转为网络形式,比起往年形式不容乐观。自己又是进又不得,退又不得,我现在26岁,家里人又一直催我找工作,早点安定下来。我却每天生活在痛苦当中,深感这三年的付出如同泡沫一般付诸东流,自己无论怎么折腾都只是在瞎蹦跶。我又该何去何从呢....真的就要这么放手吗....

      posted in C斯达克
      秋
      秋山结月
    • 多相流求解器的造波问题

      可能我这个问题比较奇葩,目标是要在multiphaseInterFoam里引入造波功能,造波功能主要基于waves2Foam。我的思路是这样的:

      首先我观察了一下intefoam跟multiphaseInterFoam对于mixture的求解区别。发现主要在于,interfoam是把两相分为主相以及背景相,求解时只求解主相,这样背景相可以简单的进行反向选取,即alpha2=1-alpha1就能获得了。

      但是相数超过2,就不存在这个背景相的问题了,在multiphaseMixture里,每一相都单独计算,最后累加出一个总相alphas出来。

      对于解速度压力那一套,两者我感觉没有太大区别。

      我打算先从两相做起,在multiphaseInterFoam里先实现一般的水-气两相造波,再考虑第三相的问题。也就是打算用multiphaseInterFoam复现interfoam的功能。

      我按照以下步骤进行:

      1.修改multiphasemixture的代码,提供一个获取各个相的相分数的接口。
      2.修改multiphaseInterFoam的代码,将这个接口提供给原先造波消波的松弛区,每一步都对波浪进行修正。
      2.修改waves2Foam的代码,对空气设定一个相分数的入射边界(实际上就是波浪相分数边界的反向选取)。
      3.在求解器中,每次循环都同时对alpha.water以及alpha.air在边界上进行更新。

      比如对于水,边界上入射的相分数可能是(假设2D算例,竖直方向共5个网格):
      (0 0 0.6 1 1)
      那么对于空气,就是
      (1 1 0.4 0 0)
      这样进行累加计算alphas的时候保证边界气液相分数累加为1,我在每一步的文件夹下观察结果,也的确没有发生越界。

      但是实际运算的时候,出现了气液相分数越界,总相分数不守恒的现象,最后造出来的波是这种形状的
      1.jpg

      因为multiphaseInterFoam我想用在化工方面比较多一点,所以我想请问一下我是遗漏了哪些部分吗?

      也请@东岳 老师指点。

      第二个问题就是,我这个问题其实也类似一个分层流动,从上到下,第一相,第二相,第三相之间,第一相跟第三相是没有任何接触的。有没有划分成两个计算域,在第二相上交换信息的办法?

      posted in OpenFOAM
      秋
      秋山结月
    • OpenFOAM中的PotentialFoam与势流理论的区别

      本人在学习OF之前,使用的商业软件基本都是依赖势流理论,引入格林函数基于边界元理论求解Laplace方程,最后用伯努利方程估算压力。

      今天发现OF中也有一个potentialFoam是基于势流理论的,但是似乎跟我了解的势流理论有很大区别,只能计算稳态问题,但是控制方程的确也是Laplace方程。看起来两种求解器不单纯是FVM与边界元这种离散方式的区别,而是还有其他的不同,请问有了解这一块的大佬能讲讲这其中的问题所在吗?

      posted in OpenFOAM
      秋
      秋山结月
    • RE: 关于重叠网格的疑问

      @wwzhao 看起来显示的只是个并行计算的信息?
      总而言之谢谢老师的指点了

      posted in OpenFOAM
      秋
      秋山结月
    • RE: 关于重叠网格的疑问

      @wwzhao 老师我想请教一个问题。我看到一些论文把重叠网格分为插值边界单元,贡献单元以及孤点单元,我想问着是否跟OpenFOAM的重叠单元的interpolated from local , mixed local/remote 以及 remoted相对应呢?

      我计算的时候遇到了remoted单元不断上升的情况(初始0,计算100s后可能上升到几千),但是我的计算确实也没有发散,能够进行下去,如果remoted对应孤点单元,是不是代表我这个计算设定是不准确的呢?

      posted in OpenFOAM
      秋
      秋山结月
    • RE: SnappyHexMesh中的使用问题

      @chaoscfd
      22233333.png
      想请问一下我是这是什么回事呢。。。
      好像snap阶段结束后..就直接真的没了...
      @东岳 顺便请教一下东岳老师

      posted in OpenFOAM
      秋
      秋山结月
    • RE: SnappyHexMesh中的使用问题

      @chaoscfd
      这里想请问您几个问题:

      1. 建立refineRegions,我这里是用distance的办法定义的refineRegion,离模型xx米内用xx等级这样,我还需要另外用box对其局部定义吗? 还是说我只需要提高distance的等级就可以了?
        2.第二就是includeAngle,我这里取得是90°,我之前尝试过取180°,的确新选取的边很多都是我想要的,但是也出现了好多我不想要的edge(比如明明就是个连续的面,结果突然出现一个edge斜着贯穿了这个面,类似一个正方型强行被分成2个三角形)。有没有什么办法,只捕捉到我想优化的edge,而抛弃一些莫名其妙捕捉到的edge呢
      posted in OpenFOAM
      秋
      秋山结月
    • RE: SnappyHexMesh中的使用问题

      @东岳 好的 明白了 不过咱的求解器是基于InterFoam自己改的...我试试算的效果如何吧
      谢谢老师

      posted in OpenFOAM
      秋
      秋山结月