首页/文章/ 详情

【Fluent案例13】大气污染物扩散数值模拟

1月前浏览727

概述


       

随着国民经济的高速发展,城市规模日益扩大,城市布局日益复杂,大量高大建筑物的出现势必对空气的流动产生显著影响,进而影响大气污染物的迁移扩散。

本案例参考Wang XIN设计的风洞实验来进行CFD模拟,计算模型为与风洞测试段尺寸相似的长方体,具体尺寸为:长3.4m,宽1.9m,高1.2m,烟囱采用内径2mm单管模拟,大气污染物采用乙烷替代,模型建筑物和排放源相对位置、坐标系布置和浓度测量位置见下图。

模型、采集位置和坐标系布局图

图中参考高度Href=0.1 m,L、W和H分别表示模型建筑物的长、宽和高,本案例L=W=H=0.1m,乙烷烟囱高度为0.5Href=0.05m,距模型建筑物迎风面距离为2Href=0.2m,取H处风速为参考风速U=2.66 m/s。

参考文献:Evaluation of CFD Simulation using RANS Turbulence Models for Building Effects on Pollutant Dispersion

1几何建模及网格划分

几何建模在ANSYS SCDM 中,网格划分在Fluent Meshing中划分多面体网格,并在烟囱以及建筑物尾迹进行网格加密,详细步骤不赘述。


2求解设置及后处理

  • 瞬态计算,Y负方向重力方向

  • 开启SST K-omega湍流模型

  • 开启组分输运模型,考虑进口扩散作用及组分扩散过程中引起的焓输运作用

  • 添加乙烷,并设置多组分材料

  • 设置烟囱入口边界条件

  • 设置风洞入口边界条件

  • 方法一:UDF速度分布

 #include "udf.h"

 #define Uref 2.66

 /* 入口风速分布 */

 DEFINE_PROFILE(in_wind_velocity,t,i)

 {

    real y, x[ND_ND]; 

    face_t f;


    begin_f_loop(f,t)

      {

         F_CENTROID(x,f,t);

         y = x[1];

         

F_PROFILE(f,t,i) = Uref*pow(y/0.1,0.256);

         

      }

    end_f_loop(f,t)

 }


 
  • 方法二:表达式

  • 其他边界条件:风洞出口采用压力出口条件、侧面和顶部采用对称条件,其他为壁面wall边界

  • 初始化,求解计算10s


残差曲线

坐标点(0、-500、1501)压力监测

截面速度分布

截面乙烷质量分数分布

来源:CFD仿真区
MeshingFluent MeshingFluent组分输运UDF湍流建筑材料ANSYS
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-10-13
最近编辑:1月前
濮小川CFD
硕士 心不唤物,物不至!
获赞 13粉丝 38文章 103课程 0
点赞
收藏
作者推荐

【CFD小贴士】CFD入门基础

