Skip to content
  • 代数多重网格法(AMG)问题请教

    Algorithm
    2
    0 赞同
    2 帖子
    3k 浏览

    自问自答一下
    C-AMG和SA-AMG的区别在于粗化方案不同以及插值矩阵的构造是不一样的,C-AMG基于细网格生成粗网格,可以理解为更精细,因此计算收敛性会更好;SA-AMG根据聚类的思想来的,所以在setup阶段会大大缩减时间,但是收敛性不如C-AMG;而两者都比较鲁棒性,因此并不是只针对某种方法有效。
    而后一个问题目前调研到的是,AMG在商软中均可以独自使用或者作为预条件使用;而Fluent、COMSOL、STAR-CCM+这三个软件的官网都有对应的介绍,基本上作为预条件都是与CG/BiCGSTAB联用,其他的还有GMRES。
    找了好久也算是把资料找到了

  • LES怎么处理壁面问题

    CFD彩虹条
    2
    0 赞同
    2 帖子
    3k 浏览
    李东岳

    这几天备课OKS,今年大幅度增加湍流的内容。前几天也在朋友圈发了个图。简单说一下。LES处理壁面几种方法:

    LES-RANS混合方程:用RANS求解壁面区域,LES求解主流区。也就是你需要调用DES、DDES等之类

    Wall-modelled LES:对壁面区域进行模化,主流区采用LES。这个进一步还区分不同的模型,比如代数类、PDE类、ODE类等。

    http://www.dyfluid.cn/theory.pdf 混合方程我在这里更新了一点内容。

    另外没理解你说的壁面雷诺数什么意思。

    捕获.JPG

    Ref: Journal of Computational Physics 305 (2016) 589–603

  • snappy加边界层

    Meshy
    4
    0 赞同
    4 帖子
    2k 浏览
    李东岳

    这个基本没有一针见血的办法,shm生成边界层已经是老大难了

  • 气泡合并模型中的系数

    Algorithm
    3
    0 赞同
    3 帖子
    3k 浏览
    WeakFormW

    更新一下,看来这里边常数都是根据模拟结果凑着调整的
    Another critical factor for the simulation of bubble column reactors is the closure model for coalescence and breakage. Various coalescence and breakage models have been examined in recent decades [39–44].

    Prince and Blanch [39] proposed a turbulent coalescence model of bubbles under the assumption of isotropic turbulence, which has been widely applied to bubble columns or stirrer vessels. In spite of its popularity, researchers concluded that the predicted coalescence rate in the Prince and Blanch model is overestimated for turbulence-induced coalescence [24,41–44]. Hence, to align the simulation results with experimental data, either coalescence or breakage rate is adjusted by multiplying with an arbitrary constant. However, there is still no general agreement in the literatures on the value of this multiplier and even contradictory conclusions were drawn. Chen et al increased the breakage rate by a factor of 10 and maintained the coalescence rate constant in their simulations [41]. Olmos et al [42] multiplied a constant of 0.075 for both the coalescence and breakage rates to fit the experimental data. As such, the justification of the coalescence or breakage model is weakened and the predictive capability of the numerical simulation is lost. Moreover, it is troublesome and time-consuming to tune the multiplier via the trial-and-error procedures. While the investigation of coalescence model will not be the focus of the present work, we improve the coalescence model with the slip velocity proposed by Bhole et al [44] and apply this model to the CFD-PBM simulations of a rectangular bubble column.气泡塔反应器模拟的另一个关键因素是聚结和破裂的闭合模型。近几十年来,人们研究了各种聚结和破裂模型 [39–44]。
    Prince 和 Blanch [39] 在各向同性湍流假设下提出了一种气泡湍流聚结模型,该模型已广泛应用于气泡塔或搅拌器容器。尽管该模型很受欢迎,但研究人员得出结论,Prince 和 Blanch 模型中预测的湍流诱导聚结的聚结率被高估了 [24,41–44]。因此,为了使模拟结果与实验数据一致,聚结率或破裂率都通过乘以任意常数来调整。然而,文献中对这个乘数的值仍然没有普遍的共识,甚至得出了相互矛盾的结论。陈等人在他们的模拟中将破裂率提高了 10 倍,并保持聚结率恒定 [41]。 Olmos 等人 [42] 将聚结和破碎率乘以 0.075 的常数以拟合实验数据。这样一来,聚结或破碎模型的合理性就被削弱了,数值模拟的预测能力也丧失了。此外,通过反复试验的过程调整乘数既麻烦又耗时。虽然聚结模型的研究不是本研究的重点,但我们改进了 Bhole 等人 [44] 提出的滑移速度聚结模型,并将该模型应用于矩形气泡塔的 CFD-PBM 模拟。

  • 0 赞同
    1 帖子
    1k 浏览
    P

    文章题目为《The inviscid impingement of a jet with arbitrary velocity profile》,主要是分析了法向二维冲击射流流场的解析解,但是个人比较疑惑的是为什么控制方程只包括了压力的泊松方程(如下)而没有包括二维的欧拉方程?
    d64e75d4-5884-42a4-b208-20c90325cc29-image.png

    全文见附件,希望有朋友可以解答。The inviscid impingement of a jet with arbitrary velocity profile.pdf

  • 0 赞同
    1 帖子
    1k 浏览
    东方白杨

    我把边界条件都改了一遍,还是有这个问题,做的是一个通道流动,用的是buoyantBoussinesqSimpleFoam,给了进口速度1.82m/s,出口大气压,流体温度20K,壁面300K,边界条件如下:
    速度:

    FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (1.82 0 0); boundaryField { inlet { type fixedValue; value uniform (1.82 0 0); } outlet { type zeroGradient; } top { type noSlip; } ground { type noSlip; } FrontandBack { type noSlip; } Trib { type noSlip; } }

    压力:

    boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0;//大气压出口 } top { type zeroGradient; } ground { type zeroGradient; } FrontandBack { type zeroGradient; } Trib { type zeroGradient; } }

    参考压力:

    FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { "(inlet|outlet|top|ground|FrontandBack|Trib)" { type fixedFluxPressure; value $internalField; } }

    一直报错如下:

    --> FOAM FATAL ERROR: Continuity error cannot be removed by adjusting the outflow. Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow. Total flux : 0.999375 Specified mass inflow : 0.0167107 Specified mass outflow : 0 Adjustable mass outflow : 0 From function bool Foam::adjustPhi(Foam::surfaceScalarField&, const volVectorField&, Foam::volScalarField&) in file cfdTools/general/adjustPhi/adjustPhi.C at line 111. FOAM exiting
  • 湍流分散力和随机游走模型区别??

    Algorithm
    6
    0 赞同
    6 帖子
    7k 浏览
    李东岳

    可以这么理解

  • 气泡受力分析

    OpenFOAM
    7
    0 赞同
    7 帖子
    6k 浏览
    Y

    @AppleKiller 哥们也做气泡拉格朗日么,用啥求解器啊,我目前用DPMFoam,添加个气泡相关方程,头痛要死

  • 0 赞同
    6 帖子
    2k 浏览
    李东岳

    要把下面那个平面也要包含进去

  • parafoam不可用

    OpenFOAM
    1
    0 赞同
    1 帖子
    1k 浏览
    A

    3193a72d-1f4f-42f5-8b05-08f530279fab-image.png
    想看拉格朗日粒子paraview不能看,parafoam就这样,我用的是ubuntu22.04.4和of11,该怎么编译parafoam?请教下

  • 大深宽比微通道结构网格绘制

    Meshy
    2
    0 赞同
    2 帖子
    1k 浏览
    李东岳

    这种没法搞,注定就是深宽比很大

  • paraview提取部分结构

    OpenFOAM
    3
    0 赞同
    3 帖子
    3k 浏览
    Z

    @Do1975 感谢大佬,我去试试

  • 请教各种系数的定义及其作用

    Algorithm
    13
    0 赞同
    13 帖子
    21k 浏览
    Y

    @lyc 我也是菜鸡,这玩意我师兄去年也想过,一直没成功,今年我做课题时候要添加自己颗粒属性方程可能接触一点。但是如何实现我也在研究中,你可以看看我之前和其他大佬的一些疑问,他们指出过在哪些文件里面进行属性添加修改,不过都不全,导致我目前现在想要修改内容也还差一点

  • 动网格

    Fluent
    5
    0 赞同
    5 帖子
    5k 浏览
    C

    @Eva-Wong 我不知道你要研究的具体是什么问题,如果严格按照贴着壁面这个物理情景,那里面涉及到的流体问题是什么呢?你的计算域应该是壁面(毛巾)和壁面(桌子)之间的充满了流体的空间吧。

  • 生物质流化床气化

    Fluent
    15
    0 赞同
    15 帖子
    16k 浏览

    @J 我们课题还是偏这个软件,谢谢您的回复

  • STAR-CCM+网格转换后,checkMesh报错

    Meshy
    2
    0 赞同
    2 帖子
    3k 浏览
    I

    1.我猜是symmetry这个边界条件,试一下改成wall,然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查,外部软件的网格节点编号和blockMesh的逻辑不一样。
    2.可能是前后的empty边界问题,你可以在starccm里生成3D网格后,重新在openfoam里extrude来生成前后empty边界
    类似这样
    constructFrom patch;
    sourceCase "../w3-d_hc1-3";
    sourcePatches (symFront);

    // If construct from patch: patch to use for back (can be same as sourcePatch)
    exposedPatchName symBack;

    // Flip surface normals before usage. Valid only for extrude from surface or
    // patch.
    flipNormals false;

    //- Linear extrusion in point-normal direction
    extrudeModel linearNormal;

    nLayers 1;

    expansionRatio 1.0;

    linearNormalCoeffs
    {
    thickness 1;
    }

    // Do front and back need to be merged? Usually only makes sense for 360
    // degree wedges.
    mergeFaces false; //true;

    // Merge small edges. Fraction of bounding box.
    mergeTol 0;

  • snappy画船体网格

    Meshy
    4
    0 赞同
    4 帖子
    2k 浏览
    李东岳

    是的

  • 0 赞同
    2 帖子
    3k 浏览
    bestucanB

    可以找已有的代码对比一下,除非特别新的东西,常用的东西一般都有模板,还有各种代码的模板。

    http://www.pudn.com/Download/item/id/3745439.html
    https://github.com/Singh-T/TDMA-2D

    这些八成是对的,把他们的代码一点一点的替换成你的代码,看看替换到哪就不对了:142:

  • 0 赞同
    2 帖子
    3k 浏览

    或者是每个cell的最大压力值也可以

  • .pvtu文件ghostLevel=1

    Algorithm
    4
    0 赞同
    4 帖子
    5k 浏览
    李东岳

    感谢分享感谢分享