OpenFOAM

4.5k 主题 26.2k 帖子
  • 5 帖子
    1k 浏览

    @bestucan 万分感谢!

  • twoPhaseEulerFoam中的IATE模型系数

    34 帖子
    8k 浏览

    @李东岳 对的,大概有十几公里,那谢谢东岳老师,我再使劲努努力钻研一下。个人觉得根据文献植入算法太难了,虽然能看懂离散和求解过程,但是想在openfoam里面实现觉得难度挺大的,当时老师您带我们植入的单相可压缩的,目前这是两相可压缩,我努力努力,再次感谢东岳老师的持续关心:xinxin:

  • 自适应网格问题

    5 帖子
    1k 浏览

    问题已解决,勾选paraView中的“Use VYkPolyhedron”选项后交错就会消失

  • openfoam的wmake的语法

    7 帖子
    2k 浏览

    @bestucan 我看了目前用到这个 -I. 的地方,包括reactingParcelFoam跟sprayFoam,都有子文件夹,比如sprayFoam下的egineFoam,应该跟compile子求解器有关

  • openfoam中的沉积

    6 帖子
    2k 浏览

    @李东岳 东岳老师是想诱惑我再报一次您的辅导班:141:

  • 3 帖子
    950 浏览

    @cresendo 我也是看了亿遍~~~~~

  • 颗粒见黏结以及颗粒形状

    3 帖子
    1k 浏览

    在opemfoam中有一个内聚力,可以使粒子之间相互黏结。它是通过粒子之间的碰撞的能量转换实现的。其中通过一个叫cohesionEnergyDensity的系数来决定它的强度。这个密度越大,黏结越明显。希望我说的这个是你想问的。

  • 文章中的网格独立性验证

    7 帖子
    3k 浏览

    @史浩 VOF和LES都是高度网格依赖类求解器。主要是LES只能捕获大于网格尺度的涡,因此,网格越小,LES越能捕获到更小的涡。类似的,VOF捕获到的界面只能大于网格尺度,如果计算域存在气泡,如果网格越小,VOF可以捕获到的气泡越小。

  • 界面温度超过饱和温度,相变异常

    7 帖子
    2k 浏览

    @李东岳 实在是有点困难~请问老师如何让两个场(如图片中的a和vDotAlpha)对应网格的值相乘,这不是内积也不是外积,用forAll循环的话,会报assignment of read-only location错误。

    forAll(a,celli) { vDotcAlphal[celli] = a[celli]*vDotcAlphal0[celli]; vDotvAlphal[celli] = a[celli]*vDotvAlphal0[celli]; }

    图片2.png

  • 传热计算中物性参数分段线性插值问题

    7 帖子
    2k 浏览

    :mihu: 想请教一下关于热物性(如$\rho,\mu,k,c_p$)对温度压力比较敏感的问题,在计算中往往会十分容易发散,这个发散的原因本质上是由于什么呢,体现在程序中是容易造成哪一块发散。而且同样都可以使用例如buoyantPimpleFoam求解器,这种与高速可压缩流体计算的区别在哪呢?在计算中是否可以采用一些特殊的方法提高程序的收敛及稳定性呢。

  • 求问如何给标量场赋值为实数?

    14 帖子
    5k 浏览

    @hongjiewang 我直接好家伙,这都两年了居然还回复我。
    由于种种原因我现在已经不怎么玩OF了,只能给你一些比较通用的建议。首先看这一堆问号大概是segment fault错误,那么我首先猜是产生了异常值导致求解器爆炸,log函数最容易出现的问题就是传入0值和负值,建议你首先检查数值问题。如果不是数值问题,那就单步调试之类的慢慢debug看看怎么回事,你的两个变量是否正确初始化了等等

  • 1 帖子
    447 浏览

    麻烦帮忙改正一下。想要规定thetaBoundary为一个常数值,然后对theta(相分数)进行对比,总是实现不了。谢谢。

    const Field<scalar>& theta();
    const Field<scalar>& thetaBoundary();
    thetaBoundary = {0.5};
    if (theta < thetaBoundary)

  • 对称结构计算结果不对称原因

    2 帖子
    753 浏览

    算出来的等值面结果如图
    2bfde5ee-55f9-497e-9109-0ab0d48a0789-image.png

  • 分区网格的一小段代码

    2 帖子
    933 浏览

    Hi,东岳大佬,请教下openfoam里可以实现分区网格计算,比如两个分区,并设置interface boundary condition来coupling这两个区域吗?如果可以,麻烦给点思路,谢谢。
    我暂时想到的只有fluid-solid interaction,即FSI toolbox in Foam-extend,改写solid domain为fluid-fluid interaction,来实现,不过有些繁琐了。

    BR,
    Shiyu

  • 在Win10-WSL2-Ubuntu20.04上

    5 帖子
    2k 浏览

    @李东岳 谢谢东岳老师,我还是在win下安装paraview吧~

  • CFDEM 非牛顿流

    9 帖子
    2k 浏览

    @lys 用的已有模型,专门针对非牛顿流体的曳力模型还没用过,希望能进一步交流啊

  • interIbFoam 湍流模型

    1 帖子
    514 浏览

    有用过Foam-extend 4.0 里的interIbFoam的吗?这个求解器是不是不能加湍流模型啊?一加湍流模型边界条件就报错

  • 3 帖子
    906 浏览

    谢谢李老师,确实是alphatWallFunction出了问题。重新修改了这个函数后,就可以了。再次感谢。:xinxin:

  • 关于网格内particle的平均速度求解

    4 帖子
    1k 浏览

    @zhe 因为最后一步你重新给p.U()速度了啊。。

  • 涡量如何在代码中引用?

    8 帖子
    2k 浏览

    @荷兰后腰阿贾克斯 嗯嗯,谢谢老哥了:xinxin3:

  • The current version of twophaseeulerfoam

    7 帖子
    1k 浏览

    Thanks a lot. I wonder if the partial elimination algorithm mentioned in Passalaqcua's paper was implemented in the current version?

  • 指定cellzone区域添加动量源项的问题

    4 帖子
    2k 浏览

    @我是河滩指定cellzone区域添加动量源项的问题 中说:

    遍历了整个网格区域进行搜索

    前辈您好,您是如何遍历整个网格区域并获取网格坐标的呢,是利用mesh.C().component()语句吗:141:

  • compressibleInterFoam计算气泡溃灭

    5 帖子
    2k 浏览

    已解决:使用foam-extend4.1版本中的compressibleInterFoam求解器,模型修改为球对称模型,计算正确。
    P.S. foam-extend4.1中的compressibleInterFoam没有耦合温度方程,计算过程中没有出现上述问题。
    参考文献:Numerical modeling of laser generated cavitation bubbles with the finite volume and volume of fluid method, using OpenFOAM

  • 构造通量(phi)方法之间的区别

    5 帖子
    1k 浏览

    @李东岳 感谢老师关心!已经在家啦,趁着有空补补基础hhh

  • μ在颗粒碰撞中是否为切向回弹系数?

    1 帖子
    543 浏览

    最近对两个概念非常混淆,即摩擦系数和切线回弹系数,这两个系数是否就是同一个概念?

    比如在以下的颗粒pairCollision参数设置中,法向碰撞的回弹系数根据alpha设定。但是切向的回弹系数,是不是就是这个mu?

    pairSpringSliderDashpotCoeffs { useEquivalentSize no; alpha 0.01; b 1.5; mu 0.1; cohesionEnergyDensity 0; collisionResolutionSteps 1; };
  • 分享:OpenFOAM不同离散格式的简单研究

  • SIMPLER in OpenFOAM

    5 帖子
    1k 浏览

    压力速度迭代感觉国际上这面玩的少了 就刚出那时候搞出来一些SIMPLEC之类的

  • 基于OpenFOAM做软件开发

    7 帖子
    2k 浏览

    @同学博 好的 十分感谢!

  • OpenFOAM中有关类模板使用的一些疑问

    1 帖子
    572 浏览

    各位朋友们好!

    我最近在看lagrangian库中的solidParticle库代码,对于里面的 类模板的使用 不太懂,solidParticleCloud.C部分代码如下:

    void Foam::solidParticleCloud::move(const dimensionedVector& g) //在solidParticleCloud类中定义一个move函数 { const volScalarField& rho = mesh_.lookupObject<const volScalarField>("rho"); const volVectorField& U = mesh_.lookupObject<const volVectorField>("U"); const volScalarField& nu = mesh_.lookupObject<const volScalarField>("nu"); interpolationCellPoint<scalar> rhoInterp(rho); interpolationCellPoint<vector> UInterp(U); interpolationCellPoint<scalar> nuInterp(nu); solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); //在solidParticleCloud类下面定义了一个trackingData类的对象td?? Cloud<solidParticle>::move(td, mesh_.time().deltaTValue()); }

    代码中的这句话:

    solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value());

    目前自己很困惑,为什么在solidParticleCloud类的函数move中可以 定义一个 solidParticle类中的 trackingData类的 对象td呢?

    自己查找了solidParticleCloud.H文件发现 solidParticleCloud的基类是 Cloud<solidParticle>类,这个Cloud<solidParticle>类 中solidParticle是实参,对Cloud <ParticleType> 类模板进行了实例化.

    我怀疑难道是因为Cloud<solidParticle>类中的 实参是solidParticle,所以其子类solidParticleCloud中函数 就可以调用/定义 solidParticle类中的trackingData类?

    请哪位高手能够指点一下呢?感激不尽:xinxin:

  • 6 帖子
    2k 浏览

    @沐风 您好! 我最近刚开始学习OF。关于您的 文献1的validation,我也想试一下,请问在什么标准求解器上修改比较好?想要应用boussinesq 假设,但看buoyantBoussinesqPimpleFoam 求解器描述为湍流的。但该算例的自然对流Ra=6*10^5,即为层流。希望您指教一下,非常感谢!!

  • 缩尺模型发散

    3 帖子
    769 浏览

    感觉是不是fluent有什么黑科技,在fluent里能算,在of里发散了。。

  • 多孔介质相变 发散

    2 帖子
    859 浏览

    补充一下~算例用的求解器是自己在原求解器基础上修改的,这个求解器已经可以确定没有语句或者逻辑上的问题了。之前用这个求解器已经正确跑过算例了,结果也是和物理相符合的。只是前两天给虚拟机扩容,操作失败,然后所有的算例文件还没来得及保存。这两天新设的这个算例总是发散,我觉得是边界条件的问题,但是又不太懂,所以还是得麻烦各位老师~

  • openfoam的壁面设置

    1 帖子
    591 浏览

    大家好,最近我在用双流体模型做一个沉积案例,其中一相由于重力会沉积到壁面,并且附着在上边,想请教一下各位大神,如何设置壁面类型或者某一相的参数,才能出现附着的效果。谢谢!!!

  • 如何插值让一个point量等于cell center

    5 帖子
    1k 浏览

    @李东岳 谢谢李老师,我试过mesh(),还是过不了。这相当于一种以cell为计算单位的简化。在一些粒子数量级较大的情况下,可以应用。

  • 30 帖子
    6k 浏览

    @veen 好的。等你的教程,嘿嘿

  • 如何找到mesh_.nCells()的源文件具体位置?

    5 帖子
    1k 浏览

    @xpqiu 谢谢邱老师的热心指点!我明白了!再次感谢:xinxin:

  • 什么是referred particles

    1 帖子
    393 浏览
    //- Interactions between real and referred (off processor) particles void realReferredInteraction();

    请问,在pairCollision中提到的referred particles 具体指什么?谢谢。

  • 8 帖子
    2k 浏览

    @jyang
    在decomposeParDict文件最后,加入:
    constraints
    {
    constraint1
    {
    type preservePatches;
    patches (region1_to_region2 region2_to_region1);
    }
    }

  • OpenFOAM如何设置specific flux

    7 帖子
    2k 浏览

    @bestucan 和论文原作者确认了,说我这样是等价的,调整了一下其他地方,目前算出来对上了。。

  • 3 帖子
    1k 浏览

    @xpqiu 对,nutW2是分配了空间了的,大小就是cellID的大小。
    造成这个的原因我怀疑是不是因为的我边界有转折导致的。。我是提取的roof这个边界上对应的第一层网格的值的。
    121.png
    因为当时我做计算的时候除以了nutW1之后计算报错了,我才发现nutW1[0]=0,不能被除。然后才改成使用了nutW2的

  • trackPart ParcelType::trackingData::在哪里被定义?

    2 帖子
    779 浏览

    OpenFOAM有个快捷命令

    src

    进入源代码所在文件夹,然后

    grep -R 'ParcelType::trackingData' ./

    就把所有文件抓住来了,如果找定义,只看.H文件

  • 流体网格的更新问题

    13 帖子
    2k 浏览

    preCICE很amazing,用好了应该比自己改求解器简单高效 https://www.bilibili.com/video/BV1iT4y1P7e7/
    (现在发帖竟然需要间隔300s:jingya: )

  • 湍流模型问题请教:

    17 帖子
    3k 浏览

    @西湖冷月 借帖求教大佬,狭窄裂缝中的湍流用哪种方法模拟好呢?

  • OpenFOAM面心数据插值到面单元节点

    5 帖子
    1k 浏览

    :xinxin2:

  • 关于绝热状态方程熵的计算过程

    1 帖子
    690 浏览

    OpenFOAM v2006中,adiabaticPerfectFluid状态方程如下:

    $\rho (p)=\rho_{ref}(\frac{p+B}{p_{ref}+B})^{\frac{1}{n_{T}}}$

    根据熵的计算公式:

    $dS=c_{p}\frac{dT}{T}-(\frac{\partial v}{\partial T})_{p}dp$

    求解器中将这一公式分为两部分进行积分计算,其中第一部分$c_{p}\frac{dT}{T}$ 的积分计算在thermophysicalModels/specie/thermo/hConst/hConstThermoI.H中实现,即$c_{p}*log(\frac{T}{T_{std}})$
    但第二部分的积分,由于公式是绝热状态的,没有温度作为参数,想请教一下官方求解器中为什么会进行以下处理:

    通过源代码公式的求导反推,官方求解器在计算第二部分积分时,在状态方程的分母上直接添加了温度T,如下式:
    CodeCogsEqn (1).gif
    公式源代码详见adiabaticPerfectFluidI.H文件的S函数:

    template<class Specie> inline Foam::scalar Foam::adiabaticperfectFluid<specie>::S ( scalar p, scalar T ) const { scalar n=1-1.0/gamma_; return -pow(p0_+B_,1.0/gamma_)*(pow((p+B_),n)-pow((Pstd+B_),n)) /(rho0_*T*n); }

    除上述绝热状态方程外,线性变化的状态方程linear也采用了同样的处理方式,想知道为什么可以这样处理?
    另外一个问题是,如果采用上述处理之后,对于气体的$c_{p}-c_{v}$,根据公式:
    CodeCogsEqn (2).gif 该如何处理?

  • Neck vortex,中文翻译

    7 帖子
    1k 浏览

    那创新性不是拉满了,厉害了!:146:

  • 最大液相分数不到1,甚至只有0.02

    2 帖子
    733 浏览

    做一个点线图,横轴是孔隙率,纵轴是最大液相分数,看看到哪里突变了。二分法找这个突变点,应该很快。

    出图的时候用同一标尺(最大值最小值相等),这样不同图之间也有可比性。看看那个刚好导致突变点的孔隙率的case,液相分数是怎么消失的,是跑出边界了还是在流场内被计算过程吃了

  • EQBMM编译小问题

    3 帖子
    873 浏览

    不知道你对QBMM这面的研究进展怎样了:135:

  • 请教一个wallFunction问题

    3 帖子
    916 浏览

    @Samuel-Tu :
    谢兄之玉言,弟牢记

    然兄亦暗示:
    有的参数需调整,以保证物理意义一致。

    另外,如兄及各位大学士不弃,望再赐一教:
    OF能自动算雷诺数吗?

    感激不尽

  • 如何把三维离散数据转化为tecplot可读的格式

    8 帖子
    2k 浏览

    @bestucan 对 用tecplot成功了 tecplot里 reverse distance的插值方法可以把无序的数据插值到有序的网格上 不过有几百组数据 还得弄个宏什么的