Characterization of size effect of natural convection in melting process of phase change material in square cavity

The accelerating effect of natural convection on the melting of phase change material(PCM) has been extensively demonstrated. However, such an influence is directly dependent on the size and shape of domain in which phase change happens, and how to quantitatively describe such an influence is still challenging. On the other hand, the simulation of natural convection process is considerably difficult, involving complex fluid flow in a region changing with time, and is typically not operable in practice. To overcome these obstacles, the present study aims to quantitatively investigate the size effect of natural convection in the melting process of PCM paraffin filled in a square latent heat storage system through experiment and simulation, and ultimately a correlation equation to represent its contribution is proposed. Firstly, the paraffin melting experiment is conducted to validate the two-dimensional finite element model based on the enthalpy method.Subsequently, a comprehensive investigation is performed numerically for various domain sizes. The results show that the melting behavior of paraffin is dominated by the thermal convection. When the melting time exceeds 50 s, a whirlpoor flow caused by natural convection appears in the upper liquid phase region close to the heating wall, and then its influencing range gradually increases to accelerate the melting of paraffin. However, its intensity gradually decreases as the distance between the melting front and the heating wall increases. Besides, it is found that the correlation between the total melting time and the domain size approximately exhibits a power law. When the domain size is less than 2 mm, the accelerating effect of natural convection becomes very weak and can be ignored in practice. Moreover, in order to simplify the complex calculation of natural convection, the equivalent thermal conductivity concept is proposed to include the contribution of natural convection to the total melting time, and an empirical correlation is given for engineering applications.     

多孔介质融化相变自然对流数值模拟

根据体积平均理论,建立多孔介质中融化相变过程自然对流模型。数值模拟了以石英砂颗粒为骨架、冰为填充相的二维融化过程。数值结果表明,在融化区上部发生较强烈的自然对流,融化界面与水平线倾角随融化时间越来越大。     

Lattice Boltzmann Simulation of Mixed Convection Heat Transfer in a Driven Cavity with Non-uniform Heating of the Bottom Wall

The goal of this article is to study numerically the mixed convection in a differentially heated lid-driven cavity with non-uniform heating of the bottom wall. The velocity field is solved by a hybrid scheme with multiple relaxation time Lattice Boltzmann(MRT-LBM) model, while the temperature field is obtained by resolution of the energy balance equation using the finite difference method(FDM). First, the model is checked and validated using data from the literature. Validation of the present results with those available in the literature shows a good agreement.A good efficiency in time simulation is confirmed. Thereafter, the model has been applied to mixed convection in a driven cavity with non-uniform heating wall at the fixed Grashof number Gr = 106. It is found that, the heat transfer is weakened as the Richardson number is augmented. For Gr = 106, we note the appearance of secondary vortices at different positions of the cavity corners.     

铸件在凝固过程中的传热和自然对流计算

论文中以流函数、涡度和温度为基本变量,建立了铸件凝固时的传热和自然对流的数学模型。控制方程采用半隐式有限差分在固定网格上求解。在模型建立过程中,着重讨论了潜热的处理和固液相区(界面)速度条件。对于潜热处理,采用了基于合金平衡相图的比热焓法;对于固液相区(界面)速度条件,提出了一种新的、基于固相率进行修正的处理方法。另外还探讨了将本文提出的解法应用于复杂几何条件的可能性。最后,通过对一个典型例子进行计算,对自然对流对铸件凝固的影响作了深入的讨论。     

格子Boltzmann方法模拟自然对流作用下融化传热过程

建立自然对流作用下融化的格子Boltzmann双分布函数模型,根据非线性对流扩散方程的格子Boltzmann模型理论提出一个新的表征融化温度场的分布函数演化方程,并通过变松弛时间方法处理固液两相变热物性传热问题.应用模型对热传导融化及自然对流融化特别固液变热物的融化过程进行模拟.模拟结果与分析解、经典的关联式结果吻合较好,模型的正确性得到了验证.模拟结果表明,自然对流对融化传热过程有着重要的影响,此外固相热传导也对融化传热、融化速率及固液两相温度分布都有一定影响.     

