Reactor field reconstruction from sparse and movable sensors using Voronoi tessellation-assisted convolutional neural networks

The aging of operational reactors leads to increased mechanical vibrations in the reactor interior.The vibration of the incore sensors near their nominal locations is a new problem for neutronic field reconstruction.Current field-reconstruction methods fail to handle spatially moving sensors.In this study,we propose a Voronoi tessellation technique in combination with convolutional neural networks to handle this challenge.Observations from movable in-core sensors were projected onto the same global field structure using Voronoi tessellation,holding the magnitude and location information of the sensors.General convolutional neural networks were used to learn maps from observations to the global field.The proposed method reconstructed multi-physics fields(including fast flux,thermal flux,and power rate)using observations from a single field(such as thermal flux).Numerical tests based on the IAEA benchmark demonstrated the potential of the proposed method in practical engineering applications,particularly within an amplitude of 5 cm around the nominal locations,which led to average relative errors below 5% and 10% in the L_(2) and L_(∞)norms,respectively.

Investigation of stress-induced progressive failure of mine pillars using a Voronoi grain-based breakable block model

The Voronoi grain-based breakable block model(VGBBM)based on the combined finite-discrete element method(FDEM)was proposed to explicitly characterize the failure mechanism and predict the deformation behavior of hard-rock mine pillars.The influence of the microscopic parameters on the macroscopic mechanical behavior was investigated using laboratory-scale models.The field-scale pillar models(width-to-height,W/H=1,2 and 3)were calibrated based on the empirically predicted stress-strain curves of Creighton mine pillars.The results indicated that as the W/H ratios increased,the VGBBM effectively predicted the transition from strain-softening to pseudo-ductile behavior in pillars,and explicitly captured the separated rock slabs and the V-shaped damage zones on both sides of pillars and conjugate shear bands in core zones of pillars.The volumetric strain field revealed significant compressional deformation in core zones of pillars.While the peak strains of W/H=1 and 2 pillars were relatively consistent,there were significant differences in the strain energy storage and release mechanism.W/H was the primary factor influencing the deformation and strain energy in the pillar core.The friction coefficient of the structural plane was also an important factor affecting the pillar strength and the weakest discontinuity angle.The fracture surface was controlled by the discontinuity angle and the friction coefficient.This study demonstrated the capability of the VGBBM in predicting the strengths and deformation behavior of hard-rock pillars in deep mine design.

Relationship between polyhedral structures formed by tangent planes of ellipsoidal particles and system sound velocity

Internal polyhedral structures of a granular system can be investigated using the Voronoi tessellations.This technique has gained increasing recognition in research of kinetic properties of granular flows.For systems with mono-sized spherical particles,Voronoi tessellations can be utilized,while radial Voronoi tessellations are necessary for analyzing systems with multi-sized spherical particles.However,research about polyhedral structures of non-spherical particle systems is limited.We utilize the discrete element method to simulate a system of ellipsoidal particles,defined by the equation(x a)2+(y1)2+(z 1/a)2=1,where a ranges from 1.1 to 2.0.The system is then dissected by using tangent planes at the contact points,and the geometric quantities of the resulting polyhedra in different shaped systems,such as surface area,volume,number of vertices,number of edges,and number of faces,are calculated.Meanwhile,the longitudinal and transverse wave velocities within the system are calculated with the time-of-flight method.The results demonstrate a strong correlation between the sound velocity of the system and the geometry of the dissected polyhedra.The sound velocity of the system increases with the increase in a,peaking at a=1.3,and then decreases as a continues to increase.The average volume,surface area,number of vertices,number of edges,and number of faces of the polyhedra decrease with the increase in sound velocity.That is,these quantities initially decrease with the increase in a,reaching minima at a=1.3,and then increase with further increase of a.The relationship between sound velocity and the geometric quantities of the dissected polyhedra can serve as a reference for acoustic material design.

A typhoon-induced storm surge numerical model with GPU acceleration based on an unstructured spherical centroidal Voronoi tessellation grid

Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.

On the calibration and verification of Voronoi-based discontinuous deformation analysis for modeling rock fracture

Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.

Role of grain boundary networks in vortex motion in superconducting films

We study the vortex dynamics of the polycrystalline superconductors in the presence of both random point defects and the generated grain boundary(GB) networks with Voronoi diagram. The synergistic effect of adjacent GBs on restricting the vortex motion in intragranular region is proposed and the corresponding intensity factor of the synergistic effect which characterizes the strength of the synergistic restriction of adjacent grain boundaries is also determined in the present work.The interconnected GBs offer easy-flow channels for vortices in addition to pinning effects on the vortices. The combined channels and the vortex flow patterns in the superconducting film are analyzed in detail from molecular dynamics simulations. Furthermore, it is discovered that the critical current increases with the decrease of magnetic field intensity,temperature, and the average grain size. The large number of vortices results in the enhanced repulsive interaction forcing the vortices to move out from the GBs. The thermal depinning from GBs leads to the lower Lorentz force range. The increase of the grain size causes the number of GBs to decrease. In summary, these effects leads the critical current to become a decreasing function of magnetic field, temperature, and grain size.

