首页/文章/ 详情

【CFD小贴士】Fluent Profile定义锂电池热源

1月前浏览325
引子

由于新能源汽车行业的快速成长,锂电池的仿真也变得日益重要起来。锂电池对热高度敏感,过高的热量很容易引发事故。电池在放电过程中,热源往往不是固定的,有时候在不同的行驶工况下,放热曲线非常复杂,这时候采用一般的函数(表达式功能),已经很难来定义这样复杂的热源。因此我们可以采用profile来定义这种复杂的热源。

问题描述

此案例为一简单演示案例,并不代表真实放热过程,只是为了演示这种点式插值型热源的定义。此案例为一18650单体电芯,热源变化不规律。

1.几何与网格

此案例采用18650电池,在SCDM中创建几何并划分结构网格


2.求解器设置

将网格文件导入Fluent,并检查网格质量

选择压力基瞬态求解器

3.打开能量方程
4.读入profile文件

此案例采用.csv格式读入,具体如下:

CSV文件格式及用法的示例说明

此文为在EXCEL编辑,保存为CSV文件,也可如上图所示在文本文件编辑,然后保存为csv文件。

需要注意的是,profile中的数据采用的是线性插值,因此计算过程中采用的时间步长需要低于profile数据时间步长,否则会造成数据信息丢失,尤其是在曲线有曲率存在的位置。
如下图所示,蓝色线条为大时间步长下得到的值,可以看到在曲率位置并未严格遵守profile指定的规律。

在 Fluent 中读取csv格式 profile 文件

5.定义电芯的热源

电芯的介质属性以及边界条件都采用默认设置

6.关闭流动方程

因为此案例没有建立流体域,因此只计算固体域的热计算

7.监测电芯的最高温度
8.初始化,并求解计算
9.结果确认
采用profile文件定义热源时,锂电池最大温度随时间的响应过程

来源:CFD仿真区
Fluent汽车新能源
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-10-13
最近编辑:1月前
濮小川CFD
硕士 心不唤物,物不至!
获赞 16粉丝 43文章 103课程 0
点赞
收藏
作者推荐

【CFD小贴士】求解壁面采用什么方法?第一层网格高度如何确定?

求解壁面采用什么方法?在近壁区域,解的梯度非常高,但近壁区域的精确计算对模拟的成功至关重要,两种近壁面湍流建模方法:A) 使用 Wall Functions(壁面函数法)B) 解析Viscous Sublayer(粘性子层解析法)1、壁面函数法壁面函数利用如图所示的可预测无量纲边界层剖面,从邻近网格单元质心的条件(速度、壁面距离)确定壁面条件(如剪切应力)这意味着单元应该位于log-layer§第一层网格应该位于 log-layer, 通常应该有一个y+值,比如30 < y+ < 300这是一个非常普遍的准则,但并不是一个绝对的规则,上限是雷诺数的函数(device Reynolds number)对于非常高的Re,如果仍然在log层,y+可能会更高(例如,一艘大型船舶,Re≈ 109,y+值大于1000是安全的)对于非常低(但仍然湍流)的Re(例如一个小泵),对数层可能只延伸到y+≈300,此时y+ = 30太粗糙,边界层上不能有足够数量的网格单元(壁面函数法不再适用)如果发生这种情况,请考虑根据粘性子层分解准则来re-meshing要完全求解重要壁面上的边界层问题,可以尝试在边界层上设置至少10层网格单元–一般来说,如果你更感兴趣的是区域中心的湍流效应,而不是壁上的力,这是合适的2、解析Viscous Sublayer(粘性子层解析法)–第一层网格单元应该在 y+≈ 1 ,并且保证足够多的边界层,增长率适中,不高于1.2-这条准则确保网格将能够充分解析粘性子层的梯度–这将大大增加网格数–一般来说,如果壁面上的力或传热是模拟的关键(气动阻力、叶轮机械叶片性能、传热),解析粘性子层就是你要采用的方法,大多数情况下推荐的湍流模型是SST k-w•使用壁面函数时,垂直于壁面的节点更少(与使用网格求解粘性子层相比)。那么第一层网格高度如何确定呢?•在前处理阶段,您需要知道网格单元的第一层(边界层)的合适大小,以便Y+位于所需的范围内•实际的流场 直到您计算出求解结果后才会知道(实际上,有时由于计算出的Y+值,返回并重新构建模型是不可避免的)•为了降低 re-mesh 的风险,我们希望在开始时通过手工计算来尝试预测单元格大小:•对于平板流,雷诺数 给定 ReL= 1.4x106•从y+的定义开始:•目标y+值和流体性质已知,需要Ut ,它被定义为:•壁面剪切应力,tw ,可以由表面摩擦系数得到,Cf:通过文献,平板表面摩擦系数的计算公式如下:•Re是已知的,因此使用这些定义来计算第一层网格的高度•我们设定的目标Y+值为 50, 因此:我们的第一层网格高度y应该大约是1毫米。如果你想要y+ = 1,只要用上面公式中的1替换50,得到y = 1.8x10-5m上述公式适用于平板流,对于内部流动, 雷诺数根据管子的水力直径计算,只要将上述公式中的摩擦系数公式替换为Cf =0.079 Red-0.25 即可。来源:CFD仿真区

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