The model study of water mass and energy exchange between the inland water body and atmosphere

Based on a one-dimensional eddy diffusion model, a model to study the water mass and energy exchange between the water body (such as lake and wetland) and the atmosphere is developed, which takes the phase change process due to the seasonal melting and freezing of water and the convection mixing process of energy caused by temperature stratification into consideration. The model uses enthalpy instead of temperature as predictive variable, which will help to deal with the phase change process and to design an efficient numerical scheme for obtaining the solution more easily. The performance of the model and the rationality of taking convection mixing into the consideration are validated by using observed data of Kinneret Lake in Israel and Lower Two Medicine Lake in Montana State in America. The comparison of model results with observed data indicates that the model presented here is capable of describing the physical process of water mass and energy between the water body (lake and wetland) and atmosphere. Comparison of the result from wetland with shallow and deep lakes under the same forcing conditions shows that the evaporation from wetland is much greater than that from lakes, which accords with the real observation fact and physical mechanism. Supported by the National Natural Science Foundation of China (Grant Nos. 40575043 and 40605024) and KZCX-sw-229     

方腔内相变材料固液相变传热研究

基于相变材料(PCM,phase change material)的相变储能设备具有储能密度高的特点。本文建立了基于相变储能元件伪焓模型的固液相变格子Boltzmann模型,研究了内部管道位置、方腔倾斜角度对PCM融化过程的影响规律。结果表明,在内管道靠近方腔上部时,由于上部界面(固液相变界面或上壁面)对自然对流阻碍作用,使PCM的融化速率减慢。但是,在此时使方腔发生倾斜,会改变管道热流体到上部界面的距离,强化PCM的热质传递过程,使融化加快。     

基于傅里叶红外光谱的SiO_2基相变储湿复合材料结构分析及性能优选预测

以SiO_2为载体材料、脂肪酸为相变材料制备具有相变调温性能与储湿调湿性能的SiO_2基相变储湿复合材料。采用等温吸放湿法和步冷曲线法测试不同脂肪酸用量的SiO_2基相变储湿复合材料相变调温性能与储湿调湿性能。利用傅里叶红外光谱测试SiO_2基相变储湿复合材料的结构组成,分析SiO_2与脂肪酸的嵌合机理。以SiO_2基相变储湿复合材料的傅里叶红外光谱特征吸收峰作为输入层,以SiO_2基相变储湿复合材料的脂肪酸用量、相对湿度52.89%下SiO_2基相变储湿复合材料吸湿平衡含湿量与放湿平衡含湿量的平均值、SiO_2基相变储湿复合材料从30~15℃降温所需的时间作为输出层,以S型激活函数作为隐含层,利用BE神经网络建立结构参数与综合相变储湿性能的SiO_2基相变储湿复合材料性能优选预测模型。结果表明,SiO_2基相变储湿复合材料中SiO_2与脂肪酸仅为物理嵌合,未发生化学作用;当脂肪酸用量0.079 mol时,所制备的SiO_2基相变储湿复合材料具有最优的综合相变储湿性能,即在相对湿度52.89%下的吸湿平衡含湿量为0.132 3 g·g~(-1)、放湿平衡含湿量0.147 5 g·g~(-1)、平衡含湿量的平均值为0.139 9 g·g~(-1),从30~15℃降温所需的时间为1 305 s;SiO_2基相变储湿复合材料的性能优选预测模型吻合性较好,具有较高的精确度,其预测值与实测值的相对误差为-2.07%和2.45%,可以用于优选预测SiO_2基相变储湿复合材料的储湿调湿性能和相变调温性能。     

Instability in Two-Sided Thermocapillary-Buoyancy Convection with Interfacial Phase Change

