Moving least squares reconstruction for sharp interface immersed boundary methods

We propose a new approach for reconstructing velocity boundary conditions in sharp-inerface immersed boundary (IB) methods based on the moving least squares (MLS) interpolation method. The MLS is employed to not only reconstruct velocity boundary conditions but also to calculate the pressure and velocity gradients in the vicinity of the immersed body, which are required in fluid structure interaction problems to obtain the force exerted by the fluid on the structure. To extend the method to arbitrarily complex geometries with nonconvex shaped boundaries, the visibility method is combined with the MLS method. The performance of the proposed curvilinear IB MLS (CURVIB-MLS) is demonstrated by systematic grid-refinement studies for two- and three-dimensional tests and compared with the standard CURVIB method employing standard wall-normal interpolation for reconstructing boundary conditions. The test problems are flow in a lid-driven cavity with a sphere, uniform flow over a sphere, flow on a NACA0018 airfoil at incidence, and vortex-induced vibration of an elastically-mounted cylinder. We show that the CURVIB-MLS formulation yields a method that is easier to implement in complex geometries and exhibits higher accuracy and rate of convergence relative to the standard CURVIB method. The MLS approach is also shown to dramatically improve the accuracy of calculating the pressure and viscous forces imparted by the flow on the body and improve the overall accuracy of FSI simulations. Finally, the CURVIB-MLS approach is able to qualitatively capture on relatively coarse grids important features of complex separated flows that the standard CURVIB method is able to capture only on finer grids.     

An efficient meshfree point collocation moving least squares method to solve the interface problems with nonhomogeneous jump conditions

We are going to study a simple and effective method for the numerical solution of the closed interface boundary value problem with both discontinuities in the solution and its derivatives. It uses a strong‐form meshfree method based on the moving least squares (MLS) approximation. In this method, for the solution of elliptic equation, the second‐order derivatives of the shape functions are needed in constructing the global stiffness matrix. It is well‐known that the calculation of full derivatives of the MLS approximation, especially in high dimensions, is quite costly. In the current work, we apply the diffuse derivatives using an efficient technique. In this technique, we calculate the higher‐order derivatives using the approximation of lower‐order derivatives, instead of calculating directly derivatives. This technique can improve the accuracy of meshfree point collocation method for interface problems with nonhomogeneous jump conditions and can efficiently estimate diffuse derivatives of second‐ and higher‐orders using only linear basis functions. To introduce the appropriate discontinuous shape functions in the vicinity of interface, we choose the visibility criterion method that modifies the support of weight function in MLS approximation and leads to an efficient computational procedure for the solution of closed interface problems. The proposed method is applied for elliptic and biharmonic interface problems. For the biharmonic equation, we use a mixed scheme, which replaces this equation by a coupled elliptic system. Also the application of the present method to elasticity equation with discontinuities in the coefficients across a closed interface has been provided. Representative numerical examples demonstrate the accuracy and robustness of the proposed methodology for the closed interface problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1031–1053, 2015     

基于激光雷达的垂直能见度反演算法及其误差评估

针对大气垂直方向上消光系数分布不均匀难以用传统方法直接测量垂直能见度的问题,提出了一种基于激光雷达探测垂直能见度的计算方法。根据大气辐射传输基本原理,借助于辐射传输方程,推导出了垂直能见度的计算公式;然后利用激光雷达原理方程和Klett算法反演出大气垂直方向上的消光系数分布,基于此提出了垂直能见度的迭代算法。最后,利用灰色模型GM(1,1)和批统计算法,对激光雷达反演得到的后向散射系数进行了评估,给出了误差置信区间为(0.760±0.339)×10^-4(srad·km)^-1。结果表明,该方法是一种特别有效的计算垂直能见度的方法,符合探测的基本需求,且误差小精度高。     

Analysis of the fringes visibility generated by a lateral cyclic shear interferometer in the retrieval of the three-dimensional surface information of an object

We report the evaluation of the topography of an object with projected fringes generated with a lateral cyclic shear interferometer (CSI) and we then compare the topography recovery obtained with the proposed method with the one obtained from a coordinate measuring machine (CMM). We also study how the fringes visibility along the z axis affects the retrieval. Finally, we discuss the advantages and drawbacks of this profilometry system.     

