对流扩散方程的QC3无网格法

对流扩散方程作为偏微分运动方程的分支,在流体力学、气体动力学等领域有着重要应用.为解决对流扩散方程难以通过解析法得到解析解的难题,采用二阶一致3点积分(Quadratically Consistent 3-Point Integration,简称QC3)提高无网格法的计算效率,通过对积分点上形函数导数的修正,改善无网格法的精度和收敛性.本文将QC3无网格法拓展到对流扩散方程问题中,时域离散采用广义特征线Galerkin法,空间离散采用QC3法.数值结果表明,应用QC3无网格法得到的对流扩散问题数值解与解析解十分接近,验证了QC3无网格法解决对流扩散问题的可行性.     

移动加热器法生长碲锌镉晶体的组分输运与界面形貌研究

移动加热器法(THM)生长碲锌镉晶体时,界面稳定性对晶体生长的质量有很大影响。本文基于多物理场有限元仿真软件Comsol建立了THM生长碲锌镉晶体的数值模拟模型,讨论了Te边界层与组分过冷区之间的关系,对不同生长阶段的物理场、Te边界层与组分过冷区进行仿真研究,最后讨论了微重力对物理场分布的影响,并对比了微重力与正常重力下的生长界面形貌。模拟结果表明,Te边界层与组分过冷区的分布趋势是一致的,在不同生长阶段,流场中次生涡旋的位置会发生移动,从而导致生长界面的形貌随着生长的进行发生变化,同时微重力条件下形成的生长界面形貌最有利于单晶生长。因此,在晶体生长的中前期,对次生涡旋位置的控制和对组分过冷的削弱,是THM生长高质量晶体的有效方案。     

Effect of wall slip on the viscoelastic particle ordering in a microfluidic channel

The formation of a line of equally spaced particles at the centerline of a microchannel, referred as “particle ordering,” is desired in several microfluidic applications. Recent experiments and simulations highlighted the capability of viscoelastic fluids to form a row of particles characterized by a preferential spacing. When dealing with non-Newtonian fluids in microfluidics, the adherence condition of the liquid at the channel wall may be violated and the liquid can slip over the surface, possibly affecting the ordering efficiency. In this work, we investigate the effect of wall slip on the ordering of particles suspended in a viscoelastic liquid by numerical simulations. The dynamics of a triplet of particles in an infinite cylindrical channel is first addressed by solving the fluid and particle governing equations. The relative velocities computed for the three-particle system are used to predict the dynamics of a train of particles flowing in a long microchannel. The distributions of the interparticle spacing evaluated at different slip coefficients, linear particle concentrations, and distances from the channel inlet show that wall slip slows down the self-assembly mechanism. For strong slipping surfaces, no significant change of the initial microstructure is observed at low particle concentrations, whereas strings of particles in contact form at higher concentrations. The detrimental effect of wall slip on viscoelastic ordering suggests care when designing microdevices, especially in case of hydrophobic surfaces that may enhance the slipping phenomenon.     

Shedding vortex simulation method based on viscous compensation technology research

Owing to the influence of the viscosity of the flow field, the strength of the shedding vortex decreases gradually in the process of backward propagation. Large-scale vortexes constantly break up, forming smaller vortexes. In engineering, when numerical simulation of vortex evolution process is carried out, a large grid is needed to be arranged in the area of outflow field far from the boundary layer in order to ensure the calculation efficiency. As a result, small scale vortexes at the far end of the flow field cannot be captured by the sparse grid in this region, resulting in the dissipation or even disappearance of vortexes. In this paper, the effect of grid scale is quantified and compared with the viscous effect through theoretical derivation. The theoretical relationship between the mesh viscosity and the original viscosity of the flow field is established, and the viscosity term in the turbulence model is modified. This method proves to be able to effectively improve the intensity of small-scale shedding vortexes at the far end of the flow field under the condition of sparse grid. The error between the simulation results and the results obtained by using fine mesh is greatly reduced, the calculation time is shortened, and the high-precision and efficient simulation of the flow field is realized.     

改进型单晶炉提拉系统偏心平衡研究

随着光伏行业的快速发展, 对硅单晶的品质和长晶装备的稳定性的要求也不断提高。直拉法是生产硅单晶的主要方法,通过提高单晶炉副室的高度以扩大单晶硅的生产规模。由于副室高度的大幅增加,且单晶炉提拉头质心相对于旋转轴心有一定距离,对单晶炉整体稳定性有较大影响,从而降低了单晶硅的生产质量。针对此问题,对单晶炉建立可靠的力学分析模型,采用数值仿真方法,对单晶炉整体进行动力学响应分析,计算得到副室高度增加后的单晶炉工作时中钨丝绳下端晶棒的运动规律以及最大摆动幅度,为改进设计提供依据。数值仿真分析表明提高单晶炉副室高度后,提拉头较大的质心偏心是单晶炉提拉系统发生摆动的主要原因。在此基础上提出在提拉头上添加质心调节装置,通过控制系统调节可保证提拉头质心位置在旋转轴线上以降低提拉系统的摆动。     