A new model of two-phase thermocapillary-buoyancy convection with phase change at gas-liquid interface in an enclosed cavity subjected to a horizontal temperature gradient is proposed,rather than the previous onesided model without phase change.We study the onset of multicellular convection and two modes of convective instability,and find four different flow regimes.Their transition map is compared with the non-phase-change condition.Our numerical results show the stabilizing effect of interfacial phase change on the thermocapillarybuoyancy convection.     

使用界面追踪对膜态沸腾的数值研究

提出一种计算带相变的自由界面的数值算法.基于分段线性界面重构(PLIC)的VOF(volume-of-fluid)方法用于追踪自由界面,并对汽液交界面上的相变导致的不连续速度场给出处理方法.此方法容易实施且被证明是有效的.流场的求解使用SIMPLE方法,表面张力使用连续表面力模型(CSF)进行计算.在三维直角坐标系下,模拟了水平壁面上的膜态沸腾,在二维适体坐标系下,模拟了竖直圆头柱体表面的自然对流膜态沸腾.计算结果与理论关系式符合较好.     

双颗粒在溶解条件下沉降的多相流动特性

对牛顿流体内溶解双颗粒在垂直管道中的沉降运动进行了直接数值模拟. 流体运动由守恒方程计算, 密度和黏性的变化考虑流场温度变化的影响, 通过积分粘性应力和压力获得颗粒的受力跟踪颗粒运动, 溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定. 通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合. 对双颗粒在等温流体无溶解条件和非等温流体溶解条件下的沉降过程进行了计算. 结果表明, 在一定雷诺数内, 热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒质量、颗粒表面形态的变化引起了颗粒的横向摆动, 并使颗粒沉降速度发生了变化.     

On the combined heat transfer in the multilayer non-gray porous fibrous insulation

A detailed numerical modeling is performed to investigate coupled heat transfer of natural convection, radiation and conduction in high-temperature multilayer thermal insulation (MTI), which consists of high-porous, non-gray semitransparent fibrous materials and reflective foils. Radiation within fibers, radiation between fibers and the reflective foils, conduction within fibers and convection from the fibers to the surrounding fluid are considered. Macroscopic (porous media) modeling is used to determine velocity, pressure and temperatures fields for fibrous insulation with a random packing geometry under natural convection, whereas the radiative transfer equation (RTE) is used to solve the radiative heat flux for non-gray materials. Key features of the macroscopic model include two-dimensional effects, non-gray radiative exchange, and the relaxation of the local thermodynamic non-equilibrium (LTNE). This model was validated by comparison with experimental data and it was used to investigate natural convection of coupled heat transfer in multilayer insulation, numerical results showed that natural convection is more likely to occur when the heated/cooled rate is low, while natural convection can be ignored in simulating steady-state coupled heat transfer in MTI.     

格子Boltzmann方法模拟融化相变过程

建立融化和自然对流耦合的格子Boltzmann双分布函数模型,采用焓法迭代求解相变非线性源项,并通过变松弛时间方法处理固液两相变热物性间题.热传导融化过程的计算值与分析解的对比分析说明该模型能够准确地模拟融化过程.自然对流条件下融化过程的模拟结果表明自然对流在一定程度上影响了融化传热及融化速率等,体现出与热传导融化不同...     

一种模拟大密度比多相流的混合算法

在扩散界面法(diffuse interface method,DIM)的基础上提出一种能够处理大密度比(large density ratio)的多相流混合算法.流场信息通过格子波尔兹曼方法(lattice Boltzmann method,LBM)获得;相界面通过直接求解Cahn-Hilliard(C-H)方程确定.为保证在大密度比情况下求解界面方程的稳定性,采用二阶迎风格式来离散方程的对流项.通过对Rayleigh-Taylor(R-T)不稳定、液体中的气泡上升及液滴撞击干燥壁面的数值模拟,验证了方法的可行性.     

Modeling of Ice-Water Phase Change in Horizontal Annulus Using Modified Enthalpy Method

