Skip to content
  • 0 赞同
    6 帖子
    8k 浏览

    @青山 我做的不涉及三相,所以只有两相。然后你这个说只能在alphat里算感觉没道理啊,那如果是层流,这个输入都没有,就不允许闪蒸了嘛,应该不会把

  • 0 赞同
    1 帖子
    3k 浏览

    最近看到一种说法,level set方程是简化版的相场方程,结合这三种方法在方程形式上非常相似,不知道从计算精度、计算量等方面来讲,这三个方法的优缺点有没有什么对比性结论?

  • 0 赞同
    8 帖子
    12k 浏览
    李东岳

    如果气含率比较大呢,比如0.3甚至更大呢,压力30 kPa和300 kPa感觉如果考虑压缩性体积差别很大,气含率这些差别应该也很大?

    我觉得你更应该从实验或者文献上获取一下信息,因为有的需要考虑,有的不需要考虑。不太普适性。

    我那个考虑了PBE,比较麻烦,你可以不考虑PBE,直接植入一个跟压力相关的线性粒径函数,比如$d=f(p)$,非常简单

    不过多谢引用我那个文章,现在引用数低了点拉低了h因子,帮我引用我很高兴 :papa:

  • 0 赞同
    9 帖子
    12k 浏览
    R

    @Tens
    您好,我参考您的代码也在compressibleInterFoam中添加了组分方程,目前出现了某种组分的质量分数超过1的情况,这个问题困扰我很久了,可以麻烦您帮我看一下代码有什么问题吗?万分感谢!

    tmp<fv::convectionScheme<scalar>> mvConvection ( fv::convectionScheme<scalar>::New//convectionScheme:Abstract base class for convection schemes. ( mesh, fields, alpha2rho2phi, mesh.divScheme("div(alpha2rho2phi,Yi)") ) ); { forAll(Y, i) { if (i != inertIndex && composition.active(i)) { volScalarField& Yi = Y[i]; fvScalarMatrix YEqn ( fvm::ddt(alpha2rho2, Yi) + mvConvection->fvmDiv(alpha2rho2phi, Yi) - fvm::laplacian(turbulence.muEff()*alpha2, Yi) == alpha2*parcels.SYi(i, Yi) + fvOptions(rho2, Yi) ); YEqn.relax(); fvOptions.constrain(YEqn); YEqn.solve(mesh.solver("Yi")); fvOptions.correct(Yi); Yi.max(0.0); Yt += Yi; } Y[inertIndex] = scalar(1) - Yt; Info<<"info::inertIndex.................="<< inertIndex <<endl; Info<<"info::Y[inertIndex].................="<< max(Y[inertIndex]) <<endl; Y[inertIndex].max(0.0);

    我在createField.H中声明了下列变量:

    volScalarField alpha2rho2 ( IOobject ( "alpha2rho2", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), alpha2*rho2 ); surfaceScalarField alpha2rho2phi ( IOobject ( "alpha2rho2phi", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), linearInterpolate(alpha2*rho2*U) & mesh.Sf() );
  • 太贵了这本书...

    CFD彩虹条
    4
    0 赞同
    4 帖子
    5k 浏览
    bestucanB

    @东岳 嗷,我以为您的70g激光打印大杀器已经服役了ಥ‿ಥ

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

    不太清楚 静候大神

  • 几本有用的CFD书

    CFD彩虹条
    1
    0 赞同
    1 帖子
    3k 浏览
    李东岳
    Finite volume methods for hyperbolic problems

    https://www.jianguoyun.com/p/DaAmei4Q9s3ZBhjBxKwB

    Thermo-fluid dynamics of two-phase flow

    https://www.jianguoyun.com/p/DaEadzUQ9s3ZBhjFxKwB

  • Jasak又发飙了

    CFD彩虹条
    2
    0 赞同
    2 帖子
    3k 浏览
    Y

    现在wsl谢谢代码调试好,ssh上超算还是很开心的

  • 0 赞同
    1 帖子
    2k 浏览
    C
    均匀各向同性湍流,其雷诺应力矩阵非对角项等于0 竖向非均匀各向异性湍流,比如槽道流(水平向均匀),只有$ \left \langle u_{1} u_{2} \right \rangle $有值,其他非对角项均为0。Pope的书中也给出了解释。 现在问题是,均匀各向异性湍流,其雷诺应力矩阵非对角项是否为0?如果参照Pope书中类似推导,在空间3个方向都是均匀的,理论上$ \left \langle u_{1} u_{2} \right \rangle $,$ \left \langle u_{1} u_{3} \right \rangle $,$ \left \langle u_{2} u_{3} \right \rangle $都为0吗?

    《Pope, S.B., Turbulent Flows. Vol. 12. 2000: Cambridge University Press.》89页
    e0bada21-7ac0-4ae4-acc3-3319a012dfc7-image.png

  • 【分享】 Residual 动态监测

    OpenFOAM
    9
    0 赞同
    9 帖子
    13k 浏览

    好人一生平安

  • 固体求解器添加两种材料

    OpenFOAM
    1
    0 赞同
    1 帖子
    2k 浏览
    洱聿

    各位大佬们好,我想在fsifoam流固耦合求解器中添加两种固体材料属性,请问有人做过相关的算例嘛?是否可以直接添加?还是需要修改求解器呢?

  • snappyHexMesh画边界层

    Meshy
    3
    0 赞同
    3 帖子
    7k 浏览
    L

    @李东岳 李老师,就是感觉有时候,不能满足通过理论计算得到的边界层信息,比如不能满足边界层厚度啊,层数也达不到,这种一般怎么处理啊?感觉workbench的mesh虽然可以,但是我希望通过这个方式画网格,这样对我后续的计算方便一点。💋

  • 0 赞同
    9 帖子
    12k 浏览
    L

    @东岳 感谢东岳老师 目前我还在完善过程中 这个模型目前还存在相当多的问题 这个模型主要是集中在工程应用上 他的大部分数据是利用单管试验来获取的 所以是一个靶向非常集中的宏观模型 对其他工艺的适用性可能非常一般 再次感谢李老师的提议 我想在更进一步完善这个模型之后 :xinxin: :xinxin: 对一些研究多孔过滤的人才能真正起到一定的意义

  • 0 赞同
    2 帖子
    4k 浏览
    疏影横斜水清浅

    @疏影横斜水清浅 自问自答吧,是组合公式
    1f3e206f-d460-45ee-9db8-ffeb3cf340df-image.png
    C(k,i)

  • 泊肃叶管道流动入口角点处理

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

    我没太看过相关的讨论。目前这种只有周期性边界才不会发生这种问题。你说的这个问题再往深了我就没关注过了,遗憾,老铁整的挺深啊:146:

  • 办公桌设计收藏

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

    https://huaban.com/boards/15971106/

    新买个房子,打算打造个看文献写代码舒适的办公环境! :xiezuoye:

    替代文字

    上面这个灯搞一个

    替代文字

    背板可以用冲孔板

    替代文字

    这个看起来不够立体,应该在屏幕上方加上一些元素

    替代文字

    搁板下面可以加个灯带

    https://www.zhihu.com/question/30550212

  • 0 赞同
    3 帖子
    5k 浏览
    田畔的风

    @李东岳 哇,感谢李老师:140:

  • meshing网格划分出现多余接触区域

    Meshy
    1
    0 赞同
    1 帖子
    4k 浏览
    A

    请问大家为什么我用meshing划分流体域网格时,明明两个面没有重叠,模型树上还会出现图片上这种接触区域呢?![6392L5}X9$J}C{GTK_ZZ7F.jpg

  • 有前辈用过sHM进行开尔文波区加密吗?

    Meshy
    4
    0 赞同
    4 帖子
    9k 浏览
    B

    请教一下你是每一层加密都对三个方向进行加密吗,shm中能能不能只对z方向进行加密呢。@一朵苍凉

  • 内部场引用编程问题

    OpenFOAM
    2
    0 赞同
    2 帖子
    3k 浏览
    C

    请教各位前辈一个简单的编程问题。我先用topoSet在网格中创建内部cellZone,名称为interior。现在想对interior每个网格体心的U赋值,代码如下:

    label patchID = mesh.cellZones().findZoneID("interior"); const labelList cells = mesh.cellZones()[patchID]; forAll(cells, i) { label celli = cells[i]; U[celli] = vector(0.1, 0.1, 0.1); }

    如果是边界,一般可以写成引用的方式:

    label patchID = mesh.boundaryMesh().findPatchID("inlet"); vectorField& patchField = U.boundaryFieldRef()[patchID];

    以上方式是可以实现赋值,想问的问题是:

    (1)内部cellZone是否可以类似写成引用的形式,再对patchField赋值?类似边界场写成以下形式会报错,但不知道怎么修改:

    label patchID = mesh.cellZones().findZoneID("interior"); vectorField& patchField = U.primitiveFieldRef()[patchID];

    报错如下:

    error: invalid initialization of reference of type 'Foam::vectorField&' {aka 'Foam::Field<Foam::Vector<double> >&'} from expression of type 'Foam::Vector<double>' 180 | vectorField& patchField = U.primitiveFieldRef()[patchID]; |

    (2)U[celli]和U.primitiveField()[celli]这两者有无区别?测试效果好像是一样的

    forAll(cells, i) { label celli = cells[i]; U[celli] = vector(0.1, 0.1, 0.1); U.primitiveField()[celli] = vector(0.1, 0.1, 0.1); }