1 何谓CFD?CFD是计算流体动力学(Computational fluid dynamics)的缩写,是预测流体流动、传热传质、化学反应及其他相关物理现象的一门学科。CFD一般要通过数值方法求解以下的控制方程组:质量守恒方程、动量守恒方程、能量守恒方程、组分守恒方程、体积力等CFD 分析一般应用在以下阶段:概念设计、产品的详细设计、发现问题、改进设计CFD分析是测试和实验的补充,减少了实验和数据采集所需的全部精力和成本2 如何开展CFD模拟?ANSYS CFD 求解器是基于有限体积法计算域离散化为一系列控制体积在这些控制体上求解质量、动量、能量、组分等的通用守恒方程偏微分方程组离散化为代数方程组用数值方法求解代数方程组以获取流场解例:将管道流动的流体区域离散为一系列的控制体* FLUENT控制体以单元cell为中心(即它们与网格直接对应),而CFX控制体以节点node为中心3 CFD模拟基本流程1 确定模拟目标我们希望得到什么样的结果(例如,压降,流量),如何使用这些结果?模拟有哪些选择?分析应该包括哪些物理模型(例如,湍流,压缩性,辐射)?需要做哪些假设和简化?能做哪些假设和简化(如对称、周期性)?需要自定义数学模型吗?–FLUENT使用UDF,CFX使用 User FORTRAN计算精度要求到什么级别?获取模拟结果的时间?CFD是否是合适的工具?2 确定计算域如何把一个完成的物理系统分割出来?计算域的起始和结束位置在这些位置能获得边界条件吗?这些边界条件类型合适吗?能把边界延伸到有合适数据的位置吗?可以简化为二维或者轴对称问题吗?3 创建几何模型如何得到流体域的几何模型?使用现有的CAD模型。从固体域中抽取出流体域?直接创建流体计算域几何模型?CAD几何可以简化吗?去除可能引起复杂网格的不必要特征(倒角、焊点等)使用对称或周期性?流场和边界条件是否都是对称或周期性的?需要切分模型以获得边界条件或者创建域吗?4 计算域网格化计算域的各个部分都需要哪种程度的网格分别率(密度)?网格必须能捕捉感兴趣的几何特征以及关注的变量梯度,如速度梯度、压力梯度、温度梯度等。可以预估出大梯度的位置吗?需要使用自适应网格来捕捉大梯度吗?哪种类型的网格是最合适的?几何的复杂度如何?四边形/六面体网格?三角形/四面体网格?多面体或马赛克 网格?需要使用非一致边界条件吗?计算机资源是否足够?需要多少个单元/节点?需要使用多少个物理模型?5 设置求解模型对给定的问题,需要定义材料属性 流体?固体?混合物?选择合适的物理模型 湍流,燃烧,多相流等。指定操作条件指定边界条件提供初始值设置求解器控制参数设置监测收敛参数* 对于复杂的问题,先求解一个简化的或2D的问题将在短时间内对模型选取和求解设置提供宝贵的经验6 求解计算通过迭代求解这些离散的守恒方程直至收敛以下情况达到收敛:两次迭代的流场结果差异小到可以忽略。监测残差趋势能帮助理解这个差异达到全局守恒:全局量的平衡感兴趣的量(如阻力、压降)达到稳定值:监测感兴趣量的变化收敛解的精度和以下因素有关:合适的物理模型,模型的精度网格密度,网格无关性数值误差7 查看模拟结果查看结果,提取有用的数据使用可视化的工具能回答以下问题:什么是全局的流动类型?是否有分离?激波、剪切层等在哪儿出现?关键的流动特征是否捕捉住了?数值报告工具能给出以下量化结果:力、动量平均换热系数面积分、体积分量通量平衡8 考虑对模型修订这些物理模型是否合适?流动是湍流的吗?流动是非稳态的吗?是否有压缩性效应?是否有三维效应?这些边界条件是否合适?计算域是否足够大?边界条件是否合适?边界值是否是合理的?网格是否是足够的?加密网格能否提高精度?网格是否有无关性?是否需要提高网格捕捉几何的细节?Fluent中常用的物理模型流动和传热动量、质量、能量方程辐射湍流雷诺平均模型 (Spalart-Allmaras, k–ε, k–ω, 雷诺应力模型)大涡模拟 (LES) 和分离涡模拟 (DES)组分输运体积反应Arrhenius 有限速率化学反应湍流快速化学反应:涡耗散, 非预混, 预混,局部预混 湍流有限速率反应:EDC, laminar flamelet, composition PDF transport,表面化学反应多相流模型:离散相模型 (DPM),VOF:混合模型,欧拉-欧拉/欧拉-粒子,液/固与空化相变模型 动网格Moving zones:Single and multiple reference frames (MRF),Mixing plane model,Sliding mesh modelMoving and deforming (dynamic) mesh (MDM)用户定义标量输运方程一些附加模型。如电池模型来源:CFD仿真区

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