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

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

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

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

配合参数对单晶炉驱动系统的影响分析

光伏产业的发展使得对硅材料的需求日益增加,同时硅单晶生产行业竞争也日趋激烈。作为生产硅单晶的重要装备,单晶炉的稳定性和可靠性关系到硅单晶生产效率的提升和成本的下降,因此其驱动系统的设计和优化成为装备制造的关键环节。本文以NVT-HG2000-V1型硅单晶生长炉的驱动系统为研究对象,用SolidWorks三维建模实现虚拟装配,采用ADAMS建立其动力学仿真模型,并对驱动系统的运动过程进行仿真模拟。采用控制变量法定量分析了铜套与升降轴的配合间隙及丝杠参数对驱动力和驱动力矩的影响规律,进而在提高硅单晶生长炉装备稳定性和可靠性方面给出合理的技术建议。结果表明,铜套与升降轴的配合间隙达到0.071 mm后能有效降低驱动系统运行所需驱动力矩,丝杠倾斜度、螺纹螺距与螺纹间摩擦系数的增大均会导致驱动系统运行所需力矩大幅增加。     

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.     

Performance evaluation of functionalized ordered mesoporous silica for VOCs adsorption by molecular dynamics simulation

Volatile organic compounds (VOCs) are growing pollutants now that cause the serious environmental pollution and threaten human health. The functionalized ordered mesoporous silica (FOMS) has attracted considerable attention in adsorbing VOCs. In this paper, the molecular dynamics simulation was used to simulate the adsorption performance of FOMS on VOCs (acetone, ethyl acetate and toluene). After simulating different pore sizes (2 nm, 3 nm and 4 nm) adsorption performances of ordered mesoporous silica (OMS) on VOCs, OMS with a pore size of 4 nm was selected to further study the influence of functional groups (vinyl, methyl, and phenyl). The following law was obtained: the saturated adsorption capacities of vinyl-functionalized OMS (V-FOMS) to acetone, ethyl acetate and toluene were 3.045 mmol.g^?1, 2.568 mmol.g^?1 and 1.976 mmol.g^?1 respectively; the saturated adsorption capacities of methyl-functionalized OMS (M-FOMS) to acetone, ethyl acetate and toluene were 2.798 mmol.g^?1, 2.312 mmol.g^?1 and 1.698 mmol.g^?1 respectively; the saturated adsorption capacities of phenyl-functionalized OMS (P-FOMS) to acetone, ethyl acetate and toluene were 2.124 mmol.g^?1, 1.941 mmol.g^?1 and 1.539 mmol.g^?1 respectively. These results show that the adsorption ability of FOMS for different adsorbates follows the sequence of acetone > ethyl acetate > toluene. Furthermore, the interaction between functional groups (vinyl, methyl and phenyl) in FOMS and VOCs was explored. It is found that the interaction between different functional groups and adsorbates is different (interaction energy effect). This interaction energy effect promotes FOMS to better adsorb VOCs. This work would provide fundamental understanding and guidance for the development of novel adsorption materials for the adsorption of VOCs.     

Direct visualization of structural defects in 2D semiconductors

Direct visualization of the structural defects in two-dimensional (2D) semiconductors at a large scale plays a significant role in understanding their electrical/optical/magnetic properties, but is challenging. Although traditional atomic resolution imaging techniques, such as transmission electron microscopy and scanning tunneling microscopy, can directly image the structural defects, they provide only local-scale information and require complex setups. Here, we develop a simple, non-invasive wet etching method to directly visualize the structural defects in 2D semiconductors at a large scale, including both point defects and grain boundaries. Utilizing this method, we extract successfully the defects density in several different types of monolayer molybdenum disulfide samples, providing key insights into the device functions. Furthermore, the etching method we developed is anisotropic and tunable, opening up opportunities to obtain exotic edge states on demand.     

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.     

Mechanical property and deformation mechanism of gold nanowire with non-uniform distribution of twinned boundaries:A molecular dynamics simulation study

The mechanical property and deformation mechanism of twinned gold nanowire with non-uniform distribution of twinned boundaries(TBs) are studied by the molecular dynamics(MD) method. It is found that the twin boundary spacing(TBS) has a great effect on the strength and plasticity of the nanowires with uniform distribution of TBs. And the strength enhances with the decrease of TBS, while its plasticity declines. For the nanowires with non-uniform distribution of TBs, the differences in distribution among different TBSs have little effect on the Young's modulus or strength, and the compromise in strength appears. But the differences have a remarkable effect on the plasticity of twinned gold nanowire. The twinned gold nanowire with higher local symmetry ratio has better plasticity. The initial dislocations always form in the largest TBS and the fracture always appears at or near the twin boundaries adjacent to the smallest TBS. Some simulation results are consistent with the experimental results.     

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

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

Pd–H and Ni–H phase diagrams using cluster variation method and Monte Carlo simulation

This work focuses on interstitial solid solutions of hydrogen in the face-centred cubic (fcc) host lattice of palladium and nickel, using a first-principles based approach. Cluster Variation Method (CVM) and Monte Carlo simulation algorithms were especially designed, allowing a coupled use of both techniques, to study hydrogen–vacancy interactions inside an fcc metallic host lattice. First-principles calculations provided the H–Vac interaction energies by structure inversion method. The phase diagrams and thermodynamic properties were computed using only theoretical inputs. The mechanisms leading to the formation of the miscibility gaps observed for both Pd–H and Ni–H systems and the hydrogen ordering on palladium interstitial lattice were reproduced without any empirical term.