移动渲染的网格细化方法

随着移动设备性能的逐渐提高,越来越多的图形渲染技术由普通PC移植到了移动平台,但移动设备的计算能力与PC相比依然存在较大差距。对于PC上能够渲染的复杂网格模型,需要进行简化后才能够在移动设备上渲染,而网格简化必然会导致模型的走样。虽然Phong着色可以实现大部分正向面的平滑效果,但是未能解决轮廓部分的不连续性。为了解决这个问题,利用一种基于顶点着色器的实时细化方法,在移动设备上实现了Phong Tessellation,达到了模型轮廓部分细化的目的。

单轴拉伸过程中双相不锈钢局部微观变形行为的细观力学有限元分析

针对双相不锈钢中奥氏体相和铁素体相分别展开了纳米压痕实验,并通过有限元反演得到两相各自的拉伸应力-应变关系,利用Voronoi Tessellation法生成代表性的微结构体积单元,对双相不锈钢的单轴拉伸行为进行了有限元仿真和模拟,研究了双相不锈钢在拉伸过程中的局部应力、应变分布和演化规律.结果表明,利用Voronoi Tessellation法建立单元模型,结合本文通过纳米压痕实验获取的两相力学性能参数,可以很好地模拟双相不锈钢的整体单拉行为,奥氏体比铁素体软,拉伸载荷下双相不锈钢的应变集中在奥氏体中,应力集中在铁素体中;局部应力应变的分布特征与两相分布特征和晶粒形状有关,最大应变值主要集中在奥氏体晶粒狭长且尖锐的区域,而最大应力则主要发生在铁素体晶粒狭长和尖锐的区域;对于奥氏体和铁素体晶粒占比相当的双相不锈钢,其虽然可以具有较为综合的宏观力学性能,但是其微观应力集中的区域和应力最大值相对较大.研究成果为进一步揭示双相不锈钢局部失效机理奠定了基础.

Effect of main inclusions on crack initiation in bearing steel in the very high cycle fatigue regime

This work aims to investigate the effect of main inclusions on crack initiation in bearing steel in the very high cycle fatigue（VHCF） regime. The size and type of inclusions in the steel were quantitatively analyzed, and VHCF tests were performed. Some fatigue cracks were found to be initiated in the gaps between inclusions（Al2 O3, Mg O-Al2 O3） and the matrix, while other cracks originated from the interior of inclusions（Ti N, Mn S）. To explain the related mechanism, the tessellated stresses between inclusions and the matrix were calculated and compared with the yield stress of the matrix. Results revealed that the inclusions could be classified into two types under VHCF; of these two, only one type could be regarded as holes. Findings in this research provide a better understanding of how inclusions affect the high cycle fatigue properties of bearing steel.

Ultrahigh cycle fatigue fracture mechanism of high-quality bearing steel obtained through different deoxidation methods

The mechanism of oxide inclusions in fatigue crack initiation in the very-high cycle fatigue(VHCF)regime was clarified by subjecting bearing steels deoxidized by Al(Al-deoxidized steel)and Si(Si-deoxidized steel)to ultrasonic tension-compression fatigue tests(stress ratio,R=−1)and analyzing the characteristics of the detected inclusions.Results show that the main types of inclusions in Si-and Al-deoxidized steels are silicate and calcium aluminate,respectively.The content of calcium aluminate inclusions larger than 15μm in Si-deoxidized steel is lower than that in Al-deoxidized steel,and the difference observed may be attributed to different inclusion generation processes during melting.Despite differences in their cleanliness and total oxygen contents,the Si-and Al-deoxidized steels show similar VHCF lives.The factors causing fatigue failure in these steels reveal distinct differences.Calcium aluminate inclusions are responsible for the cracks in Al-deoxidized steel.By comparison,most fatigue cracks in Si-deoxidized steel are triggered by the inhomogeneity of a steel matrix,which indicates that the damage mechanisms of the steel matrix can be a critical issue for this type of steel.A minor portion of the cracks in Si-deoxidized steel could be attributed to different types of inclusions.The mechanisms of fatigue fracture caused by calcium aluminate and silicate inclusions were further analyzed.Calcium aluminate inclusions first separate from the steel matrix and then trigger crack generation.Silicate inclusions and the steel matrix are closely combined in a fatigue process;thus,these inclusions have mild effects on the fatigue life of bearing steels.Si/Mn deoxidation is an effective method to produce high-quality bearing steel with a long fatigue life and good liquid steel fluidity.

