Skip to content

OpenFOAM

OpenFOAM交流区

4.5k 主题 26.5k 帖子
  • OpenFOAM版本区别

    3
    3 帖子
    1k 浏览
    Y

    @李东岳 好的,感谢李老师回复,我们服务器系统是centos6.5,好像不太支持新版本:136:

  • 想知道能量方程粘性项作用该怎么写

    7
    7 帖子
    1k 浏览
    L

    @lyt 难道是后面一定要加两个()(),才行?这是什么逻辑

  • 计算中断的问题

    3
    3 帖子
    958 浏览

    @caijinjin 是计算停在这里还是终端会停止,方便的话可以交流一下吗,我微信15098932699

  • 《无痛苦N-S方程笔记》纠错

    2
    2 帖子
    816 浏览
    bestucanB

    按说,那个 pdf 的大部分勘误记录都会跳转到这里
    https://cfd-china.com/topic/2338/无痛苦n-s方程笔记-勘误/38

    “错”和“误”不一样,(所以某些理论上无法消除、工具条件限制带来的与事实不符的差距叫“误差”不叫“错差”)(我对你措辞没意见,只是看你挺有互助热情,也“互助”你一下,说给你听听,参考参考。)

    可以直接说出笔误位置,细微的笔误难以发现(容易发现也不会误了)。比如:

    "接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x1, y1, z1) 点的时候"
    这句第二个坐标(或者第二个(x1, y1, z1))标号应该是(x2, y2, z2)

    我看你之前也发过类似的帖子,可以看完了一遍书,把发现的笔误统一发一下,也方便修改。

    欢迎继续:146:

  • 6 帖子
    1k 浏览
    O

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

  • 2 帖子
    776 浏览
    bestucanB

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

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

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

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

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

  • 如何过滤准DNS结果

    2
    2 帖子
    706 浏览
    Z

    顶一下

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

    2
    2 帖子
    573 浏览
    bestucanB

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

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

    6
    6 帖子
    2k 浏览
    strive_DUTS

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

  • 4 帖子
    885 浏览
    李东岳

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

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

    已移动
    3
    3 帖子
    870 浏览
    L

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

  • 求多相流经典实验论文

    1
    1 帖子
    429 浏览
    W

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

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

    7
    7 帖子
    2k 浏览
    Y

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

  • 速度标准差

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

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

  • 速度残差

    3
    3 帖子
    862 浏览
    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’?
  • SST k-ω湍流模型的初始条件咋取呀?

    3
    3 帖子
    1k 浏览
    李东岳

    基本上,不发散能对的上就行 :mihu:

  • OpenFOAM中刚盖假定

    2
    2 帖子
    569 浏览
    李东岳

    刚盖是不动的盖子? 速度为0不就可以了么

  • paraview批量保存descriptiveStatistics的数据

    3
    3 帖子
    860 浏览

    @氕氘-川 批量保存的脚本老报错: 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 帖子
    1k 浏览

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

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

    4
    4 帖子
    870 浏览
    T

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