FLUENT根据自定义条件停止计算_Fluent_UDF

FLUENT根据自定义条件停止计算

12小时前浏览4

正文共： 1249字 4图预计阅读时间： 4分钟

1 前言

有时候，我们在做仿真计算时，可能想当关注的某个变量达到预设值时，就自动停止计算而不用人盯着。今天我们用以前的一个高压氢气充注案例来演示一下上述需求的实现方法。

2 方法描述

第一步，假设我们想在充注口压力达到22atm时就停止计算，编写如下的UDF，在每次迭代结束后计算充注口压力，判断是否达到目标值，如果达到了，就将一个变量“interruptflag”置1。

#include "udf.h"#define INLET_ZONE_ID 23#define pressure_threshold 22e5real current_pressure = 0.0;//calculating inlet pressureDEFINE_EXECUTE_AT_END(check_pressure_stop){    real sum_p = 0.0, sum_area = 0.0;    Domain *domain=Get_Domain(1);    Thread *thread = Lookup_Thread(domain, INLET_ZONE_ID);    real area[ND_ND];    face_t f;    begin_f_loop(f, thread)    if PRINCIPAL_FACE_P(f, thread)    {            F_AREA(area, f, thread);            sum_p += F_P(f, thread) * NV_MAG(area);            sum_area += NV_MAG(area);    }    end_f_loop(f, thread)    sum_p = PRF_GRSUM1(sum_p);    sum_area = PRF_GRSUM1(sum_area);    current_pressure = sum_p / sum_area;  // calculating average pressure    node_to_host_real_1(current_pressure);    #if !RP_NODE    Message("the inlet pressure is %f\n",current_pressure);    if(current_pressure>=pressure_threshold)    {        RP_Set_Integer("interruptflag",1);        Message("the inlet pressure threshold has reached \n");  }    else    RP_Set_Integer("interruptflag",0);    #endif}

第二步，在FLUENT TUI命令声明上述的标志变量“interruptflag”。

(rp-var-define 'interruptflag 0 'interger #f)

第三步，在Execute commands创建如下的指令语句，每个时间步长都执行，当标志变量interruptflag大于0时，停止计算

(if (> (%rpgetvar 'interruptflag) 0)(set! mstop? #t))

测试一下结果，可以看到，当入口压力达到22atm时，计算自动停止，并且就算后续手动再计算也会停止。

来源：仿真与工程

FLUENT大型保温箱结构热桥模拟

正文共： 789字 8图预计阅读时间： 2分钟1 前言本案例来自于现场实践，结构设计是大型保温箱的核心内容，保温材料自身是难以作为承载结构的。有如下图的一种做法，外部制作金属结构，具有足够的承载能力，然后在金属结构内表面安装保温材料，但是保温材料怎么固定在结构上就很重要了。下图的做法是，通过金属销钉将保温材料压在结构内部。但是这必然有个问题，由于金属销钉穿过了保温厚度，因此将产生热桥，对于保温是不利的。今天，我们来讨论一下这个问题。2 建模与网格创建上图所示的二维保温箱截面模型，计算域包含保温材料和承载结构。划分四边形结构化网格，节点数11167，最小正交质量1.0。3 边界条件与求解设置由于保温材料和承载结构都是固体，因此只需要求解固体导热方程即可。按保温材料和金属材料的导热系数施加，内部高温壁面设置为700℃恒温面，结构外部壁面设置为对流边界，对流传热系数10W/m2℃，环境温度26.85℃。注意，能量方程的收敛标准设置为1e-8，采用默认的无法满足能量平衡，读者务必要report一下质量流、能量流这些，用以判断是否达到收敛。4 计算结果我们看一下温度场，可以看到金属销钉处有明显的热桥，热量通过该位置传到外部结构。也可以从外部壁面的温度曲线看出热桥效应，外壁平均温度54.82℃，平均散热密度279.67W/m2，单独考虑保温材料，保温外壁的稳态温度应为48.56℃，散热密度217.14W/m2，因此热桥的影响范围相当广，增加了约22%的散热，如果金属销钉越多，影响越严重。接下来，考虑在外壳的蒙板处贴一层保温材料，结构梁不贴。同样的边界条件，温度分布如下，可以看出蒙板处的温度基本无问题，但是结构梁的温度上升显著，这个很好理解，热肯定往热阻小的地方跑。此时平均散热密度246.54W/m2，相比纯保温材料增加了13.5%的散热，可见结构热桥的影响还是非常大的。可以预测，将结构梁也贴保温后，效果会进一步改善。来源：仿真与工程