用于反演半导体激光云高仪垂直能见度的算法

根据云高仪工作原理,从米散射激光雷达方程出发,推导出一种低层大气中气溶胶消光系数边界值的获取方法,并在此基础上利用Fernald方法反演得到了大气消光系数的垂直分布;通过垂直路径上的消光系数反演得到垂直路径的平均能见度,结果表明,该平均能见度更接近垂直方向真实的能见度。用该算法处理了半导体激光云高仪实际测量的数据,并与Vaisala云高仪测量的垂直能见度进行了对比,结果表明:该算法在反演大气垂直能见度时具有较高的准确性和可靠性。     

1．06μm激光大气后向散射强度仿真计算研究

针对激光大气后向散射信号对激光导引头制导性能的影响,分析了激光大气后向散射的形成,利用MATLAB语言对1．06μm激光大气后向散射信号强度进行了仿真计算,同时设计开发了一套模拟计算软件,并通过试验验证了软件计算激光大气后向散射信号强度的准确性。     

基于多路可视图的健康与心梗患者心电图信号复杂网络识别

可视图(visibility graph, VG)算法已被证明是将时间序列转换为复杂网络的简单且高效的方法,其构成的复杂网络在拓扑结构中继承了原始时间序列的动力学特性.目前,单维时间序列的可视图分析已趋于成熟,但应用于复杂系统时,单变量往往无法描述系统的全局特征.本文提出一种新的多元时间序列分析方法,将心梗和健康人的12导联心电图(electrocardiograph, ECG)信号转换为多路可视图,以每个导联为一个节点,两个导联构成可视图的层间互信息为连边权重,将其映射到复杂网络.由于不同人群的全连通网络表现为完全相同的拓扑结构,无法唯一表征不同个体的动力学特征,根据层间互信息大小重构网络,提取权重度和加权聚类系数,实现对不同人群12导联ECG信号的识别.为判断序列长度对识别效果的影响,引入多尺度权重度分布熵.由于健康受试者拥有更高的平均权重度和平均加权聚类系数,其映射网络表现为更加规则的结构、更高的复杂性和连接性,可以与心梗患者进行区分,两个参数的识别准确率均达到93.3%.     

X‐ray speckle visibility spectroscopy in the single‐photon limit

The technique of speckle visibility spectroscopy has been employed for the measurement of dynamics using coherent X‐ray scattering. It is shown that the X‐ray contrast within a single exposure can be related to the relaxation time of the intermediate scattering function, and this methodology is applied to the diffusion of 72 nm‐radius latex spheres in glycerol. Data were collected with exposure times as short as 2 ms by employing a resonant shutter. The weak scattering present for short exposures necessitated an analysis formalism based on the spatial correlation function of individual photon charge droplets on an area detector, rather than the usual methods employed for intensity correlations. It is demonstrated that this method gives good agreement between theory and experiment and thus holds promise for extending area‐detector‐based coherent scattering methods to the study of faster dynamics than previously obtainable.     

两相流有限穿越可视图演化动力学研究

针对小管径两相流流动特性, 全新优化设计弧形对壁式电导传感器. 通过动态实验在获取传感器测量信号的基础上, 采用有限穿越可视图理论构建对应于不同流型的两相流复杂网络. 通过分析发现, 有限穿越可视图网络异速生长指数和网络平均度值的联合分布可实现对小管径两相流的流型辨识; 有限穿越可视图度分布曲线峰值可有效刻画与泡径大小分布相关的流动物理结构细节特征; 网络平均度值可表征流动结构的宏观特性; 网络异速生长指数对流体动力学复杂性十分敏感, 可揭示不同流型演化过程中的细节演化动力学特性. 两相流测量信号的有限穿越可视图分析为揭示两相流流型的形成及演化动力学机理提供了新途径. 关键词： 两相流 复杂网络 有限穿越可视图 网络异速生长指数     

CCD激光测距实验

CCD激光测距实验是集光、机、电为一体的综合性实验，学生可以用简单的几何光学知识、采用激光遥感方式实现距离的非接触精确测量．本文简述了实验原理，给出了计算机自动处理和显示的系统硬件和软件环境，以及系统的校准方法，最后介绍了模拟实验结果。