Skip to content

OpenFOAM

OpenFOAM交流区

4.5k 主题 26.5k 帖子
  • 怎么展示一个切面上的拉格朗日粒子

    6
    6 帖子
    890 浏览
    J

    @星星星星晴 感谢大佬指点!

  • sixDoFMotion 并行计算出错问题

    2
    2 帖子
    474 浏览
    李东岳

    看起来像个bug。你可以试试双核的,然后手动把traslationRotation写进processor里面,然后看看报错不。不报错的话可能是bug

  • 用MPPICInterFoam做泥沙冲淤的问题

    3
    3 帖子
    599 浏览
    T

    @李东岳 谢谢李老师,我再看看

  • 关于2D激波管验证

    8
    8 帖子
    1k 浏览
    ZY-GONGZ

    @李东岳 :135: 不好意思,李老师,再打扰您一下,我还有一个问题,我目前使用的是v10版本的multiphaseEulerFoam,我在dragModel里可以通过interface_.dispersed().d()这种来调用直径参数,现在我想在求解器的.c程序里调用它,应该用什么方式呢
    我注意到您再小代码里提到过以下代码

    meanDiameter { type coded; libs ("libutilityFunctionObjects.so"); name error; codeExecute #{ const volScalarField& d = mesh().lookupObject<volScalarField>("d.alpha.oil"); scalar d32 = d.weightedAverage(mesh().V()).value(); if (Pstream::master()) { std::ofstream file; file.open ("d32", std::ofstream::out | std::ofstream::app); file << mesh().time().timeName() << " " << d32 << "\n"; file.close(); } #}; }

    这段我是否能用呢

  • interPlicFoam

    90
    90 帖子
    86k 浏览
    J

    @队长别开枪 好的 谢谢前辈!

  • 请教simpleFoam计算时smoothsolver无后续的问题

    9
    9 帖子
    1k 浏览

    @李东岳 好的,谢谢老师,有错误信息就好办,看来还是得转战linux:135:

  • 边界赋值问题

    2
    2 帖子
    505 浏览
    李东岳

    下面这个代码是对壁面一层网格进行赋值的。你可以参考一下。mag可以。

    //- Initialise the near-wall P field to zero forAll(patches, patchi) { const fvPatch& currPatch = patches[patchi]; if (isA<wallFvPatch>(currPatch)) { forAll(currPatch, facei) { label faceCelli = currPatch.faceCells()[facei]; epsilon[faceCelli] = 0.0; G[faceCelli] = 0.0; } }
  • 并行计算时间步储存文件有小数点

    1
    1 帖子
    281 浏览
    yfleungY

    a9956a74-d12c-48a9-8f31-5da7020a5fc6-image.png 各位大佬下午好,我最近在跑一个微流道的圆柱绕流,网格尺寸非常小,然后计算的时间步长选取的也比较小,一个时间步长是5e-5,存储时间为0.5秒存储一次,当时在存储的时候会有小数点后几位出现,各位大佬有知道这是什么原因导致的吗:chouchou:

  • openfoam中黏弹流算例

    3
    3 帖子
    520 浏览
    yfleungY

    @李东岳 好的,谢谢:chouchou: 我已经找到对应的黏弹性算例了。

  • 22 帖子
    6k 浏览

    你好,想问一下你问的这个问题解决了吗,我也遇到了一样的问题。如果你已经解决了,能麻烦你提点指教一下吗

  • 摩擦速度比论文偏小

    5
    5 帖子
    2k 浏览
    H

    请问您最后是怎样计算磨擦速度的呢,我的算例计算磨擦速度也偏小一点。

  • 关于网格高度变化的动网格设置

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

    @田畔的风 好的,谢谢。我是看到你在后面又把空气密度除了一次,所以存在疑惑。如果你取空气密度为1,那么在数值上应该是没问题的。

  • 壁面剪切力的计算

    5
    5 帖子
    978 浏览
    疏影横斜水清浅

    @李东岳 老师,求出的数值结果是6e9dde13-6ec5-40ac-98f1-7d6df59514f0-image.png ,还是2c1168c0-20c4-45ef-8e51-45a5f53de754-image.png

  • icem网格导入openfoam错误

    4
    4 帖子
    300 浏览
    F

    @李东岳 好的,谢谢老师指点:chouchou:

  • LES-RANS对比,freejet

    1
    1 帖子
    377 浏览
    李东岳

    前几天好像论坛有个老铁用LES算freejet,图看起来不太多。刚才我发现了一个图,这个图看起来挺符合物理的

    776ed1a2-8058-43fb-91e1-b5f3013a7b1a-image.png

    https://github.com/unicfdlab/QGDsolver

    仔细看了下,好像是超音速

  • 使用mapFields之后报错,出现No MRF models present

    2
    2 帖子
    452 浏览
    李东岳

    No MRF models present不是问题

    你这个看起来是壁面函数修正的时候,有0存在。我看你是高超音速。是不是你算的是无粘流。

  • 8 帖子
    2k 浏览
    李东岳

    一楼那个图,如果不想让速度分散的话,感觉只能在上下两层加一层多孔介质,不知道其他大佬怎么看?

  • SOWFA问题请教

    4
    4 帖子
    2k 浏览
    李东岳

    @zzcfd

    这个边界条件我也没有用过。不过我刚才看了一下源代码。你看一下mapped和TimeVaryingMappedFixedValue的区别:

    //- Name of the field data table, defaults to the name of the field word fieldTableName_; //- If true adjust the mapped field to maintain average value Switch setAverage_; //- Fraction of perturbation (fraction of bounding box) to add scalar perturb_; //- Interpolation scheme to use word mapMethod_; //- 2D interpolation (for 'planarInterpolation' mapMethod) autoPtr<pointToPointPlanarInterpolation> mapperPtr_; //- List of boundaryData time directories instantList sampleTimes_; //- Current starting index in sampleTimes label startSampleTime_; //- Interpolated values from startSampleTime Field<Type> startSampledValues_; //- If setAverage: starting average value Type startAverage_; //- Current end index in sampleTimes label endSampleTime_; //- Interpolated values from endSampleTime Field<Type> endSampledValues_; //- If setAverage: end average value Type endAverage_; //- Time varying offset values to interpolated data autoPtr<Function1<Type>> offset_; Example of the boundary condition specification: \verbatim <patchName> { type mapped; fieldName T; setAverage no; average 0; interpolationScheme cell; value uniform 0; } \endverbatim

    TimeVaryingMappedFixedValue看起来就是附加了开始sample的时间以及结束sample时间的mapped边界条件。在时间之外,就是fixedValue。你核实下。

  • OpenFOAM并行运算。

    2
    2 帖子
    537 浏览
    李东岳

    @zz给我趴下 在 OpenFOAM并行运算。 中说:

    前几天我服务器里在跑着算例,但办公室停电了。来电之后我想接着跑算例,结果发现算例跑不了了

    这是一个很诡异的问题。毫无理由啊

  • 2 帖子
    491 浏览
    田畔的风

    算是解决了吧,其实大可不必费力去猜,先写个Info看看这玩意执行的时间点合不合适,然后直接访问fvMesh对象,拿个引用出来就完事了:mihu:

    codeCorrect #{ volVectorField &U = mesh().lookupObjectRef<volVectorField>("U"); #}