环形浅液层内热流体波的可视化实验研究

对外壁加热、内壁冷却、厚度为1.0 mm的环形硅油浅液层内的热毛细对流及其稳定性进行了可视化实验研究,观察到了旋转的螺纹状的热流体波,确定了发生热流体波的临界温差值△Tc=4.8 K,证实了热流体波从两组逐渐演变为一组的发展过程,验证了已有的数值模拟结果.     

Identification of Regions with Inhomogeneous Particle Behavior in Dilute Gas‐solid Flows

The behavior of the particulate phase in a highly turbulent gas flow has been investigated in a vertical channel. Variations of the flow configuration (1. Flow past a cylinder, 2. flow past a wall‐mounted obstacle and 3. flow around a horizontally injected jet) have been subject to both experiments and numerical simulations. The velocity vector field of the solid phase has been measured by digital particle image velocimetry (DPIV). The measurements have been focused on particle‐obstacle collisions and crossflow in the vicinity of the jet nozzle using the lately developed twinpeak detection method. By application of this method regions of highly inhomogeneous particle behavior could be detected mainly upstream of the flow perturbation. Numerical results have been obtained by an Eulerian‐Lagrangian method on boundary‐fitted grids. Particle‐particle interactions as well as interphase exchange of momentum have been taken into account. The simulation results showed to be well in accordance with the experimental results.     

Potential of digital holography in particle measurement

This paper describes the potential power of digital holography in particle measurement and its expected development in the near future. In digital holography, image reconstruction is carried out numerically on a computer using observed hologram patterns and some quantitative information can be derived from the reconstructed images. In this paper, the basic concept and procedure of digital in-line holography are shown mainly for particle depth measurement and the performance test results obtained in numerical simulations and experiments are demonstrated to examine the potential of the present method.     

颗粒碰撞的直接模拟算法

为研究细微颗粒间平均碰撞率,构建了直接模拟三维空间内颗粒碰撞的数值方法.对剪切流内夹带的无惯性颗粒和自由运动颗粒的碰撞过程分别进行了数值模拟,得到的颗粒平均碰撞率的模拟值与理论值的相对误差小于1%.在碰撞颗粒对的搜索过程中,引入网络排除法等优化措施,使计算量降低两个量级以上,而不增大相对误差.采用不同的Δt模拟同一颗粒群的运动,发现较小的时间步长可以避免因为截断误差增大的相对误差.在不遗失颗粒碰撞的条件下,较密的计算网格可明显减少计算量.     

Material point method applied to multiphase flows

The particle-in-cell method (PIC), especially the latest version of it, the material point method (MPM), has shown significant advantage over the pure Lagrangian method or the pure Eulerian method in numerical simulations of problems involving large deformations. It avoids the mesh distortion and tangling issues associated with Lagrangian methods and the advection errors associated with Eulerian methods. Its application to multiphase flows or multi-material deformations, however, encounters a numerical difficulty of satisfying continuity requirement due to the inconsistence of the interpolation schemes used for different phases. It is shown in Section 3 that current methods of enforcing this requirement either leads to erroneous results or can cause significant accumulation of errors. In the present paper, a different numerical method is introduced to ensure that the continuity requirement is satisfied with an error consistent with the discretization error and will not grow beyond that during the time advancement in the calculation. This method is independent of physical models. Its numerical implementation is quite similar to the common method used in Eulerian calculations of multiphase flows. Examples calculated using this method are presented.     

单基阵三维纯方位水下信标声学定位方法

针对水下静止信标的纯方位定位问题,提出一种基于单基阵平台纯方位信息的定位方法。通过水下无人运动平台两次量测的空间角度信息建立非线性方程,利用拟牛顿迭代方法实现稳定数值求解。同时,对量测信息进行预处理和定位后置处理,以进一步提高算法的稳定性和精度。利用计算机模拟水下声传播环境进行定位仿真,结果表明,基于常规波束形成的方位估计结果,在方位量测误差标准差为1°时,所提方法的定位结果相对误差可达3.5%以内。同时,研究了不同方位量测误差条件下、不同海深测量误差条件、观测平台和目标不同态势下定位算法的定位性能。结果表明,此方法可用于对水下静止信标的声学定位,且鲁棒性好。     

A method for calibrating the position error of phase retrieval with transverse translation diversity in optical wavefront measurement

We present a method for calibrating the position error of phase retrieval with transverse translation diversity in optical wavefront in situation where the position error sharply influences algorithm precision. This method involves testifying that the essence of the position error in phase retrieval is the translation in frequency domain and the minimum of iterative error against position error reaches when there is no position error. Then the least square method is used for fitting the relationship between the iterative error and the assumed position error in polynomial function. The computer simulations used to prove the validity of this method are described. The results indicated that this method can calibrate the position error to approximately 1/500 mm which nearly has no influence on the phase retrieval.     

