• open foam的参考压力

    OpenFOAM
    0 赞同
    1 帖子
    5 浏览

    open foam中使用interphasechangefoam进行空化流仿真时,因为出现负压,所以想设置参考压力,使得负压变成正值
    pRefCell 0;
    pRefValue 设置为出口压强;
    但是这样也没有起到效果,有没有大佬知道怎么设置呢

  • 0 赞同
    3 帖子
    250 浏览

    请问题主解决此问题了吗?

  • 关于paraview画图

    OpenFOAM
    0 赞同
    16 帖子
    1k 浏览

    一眼望去,全是进步空间

    哈哈,幽默,有意思

  • OpenFOAM壁面y+验证:低雷诺数模型

    OpenFOAM
    0 赞同
    4 帖子
    397 浏览

    加密网格不能使用壁面函数,同时要使用低雷诺数湍流模型

  • 0 赞同
    2 帖子
    25 浏览

    继续更新,后面发生了很奇怪的事情,我试着重新编译,但是报错了,报错提示如下:看他的意时缺少了一个分号,但是我觉得这不合理,我仔细检查了语法,没有错误。

    image.png

    之后我重新尝试把代码复制到另外一个文件夹:但是却报了一个新的错误,不知道咋办了。有大佬给点建议吗,这可能是什么原因呢?我试着编译一些其他的代码,编译能成功。

    image.png

  • open foam空化流出现负压

    OpenFOAM
    0 赞同
    1 帖子
    22 浏览

    在open foam用interphasechangefoam进行两相流空化计算,空化模型SS模型,入口为速度入口,出口为压力出口,出口压力设置为122300Pa,空化压强2300,计算出的p_rgh和p都有负压区,该怎么解决呢

  • paraView调整显示比例

    OpenFOAM
    0 赞同
    6 帖子
    2k 浏览

    @Samuel-Tu 请问这个功能在paraview的哪里可以找到呢

  • 关于LES的dynamicSmagorinsky一些问题

    OpenFOAM
    0 赞同
    12 帖子
    230 浏览

    @李东岳 f084b820-6ac4-4a71-aa27-4cc598955778-image.png 李老师,这是我文献中看到的,这种是咋实现的,没什么思路现在

  • 关于晃荡惯性力的植入

    OpenFOAM
    0 赞同
    10 帖子
    117 浏览
    应该是重力作用大于惯性力的作用, 你得植入效果已经出来了,你核实一下这个乘的rho应该是水的固定值密度,还是水跟空气混合的密度
  • 在安装CFDEM时,安装不了LPP

    OpenFOAM
    0 赞同
    3 帖子
    86 浏览

    请问,就是改成下面这个解决的吗

  • 0 赞同
    8 帖子
    5k 浏览

    您好,安装编译所需文件时报错:Package libvtk6-dev is not available, but is referred to by another package.
    This may mean that the package is missing, has been obsoleted, or
    is only available from another source

    E: Package 'libvtk6-dev' has no installation candidate
    E: Unable to locate package python-numpy
    怎么解决呢

  • bound 函数的实现方法

    OpenFOAM
    0 赞同
    6 帖子
    48 浏览

    bound 函数的实现到底为什么要这样。

    这个看起来是个很简单的数值操作。不是从物流守恒定律搞出来的。这就是Henry自己搞的。他很多东西都自己搞。

    average这个函数我最大的印象就是dynSmagorinsky模型,Henry搞这么个东西进去,整个场的动态Cs变成了一个标量,后来让alberto怼有bug,Henry直接给dynSmagorinsky模型给删了。到现在OpenFOAM也没有动态Smagorinsky模型。 https://bugs.openfoam.org/view.php?id=99

  • Smagorinsky模型系数问题

    OpenFOAM
    0 赞同
    42 帖子
    683 浏览

    最近没细看这个玩别的去了。我参考Spectral analysis of large eddy simulations of the convective boundary layer昨天晚上推了个八九不离十。这个ck ce是符合湍流能谱的一个常数。这是前人已经做过的工作。没法写英文了。我以为别人没玩过这个。哈哈

  • 0 赞同
    5 帖子
    50 浏览

    仿真的结果、数据全都保存于虚拟机之外的Windows文件夹中,而不是储存在虚拟机中,因此只给虚拟机20~30G,它基本上都用不完。

    @Prometheus10 这太流弊了,CFD中文网yyds,82年的salute送上!

    :high: :high: :high: :high: :high:

  • 0 赞同
    11 帖子
    134 浏览

    @coolhhh
    从pEqn.H的代码来看,并没有这种water / gas之分
    @李东岳

    没试过单核跑,单核算得会非常慢,而且估计得算几步才能看结果,我试一下。 correctPhi是必须要是yes的,改成no的话会造成自适应加密出问题,流场所有地方的alpha.water会变得不严格是1,比1小一些,然后流场的所有地方都会加密。
  • 如何人为将曲率给定为一个定值

    OpenFOAM
    0 赞同
    8 帖子
    108 浏览

    @李东岳 好的,东岳老师,我试试按照您说得这样算一下。

    另外我应该补充一下,100μm直径气泡的实验测得的浮升终速度大约是0.005~0.008m/s.

    e584ada7-1a05-4bb6-8cd8-f1e1997c43cc-image.png

    参考文献 Pawliszak, P., Ulaganathan, V., Bradshaw-Hajek, B.H., Manica, R., Beattie, D.A., Krasowska, M., 2019. Mobile or Immobile? Rise Velocity of Air Bubbles in High-Purity Water. J. Phys. Chem. C 123, 15131–15138.
  • 0 赞同
    2 帖子
    23 浏览

    可以没问题。估计要跑好久。

  • 求一篇关于 VOF 方法 MULES 的文章

    Algorithm
    0 赞同
    9 帖子
    4k 浏览

    @thegame 感谢老铁,将在下一版本更新

  • 0 赞同
    45 帖子
    9k 浏览

    @xiezhuoyu Nice! 挺有意思

    Description 28 Two-scheme Courant number based blending differencing scheme. 29 30 Similar to localBlended but uses a blending factor computed from the 31 face-based Courant number and the lower and upper Courant number limits 32 supplied: 33 \f[ 34 weight = 1 - max(min((Co - Co1)/(Co2 - Co1), 1), 0) 35 \f] 36 where 37 \vartable 38 Co1 | Courant number below which scheme1 is used 39 Co2 | Courant number above which scheme2 is used 40 \endvartable 41 42 The weight applies to the first scheme and 1-weight to the second scheme. 43 44 Example of the CoBlended scheme specification using LUST for Courant numbers 45 less than 1 and linearUpwind for Courant numbers greater than 10: 46 \verbatim 47 divSchemes 48 { 49 . 50 . 51 div(phi,U) Gauss CoBlended 1 LUST grad(U) 10 linearUpwind grad(U); 52 . 53 . 54 } 55 \endverbatim 56 57 SourceFiles 58 CoBlended.C
  • 刚入门fluent问一个简单的问题

    Fluent
    0 赞同
    4 帖子
    64 浏览

    @李东岳 感谢东岳老师的回复!

  • 初始速度场给定

    OpenFOAM
    0 赞同
    10 帖子
    1k 浏览

    @fubianhanshu 你好可以加个联系方式吗,想请教一些问题

  • 0 赞同
    20 帖子
    487 浏览

    @Prometheus10

    期待期待

  • 0 赞同
    8 帖子
    2k 浏览

    您好Ivan,想请问油水分离的模拟中,水连续相,微小油滴作为分散相,这种两相液体如何设置,微小油滴可以像拉格朗日粒子那样设置一定的大小和占比,然后自动生产嘛?谢谢您

  • 0 赞同
    9 帖子
    65 浏览

    @李东岳
    懂了,谢谢李老师!

  • 质量入口边界

    OpenFOAM
    0 赞同
    3 帖子
    64 浏览

    @李东岳质量入口边界 中说:

    你考虑的是多相流还是单相流

    李老师您好,感谢您的回复我考虑的是多相流,但是入口处只有气体,也只有气体进入

  • 运行结束后输出密度rho

    OpenFOAM
    0 赞同
    11 帖子
    111 浏览

    @wangfei9088 了解了,感谢:142:

  • 重新编译热物理库

    OpenFOAM
    0 赞同
    6 帖子
    116 浏览

    请问您成功了吗,我也有这方面的需要

  • 0 赞同
    82 帖子
    4k 浏览

    看到 cfd online 的后处理代码:compressibleyPlusLES,感觉修改下作为后处理工具,不失为一个更简单直接的方法。

    volScalarField yPlusTemp ( IOobject ( "yPlusTemp", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedScalar("yPlusTemp", dimless, 0.0) ); volVectorField uPlus ( IOobject ( "uPlus", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedVector("uPlus", dimless, vector::zero) ); volScalarField uTau ( IOobject ( "uTau", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedScalar("uTau", dimVelocity, 0.0) ); volScalarField yPlus ( IOobject ( "yPlus", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedScalar("yPlus", dimless, 0.0) ); Info<< "Reading field UMean\n" << endl; volVectorField UMean ( IOobject ( "UMean", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh ); Info<< "Reading field U\n" << endl; volVectorField U ( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh ); Info<< "Reading field rho\n" << endl; volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ), mesh ); # include "compressibleCreatePhi.H" Info<< "successs1\n" << endl; autoPtr<basicThermo> pThermo ( basicThermo::New(mesh) ); basicThermo& thermo = pThermo(); Info<< "successs2\n" << endl; autoPtr<compressible::LESModel> sgsModel ( compressible::LESModel::New(rho, UMean, phi, thermo) ); /*---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License This file is part of OpenFOAM. OpenFOAM is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. OpenFOAM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>. Application yPlusLES Description Calculates and reports yPlus and uPlus for all wall patches, for the specified times, for compressible LES. \*---------------------------------------------------------------------------*/ #include "fvCFD.H" //#include "incompressible/singlePhaseTransportModel/singlePhaseTransportModel.H" //aali #include "compressible/LES/LESModel/LESModel.H" #include "basicThermo.H" //ali #include "wallFvPatch.H" #include "LESModel.H" #include "nearWallDist.H" #include "wallDist.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // int main(int argc, char *argv[]) { timeSelector::addOptions(); #include "setRootCase.H" # include "createTime.H" instantList timeDirs = timeSelector::select0(runTime, args); # include "createMesh.H" forAll(timeDirs, timeI) { runTime.setTime(timeDirs[timeI], timeI); Info<< "Time = " << runTime.timeName() << endl; fvMesh::readUpdateState state = mesh.readUpdate(); // Wall distance wallDist y(mesh, true); if (timeI == 0 || state != fvMesh::UNCHANGED) { Info<< "Calculating wall distance\n" << endl; Info<< "Writing wall distance to field " << y.name() << nl << endl; y.write(); } #include "createFields.H" volScalarField::GeometricBoundaryField d = nearWallDist(mesh).y(); volScalarField muEff(sgsModel->muEff()); const fvPatchList& patches = mesh.boundary(); dimensionedScalar uTauAvg("uTauAvg", dimVelocity, 0); const volScalarField muLam(sgsModel->mu()); scalar nPatch = 0; Info<< "Summary: " << nl << endl; forAll(patches, patchi) { const fvPatch& currPatch = patches[patchi]; if (typeid(currPatch) == typeid(wallFvPatch))//isA<wallFvPatch>(currPatch)) { yPlusTemp.boundaryField()[patchi] = d[patchi] *sqrt ( muEff.boundaryField()[patchi]/rho.boundaryField()[patchi] *mag(UMean.boundaryField()[patchi].snGrad()) ) /muLam.boundaryField()[patchi]/rho.boundaryField()[patchi]; const scalarField& YpTemp = yPlusTemp.boundaryField()[patchi]; //effective viscosity used in the calculation of u* uTau.boundaryField()[patchi] = sqrt ( sgsModel->muEff() *mag(UMean.boundaryField()[patchi].snGrad()) /sgsModel->rho() ); const fvPatchScalarField& uTauWall = uTau.boundaryField()[patchi]; dimensionedScalar uTauTmp("uTauTmp", dimVelocity, average(uTauWall)); uTauAvg += uTauTmp; nPatch ++; Info<< "Patch " << patchi << " named " << currPatch.name() << " y+ : min: " << min(YpTemp) << " max: " << max(YpTemp) << " average: " << average(YpTemp) << " avgUGradWall: " << average(mag(U.boundaryField()[patchi].snGrad())) << nl << endl; } } uTauAvg /= nPatch; Info << " avg. friction velocity uTau is: " << uTauAvg.value() << " (averaged over " << nPatch << " wall(s))" << nl <<endl; //yPlus is = u* x y/nu yPlus = y.y() * uTauAvg / ((sgsModel->mu()) / (sgsModel->rho())); uPlus = UMean / uTauAvg; Info << "Writing yPlus and uPlus to corresponding fields." << nl <<endl; yPlus.write(); uPlus.write(); // Info<< "Writing yPlus to field " // << yPlus.name() << nl << endl; // // yPlus.write(); } Info<< "End\n" << endl; return 0; } // ************************************************************************* //
  • openfoam添加湍流入口方法?

    OpenFOAM
    0 赞同
    43 帖子
    11k 浏览

    @dyb 你好我可以请问一下这个边界的几个参数的设置方法么,就是LField、RField和refField,我拜读过你的硕士论文但是里面我看没有做细致介绍,所以可以请教一下么,我也是搞风电现在也需要这个和les做耦合

  • 0 赞同
    7 帖子
    124 浏览

    @队长别开枪
    内柱:Cylinder z1(0 0 -20) z2 (0 0 15) 半径 6.25
    外环:z1(0 0 -7.25) z2(0 0 5.25) 内半径 7.5 外半径 12.5
    外环原先是关于水面对称,下降一米让其上浮

    decayHeave1V0.rar

    麻烦大佬给看看

    据另一位做六自由度双浮箱的同学说,他也遇到了类似的问题 只要overset之间相互接触就会发散,不知问题到底出在哪

  • openfoam中添加新变量并输出

    OpenFOAM
    0 赞同
    3 帖子
    65 浏览

    我也不是很懂,从代码角度来看,貌似只有进入了if判断语句中才会进行变量的输出?是不是得把IO操作放在外层

  • 湍流模型中的湍动能输运方程源项

    OpenFOAM
    0 赞同
    5 帖子
    52 浏览

    @李东岳 好的,谢谢老师

  • 0 赞同
    5 帖子
    52 浏览

    空的是指真空还是空气?
    空气就两相流(入口是water,计算域本来是air),真空就知识盲区了。

  • 0 赞同
    9 帖子
    1k 浏览

    @wsy11 那就是你的blockmesh 设置 把底面 ground 在设置低一些

  • 0 赞同
    4 帖子
    98 浏览

    多谢两位老师

  • scalarField和volScalarField

    OpenFOAM
    0 赞同
    6 帖子
    84 浏览

    @kcol 相当于volScalarField,但是只定义内部网格,在边界网格上没有定义

  • 0 赞同
    7 帖子
    177 浏览

    @yhdthu 好的 谢谢大佬

  • fluent湍流粘度超限,多相流发散

    Fluent
    0 赞同
    1 帖子
    60 浏览

    IMG_20231127_193223.jpg
    在计算单相流作为初场,可以收敛,但是湍流粘度比超限,继续设置为多相流空化计算的时候,开始发散,不知道有哪些解决办法,已经重新画过网格。
    边界条件是速度入口,压力出口,湍流sst k-w,稳态计算

  • FLUENT中湍流粘性比超标

    Fluent
    0 赞同
    8 帖子
    3k 浏览

    @李东岳 老师,最近我也遇到了这个问题,单相流计算的时候会提醒,但是结果也可以收敛,但是计算完初场的时候,改为多相流空化计算,马上就发散了,不知道怎么可以解决这个问题,已经重新画了好几次网格了

  • Wray-Agarwal湍流模型

    已移动 OpenFOAM
    0 赞同
    16 帖子
    2k 浏览
    一个流弊的低雷诺数SA模型在CFD中文网于2023年11月27日诞生!
  • 不同湍流强度下的来流

    OpenFOAM
    0 赞同
    6 帖子
    96 浏览

    如果要改湍流强度,你直接用codedFixedValue去重写一下入口InletK会更好操作一些。

  • 编译overWaveDyMFoam后生成波浪出错

    OpenFOAM
    0 赞同
    2 帖子
    52 浏览

    我没有这两个求解器啊老铁 ..

  • 0 赞同
    9 帖子
    125 浏览

    @coolhhh 对,我使用这个命令:alias Anaconda3='PATH=/home/install/Anaconda3/bin:$PATH'时,还没有把anaconda的环境给注释掉,OpenFOAM也可以正常使用,但是并行命令不能使用。
    我大概懂大佬你的意思了,如果使用了alias命令,只要没有输入Anaconda3应该就不会激活anaconda环境。可能是因为我这边前人在给anaconda设置环境的时候,有点繁琐(具体我也没弄清楚),所以我可能只使用一条alias命令,不能完全把相应的anaconda环境给包含进去。所以才导致使用alias Anaconda3='PATH=/home/install/Anaconda3/bin:$PATH'命令还不能消除使用并行命令带来的报错。再次感谢您的回答,谢谢!

  • OpenFOAM 浮体铰接可以计算吗

    OpenFOAM
    0 赞同
    6 帖子
    111 浏览

    @hachikoi 就是说目前还没有实现需要自己修改源码是吗

  • 0 赞同
    1 帖子
    88 浏览

    分享几个自己了解的开源软件:

    OpenFOAM,不用说了。

    T-Flows,也是一款基于FVM的求解器,没咋用过,听说是小型OF,没有OF那么复杂的类和求解器类型。链接文本

    REEF3D,一款船海领域的开源CFD软件。链接文本

  • 《无痛苦N-S方程笔记》勘误

    C斯达克
    0 赞同
    52 帖子
    13k 浏览

    @疏影横斜水清浅 多谢老铁,将在下一版本更新

    @知行 老铁我这之前忘了回你了,才看到,我处理后将在下一版本更新

  • 0 赞同
    16 帖子
    3k 浏览

    老哥们,你们是如何解决这个问题
    我在使用两套overset做浮体自衰减的时候也遇到类似的问题

  • 0 赞同
    8 帖子
    346 浏览

    @LuoSsnappyHexMesh对复杂地形添加边界层报“浮点数例外 (核心已转储)”错误 中说:

    教程中似乎有个建筑风的案例:
    /Tutorials/inco

    好的谢谢 我试一下

  • SHM尖锐直角边界层添加

    OpenFOAM
    0 赞同
    4 帖子
    92 浏览

    @xpqiu 感谢大佬回复,问题已解决 请问大佬有遇到过多重overset计算错误的问题吗

CFD中文网2016-2023 | 京ICP备15017992号-2