Skip to content

OpenFOAM

OpenFOAM交流区

4.5k 主题 26.5k 帖子
  • 1 帖子
    293 浏览
    学流体的小明

    话说这个latex代码怎么转义不了
    使用LES计算气泡槽道流,因为是多相流,某一点可能会有气泡经过,解析的雷诺应力$u\prime_i u\prime_j$会受到很大影响,所以我想对 $u\prime_i u\prime_j$ 进行加权处理。
    具体方法是:生成张量场,计算每一个cell每一个瞬时的$\alpha u\prime_i u\prime_j$,$(1-\alpha) u\prime_i u\prime_j$,$\rho u\prime_i u\prime_j$,然后使用fieldAverage进行系综平均。
    我想在controlDict中生成场$\rho u\prime_i u\prime_j$,然后每个时刻计算一下这个场,怎样通过functions{ }操作呢?因为是要每一个时间步都计算以获得正确的值,所以不太想按一定时间间隔输出数据之后再进行后处理。
    当然也可以直接在solver的createFields.H中创建,我就是想找一个更通用的不用修改求解器的方法,直接很简便地在controlDict - functions{ }中实现。
    毕竟functions{ }可以用writeObjects输出

    /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Version: v2012 \\ / A nd | Website: www.openfoam.com \\/ M anipulation | ------------------------------------------------------------------------------- Description Writes out specified objects, e.g. fields, stored on the case database. \*---------------------------------------------------------------------------*/ type writeObjects; libs (utilityFunctionObjects); objects (<object names>); writeControl writeTime; // ************************************************************************* //
  • 6 帖子
    916 浏览
    洱聿

    @洱聿 目前主动停止算例运行,会出error文件(但是之前以为是主动停止报的错误,就一直没有注意过),显示错误如下:

    [manage01:228000] [[60655,0],0] tcp_peer_recv_connect_ack: received different version from [[60655,1],1]: 1.8.8 instead of 1.10.2

    目前集群使用的mpi版本是1.10.2,是因为mpi版本太高,所以foam extend4.0不适用嘛?有fe40专门对应的openmpi版本嘛?

  • 如何获得固体边界上水相部分某一场的数值

    2
    2 帖子
    257 浏览
    李东岳

    我觉得你可以写个判断语句。对整个边界patch做循环,比如边界上alpha>0.5的就是水箱,然后提取这个压力。

  • cyclicAMI边界

    2
    2 帖子
    412 浏览
    李东岳

    把你的boundary文件贴一下看看

  • 时间步长对计算结果的影响

    2
    2 帖子
    501 浏览
    李东岳

    带化学反应的时候存在刚性方程组。所以还不能像是普通流场那样放大时间步。

  • 关于icoReactingMultiphaseInterFoam中的组分方程

    3
    3 帖子
    507 浏览
    李东岳

    这在以前的openfoam里面是出现过这种问题的。我详细的跟踪过twoPhaseEulerFoam的算法演变。你这个求解器应该是ESI那面的吧?可以问问他们什么时候能把alpha、rho搞进去

  • 2 帖子
    437 浏览
    李东岳

    这个肯定是没法提说不出口的。咱们要知道,官方开发openfoam40年了。他们对自己的代码/算法非常自信。除非明显的植入错误。他们不会做类似处理。我比较熟悉基金会那面。类似这种需求,做一个统一化的架构整合,属于技术支持,提出来他们会需要你提供fund,然后他们会处理。

    ESI那面我联系不多。不过这种大概率来看。自己整合一个然后提交给他们会比较好。

  • 关于在compressibleInterFoam中添加组分方程

    4
    4 帖子
    722 浏览
    R

    您好,请问一下您实现在compressibleInterFoam中添加组分方程了嘛:chitang:

  • 2000万网格OpenFOAM算例

    1
    1 帖子
    438 浏览
    李东岳

    2000w.tar.xz

    加特林,2000步,906秒

    20000w.tar.xz

  • 后处理得到momentum thickness

    7
    7 帖子
    1k 浏览
    F

    @luofq-sysu 第一个问题:是的,我取的是流向速度时均值。
    第二个问题,是的,是不同流向位置。我当时做的是一个平板流动,平板上流向的网格113个。150个不是一个必须的值,看你算例的情况跟需求。我的是二维流动,因此没有展向的问题。

  • 如何处理时间步

    3
    3 帖子
    488 浏览
    L

    @李东岳
    好,谢谢李老师。。。
    等到了李老师回复

  • 潜坝冲刷试验模拟求助

    3
    3 帖子
    527 浏览
    A

    @李东岳 这个还没有用过,我尝试一下,谢谢李老师。

  • XiFoam中是否可以增加详细反应机理进行计算

    1
    1 帖子
    279 浏览

    诸位大佬,XiFoam中所使用的燃烧模型,是基于反应变量的湍流燃烧模型。
    b=1:未燃混合气b=0:已燃混合气
    是否可以考虑在其中增加详细的化学反应机理,进行计算。

  • PimpleFoam和PisoFoam对于槽道湍流计算的影响。

    3
    3 帖子
    564 浏览
    2

    @学流体的小明 我的那个倒也没有发散,还在正常的计算,就是感觉要算的时间比pisoFoam要久。

  • 2 帖子
    432 浏览
    bestucanB

    一个是锥子喷射,一个是锥子孔喷射。

    可能是:

    锥子喷射能保证喷着时是锥形分布。

    锥子孔空射,是保证出口形状是锥子孔,至于能不能喷射成锥形,得看喷射压力和喷射流量了。

    大口径小流量,是会单边喷射的,就像洗澡的花洒没水了。。。。:mianmo:

  • reactingFoam修改化学反应机理后遇到了问题

    5
    5 帖子
    1k 浏览
    Y

    又又又遇到了新问题,一步都没算下去,请问这种问题一般出在哪里:xinlei:
    Create time

    Create mesh for time = 0

    PIMPLE: No convergence criteria found

    PIMPLE: Operating solver in transient mode with 1 outer corrector
    PIMPLE: Operating solver in PISO mode

    Using LTS
    Reading thermophysical properties

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

    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 thermophysical transport model

    Selecting thermophysical transport type RAS
    Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
    Creating reaction model

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

    standardChemistryModel: Number of species = 544 and reactions = 2446
    Selecting ODE solver seulex
    Creating field dpdt

    Creating field kinetic energy K

    No MRF models present

    No fvModels present
    No fvConstraints present

    Starting time loop

    Time scales min/max:
    Flow = 9.728569846e-08, 0.0001
    Temperature = 4.494232837e+307, 4.494232837e+307
    Composition = 4.494232837e+307, 4.494232837e+307
    Overall = 9.728569846e-08, 0.0001
    Time = 1e-07

    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.0001092756986, No Iterations 1
    DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.006862564884, No Iterations 1
    DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.0001167898495, No Iterations 1
    #0 Foam::error::printStack(Foam::Ostream&) at ??:?
    #1 Foam::sigFpe::sigHandler(int) at ??:?
    #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3 double Foam::standardChemistryModel<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGasFoam::specie >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) at ??:?
    #4 Foam::combustionModels::EDC::correct() at ??:?
    #5 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
    #6 ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #8 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
    Floating point exception (core dumped)

  • reactingFoam报错

    5
    5 帖子
    636 浏览
    vbcwlV

    @李东岳 谢谢东岳老师

  • 14 帖子
    2k 浏览
    X

    @lrl3512 要不您加一下我qq 2559727310

  • 气泡槽道流的计算,求指点

    11
    11 帖子
    2k 浏览
    学流体的小明

    更新一下进度
    暂时不用可压缩求解器计算了,因为发现当前的网格算不到声波那种程度。不过我倒是实现了向温度方程内添加源项这件事,就是不太确定公式以及编程对不对。

  • 在fsifoam中运行算例报错

    8
    8 帖子
    941 浏览
    洱聿

    目前算例移植到计算机集群上计算,模型是水面上漂浮的平板。大约计算到0.8s时出现这个报错,查找原因是内存分配问题;
    8fda3931-8566-4c1d-96e3-08c9fb3e932f-image.png
    有查看了相关的log文件,文件中给出的报错是如下所示
    d5c39409-328d-471d-be7b-801f479d1faa-image.png
    为什么算例运行中会出现nu不被定义的报错呢?是因为计算过程中,网格出现问题嘛?