Nonequilibrium Green Functions Approach to Strongly Correlated Fermions in Lattice Systems

Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

进口浓度对水力旋流器颗粒分级的影响

通过对水力旋流器内液固多相流动的固体颗粒运动进行理论分析、数值模拟和实验测试,探讨了进口颗粒浓度对采用水力旋流器进行高炉污泥颗粒分级的影响。液固多相流动的数值模拟证实,固体颗粒在水力旋流器内的径向沉降速度近似与颗粒粒径的平方成正比。实验给出了不同进口浓度的颗粒分级效率曲线。本文的实验结果、数值模拟结果以及颗粒离心沉降的理论分析都表明,进口浓度对水力旋流器内高炉污泥的颗粒分级影响不大,该结果诠释了在高炉污泥脱锌操作所考虑的颗粒浓度范围内为什么低浓度下水力旋流器颗粒分级的数值模拟结果与较高浓度下的实验结果基本一致。     

Laser-based volumetric colour-coded three-dimensional particle velocimetry

We present a method of three-dimensional particle velocimetry with a single digital colour camera using multiple colour illumination to encode image depth over a large volume. A copper vapour laser operating at 511 nm is used to pump an optical fibre producing a multiple-wavelength beam via multiple order stimulated Raman scattering. The beam is dispersed and formed into a stack of thin sheets to illuminate a volume of space. The spatial co-ordinates of particles imaged within the illuminated volume are obtained from their imaged x,y positions with depth discerned from particle hue (set by the wavelength of illumination). The method exhibits an RMS depth error of 3% in relation to the thickness of the illuminated region. This paper reports a proof-of-principle of three-dimensional particle imaging using a multi-wavelength laser source with a view to 3D-3C particle velocimetry.     

Application of the dissipative particle dynamics method to magnetic colloidal dispersions

The validity of the application of the dissipative particle dynamics (DPD) method to ferromagnetic colloidal dispersions has been investigated by conducting DPD simulations for a two–dimensional system. First, the interaction between dissipative and magnetic particles has been idealized as some model potentials, and DPD simulations have been carried out using such model potentials for a two magnetic particle system. In these simulations, attention has been focused on the collision time for the two particles approaching each other and touching from an initially separated position, and such collision time has been evaluated for various cases of mass and diameter of dissipative particles and model parameters, which are included in defining the equation of motion of dissipative particles. Next, a multi–particle system of magnetic particles has been treated, and particle aggregates have been evaluated, together with the pair correlation function along an applied magnetic field direction. Such characteristics of aggregate structures have been compared with the results of Monte Carlo and Brownian dynamics simulations in order to clarify the validity of the application of the DPD method to particle dispersion systems. The present simulation results have clearly shown that DPD simulations with the model interaction potential presented here give rise to physically reasonable aggregate structures under circumstances of strong magnetic particle–particle interactions as well as a strong external magnetic field, since these aggregate structures are in good agreement with those of Monte Carlo and Brownian dynamics simulations.     

光子计数深度获取系统误差机理研究和矫正

采用两路光子时间到达点构建光纤式光子计数伪随机码深度获取系统.为了研究降低深度误差的方法,以高斯函数为激光回波脉冲,计算瞬时概率密度函数,引入"时间行走"效应数学模型,推导深度误差克拉美罗下限.随着激光回波能量的增大,深度误差先降低再增大,并且码长越长,深度误差越小.采用理论推导的累积分布函数,生成光子时间到达点,蒙特卡洛仿真伪随机序列光子探测过程,结果大于理论数值模拟,符合克拉美罗下界原理.17组标定实验表明:由于目标表面特性的不同而导致探测到激光回波中光子数的浮动,该浮动引发光子"时间行走"效应,并带来伪随机码深度获取系统的互相关函数的整体偏移.采用数值拟合方程拟合不同光子计数比例值下的深度误差,测量得到的光子计数比例值,代入拟合的矫正方程,矫正后的深度均方误差下降至1cm.     

修正内嵌边界法对颗粒运动的直接模拟

本文提出了一种基于内嵌边界法(IBM)的思想,在描述颗粒只有较少网格时能较为准确地计算气固多相流的模型.模型依据边界层理论,通过定义颗粒表面的速度分布函数来修正由于计算网格尺度大于边界层厚度所带来的误差.本文通过计算单颗粒绕流和单颗粒沉降的最终沉降速度并与实验进行比较验证了本文模型的准确性.     

A study of process parameters during pulsed Nd:YAG laser notching of C70S6 fracture splitting connecting rods