Advances in Studies and Applications of Centroidal Voronoi Tessellations

Centroidal Voronoi tessellations(CVTs) have become a useful tool in many applications ranging from geometric modeling,image and data analysis,and numerical partial differential equations,to problems in physics,astrophysics,chemistry,and biology. In this paper,we briefly review the CVT concept and a few of its generalizations and well-known properties.We then present an overview of recent advances in both mathematical and computational studies and in practical applications of CVTs.Whenever possible,we point out some outstanding issues that still need investigating.

Electron Diffraction Patterns of Plesiotwinning of Fe-Illite

A new selected area electron diffraction (SAED) pattern was found from Fe-illite. The lattice of the ED pattern can be described as a compound tessellation of {3,6}[79{3,6}], which is based on a large coincidence-site lattice (CSL) with oriented crystal associations rotated by a non-crystallographic angle 26°. This phenomenon corresponds to the plesiotwinning of Fe-illite and has been never reported in the nature so far.

10-year fundus tessellation progression and retinal vein occlusion

AIM： To access the 10-year fundus tessellation progression in patients with retinal vein occlusion. METHODS： The Beijing Eye Study 2001/2011 is a populationbased longitudinal study. The study participants underwent a detailed physical and ophthalmic examination. Degree of fundus tessellation was graded by using fundus photographs of the macula and optic disc. Progression of fundus tessellation was calculated by fundus tessellation degree of 2011 minus degree of 2001. Fundus photographs were used for assessment of retinal vein occlusion. RESULTS： The Beijing Eye Study included 4403 subjects in 2001, 3468 subjects was repeated in 2011. Assessment of retinal vein obstruction and fundus tessellation progression were available for 2462 subjects（71.0%）, with 66 subjects fulfilled the diagnosis of retinal vein occlusion. Of the 66 participants, 59 participants with unilateral branch retinal vein occlusion, 5 participants with unilateral central retinal vein occlusion, 1 participant with bilateral branch retinal vein occlusion, and 1 participant with branch retinal vein occlusion in one eye and central retinal vein occlusion in the other eye. Mean degree of peripapillary fundus tessellation progression were significantly higher in the whole retinal vein occlusion group（0.33±0.39, P〈0.001）, central retinal vein occlusion group（0.71±0.8, P=0.025） and branch retinal vein occlusion group（0.29±0.34, P=0.006） than the control group（0.20±0.26）. After adjustment for age, prevalence of tilted disc, change of best corrected visual acuity, axial length, progression of peripapillary fundus tessellation was associated with the presence of retinal vein occlusion（P=0.004; regression coefficient B, 0.094; 95%CI, 0.029, 0.158; standardized coefficient B, 0.056）. As a corollary, after adjusting for smoking duration, systolic blood pressure, anterior corneal curvature, prevalence of RVO was associated with more peripapillary fundus tessellation progression（P〈0.001; regression coefficient B： 1.257; OR： 3.517; 95%CI： 1.777, 6.958）. CONCLUSION： Peripapillary fundus tessellation progresses faster in individuals with retinal vein occlusion. This may reflect the thinning and hypoperfusion of choroid in patients with retinal vein occlusion.

3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing

A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Voronoi Based Discrete Least Squares Meshless Method for Heat Conduction Simulation in Highly Irregular Geometries

A new technique is used in Discrete Least Square Meshfree（DLSM） method to remove the common existing deficiencies of meshfree methods in handling of the problems containing cracks or concave boundaries. An enhanced Discrete Least Squares Meshless method named as VDLSM（Voronoi based Discrete Least Squares Meshless） is developed in order to solve the steady-state heat conduction problem in irregular solid domains including concave boundaries or cracks. Existing meshless methods cannot estimate precisely the required unknowns in the vicinity of the above mentioned boundaries. Conducted researches are limited to domains with regular convex boundaries. To this end, the advantages of the Voronoi tessellation algorithm are implemented. The support domains of the sampling points are determined using a Voronoi tessellation algorithm. For the weight functions, a cubic spline polynomial is used based on a normalized distance variable which can provide a high degree of smoothness near those mentioned above discontinuities. Finally, Moving Least Squares（MLS） shape functions are constructed using a varitional method. This straight-forward scheme can properly estimate the unknowns（in this particular study, the temperatures at the nodal points） near and on the crack faces, crack tip or concave boundaries without need to extra backward corrective procedures, i.e. the iterative calculations for modifying the shape functions of the nodes located near or on these types of the complex boundaries. The accuracy and efficiency of the presented method are investigated by analyzing four particular examples. Obtained results from VDLSM are compared with the available analytical results or with the results of the well-known Finite Elements Method（FEM） when an analytical solution is not available. By comparisons, it is revealed that the proposed technique gives high accuracy for the solution of the steady-state heat conduction problems within cracked domains or domains with concave boundaries and at the same time possesses a high convergence rate which its accuracy is not sensitive to the arrangement of the nodal points. The novelty of this paper is the use of Voronoi concept in determining the weight functions used in the formulation of the MLS type shape functions.