Experimental and numerical analysis of cylindrical turbulators in a double pipe heat exchanger under turbulent flow conditions

In the present work, the use of cylindrical turbulators in a double pipe heat exchanger has been investigated. Cylindrical fin type of turbulators has been placed circumferentially separated by 90° on the outer side of an inner pipe at a regular pitch. Experimental studies were undertaken for different air flow rates in a turbulent regime whose Reynolds number range between 2500 and 10000. Heat transfer characteristics like Nu and friction factor have been experimentally determined. Parametric studies were conducted by changing the pitch and also the orientation of the turbulators. Nu and friction factor were found to increase as the pitch is reduced. A model with alternatively changed orientation outperformed others by exhibiting highest Nu and reduced friction factor.     

Thermoelastic solutions for thermal distributions moving over thin slim rod under memory-dependent three-phase lag magneto-thermoelasticity

Abstract

Enlightened by the Caputo fractional derivative, the present study deals with a novel mathematical model of generalized thermoelasticity to investigate the transient phenomena due to the influence of magnetic field and moving heat source in a rod in the context of three-phase lag (TPL) theory of thermoelasticity. Both ends of the rod are fixed and heat insulated. Employing Laplace transform as a tool, the problem has been transformed into the space-domain and solved analytically. Finally, solutions in the real-time domain are obtained by applying the inverse Laplace transform. Numerical calculation for stress, displacement, and temperature within the rod is carried out and displayed graphically. The effect of moving heat source speed on temperature, stress, and temperature is studied. It is found from the distributions that the temperature, thermally induced displacement and stress of the rod are found to decrease at large source speed. For the better understanding of the effect of moving heat source on all the distributions, three animations are added.     

An improved depth-averaged nonhydrostatic shallow water model with quadratic pressure approximation

Phase-resolved information is necessary for many coastal wave problems, for example, for the wave conditions in the vicinity of harbor structures. Two-dimensional (2D) depth-averaging shallow water models are commonly used to obtain a phase-resolved solution near the coast. These models are in general more computationally effective compared with computational fluid dynamics software and will be even more capable if equipped with a parallelized code. In the current article, a 2D wave model solving the depth-averaged continuity equation and the Euler equations is implemented in the open-source hydrodynamic code REEF3D. The model is based on a nonhydrostatic extension and a quadratic vertical pressure profile assumption, which provides a better approximation of the frequency dispersion. It is the first model of its kind to employ high-order discretization schemes and to be fully parallelized following the domain decomposition strategy. Wave generation and absorption are achieved with a relaxation method. The simulations of nonlinear long wave propagations and transformations over nonconstant bathymetries are presented. The results are compared with benchmark wave propagation cases. A large-scale wave propagation simulation over realistic irregular topography is shown to demonstrate the model's capability of solving operational large-scale problems.     

基于平面度的燃料电池电化学性能模拟

 固体氧化物燃料电池的翘曲会影响电极-盖板界面的接触情况，从而影响电化学性能，对相关制造工艺提出了很大的挑战．为了分析燃料电池平面度对放电过程的影响，揭示其潜在的风险，我们建立了两个基于有限元法的仿真模型，对考虑平面度缺陷的燃料电池封装和放电进行分析．在对固体氧化物燃料电池进行平面度测量的基础上，首先建立了具有真实燃料电池翘曲特性的几何模型，分析封装过程中接触压力的分布情况．然后将接触电阻的仿真结果导入到三维多物理场耦合模型中，模拟具有平面度缺陷的燃料电池电化学性能．计算结果展示了燃料电池两侧封装过程中接触压力的分布情况．通过对比有接触电阻和无接触电阻的燃料电池电流密度，分析了电池与盖板的接触对放电过程的影响．结果表明，燃料电池的凹陷面较难达到满意的接触状态，需要比凸起面更大的封装压力．燃料电池表面接触电阻的变化将引起电流传导路径的变化，产生局部高电流或低电流．这项工作强调了在燃料电池中保持均匀分布的接触电阻的重要性，为今后的优化工作奠定了基础．     

A fractional order HIV-TB coinfection model with nonsingular Mittag-Leffler Law

The biological models for the study of human immunodeficiency virus (HIV) and its advanced stage acquired immune deficiency syndrome (AIDS) have been widely studied in last two decades. HIV virus can be transmitted by different means including blood, semen, preseminal fluid, rectal fluid, breast milk, and many more. Therefore, initiating HIV treatment with the TB treatment development has some advantages including less HIV-related losses and an inferior risk of HIV spread also having difficulties including incidence of immune reconstitution inflammatory syndrome (IRIS) because of a large pill encumbrance. It has been analyzed that patients with HIV have more weaker immune system and are susceptible to infections, for example, tuberculosis (TB). Keeping the importance of the HIV models, we are interested to consider an analysis of HIV-TB coinfected model in the Atangana-Baleanu fractional differential form. The model is studied for the existence, uniqueness of solution, Hyers-Ulam (HU) stability and numerical simulations with assumption of specific parameters.