This paper reports on the laser notching technology in C70S6 steel for fracture splitting connecting rod using a Nd:YAG pulsed laser. The effects of process parameters on starting notch (SN) dimension and morphology were investigated by both finite element method (FEM) simulations and physical experiments for various process parameters. Optical microscopy and the scanning electron microscope (SEM) were used for the measurement of SN dimension and the observation of SN morphology in the experiment. The results were compared with the predictions. It was found that the FEM simulations results showed good consistency with the experiments, which indicates that the finite element model is feasible and reliable. Based on the principal findings from the two methods, optimum ranges of process parameters for different fracture splitting connecting rods were predicted, which are a flexibly adjusting notch depth, a curvature radius less than 0.08 mm and an opening angle within the range 18–26°. The results indicate that the predicting ranges are suitable for making good SNs, which has also been proved by the fracture splitting experiments.     

窄通道内热薄燃料表面火焰传播特性研究

利用实验和数值模拟对微重力和常重力条件下高度为14mm和10mm的窄通道内热薄纸张表面火焰传播特性进行了研究。不同重力条件下窄通道内火焰传播速度随气流速度变化的规律符合得较好,说明地面窄通道实验能够模拟微重力条件下材料表面火焰传播的主要特征。地面窄通道中浮力流动速度的最大值约为5cm／s,与常规实验通道（高度较大）相比...     

Nd: YAG激光烧蚀裂解加工技术模拟分析与实验研究

激光裂解技术能够极大改善发动机缸体主轴承座的加工质量并显著提高加工效率. 为探寻Nd:YAG激光烧蚀球墨铸铁材料裂解槽的裂解性能, 本文基于有限元法成功构建了发动机缸体主轴承座激光裂解加工过程仿真模型, 针对QT500-7球墨铸铁主轴承座的裂解参数进行了仿真分析. 研究结果表明: 在影响裂解质量的三个裂解槽几何参数中, 槽深较张角及曲率半径对裂解载荷的影响效应更为明显; 裂解载荷随槽深的增加而迅速降低, 随槽张角和曲率半径的增加而升高; QT500-7球墨铸铁发动机缸体主轴承座激光裂解加工优化参数应为裂解槽深选为0.5 mm, 裂解槽张角选为60^o, 裂解槽半径选为0.2 mm. 有限元模拟分析结果得到了单向拉伸实验结果的验证. 本工作通过ABAQUS仿真模拟及大量裂解载荷试验确立了裂解槽几何形状的优化参数, 为显著降低裂解载荷和优化裂解工艺提供了数值参考, 有利于实现发动机缸体加工的快速发展, 从而促进汽车工业实现绿色制造.     

Error estimates for the fast multipole method

Error estimates for algorithms based on truncations for evaluating electrostatic interactions in molecular dynamics applications are very important for several reasons. For example, the estimates are necessary to establish the validity of the simulations and can be used to estimate various simulation parameters. Very precise estimates have been found for the Ewald method and the related particle mesh Ewald method. However, for the very popular fast multipole method such a precise estimate is not available. In this paper, we illustrate the rather complicated error behavior of the fast multipole method and we use statistical methods to derive an estimate for the root mean square error on the forces. Furthermore, the expected maximum error on the force acting on a single particle is studied. The estimates are tested against errors obtained from simulations and are found to be very precise.     

Modified periodic-shell boundary condition for dissipative particle dynamics simulations

We have developed the modified periodic-shell boundary condition (BC) for dissipative particle dynamics (DPD) simulations, which enables us to simulate an outer flow problem around an obstacle using a small simulation region. In order to clarify the validity of this BC, we have treated a uniform flow past a circular cylinder. The present BC has been compared with the ordinary BC such as the uniform flow condition. Also, the present results have been compared with those of the numerical results of the Navier–Stokes equation. The ordinary uniform BC is seen to give rise to significantly distorted flow fields and also to significant disappearance of dissipative particles from the simulation region. In contrast, for the present modified periodic-shell BC, the number density of dissipative particles is kept almost constant during a simulation run, and the flow field is in reasonable agreement with the result, which has been obtained by numerical simulations of the Navier–Stokes equation.     

Numerical modeling of thermal force in a plasma for test-ion transport simulation based on Monte Carlo Binary Collision Model

We present a new numerical model of the thermal force in a plasma, based on the Monte Carlo Binary Collision Model (BCM) [T. Takizuka, H. Abe, J. Comput.Phys. 25 (1977) 205]. This model can be applied for the transport simulation of test ions. The model consists of two steps: (i) choosing a background plasma ion velocity from a distorted Maxwell distribution under the temperature gradient, and (ii) calculating a Coulomb collision between a test particle and the above chosen ion by using the BCM. For the step (i), we developed a velocity sampling method from a distorted Maxwellian, which enables the BCM to bring the thermal force on a test particle in the step (ii).A systematic series of simulations has been performed under various conditions to examine the model. The results of these simulations have been compared with the theoretical values, and it is shown that our model simulates the thermal force correctly for important characteristic features; dependences on the temperature gradient, the test particle velocity, and the background plasma density.