双重点最小二乘配点无网格法

配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。     

材料三维微结构表征及其晶体塑性有限元模拟

材料的力学性能，尤其是在有限变形下所呈现的宏观各向异性，是材料结构设计和服役寿命考虑的关键因素。由于宏观模型不能较好地反映材料微观结构（晶粒的形貌和取向等）对宏观塑性各向异性的影响，因此，本文建立了能实际反映晶粒形貌的三维Voronoi模型，并基于晶体塑性理论对铝合金在有限变形下的响应进行计算。首先，建立反映材料微结构的代表性体积单元RVE模型进行计算，并与实验结果进行对比验证。然后，以单向拉伸为例，分析了有限变形过程中试件的晶粒形貌和取向分布等微观因素对宏观各向异性演化的影响，并从材料和结构两个层面讨论了微观结构对宏观力学性能的影响。结果表明，本文模型能够反映微观结构对宏观力学性能的影响，为实际生产制造领域构件的力学性能提供可靠的预测。     

Lattice Boltzmann simulation on thermal performance of composite phase change material based on Voronoi models

Phase change materials(PCMs) are important for sustaining energy development. For the thermal performance enhancement, the composite PCM with metal foam reconstructed by the Voronoi method is investigated in this work. The lattice Boltzmann method(LBM) is used to analyze the melting process on a pore scale. The melting interface evolution and temperature contour of the composite PCM are explored and compared with those of pure PCM. Moreover, structure parameters including the pore density, porosity and irregularity are investigated comprehensively, indicating that the additive of metal foam strengthens the melting performance of PCM obviously. Compared with pure PCM, the composite PCM has quick rates of the melting front evolution and heat transfer. The heat conduction plays a great role in the whole melting process since the convection is weakened for the composite PCM. To improve the melting efficiency, a larger pore density and smaller irregularity are recommended in general. More significantly, a suitable porosity is determined based on the requirement for the balance between the melting rate and heat storage capacity in practical engineering.     

基于Voronoi切割算法的碎冰区构造及离散元分析

离散元方法广泛应用于海冰，特别是碎冰区的动力过程及其对海洋结构作用过程的数值模拟。为构造碎冰区中的冰块几何特性，基于二维Voronoi图方法对计算域进行随机切割以生成碎冰区中冰块的几何形态，并采用球体单元对每个碎冰块单元进行填充，从而确定碎冰区的初始分布场。在采用Voronoi图进行碎冰区构造时，可对冰块尺寸、几何形态和密集度等海冰参数进行设定。为确定冰块的不同几何规则度，综合采用排斥法和扰动法以定量地控制碎冰块几何形态从完全随机分布到规则分布的连续变换。为分析不同几何规则度下碎冰块的几何特性概率分布规律，对计算域内冰块的面积和边数等参数进行统计分析，从而可更合理地参数化控制初始冰场中碎冰块的几何特性。在此基础上，本文基于粘接-破碎的球体离散元方法对不同冰况下锥体结构的冰荷载进行了数值计算，讨论分析了碎冰区的海冰密集度、冰块面积和几何规则度对冰载荷的影响。     

Dislocation density crystalline plasticity modeling of lath martensitic microstructures in steel alloys

A three-dimensional multiple-slip dislocation density-based crystalline formulation, specialized finite-element formulations and Voronoi tessellations adapted to martensitic orientations were used to investigate large strain inelastic deformation modes and dislocation density evolution in martensitic microstructures. The formulation is based on accounting for variant morphologies and orientations, retained austenite and initial dislocation densities that are uniquely inherent to martensitic microstructures. The effects of parent austenite orientation and retained austenite were also investigated for heterogeneous fcc/bcc crystalline structures. Furthermore, the formulation was used to investigate microstructures mapped directly from SEM/EBSD images of martensitic steel alloys. The analysis indicates that variant morphology and orientations have a direct effect on dislocation density accumulation and inelastic localization in martensitic microstructures, and that lath directions, orientations and arrangements are critical characteristics of high strength martensitic deformation and behavior.     

Coupled Models and Parallel Simulations for Three-Dimensional Full-Stokes Ice Sheet Modeling

A three-dimensional full-Stokes computational model is considered for determining the dynamics,temperature,and thickness of ice sheets.The goveming thermomechanical equations consist of the three-dimensional full-Stokes system with nonlinear rheology for the momentum,an advective-diffusion energy equation for temperature evolution,and a mass conservation equation for ice-thickness changes.Here,we discuss the variable resolution meshes,the finite element discretizations,and the parallel algorithms employed by the model components.The solvers are integrated through a well-designed coupler for the exchange of parametric data between components.The discretization utilizes high-quality,variable-resolution centroidal Voronoi Delaunay triangulation meshing and existing parallel solvers.We demonstrate the gridding technology,discretization schemes,and the efficiency and scalability of the parallel solvers through computational experiments using both simplified geometries arising from benchmark test problems and a realistic Greenland ice sheet geometry.     

Numerical simulations of the steady Navier–Stokes equations using adaptive meshing schemes

In this paper, we consider an adaptive meshing scheme for solution of the steady incompressible Navier–Stokes equations by finite element discretization. The mesh refinement and optimization are performed based on an algorithm that combines the so‐called conforming centroidal Voronoi Delaunay triangulations (CfCVDTs) and residual‐type local a posteriori error estimators. Numerical experiments in the two‐dimensional space for various examples are presented with quadratic finite elements used for the velocity field and linear finite elements for the pressure. The results show that our meshing scheme can equally distribute the errors over all elements in some optimal way and keep the triangles very well shaped as well at all levels of refinement. In addition, the convergence rates achieved are close to the best obtainable. Extension of this approach to three‐dimensional cases is also discussed and the main challenge is the efficient implementation of three‐dimensional CfCVDT generation that is still under development. Copyright © 2007 John Wiley & Sons, Ltd.     

基于Voronoi图和动态自适应蚁群算法的UAV航迹规划

以UAV航迹规划为应用背景,提出了一种基于Voronoi图和动态自适应蚁群算法的航迹规划方法;为了提高航迹规划问题最优解的质量及全局求解能力,克服传统蚁群算法收敛速度慢、容易陷入局部最优等缺点,提出了一种动态自适应蚁群算法;采用动态自适应航迹点选择策略并将信息素更新规则和挥发系数进行动态自适应调整变化来对蚁群算法进行了改进,提高了算法的求解效率;根据战场已知威胁源生成Voronoi加权图,并与所提的动态自适应蚁群算法相结合求解规划空间中的最优航迹;考虑到UAV的物理约束限制,对生成的可行航迹进行平滑优化;仿真结果表明,该方法能够为UAV规划出一条满足要求的可飞航迹,验证了所提方法在解决航迹规划问题时是可行、有效的;     

Approximations of a Ginzburg-Landau model for superconducting hollow spheres based on spherical centroidal Voronoi tessellations

In this paper the numerical approximations of the Ginzburg- Landau model for a superconducting hollow spheres are constructed using a gauge invariant discretization on spherical centroidal Voronoi tessellations. A reduced model equation is used on the surface of the sphere which is valid in the thin spherical shell limit. We present the numerical algorithms and their theoretical convergence as well as interesting numerical results on the vortex configurations. Properties of the spherical centroidal Voronoi tessellations are also utilized to provide a high resolution scheme for computing the supercurrent and the induced magnetic field.