低温回热制冷机内交变流动与换热的格子-Boltzmann方法模拟

采用格子-Boltzmann方法(LBM)分别通过对二维通道及多孔介质内交变流动与换热进行了数值模拟,以研究低温回热制冷机中脉冲管内及回热器内交变流动与换热规律.结果表明:二维通道交变流动结果与解析解吻合得很好.当交变流动数较大时,出现速度环形效应,温度表现为周期波动变化.多孔介质内各点速度总体均呈现交替波动变化,但振...     

基于格子波尔兹曼方法的回热器数值模拟

利用格子玻尔兹曼方法,直接对蚀刻薄片和层叠丝网回热器的微观结构流场进行了模拟.得到了两种回热器填料的微观流场和两端的压差.模拟结果显示,当回热器的直径、水力直径和填充率相近情况下,不同流速下蚀刻薄片卷裹式回热器的稳态阻力系数均比层叠丝网回热器小.稳态阻力系数的模拟变化趋势与实验一致.     

基于格子Boltzmann方法的疏水表面润湿性数值模拟

润湿性对固体表面上液体的各种动力学行为具有重要影响,疏水表面的特殊润湿性是其在减阻、降噪、防污等领域有着广泛应用前景的根本原因。基于Shan-Chen模型的格子Boltzmann方法对疏水表面润湿性进行数值模拟,获得了材料属性和微形貌对疏水表面润湿性的影响规律。研究表明,要使疏水表面处于Cassie-Baxter润湿状态,微形貌高度必须大于某一临界值,而当疏水表面一旦处于Cassie-Baxter润湿状态后,继续增加微形貌高度也不会提高其疏水性能;疏水表面的表观接触角随气液界面分数先增大后减小,且存在一个最佳的气液界面分数使表观接触角达到最大。     

多孔介质模型镁合金流变铸轧LBM模拟

根据REV尺度多孔介质格子-Boltzmann方法基本理论,建立求解高固相率半固态浆料渗流过程的热流耦合数值模型.应用模型对变孔隙率Poiseuille流动及纯扩散凝固进行模拟.LBM模拟结果与有限差分模拟吻合.分析了压差作用下浆料在平板间渗流冷却与热流耦合过程.模拟结果表明浆料渗流速度随温度降低而减慢,相同时刻下凝固潜热较大的浆料渗流Re数高.利用提出的模型初步探索了AZ91D镁合金半固态流变铸轧过程的热流场分布.     

改善瓶用聚酯结晶性能的方法

本提出了几种通过共聚改性提高聚酯冷结晶温度的方法，同时探索了缩聚工艺条件对聚酯结晶性能的影响，实验结果表明，通过降低缩聚初期的温度，可显减慢聚酯的结晶速度。     

微通道气体流动的格子波尔兹曼模拟

通过隐式格子波尔兹曼方程,并采用壁面平衡边界条件以及二阶关系,模拟了微通道气体流动中的非线性压力和壁面滑移速度,模拟结果与Arkilic的解析结果十分吻合,验证了格子波尔兹曼方法在滑移流区的有效性.     

基于格子玻尔兹曼方法的低We数液滴铺展凝固模拟研究

为了研究3D喷墨打印中液滴在壁面的铺展、凝固现象与机理,基于格子玻尔兹曼方法(Lattice Boltzmann method,LBM),建立了三维多组分相变模型,模拟计算了单液滴在低We数条件下与低温基板碰撞后的演变过程。在模拟过程中考虑了壁温、壁面润湿性等因素对于液滴的铺展、凝固的影响。模拟结果表明,在非润湿性壁面,液滴铺展产生震荡阻尼现象,通过改变壁温控制液滴凝固速度可以达到阻碍或者促进液滴铺展;而在润湿性壁面,凝固会阻碍液滴铺展,液滴最终铺展因子随接触角减小而降低,并且与壁温成正比关系。此外,壁温在低于一定范围后才会对液滴形貌(铺展因子、接触角)造成明显影响。     

霜层生长过程中的导热模型

在逾渗理论的基础上,对霜层的导热系数进行分析.现有研究表明由于冰晶体的不断生长,会形成不同的集团,当在霜层中形成大的连通集团的时候,霜层导热系数的上升速度会突然加快.在此基础上提出了基于逾渗理论的导热模型,并采用数值模拟的方法验证了上述过程.以临界点为分界线,分别导出了霜层生长过程中,不同孔隙率下的导热系数的通用数学表达式.与其他导热模型的结果比较表明,本模型更加符合实际,并且有较大的使用范围.     

聚氧化乙烯等温结晶过程的计算机模拟

用Monte Carlo方法模拟了聚氧化乙烯(PEO)在预先成核条件下的等温结晶过程,并用Avrami方程进行了结晶动力学处理.结果显示,Avrami方程能够成功地描述被模拟的PEO等温结晶的初期过程.随着结晶温度的升高,结晶速率和球晶的线生长速率减小,而Avrami指数值基本不变,都接近于3.采用热台偏光显微镜(HSPOM)实验证实了模拟实验模型和结果的正确性.     

气相生长硒化镉单晶体的生长速度研究

以晶体生长理论为基础，计算了气相提拉生长CdSe晶体时的生长速率。结果表明用提拉法气相生长CdSe晶体时，晶体的气相生长速率将会随着时间的延长以指数关系快速的趋近于提拉速度，以后不再变化。由此，在选定温场的前提下，对CdSe单晶体的气相生长速度进行了优化，确定了3mm／d的生长速度，得到了平界面生长的尺寸为Ф20mm×30mm，电阻率高达10^9Ω·cm，且未观察到深能级陷阱的优质CdSe大单晶体.     

