Skip to content

OpenFOAM

OpenFOAM交流区

5.3k 主题 31.5k 帖子
  • 6 帖子
    4k 浏览
    O

    @chaoscfd 我用了 但还是mesh.findcell最快

  • 2 帖子
    2k 浏览
    bestucanB

    我猜是这样的,
    process 里必须有动网格,这样计算的过程只涉及网格如何动,

    如果 process 里没有动网格,那么当动网格运动到该区域,还要造动网格。

    至于 process 里没有动网格也能算,那是动网格没有运动到该 process 覆盖的区域。

    可以测试,在计算区域内设置动网格,但是不让它动,比如平静的水面。应该可以一直算。

    paraview 里可以看每个 process 的区域,算不动的那步估计就是跨区域运动了

  • 如何过滤准DNS结果

    2
    2 帖子
    2k 浏览
    Z

    顶一下

  • 请问小红帽系统版本不一致还能安装OF吗?

    2
    2 帖子
    2k 浏览
    bestucanB

    这个是安装 docker 的问题,不是 OF 的问题 :wocao:
    去 docker 相关社区问更容易得到答案:chigua:

  • OpenFOAM后处理怎么提取某点随时间变化的数据

    6
    6 帖子
    5k 浏览
    strive_DUTS

    @tens 再次感谢您!重新尝试了一下,可以输出了!

  • 4 帖子
    3k 浏览
    李东岳

    应该影响不大。这种颗粒应该不多,不会所有的颗粒都恰好那样。所以对alpha的影响也很小。除非颗粒直径非常大。不过直径非常大的话,那就会引起收敛性的问题。

  • 关于相变求解器interphasechangeFoam源项的疑问

    已移动
    3
    3 帖子
    3k 浏览
    L

    b75ecf32b61ea077acecbb62002a6c8.png
    这是我后来自己的解释

  • 求多相流经典实验论文

    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 帖子
    2k 浏览

    @氕氘-川 批量保存的脚本老报错: 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左右,计算异常缓慢,并且继续运行算例残差也下降不了。请问这个可能是什么原因导致的呢?