Fluent

872 主题 4.1k 帖子
  • 关于Fluent中壁面润滑力的UDF

    2 帖子
    2k 浏览

    这个问题 我自己想明白了,这个方向向量不需要指定,该udf使用DEFINE_SOURCE,在指定source时同时指定udf作用的相和方向即可。

  • 关于fluent中壁面吸收组分的问题

    4 帖子
    2k 浏览

    哪位有过类似项目案例经验的朋友,可否赐教一二

  • 关于VOF液面上升的问题

    4 帖子
    2k 浏览

    @东岳 是呀,想模拟一个静水面,如何才能让他静止不动?(不让下方压力入口顶着进入)

  • 请问UDF给定速度型入口问题

    5 帖子
    3k 浏览

    @逆变张量 是的速度入口不合适,给的总压

  • 关于Fluent瞬态计算残差的设置

    5 帖子
    5k 浏览

    @l-j刘侃
    感谢您的回复,您可以说一下他用的是absolute还是relative么? 我以前都用absolute,然后判据给的是-4。他这个是不是用的是relative格式的准则呢?

  • 3 帖子
    2k 浏览

    首先是CPU 但是CPU的优化性和效果是非常复杂的问题 一般来说 目前的主流 最好最大主频不要低于2.8 同时并行核心数不要低于32 从我的经验48核心比较好
    其次是内存 内存对一些内存要求多的算例的提升效果比较明显 能够占到30%左右
    硬盘一般跟软件面板操作相关 有的时候也是一些时间
    显卡的话 如果你是Fluent可以适用CUDA的方法 提升效果比较明显 两张4W的卡 大概coupled的下 提升30%左右
    最后就是操作系统 你的多核心下 有些操作系统能够提升速度几百倍(显然有些系统对于40以上的核心效果很差) 一般Linux下比同等的windows server要快5%

  • 关于东岳老师朋友圈文献的求助~

    5 帖子
    2k 浏览

    @walten :ok2: :ok2: :ok2: :ok2: :ok2:

  • fluent中标准壁面函数的代码是怎样的呢?

    3 帖子
    2k 浏览

    @云风 主要是对方程和C编程理解较浅,想找个例子对照着改一下,代码中与fluent相关的可能就是fluent中变量的形式,比如说x,y,z是x[0],x[1],x[2]这样的

  • 如何选择scalable wall functions呢?

    2 帖子
    2k 浏览

    scalable wall functions应该跟 Standard Wall Functions差不多,只是多了个,无论边界层网格怎么细化,它都能得出一致的结果

  • tecplot如何提取某一时刻的自由面

    4 帖子
    2k 浏览

    某一时刻是不是你需要有那个时刻的数据啊

  • 利用fluent动网格方法更改网格可行吗

    2 帖子
    1k 浏览

    实际结构垂直距离差了5m,

    要看你整个计算域多大,如果整个计算域50m,偏差10%应该是可以的。不过动网格设置从算法上要改动,为何不直接重做一个网格呢

  • 关于udf载入出错原因

    3 帖子
    2k 浏览

    应该是,udf.h头文件无法读取

  • 关于在哪里设置温度和压强的udf

    4 帖子
    3k 浏览
    adjust init on_demand
  • 才发现Fluent里面密度基全是耦合求解器

    7 帖子
    3k 浏览

    @东岳 对,现在打开了

  • 3 帖子
    2k 浏览

    (2)
    壁面函数的核心是一个简化模型。未必非常准确,需要结合case对比
    (3) 标准做法下,计算结果具有网格无关性,即不影响。 对的

  • FLUENT可以创建新的zone吗?

    3 帖子
    2k 浏览

    可以尝试一下标记一下这个区域,然后separate一下cell zone,但是多数情况下分割出来的区域边界锯齿状比较严重,尤其是非结构网格。

  • fluent模拟异质催化反应

    4 帖子
    2k 浏览

    @zousiyu 感谢,我不做这方面,只是对这个问题进行一个有益的探索

  • CFD计算中的定常与非定常的问题

    3 帖子
    2k 浏览

    @ibelief 问题已经解决了差不多了,应该是由于在进行定常计算作为初始化的时候,这个初始值就是有问题的,因为本人算的是管内的流体的涡旋运动,之前不知道为啥使用simple或者simplec的时候,会出现气相,虽然不是持续存在的,之后改为coupled算法后,定常计算中就在没有出现气相,所以我估计可能跟算法有一定的关系,也许在我的例子里使用coupled算法更加合理些

  • 3 帖子
    2k 浏览

    在fluent中VOF在缺省情况下是只能设置一项为可压缩性,如楼上所言,气体可以设置为可压缩气体,液体的压缩性可以使用Tait方程通过UDF插入fluent中,Ps:这个UDF很简单,可以参考help文件中的UDF案例,里面有一个是关于液相可压缩性的例子

  • 3 帖子
    6k 浏览

    感谢回复!帖子很长,还要消化一下。。但总体感觉当时花在GPU上的预算换在CPU上可能会更合适一点:zoule:

  • 3 帖子
    2k 浏览

    @ibelief 多谢,但是这个是含有epsilon的湍流模型中浮力项对湍动能的影响,那如果是LES方法怎么考虑浮力呢?

  • 2 帖子
    2k 浏览

    没有遇到过,您探索一下,然后总结一下岂不是更好!
    功德无量!

  • 关于三维燃烧室数值模拟中喷嘴的设置

    2 帖子
    2k 浏览

    (1) 计算中inject 的计算时间的长度等
    (2) fluent 手册里面有详细介绍
    (3)算一下速度,根据某经验公式进行计算,fluent 手册中也有提到
    (4) report 出口情况

  • fluent 多节点并行

    2 帖子
    2k 浏览

    写jou.
    设置并行

    比如CSDN 上提供的解决方案

    /define/user-defined/compiled-functions/compile
            libraryName
            yes
            sourcefilename.c

    /define/user-defined/compiled-functions/load
           libraryName

    ;#编译并载入UDF

    /file read-case
           yourcasename.cas
           /file read-data
           yourdatname.dat

    ;#读入调试好的case、data文件

    /define/user-defined/function-hooks/adjust "youradjustname::libraryName"""

    ;#hook像诸如init adjust这些需要hook的udf,不知道为何除了udf里定义的DEFINE_ADJUST(youradjustname,d)外,还会有一个名字""为空的,并且三个"要连写,两者之间不能有空格,不知道代表什么

    /solve/patch

    phase-2

    12

    mp

    1

    ;#我计算的是一个液面晃荡的问题,所以是两相流,要Patch一下,选择phase-2作为我们的domain,然后我要patch的zone的编号是12,然后要空格一下进入下一步,然后选择mp即体积分数,令其为1
           /file/auto-save/root-name yourrootname
           /file/auto-save/data-frequency 10000

    /file/auto-save/case-frequency 10000

    ;#自动保存data,case,每10000个时间步或是10000迭代步
           /solve/iterate
           200000

    ;#如果是稳态问题,迭代200000次

    /solve/set/time-step 0.0001
           /solve/dual-time-iterate 
           50000
           200

    ;#如果是非稳态问题,时间步为0.0001,时间步数50000,每步迭代200次
           /file/write-case-data
           yourfinalname.out
           exit
           yes

    ;#最后导出最后的case与data文件

    然后最后提交任务可以用lsf、pbs啥的,但是我不会,不过直接用命令提交貌似也麻烦不了多少。

    新建一个文件new

    然后文件里面写

    node9

    node9

    node9

    ....

    node9

    这样共n列代表你使用节点9,里面的n个核。

    nohup fluent 3d -t12 -cnf=new -g -i yourjouname.jou >&youroutname.out&

    其实-t12表示一共12个核,-cnf=new使用node9的12核,然后运行过程输出到youroutname.out中,这是使用一个节点的命令,使用多个节点只需对命令作些许修改如下 只需加入-pethernet 和 -ssh

    nohup fluent 3ddp -pethernet -cnf=myhosts -t36 -ssh -g -i f30.jou>&out.out &

    作者:aleigus
    来源:CSDN
    原文:https://blog.csdn.net/aleigus/article/details/77161181
    版权声明:本文为博主原创文章,转载请附上博文链接!

  • 关于fluent中定义溶液的问题

    2 帖子
    2k 浏览

    UDS 是一个方法

  • 这是上述问题的图片

    2 帖子
    2k 浏览

    后来解决的怎么样?

  • 求助冷凝udf

    2 帖子
    2k 浏览

    #include "udf.h"
    /Constants used in psat_h2o to calculate saturation pressure/
    #define PSAT_A 0.01
    #define PSAT_TP 338.15
    #define C_LOOP 8
    #define H2O_PC 22.089E6
    #define H2O_TC 647.286
    /user inputs/
    #define MAX_SPE_EQNS_PRIM 2 /total number of species in primary phase/
    #define index_evap_primary 0 /evaporating species index in primary phase/
    #define prim_index 0 /index of primary phase/
    #define P_OPER 101325 /operating pressure equal to GUI value/
    /end of user inputs/
    //
    / UDF for specifying an interfacial area density /
    /
    /
    double psat_h2o(double tsat)
    /* /
    / Computes saturation pressure of water vapor /
    / as function of temperature /
    / Equation is taken from THERMODYNAMIC PROPERTIES IN SI, /
    / by Reynolds, 1979 /
    / Returns pressure in PASCALS, given temperature in KELVIN /
    {
    int i;
    double var1,sum1,ans1,psat;
    double constants[8]={-7.4192420, 2.97221E-1, -1.155286E-1,
    8.68563E-3, 1.094098E-3, -4.39993E-3, 2.520658E-3, -5.218684E-4};
    / var1 is an expression that is used in the summation loop /
    var1 = PSAT_A(tsat-PSAT_TP);
    /* Compute summation loop /
    i = 0;
    sum1 = 0.0;
    while (i < C_LOOP){
    sum1+=constants[i]pow(var1,i);
    ++i;
    }
    ans1 == sum1(H2O_TC/tsat-1.0);
    / compute exponential to determine result /
    / psat has units of Pascals /
    psat = H2O_PCexp(ans1);
    return psat;
    }
    DEFINE_HET_RXN_RATE(user_evap_condens_react, c, t, hr, mw, yi, rr, rr_t)
    {
    Thread **pt = THREAD_SUB_THREADS(t);
    Thread *tp = pt[0];
    Thread *ts = pt[1];
    int i;
    real concentration_evap_primary, accum = 0., mole_frac_evap_prim,
    concentration_sat ;
    real T_prim = C_T(c,tp); /primary phase (gas) temperature/
    real T_sec = C_T(c,ts); /secondary phase (droplet) temperature/
    real diam = C_PHASE_DIAMETER(c,ts); /secondary phase diameter/
    real D_evap_prim = C_DIFF_EFF(c,tp,index_evap_primary)

    0.7C_MU_T(c,tp)/C_R(c,tp);
    /primary phase species turbulent diffusivity/
    real Re, Sc, Nu, urel, urelx,urely,urelz=0., mass_coeff, area_density,
    flux_evap ;
    if(Data_Valid_P())
    {
    urelx = C_U(c,tp) - C_U(c,ts);
    urely = C_V(c,tp) - C_V(c,ts);
    #if RP_3D
    urelz = C_W(c,tp) - C_W(c,ts);
    #endif
    urel = sqrt(urelxurelx + urelyurely + urelzurelz);
    /relative velocity/
    Re = urel * diam * C_R(c,tp) / C_MU_L(c,tp);
    Sc = C_MU_L(c,tp) / C_R(c,tp) / D_evap_prim ;
    Nu = 2. + 0.6 * pow(Re, 0.5)* pow(Sc, 0.333);
    mass_coeff = Nu * D_evap_prim / diam ;
    for (i=0; i < MAX_SPE_EQNS_PRIM ; i++)
    {
    accum = accum + C_YI(c,tp,i)/mw[i][prim_index];
    }
    mole_frac_evap_prim = C_YI(c,tp,index_evap_primary)
    / mw[index_evap_primary][prim_index] / accum;
    concentration_evap_primary = mole_frac_evap_prim * P_OPER
    / UNIVERSAL_GAS_CONSTANT / T_prim ;
    concentration_sat = psat_h2o(T_sec)/UNIVERSAL_GAS_CONSTANT/T_sec ;
    area_density = 6. * C_VOF(c,ts) / diam ;
    flux_evap = mass_coeff *
    (concentration_sat - concentration_evap_primary) ;
    *rr = area_density * flux_evap ;
    }
    }
  • 2 帖子
    2k 浏览

    图在哪?
    感觉描述的很有意思
    再试试??

  • 环境温度的用户自定义函数从哪里设置呀

    2 帖子
    2k 浏览

    请自己看手册。。。
    祝早日脱离苦海

  • 关于加载用于自定义函数

    2 帖子
    2k 浏览

    可以加载多个

  • 2 帖子
    1k 浏览

    我真的发现,凡是写“求大神指导的”,下面的回复就是这样的

  • 2 帖子
    2k 浏览

    编译方式不一样?

  • 请教大家一个有关反应速率的udf急。

    2 帖子
    2k 浏览

    #include "udf.h"
    DEFINE_EC_RATE(user_ec_rate, f, fthread, r, V, current, didV, Eeq)
    {
    double alpha_a = r->alpha_a, alpha_c = r->alpha_c;
    double io = r->io
    double T = F_T(f,fthread);
    double arg1, arg2;
    cxboolean tafelmethod = r->tafelmethod;
    int i;
    double eta;
    if (tafelmethod)
    {
    alpha_a = 2.303 * UNIVERSAL_GAS_CONSTANT * 298.15 /(alpha_a * FARADAY_CONSTANT);
    alpha_c = 2.303 * UNIVERSAL_GAS_CONSTANT * 298.15 /(alpha_c * FARADAY_CONSTANT);
    }
    Eeq = r->Eeq;
    eta = V - Eeq;
    for(i = 0; i<r->n_reactants; i++)
    if( ABS( r->exp_reactant[i] ) > SMALL_S )
    {
    int ni = r->reactant[i];
    io = pow((F_YI(f,fthread,ni)/MAX(r->yi_ref[ni],SMALL) + 1.0e-20), r->exp_reactant[i]);
    }
    for(i = 0; i<r->n_products; i++)
    if( ABS( r->exp_product[i] ) > SMALL_S )
    {
    int ni = r->product[i];
    io = pow((F_YI(f,fthread,ni)/MAX(r->yi_ref[ni],SMALL) + 1.0e-20), r->exp_product[i]);
    }
    arg1 = FARADAY_CONSTANT / (UNIVERSAL_GAS_CONSTANTT);
    arg2 = arg1eta;
    current = io( exp( arg2alpha_a ) - exp( -arg2alpha_c ) );
    didV = io( arg1alpha_aexp( arg2alpha_a ) + arg1alpha_cexp( -arg2alpha_c ) );
    /* If multiple electrochemical reactions are used, you can define rate for each reaction
    using the following if-statement /
    /
    if (STREQ(r->name, "reaction-1"))
    {
    ...
    }
    else if (STREQ(r->name, "reaction-2"))
    {
    ...
    }
    */
    }

  • 二维双圆柱串联

    2 帖子
    2k 浏览

    Cd 变化,物理上就是这样吧

  • 2 帖子
    2k 浏览

    getdata取点;
    给原作者发个邮件问一下;
    这个图是excel 做的吧

    你可以用excel 也可以用origin

  • 有关ANSYS Fluent和PV panels simulation的问题

    2 帖子
    2k 浏览

    大气温度可以通过udf patch 进去,太阳辐照我推测可以采用边界条件设置,但是Fluent guide里面完全没有提到关于weather file 的东西,请您给些提示吧

  • openfoam,导入,ansys

    2 帖子
    2k 浏览

    进展怎么样?

  • 咨询一个关于fluent并行计算的问题

    2 帖子
    1k 浏览

    MPI 安装都一样吗?

  • 2 帖子
    1k 浏览
  • 2 帖子
    1k 浏览

    你查一下这方面的实验论文吧?
    壁面速度低,容易造成管壁聚集

  • 特征长度如何选取

    2 帖子
    2k 浏览

    雷诺数计算和其定义有关。比如槽道流动中特征长度采用的是槽径,平板流动就是板长,而颗粒雷诺数就是采用颗粒粒径。
    你搜过其他论文里面的选取原则吗?

  • Fluent周期网格模拟有什么限制吗?

    2 帖子
    1k 浏览

    可以吧,fluent给定的周期边界的案例就是类似的

  • fluent 大渦模拟温度变化问题

    2 帖子
    1k 浏览

    我搜索了一下“大神”这个关键词,在本站目前有355个帖子,下面几乎都是和我这条回复一样。

  • 2 帖子
    1k 浏览

    RSM 怎么样? 对雷诺数没有的限制

  • fluent theory guide翻译

    2 帖子
    2k 浏览

    @东岳 在 关于correctPhi.H这个函数 中说:

    哦对了还有fluent理论指南。最近琐事太多了:sad: 。都已经翻译完了。就是格式太难统一了。因为是各自翻译的。用的是word。处理格式比较费时。

    这个可能是最近的进展了

  • Mixture model用UDF定义材料属性?

    2 帖子
    1k 浏览

    DEFINE_PROFILE 和 DEFINE_SOURCE 应该可以实现

  • [求教]关于复合运动的overset设置问题

    2 帖子
    2k 浏览

    你最后是如何实现的呢?可以分享一下吗?

  • 密度基非稳态点火问题

    2 帖子
    2k 浏览

    你物理模型上重新设置一下

  • 导出壁面热流密度profile文件问题

  • 2 帖子
    2k 浏览

    我也没有做个类似的问题,我自己的射流模拟中,这样的压力波动是非常正常的。

    物理过程中是不是也真实存在压力波动呢?