每页:
搜索

All posts by Nancy Bannach

多孔介质中的热平衡与热非平衡传热

2020年 3月 3日

由于具有适用性强、低成本和特殊的热性能等特点,多孔材料的应用范围越来越广泛。例如,因具有优异的机械和热性能,泡沫材料越来越多地用于不同的航空应用。在电动汽车所用的电池中也发现了多孔结构。我们甚至在自然界中发现了无数多孔材料,例如土壤,岩石和木材。当使用它们时,我们会利用它们的热性能。多孔材料的许多工业应用均要求它们具有优良的热性能。 微观层面的热传递 让我们从微观层面上仔细研究多孔结构中的热传递。正如我们之前的博客文章中所讨论的,我们使用这些发现来验证和理解宏观层面的流动方程。在该示例中,流动是等温的,因此我们不研究孔隙几何结构在热传输中的意义。由于流体的热性能可能与固体的性质显著不同,因此,它们之间的相互作用对于理解热传递如何在多孔介质中工作至关重要。   一个冷却的多孔结构的温度演变过程。初始局部不平衡会随时间达到热平衡。 使用与之前的博客文章中相同的示例,并注入比多孔基质热得多的流体。我们观察到,多孔基质 T\textrm{s} 和流体的温度 T\textrm{f} 最初不相同,并且随着时间变化逐渐达到平衡。当然,这取决于边界条件以及流体和固体的热性质。在许多应用中,该假设 T\textrm{s}=T\textrm{f} 是有效的,我们称其为(局部)热平衡;而在其他应用中,该假设是无效的,我们称其为(局部)热非平衡T\textrm{s}\neq T\textrm{f} “局部”是指温度T\textrm{f} 和 T\textrm{s}的逐点比较。 热平衡下的热传递 在局部热平衡假设下,我们只需要一个方程来描述整个(固体和流体)多孔结构的平均温度。基于能量守恒以及应用混合规则,热传递方程式可以表达为 (1) \left(\rho Cp\right)\textrm{eff}\frac{\partial T}{\partial t}+\rho\textrm{f} C{p,\textrm{f}} \mathbf{u}\cdot\nabla T+\nabla\cdot(-k\textrm{eff}\nabla T)=Q 显然,这与众所周知的传热方程式非常相似。流体和多孔介质的热特性被组合为有效特性,即有效体积热容和有效导热率。 \left(\rho Cp\right)\textrm{eff}=\theta\textrm{p}\rho\textrm{s}C{p,\textrm{s}}+\theta\textrm{f}\rho\textrm{f}C{p,\textrm{f}} 式中,指数\textrm{f}和 \textrm{s}分别代表流体和固体, \rho 是密度 Cp 为恒压下的热容,\theta\textrm{s} 为固体体积分数。假设为完全饱和的多孔介质,孔隙率将对应于流体体积分数 \theta\textrm{f} =1-\theta\textrm{s}。 对于热传导,有效的热导率 k\textrm{eff} 取决于多孔介质的结构以及固体和流体的热导率。该软件提供了三个选项来计算有效导热率k\textrm{eff} : 体积平均值,代表与热通量平行的固体和流体条纹 k\textrm{eff}=\theta\textrm{s} k\textrm{s} + \theta\textrm{f} k\textrm{f} 倒数平均值,垂直于热通量的固体和流体条纹 \frac{1}{k\textrm{eff}}=\frac{\theta\textrm{s}}{k\textrm{s}} + \frac{ \theta\textrm{f}}{k\textrm{f}} 幂定律,对于固体和流体具有相似热导率的随机几何k\textrm{eff}=k\textrm{s}^{\theta\textrm{p} }\cdot k\textrm{f}^{\theta\textrm{f}} 下面,我们通过使用多孔材料的人工示例来说明这三种平均技术,并将不同选项的结果与计算值进行比较。 使用不同有效导热系数选项计算的平均温度比较。从左至右:固体(灰色)和液体(蓝色)材料分别以并行,平行和格纹的方式排列。热通量是由上下边界之间的规定温差产生的。 上述图示表明,结构越精细,由倒数平均值和幂定律计算的近似值越好。真正的有效导热系数介于体积平均值和倒数平均值之间,根据 混合物的规则,分别对应上限和下限。如果对流是主要作用,那么混合规则对热导率的作用就不那么重要了。 多孔材料也可以由几种固体和不流动的流体组成。例如,由不同矿物质和截留的液体组成的岩石。也可以在模型中考虑这一点,有效材料属性也相应地被计算。例如,可根据k\textrm{s}=\displaystyle\sum{i=1}\theta{\textrm{p}i}k{\textrm{s}i}来计算由 i 种不同材料组成的多孔基质的体积平均导热率。 热分散 热分散是与多孔微结构有关的另一个重要作用。通常,对于以对流为主的状态,流体在孔隙尺度上遵循旋涡状路径,从而增强了固相和液相之间的热交换。宏观上,这是通过对热传递方程式 (等式1), k\textrm{disp}=\rho\textrm{f} C{p,\textrm{f}} D{ij}有贡献的附加热导率来描述的, 其中D_{ij} 是由于快速速度场而导致的分散张量。 我们将计算结果与之前博客文章中示例的平均温度结果相比较。下图显示了通过微尺度方法计算出的平均温度,以及从具有和不具有热分散的平均宏观方程式中获得的值。 微观和宏观方法的平均温度比较。包括热分散时,会更好地匹配。 热非平衡下的热传递 如本篇博客文章开头所述,局部热平衡并不总是能达到的。特别的,对于快速的非等温流动,较短的时间尺度,或在强烈依赖于其他影响(例如相变)的情况下,固体和流体温度之间的差异可能很大。此时 等式 1 并不完全有效,必须分别考虑各相的能量平衡,并且必须以显式方式考虑两相之间的热交换。这是通过两个温度模型完成的。局部热非平衡方法 方法 […]

Up 和 Down 算子助力薄结构的分析

2016年 1月 12日

模拟含有薄结构的复杂几何时,计算量会相当大,因为为了解析薄结构需要大量的网格单元。COMSOL Multiphysics 提供了专门的特征来模拟薄结构,从而对这类模型实现高效求解,同时保持较高的精度。为建立薄结构并执行后处理,COMSOL Multiphysics 还提供了专用算子,帮助您考虑获得精确结果所需的所有相关参数。

蒸发冷却模拟简介

2014年 12月 8日

当您想到蒸发时,很可能想的是办公桌上正散发着咖啡或茶香的杯子。实际上,蒸发也是从气象学到食品加工等许多工业和科学应用的环节之一。本篇博文是蒸发冷却模拟系列博文中的第一篇。这里,我们将以您的咖啡杯为例来介绍其中的基本概念。


第一页
上一页
1–3 of 3
下一页
最后一页
浏览 COMSOL 博客