Frost growth around a cylinder in a wet air stream

This paper describes a two-dimensional model which permits the evaluation of the local properties during the frost formation. To achieve this objective it is necessary to know the local coefficients of heat and mass transfer which were determined by solving the flow, temperature and humidity fields. To solve the flow field, the governing equations were developed in terms of the stream and vorticity functions. A numerical solution is obtained for low Reynolds number up to 400. With the flow field solved, the temperature and humidity fields are determined and hence the local heat and mass coefficients. These coefficients were then used in the solution of a two stage model for frost formation permitting the prediction of the local frost properties around the cylinder such as density, thickness and temperature.     

Two-phase lattice Boltzmann simulation of natural convection in a Cu-water nanofluid-filled porous cavity: Effects of thermal boundary conditions on heat transfer and entropy generation

The present work investigates the effect of four different thermal boundary conditions on natural convection in a fluid-saturated square porous cavity to make a judicious choice of optimal boundary condition on the basis of entropy generation, heat transfer and degree of temperature uniformity. Four different heating conditions- uniform, sinusoidal and two different linear temperature distributions are applied on the left vertical wall of the cavity respectively, while maintaining the right vertical wall uniformly cooled and the horizontal walls thermally insulated. The two-phase thermal lattice Boltzmann (TLBM) model for nanofluid is extended to simulate nanofluid flow through a porous medium by incorporating the Brinkman–Forchheimer-extended Darcy model. The close agreement between present LBM solutions with the existing published results lends validity to the present findings. The current results indicate that the uniform and bottom to top linear heating are found to be efficient heating strategies depending on Rayleigh number (10³?≤?Ra?≤?10⁵) and Darcy number (10^?1?≤?Da?≤?10^?6). It is observed that the nanofluid improves the energy efficiency by reducing the total entropy generation and enhancing the heat transfer rate although its augmentation depends on the optimal volume fraction of nanoparticles.     

Numerical simulation of dendritic crystal growth in a channel

A quasi-stationary mathematical model of dendritic growth from an undercooled melt in a channel is considered. This is based on the integral representation formula of the problem which assumes that the crystals grow steadily so that the temperature field satisfies a convection-diffusion equation in the moving frame of reference. Numerically, the above phenomenon is very unstable and any small perturbation in the interface shape may give rise to the growth of an artificial (numerical) branch in the dendrite. Although important advances have been made in the numerical simulation of this process in the last decade, more stable numerical schemes are required in order to understand, in detail, the dynamics of pattern formation. The proposed numerical approach is based on a BEM formulation where cubic B-splines are used to represent the contour geometry, as a result of the requirement of C²-continuity for the evaluation of the surface curvature. A standard quadratic interpolation is used for the densities within each boundary element, with an appropriate smoothing scheme to avoid the zig-zag instability which develops at the interface after a large number of time steps. The results obtained with the proposed numerical formulation are compared with experimental, theoretical and other numerical results.     

泥石流的增强碰撞算子格子Boltzmann模型

利用增强碰撞算子格子的Boltzmann方法,分析泥石流作为非Newton体的特点,根据泥石流的Bingham体模型,构造了一个特殊的增强碰撞算子的格子Boltzmann模型,并由此成功地模拟了泥石流入汇主河的运动过程,为预测和防治泥石流入汇主河所造成的灾害提供了一些理论依据。     

Detection of surface lattice defects using line scan method

The presence of different kinds of surface lattice defects such as missing atom, interstitial atom, line defects, in graphite single crystal have been identified by using scanning tunneling microscope. These defects cause displacement of atoms from their mean position and lattice strain is introduced. By measuring the displacement of atoms from their mean position. lattice strain has been calculated. It is found that among single point defects, vacancies cause maximum lattice strain. Paper presented at the poster session of MRSI AGM V1, Kharagpur, 1995     

Lattice Boltzmann simulation of flow past a non-spherical particle

Lattice Boltzmann method was used to predict the fluid-particle interaction for arbitrary shaped particles. In order to validate the reliability of the present approach, simulation of flow past a single stationary spherical, cylindrical or cubic particle is conducted in a wide range of Reynolds number (0.1 < Re_p < 3000). The results indicate that the drag coefficient is closely related to the particle shape, especially at high Reynolds numbers. The voxel resolution of spherical particle plays a key role in accurately predicting the drag coefficient at high Reynolds numbers. For non-spherical particles, the drag coefficient is more influenced by the particle morphology at moderate or high Reynolds numbers than at low ones. The inclination angle has an important impact on the pressure drag force due to the change of projected area. The simulated drag coefficient agrees well with the experimental data or empirical correlation for both spherical and non-spherical particles.     

Particulate flow simulation using Lattice Boltzmann Method: A rheological study

This work deals with calculating rheological properties of a suspension of particles in a fluid. A suspension of mono- and poly-disperse circular particles in shear flow is studied using two different methods for application of shear force: (a) by placing parallel walls at the top and bottom of the domain which are moving in opposite directions with the same velocity, and (b) using the Lees–Edwards boundary condition. The system which starts moving from rest, is allowed to reach a statistically steady state. Rheological properties namely, bulks shear stress, effective viscosity and normal stress difference of the suspension at different particle-based Reynolds numbers and different mean particle area fractions are calculated. Furthermore, the effect of size distribution on the relative effective viscosity of the suspension is investigated. Comparison of the present results with empirical formulations found in the literature shows reasonable agreement.