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

    /---------------------------------------------------------------------------\

    \ / F ield OpenFOAM: The Open Source CFD Toolbox \ / O peration Website: https://openfoam.org \ / A nd Version: 6 \/ M anipulation

    *---------------------------------------------------------------------------*/
    Build : 6-fa1285188035
    Exec : reactingFoam
    Date : Dec 24 2024
    Time : 10:45:39
    Host : "DESKTOP-MGCJK2F"
    PID : 3256
    I/O : uncollated
    Case : /home/mywsl/OpenFOAM/mywsl-6/run/turbulent
    nProcs : 1
    sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
    fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
    allowSystemOperations : Allowing user-supplied system call operations

    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Create time

    Create mesh for time = 0

    PIMPLE: No convergence criteria found

    PIMPLE: Operating solver in PISO mode

    Reading thermophysical properties

    Selecting thermodynamics package
    {
    type hePsiThermo;
    mixture reactingMixture;
    transport sutherland;
    thermo janaf;
    energy sensibleEnthalpy;
    equationOfState perfectGas;
    specie specie;
    }

    Selecting chemistryReader foamChemistryReader
    elements not defined in "/home/mywsl/OpenFOAM/mywsl-6/run/turbulent/constant/chem_Li_2007"
    Reading field U

    Reading/calculating face flux field phi

    Creating turbulence model.

    Selecting turbulence model type RAS
    Selecting RAS turbulence model kEpsilon
    kEpsilonCoeffs
    {
    Prt 0.85;
    Cmu 0.09;
    C1 1.44;
    C2 1.92;
    C3 0;
    sigmak 1;
    sigmaEps 1.3;
    }

    Creating reaction model

    Selecting combustion model EDC
    Selecting chemistry solver
    {
    solver ode;
    method standard;
    }

    StandardChemistryModel: Number of species = 9 and reactions = 21
    Selecting ODE solver seulex
    using integrated reaction rate
    Creating field dpdt

    Creating field kinetic energy K

    No MRF models present

    No finite volume options present
    Courant Number mean: 0 max: 0

    Starting time loop

    Courant Number mean: 0 max: 0
    deltaT = 1.111111111e-05
    Time = 1.11111e-05

    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 7.237527843e-08, No Iterations 1
    DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 6.496027629e-08, No Iterations 1
    DILUPBiCGStab: Solving for H, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
    DILUPBiCGStab: Solving for H2, Initial residual = 0.9999999548, Final residual = 1.249524636e-08, No Iterations 1
    DILUPBiCGStab: Solving for O, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
    DILUPBiCGStab: Solving for OH, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
    DILUPBiCGStab: Solving for H2O, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
    DILUPBiCGStab: Solving for O2, Initial residual = 0.9999997617, Final residual = 1.249524395e-08, No Iterations 1
    DILUPBiCGStab: Solving for HO2, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
    DILUPBiCGStab: Solving for H2O2, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
    DILUPBiCGStab: Solving for h, Initial residual = 0.001698464143, Final residual = 6.399191709e-07, No Iterations 1
    min/max(T) = 300, 1800.356388
    DICPCG: Solving for p, Initial residual = 0.8154141575, Final residual = 0.05901504271, No Iterations 4
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 0.1648808522, global = -0.04664958196, cumulative = -0.04664958196
    DICPCG: Solving for p, Initial residual = 0.5937630302, Final residual = 9.1448632e-07, No Iterations 49
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 1.031207163e-06, global = 1.108666077e-07, cumulative = -0.0466494711
    DILUPBiCGStab: Solving for epsilon, Initial residual = 0.7898692099, Final residual = 7.356031497e-07, No Iterations 13
    bounding epsilon, min: -57.68897069 max: 199.8018274 average: 2.853407683
    DILUPBiCGStab: Solving for k, Initial residual = 1, Final residual = 7.436499697e-07, No Iterations 14
    bounding k, min: -53.76693743 max: 137.3801042 average: 1.662759794
    ExecutionTime = 7.11 s ClockTime = 7 s

    Courant Number mean: 0.3263677773 max: 38.86161238
    deltaT = 5.71632726e-08
    Time = 1.11683e-05

    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCGStab: Solving for Ux, Initial residual = 0.01088173554, Final residual = 9.291750018e-07, No Iterations 2
    DILUPBiCGStab: Solving for Uy, Initial residual = 0.126974428, Final residual = 5.413457541e-07, No Iterations 3
    DILUPBiCGStab: Solving for H, Initial residual = 0.05031914993, Final residual = 2.526806224e-07, No Iterations 3
    DILUPBiCGStab: Solving for H2, Initial residual = 0.04988384318, Final residual = 2.536567037e-07, No Iterations 3
    DILUPBiCGStab: Solving for O, Initial residual = 0.05023601744, Final residual = 2.528851275e-07, No Iterations 3
    DILUPBiCGStab: Solving for OH, Initial residual = 0.0498739039, Final residual = 2.53709804e-07, No Iterations 3
    DILUPBiCGStab: Solving for H2O, Initial residual = 0.04983585131, Final residual = 2.537599816e-07, No Iterations 3
    DILUPBiCGStab: Solving for O2, Initial residual = 0.04987406878, Final residual = 2.536836372e-07, No Iterations 3
    DILUPBiCGStab: Solving for HO2, Initial residual = 0.04887920295, Final residual = 2.557918423e-07, No Iterations 3
    DILUPBiCGStab: Solving for H2O2, Initial residual = 0.04849527326, Final residual = 2.565684649e-07, No Iterations 3
    DILUPBiCGStab: Solving for h, Initial residual = 0.05497148043, Final residual = 4.9437153e-07, No Iterations 3
    min/max(T) = 299.9248087, 1796.462561
    DICPCG: Solving for p, Initial residual = 0.01518365283, Final residual = 3.016271454e-08, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 6.870692989e-09, global = -5.333936149e-09, cumulative = -0.04664947643
    DICPCG: Solving for p, Initial residual = 0.0003488105302, Final residual = 1.108179983e-09, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 2.557155586e-10, global = 8.320508309e-11, cumulative = -0.04664947635
    DILUPBiCGStab: Solving for epsilon, Initial residual = 0.1451083224, Final residual = 2.148726213e-07, No Iterations 2
    DILUPBiCGStab: Solving for k, Initial residual = 0.8290663092, Final residual = 5.98405794e-08, No Iterations 4
    ExecutionTime = 8.97 s ClockTime = 9 s

    Courant Number mean: 0.0013035774 max: 0.1228067562
    deltaT = 6.644108124e-08
    Time = 1.12347e-05

    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCGStab: Solving for Ux, Initial residual = 0.1832008417, Final residual = 3.515359245e-07, No Iterations 8
    DILUPBiCGStab: Solving for Uy, Initial residual = 0.8416763107, Final residual = 7.756203776e-07, No Iterations 10
    DILUPBiCGStab: Solving for H, Initial residual = 0.4428609204, Final residual = 9.062138129e-07, No Iterations 9
    DILUPBiCGStab: Solving for H2, Initial residual = 0.4410214426, Final residual = 9.736737868e-07, No Iterations 9
    DILUPBiCGStab: Solving for O, Initial residual = 0.4424908421, Final residual = 9.30398611e-07, No Iterations 9
    DILUPBiCGStab: Solving for OH, Initial residual = 0.4409439717, Final residual = 9.745009796e-07, No Iterations 9
    DILUPBiCGStab: Solving for H2O, Initial residual = 0.4408222366, Final residual = 9.738422502e-07, No Iterations 9
    DILUPBiCGStab: Solving for O2, Initial residual = 0.4409741964, Final residual = 9.739327596e-07, No Iterations 9
    DILUPBiCGStab: Solving for HO2, Initial residual = 0.4368143603, Final residual = 7.794393661e-07, No Iterations 9
    DILUPBiCGStab: Solving for H2O2, Initial residual = 0.4352225269, Final residual = 9.850169941e-07, No Iterations 9
    DILUPBiCGStab: Solving for h, Initial residual = 0.4603404275, Final residual = 8.710686777e-07, No Iterations 10
    min/max(T) = 299.9999982, 1791.283236
    DICPCG: Solving for p, Initial residual = 0.04989128258, Final residual = 1.728078903e-08, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 3.9720677e-09, global = -2.805246459e-09, cumulative = -0.04664947915
    DICPCG: Solving for p, Initial residual = 0.00040507004, Final residual = 2.635572617e-10, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 6.291142477e-11, global = 5.449968679e-11, cumulative = -0.0466494791
    DILUPBiCGStab: Solving for epsilon, Initial residual = 0.2709179929, Final residual = 2.790707667e-07, No Iterations 9
    DILUPBiCGStab: Solving for k, Initial residual = 0.6253639963, Final residual = 3.507848237e-07, No Iterations 10
    ExecutionTime = 10.81 s ClockTime = 11 s

    Courant Number mean: 0.0008942027266 max: 0.05529325678
    deltaT = 7.968158396e-08
    Time = 1.13144e-05

    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCGStab: Solving for Ux, Initial residual = 0.06593954665, Final residual = 3.351747253e-07, No Iterations 8
    DILUPBiCGStab: Solving for Uy, Initial residual = 0.4102152305, Final residual = 9.318726704e-07, No Iterations 10
    DILUPBiCGStab: Solving for H, Initial residual = 0.2892697203, Final residual = 5.28621965e-07, No Iterations 9
    DILUPBiCGStab: Solving for H2, Initial residual = 0.2624968155, Final residual = 4.376897291e-07, No Iterations 9
    DILUPBiCGStab: Solving for O, Initial residual = 0.2808888368, Final residual = 5.574822147e-07, No Iterations 9
    DILUPBiCGStab: Solving for OH, Initial residual = 0.2591248329, Final residual = 4.295608139e-07, No Iterations 9
    DILUPBiCGStab: Solving for H2O, Initial residual = 0.2597666461, Final residual = 4.522225416e-07, No Iterations 9
    DILUPBiCGStab: Solving for O2, Initial residual = 0.2612620687, Final residual = 4.395858975e-07, No Iterations 9
    DILUPBiCGStab: Solving for HO2, Initial residual = 0.2220485989, Final residual = 5.99379148e-07, No Iterations 9
    DILUPBiCGStab: Solving for H2O2, Initial residual = 0.1820866327, Final residual = 6.019580413e-07, No Iterations 10
    DILUPBiCGStab: Solving for h, Initial residual = 0.09990726406, Final residual = 5.829150064e-07, No Iterations 10
    min/max(T) = 299.9999961, 1784.763549
    DICPCG: Solving for p, Initial residual = 0.02262674874, Final residual = 3.398546862e-08, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 8.084807953e-09, global = -4.286910994e-09, cumulative = -0.04664948339
    DICPCG: Solving for p, Initial residual = 0.0001704309813, Final residual = 3.621876029e-10, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 8.790648795e-11, global = 1.74305778e-11, cumulative = -0.04664948337
    DILUPBiCGStab: Solving for epsilon, Initial residual = 0.0897510291, Final residual = 6.434387951e-07, No Iterations 3
    DILUPBiCGStab: Solving for k, Initial residual = 0.2857450205, Final residual = 8.582666035e-07, No Iterations 9
    ExecutionTime = 12.45 s ClockTime = 12 s

    Courant Number mean: 0.001010967894 max: 0.05839075857
    deltaT = 9.556638249e-08
    Time = 1.141e-05

    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCGStab: Solving for Ux, Initial residual = 0.01353903567, Final residual = 4.237016717e-07, No Iterations 6
    DILUPBiCGStab: Solving for Uy, Initial residual = 0.1965833104, Final residual = 4.895900842e-07, No Iterations 7

  • 关于求导精度是否会降阶

    Algorithm
    3
    0 赞同
    3 帖子
    7k 浏览
    Elibathe康E

    用有限差分法得到的二阶精度流场变量和用有限体积法得到的二阶精度流场变量,再分别求涡量场,用两种方法得到的涡量场是否都具有二阶精度或者一阶精度?我不太清楚这其中的原理。你说不会降阶,能否具体说一下原因?

  • IOdictionary.set 语句报错求助

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

    已解决,解决方案是scalarDict.set("vv", "vv [0 0 0 0 0 0 0] -0.1");

  • 0 赞同
    2 帖子
    4k 浏览
    I

    不太懂。请问是2D 还是3D ? 使用了对称轴?
    推测原因:
    1.上下边界距离引起的压差?
    2.湍流模型不行?

  • 关于turbinefoam的网格大小问题

    OpenFOAM
    4
    0 赞同
    4 帖子
    4k 浏览
    C

    @zhou 谢谢,不过这个方法还是挺难搞得,真难搞这东西

  • James Gerald McDonald

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

    嗯 AP 很多AP都很强力 强力AP!

  • fluent一系列云图展示问题

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

    @纵浪大化 哦哦,谢谢您,我知道了!

  • 有关高斯积分积分点的确定

    Algorithm
    3
    0 赞同
    3 帖子
    8k 浏览
    C

    @kasitefy
    非常对。

    另外一个问题是这种加和的数值积分忒慢了

  • 0 赞同
    6 帖子
    10k 浏览
    W

    @李东岳 哈哈,其实是从nodebb的帮助中抄来的:cheeky:

  • 如何理解PDE的弱解?

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

    @random_ran

    不满足微分方程,但是满足积分方程的解就是弱解吧。本来俩方程应该等价的。

    其实我觉得关键单不在于什么是弱解,而是这帮数学家为啥要引入弱解这个概念,我的理解是丫们发现物理上一些有意义的解不严格满足他们推导出来的微分方程,但是似乎又差得不多,然后整吧整吧搞出一个弱解的概念来。

    有点感觉想和稀泥。。。。

  • 服务器和工作站的区别?

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

    我最近也想配一台工作站,导师也不太懂就交给我自己来配了,我之前小打小闹配的第一台机器用的2950x觉得amd挺好的,看到新出的3990x觉得太香了,有没有大佬可以讨论一下具体怎么配置好呀,3990x怎么才能双路并行呀

  • 英语记录

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

    I want to keep my bike close by, so I can fix it

    I want to take off (remove) the bike wheel

    Allen key

    It doesn't fit

    Move over, scoot over, can I get passed, please

    She likes more and more like me everyday

    I rarely (don't often) watch TV

    Why did you set the temperature to 20 degrees

    Dont rush me

    I want to get rid of the extra oil

    I can take some, I will finish what's left

    She has a loose tooth, and its goin to fall out

    She is not fully awake, she is a bit groggy

    She is the one who loves me the most, she loves me the most

    Put on the striped onesie

    Move forwards and backwards

    Her ton nail is broken

    This is her first time seeing it

    This drink is non-alcoholic, you can have a sip

    Take a big step forward

    A loose thread

    Lift my leg and put it down

    Stick out your tongue

    She is lying on her side

    I did it all in one go

    Do you mind if I add her to our chat?

    Please comment on the sentences above

    Look at the lines on your arm. You leaned against something, that left lines on your arm.

    Pull the handle to open the wardrobe

    She just wants to be held up one your shoulder, she just wants to be held close to your chest

    I am afraid of the dark

    When she is sleeping, how does she know who is beside her

    She does not care about anything, she is not easily bothered by things, she has high tolerance. I have a high pain tolerance.

    If someone's ideas can be easily influenced by others,

    She does not want to lean back

    Give yourself a round of applause

    Candy wrapper

    I am meshing the pataya???

    Are u scared of mosquitoes biting you?

    Do you hear me calling you?

    Another spoonful

    How about I weigh you?

  • 关于多块网格并行剖分的负载平衡问题

    Meshy
    4
    0 赞同
    4 帖子
    7k 浏览
    小狗狗

    @李东岳 同一个计算域的网格,有单元的节点信息,单元相连的单元信息,面信息。但是两个不同计算域网格之间没有这种拓扑信息。这种情况下如何并行剖分网格同时保证负载平衡,减少通信量

  • CFD青年成长支持计划(2021)

    公告
    137
    0 赞同
    137 帖子
    165k 浏览
    F

    想问一下李老师这个计划还有吗!!很想要这本湍流模型:mianmo:

  • 0 赞同
    2 帖子
    2k 浏览
    M

    你好,请问解决了吗?

  • 0 赞同
    11 帖子
    13k 浏览
    I

    @Shihang-Chen
    你好,我在做类似的代码。想请教伪时间步植入后,计算结果中出现松弛要怎么处理(计算的周期和实际周期不一致)。我是在动量方程中加入的隐式伪时间步。
    512ceafd-f193-4d59-a045-c2859969922c-9e6a4889e6591786ea13032244e9111.png

  • 0 赞同
    1 帖子
    3k 浏览
    Cp_ZhaoC

    这几天在看有限体积法的书,大佬的论文、《数值传热学》和 《The FVM inCFD OpenFOAM..Matlab》都看了,逛论坛发现都在推荐《Computational Methods for Fluid Dynamics》这本;去SpringerLink一搜正好今年出了新版,网上资源好像还不多,在这里分享一下: 点击下载

  • CFD中文网 定位与规范

    已固定 公告
    101
    0 赞同
    101 帖子
    651k 浏览
    H

    我写 固我在,嘿嘿

  • 0 赞同
    10 帖子
    8k 浏览
    L

    统一回复一下:

    感谢给的热心解答。

    用线提取密度的垂向变化(内波),表面波的probe暂时不适用这个问题。

    有效解决方案:
    方案A:用toposet选取一定区域的网格,输出rho,这就有了xyzrho的数据。然后根据需要的xz坐标,做插值即可。
    方案B:选取全局的网格,输出cuttingPlane上的数据。
    方案C:输出特定等直面的数据。

    方案A灵活,但操作复杂;方案B数据最全,但存储空间比较大;方案C直接,对三维工况不友好。

    这都是后处理的技巧,能获得需要的数据。

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

    不太清楚 静候大神