Phase change in ice-water systems in the geometry of horizontal cylindrical annulus with constant inner wall temperature and adiabatic outer wall is modeled with an enthalpy-based mixture model. Solidification and melting phenomena under different temperature conditions are analyzed through a sequence of numerical calculations. In the case of freezing of water, the importance of convection and conduction as well as the influence of cold pipe temperature on time for the complete solidification is examined. As for the case of melting of ice, the influence of the inner pipe wall temperature on the shape of the ice-water interface, the flow and temperature fields in the liquid, the heat transfer coefficients and the rate of melting are analyzed. The results of numerical calculations point to good qualitative agreement with the available experimental and other numerical results.     

考虑Soret和Dufour效应的方腔内双扩散自然对流格子Boltzmann模拟

采用格子Boltzmann方法,考虑Soret和Dufour效应,对内置高浓度发热圆的方腔内部双扩散自然对流现象进行数值模拟.高浓度发热圆位于方腔中心,四周壁面均为低温低浓度.在该模型中,用三个独立的LBGK方程分别模拟速度场、温度场和浓度场,并通过Boussinesq近似将它们耦合起来.分析Soret数和Dufour数对方腔内部双扩散自然对流的影响,得到流线图、等温线图、等浓度线图、发热圆表面平均Nusselt数和平均Sherwood数.结果表明：Soret和Dufour效应对方腔内双扩散自然对流影响明显,不能忽略.     

潜热型热功能流体强化换热分析

本文建立了分析带有相变微胶囊的潜热型热功能流体的流动和换热过程数理模型,应用有限差分法和移动热源法进行联合求解。计算结果表明,相变微胶囊的加入:较好地提高了流体的换热性能。获得了相变颗粒大小和体积分数对强化换热影响等结果。计算结果对该类流体的设计和应用提供了理论依据。     

溶解与热对流对固体颗粒运动影响的直接数值模拟

对牛顿流体内溶解与热对流对单颗粒在垂直管道中的沉降运动进行了直接数值模拟.流体运动由守恒方程计算,密度和黏性的变化考虑流场温度变化的影响,通过积分黏性应力和压力获得颗粒的受力跟踪颗粒运动,溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定.通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合.分别模拟了颗粒在等温流体、热流体、冷流体及颗粒溶解四种情况下的沉降过程.结果表明,在一定雷诺数内,热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒表面形态的变化引起了颗粒的横向摆动,并使颗粒沉降速度发生了变化. 关键词： 溶解 热对流 颗粒两相流 直接数值模拟     

超短激光打孔中快速相变的格子玻尔兹曼模拟

采用双重分布函数的格子玻尔兹曼模型，对单脉冲激光金属打孔过程中的快速相变传热进行研究.模型考虑了金属材料熔化后熔体的流动换热，并采用浸没移动边界方案对过程中的固液界面进行追踪.采用纯导热模型和考虑对流的换热模型计算，将结果和试验进行对比，结果表明：在激光打孔过程中熔体的流动对相变传热产生较大影响，采用考虑流动换热模型的结果与实验更接近.进而对熔化速度、熔化深度以及温度场的变化进行分析，并探讨不同激光工艺参数对相变过程的影响.模拟发现一个脉冲结束后，激光的脉宽越大，孔深越小，孔径越大，且最高温度较短脉冲激光越低.     

Two-dimensional MRT LB model for compressible and incompressible flows

In the paper we extend the Multiple-Relaxation-Time （MRT） Lattice Boltzmann （LB） model pro- posed in [Europhys. Lctt., 2010, 90： 54003] so that it is suitable also for incompressible flows. To decrease tile artificial oscillations, the convection term is discretized by the flux linfiter scheme with splitting technique. A new model is validated by some well-known benchmark tests, including Rie- mann problem and Couette flow, and satisfying agreements are obtained between the sinmlation results and ana.lytical ones. In order to show the merit of LB model over traditional methods, the non-equilibrium characteristics of system are solved. The simulation results are consistent with the physical analysis.