Skip to content

OpenFOAM

OpenFOAM交流区

5.3k 主题 31.6k 帖子
  • 求多相流经典实验论文

    1
    1 帖子
    1k 浏览
    W

    要做简单例子的定量对比,来拿验证方法可靠性的,涉及到水破碎分离的,列如dambreak就很经典 但是那个实验没有涉及到破碎

  • 使用SSTkomega和laminar算出来的雷诺应力问题

    7
    7 帖子
    6k 浏览
    Y

    @fu 在使用laminar(Quasi-DNS)计算的时候,壁面处的网格Y+<1 ,且整个区域网格也要相对较小;当使用湍流模型的时候,壁面处理跟laminar应该不同,使用壁面函数之类的;如果题主使用一套网格应用于两种流动模型,是否会导致计算结果的偏差呢?还请大佬指点

  • 速度标准差

    3
    3 帖子
    2k 浏览
    疏影横斜水清浅

    @李东岳 好的,谢谢老师!具体的公式我也在查,文献中没有直接给出

  • 速度残差

    3
    3 帖子
    2k 浏览
    H

    @李东岳 谢谢老师回复~
    老师可以帮我看一下这个错误怎么修改吗,按照您之前的一个帖子,定义了一个随时间变化的热流边界~

    HEATED { type externalWallHeatFluxTemperature; mode flux; q coded; // Heat flux code #{ scalar q=*this; scalar t = this->db().time().value(); if(t >= 0 && t<=1.5) { q = 180000; } else if(t<=3) { q = 200000; } else if(t<=4.5) { q = 250000; } else if(t<=6) { q = 300000; } else { q = 0; } q.writeEntry("", os); #}; value uniform 300; kappaMethod lookup; kappa "K"; } /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:32:23: error: cannot convert ‘const Foam::PatchFunction1Types::qPatchFunction1ScalarField’ to ‘Foam::scalar’ {aka ‘double’} in initialization /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:33:31: error: ‘const class Foam::PatchFunction1Types::qPatchFunction1ScalarField’ has no member named ‘db’ /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:55:9: error: request for member ‘writeEntry’ in ‘q’, which is of non-class type ‘Foam::scalar’ {aka ‘double’} /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:55:24: error: ‘os’ was not declared in this scope; did you mean ‘cos’?
  • paraview批量保存descriptiveStatistics的数据

    3
    3 帖子
    3k 浏览

    @氕氘-川 批量保存的脚本老报错: ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
    vtkCSVWriter (000001F6DAFC08C0): Unable to open file: /-0.1150.csv
    ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
    vtkCSVWriter (000001F6DAFCC8C0): Unable to open file: /-0.0650.csv
    ...

    这是我的脚本:
    import numpy as np
    from paraview import numpy_support as ns
    import paraview as pv
    import paraview.simple as ps

    afoam = FindSource('a.foam')

    for y in np.arange(-0.115,0.115,0.05):

    # Make a cut of the wing to have a local 2D profile Slice1 = Slice( Input=afoam ) Slice1.SliceType = "Plane" Slice1.Triangulatetheslice = 0 Slice1.SliceOffsetValues = [0.0] Slice1.SliceType.Origin = [0.0, y, 0.0] Slice1.SliceType.Normal = [0.0, 1.0, 0.0] descriptiveStatistics1 = DescriptiveStatistics(Slice1) descriptiveStatistics1.VariablesofInterest = ['T'] descriptiveStatistics1.AttributeMode = 0 Render() writer = CreateWriter(str(y)+".csv",OutputPort(descriptiveStatistics1, 1)) writer.FieldAssociation = "Cells" writer.UpdatePipeline() del writer #SaveData(str(y)+".csv", proxy=OutputPort(descriptiveStatistics1, 1), FieldAssociation='Cells') descriptiveStatistics1.UpdatePipeline()
  • 6 帖子
    4k 浏览

    @wsxfyy 你可以多试试,有的可以改,有的会报错

  • 请问各位有人遇到过这个问题吗,是什么原因

    4
    4 帖子
    2k 浏览
    T

    @李东岳 多谢东岳老师,重新划分网格后可以运行了

  • 1 帖子
    1k 浏览
    L

    0b9f25fa49feec7e6e3148b3f3df445.png
    这个意思是说我Foam::twoPhaseMixtureThermo::New没定义?
    2099f75c499411d7558623e7b73b800.png
    这个应该算定义吗?能通过报错中的(.text.startup+0X35ed)找到为啥错吗。。
    有遇到过相同问题的朋友吗

  • 2 帖子
    2k 浏览
    Number44N

    试试zeroGradient

  • codedFixedValue边界条件

    3
    3 帖子
    3k 浏览
    D

    @李东岳 好滴好滴,非常感谢,我这就试试:xinxin:

  • coalChemistryFoam求解器添加新物质

    1
    1 帖子
    1k 浏览

    各位前辈中午好,coalChemistryFoam自带case是求解两种物质:煤和石灰石,我想添加新的需要反应的颗粒,请问需要怎么修改求解器呢?谢谢大家。

  • Simple算法中通量phi的计算

    16
    16 帖子
    15k 浏览
    Q

    @李东岳 在 Simple算法中通量phi的计算 中说:

    也就是说在边界处phi是不正确的。不知道是我没有找到phi边界修正,

    边界处的phi不需要计算,是用户给的。加入速度1m/s,乘以网格边界面积,就是phi啊,并且不变。

    如果是pressureInletVelocity 边界呢? 代码里可以看出pressureInletVelocity 需要用边界上的phi来更新出新的边界速度。如果是这种边界边界上的phi又应该如何确认呢?

  • pimpleFoam当运行到nOuterCorrectors=3时报错

    7
    7 帖子
    5k 浏览
    S

    @李东岳 李老师,您好,我按照您的建议进行了相关的测试。结果发现当更换为GAMG时计算比PCG慢非常多,如果我在计算一开始取残差与相对残差值较大时,完全不收敛的情况基本上到500步左右就改善了。但是我又发现一个特别奇怪的现象:
    首先通过pimpleFoam (nOuterCorrectors=1) 计算到一定步时计算结果收敛很好,残差基本下降至1e-5左右;但是当我在运行时更改nOuterCorrectors=2时,初始残差立即上升至0.4左右,计算异常缓慢,并且继续运行算例残差也下降不了。请问这个可能是什么原因导致的呢?

  • 关于dyfluid.com上reactingTwoPhaseEulerFoam的验证算例

    4
    4 帖子
    4k 浏览
    李东岳

    替代文字

    你说上面这个图。这个算力对应的不是这个图。所以不一样。不过如果你增加近气速度的话,应该也是可以的。你加大10倍试试。

  • 3 帖子
    3k 浏览
    Z

    您好,请问这个问题您是否有好的解决办法了吗?
    强制设为0可能不太对,反映出来的可能不是压实效果。。

  • 关于量纲对齐

    5
    5 帖子
    4k 浏览

    @李东岳 东岳老师,这个主要是因为我的是粘弹性流体,在动量方程里边加入了弹性相,不只有雷诺数,还有魏森贝格数,可能是这两个无量纲数导致的:135:

  • owner neighbour 网格寻址

    3
    3 帖子
    4k 浏览
    李东岳

    4年过去了,当时发完这个贴就出去玩了。:jingya:

  • 4 帖子
    3k 浏览
    李东岳

    @ruibin_tongji 你看你的log是不是有些迭代次数太多了

  • 6 帖子
    4k 浏览
    Z

    @是嘿嘿牛啊 密度是被动地,是根据流体种类、温度和压力变化的,你可以设置上下两层温差或不同的流体种类

  • 湍流模型增加一laplacian项问题请教

    3
    3 帖子
    2k 浏览
    W

    @李东岳 噢噢,谢谢东岳老师。