Modified virtual internal bond model based on deformable Voronoi particles

In last time,the series of virtual internal bond model was proposed for solving rock mechanics problems.In these models,the rock continuum is considered as a structure of discrete particles connected by normal and shear springs(bonds).It is well announced that the normal springs structure corresponds to a linear elastic solid with a fixed Poisson ratio,namely,0.25 for threedimensional cases.So the shear springs used to represent the diversity of the Poisson ratio.However,the shearing force calculation is not rotationally invariant and it produce difficulties in application of these models for rock mechanics problems with sufficient displacements.In this letter,we proposed the approach to support the diversity of the Poisson ratio that based on usage of deformable Voronoi cells as set of particles.The edges of dual Delaunay tetrahedralization are considered as structure of normal springs(bonds).The movements of particle’s centers lead to deformation of tetrahedrals and as result to deformation of Voronoi cells.For each bond,there are the corresponded dual face of some Voronoi cell.We can consider the normal bond as some beam and in this case,the appropriate face of Voronoi cell will be a cross section of this beam.If during deformation the Voronoi face was expand,then,according Poisson effect,the length of bond should be decrees.The above mechanism was numerically investigated and we shown that it is acceptable for simulation of elastic behavior in 0.1–0.3 interval of Poisson ratio.Unexpected surprise is that proposed approach give possibility to simulate auxetic materials with negative Poisson’s ratio in interval from–0.5 to–0.1.

An accelerated scheme with high quality mesh based on Lloyd iteration

High quality mesh plays an important role for finite element methods in science computation and numerical simulation.Whether the mesh quality is good or not,to some extent,it determines the calculation results of the accuracy and efficiency.Different from classic Lloyd iteration algorithm which is convergent slowly,a novel accelerated scheme was presented,which consists of two core parts:mesh points replacement and local edges Delaunay swapping.By using it,almost all the equilateral triangular meshes can be generated based on centroidal Voronoi tessellation(CVT).Numerical tests show that it is significantly effective with time consuming decreasing by 40%.Compared with other two types of regular mesh generation methods,CVT mesh demonstrates that higher geometric average quality increases over 0.99.

Edge-Weighted Centroidal Voronoi Tessellations

Most existing applications of centroidal Voronoi tessellations(CVTs) lack consideration of the length of the cluster boundaries.In this paper we propose a new model and algorithms to produce segmentations which would minimize the total energy—a sum of the classic CVT energy and the weighted length of cluster boundaries.To distinguish it with the classic CVTs,we call it an Edge-Weighted CVT(EWCVT).The concept of EWCVT is expected to build a mathematical base for all CVT related data classifications with requirement of smoothness of the cluster boundaries.The EWCVT method is easy in implementation,fast in computation,and natural for any number of clusters.

LED Adaptive Deployment Optimization in Indoor VLC Networks

Driven by the continuous penetration of high data rate services and applications,a large amount of unregulated visible light spectrum is used for communication to fully meet the needs of 6th generation(6G)mobile technologies.Visible light communication(VLC)faces many challenges as a solution that complements existing radio frequency(RF)networks.This paper studies the optimal configuration of LEDs in indoor environments under the constraints of illumination and quality of experience(QoE).Based on the Voronoi tessellation(VT)and centroidal Voronoi tessellation(CVT)theory,combined with the Lloyd’s algorithm,we propose two approaches for optimizing LED deployments to meet the illumination and QoE requirements of all users.Focusing on(i)the minimization of the number of LEDs to be installed in order to meet illumination and average QoE constraints,and(ii)the maximization of the average QoE of users to be served with a fixed number of LEDs.Monte Carlo simulations are carried out for different user distribution compared with hexagonal,square and VT deployment.The simulation results illustrate that under the same conditions,the proposed deployment approach can provide less LEDs and achieve better QoE performance.

Inferring seismic hazard in Sichuan-Yunnan region of China based on the modern earthquake catalogue (1980-2019)

Based on the modern earthquake catalogue,the incomplete centroidal voronoi tessellation(ICVT)method was used in this study to estimate the seismic hazard in Sichuan-Yunnan region of China.We calculated spatial distributions of the total seismic hazard and background seismic hazard in this area.The Bayesian delaunay tessellation smoothing method put forward by Ogata was used to calculate the spatial distributions of b-value.The results show that seismic hazards in Sichuan-Yunan region are high,and areas with relatively high hazard values are distributed along the main faults,while seismic hazards in Sichuan basin are relatively low.