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.

Convergent Adaptive Finite Element Method Based on Centroidal Voronoi Tessellations and Superconvergence

We present a novel adaptive finite element method(AFEM)for elliptic equations which is based upon the Centroidal Voronoi Tessellation(CVT)and superconvergent gradient recovery.The constructions of CVT and its dual Centroidal Voronoi Delaunay Triangulation(CVDT)are facilitated by a localized Lloyd iteration to produce almost equilateral two dimensional meshes.Working with finite element solutions on such high quality triangulations,superconvergent recovery methods become particularly effective so that asymptotically exact a posteriori error estimations can be obtained.Through a seamless integration of these techniques,a convergent adaptive procedure is developed.As demonstrated by the numerical examples,the new AFEM is capable of solving a variety of model problems and has great potential in practical applications.

Truncated Newton-Based Multigrid Algorithm for Centroidal Voronoi Diagram Calculation

In a variety of modern applications there arises a need to tessellate the domain into representative regions,called Voronoi cells.A particular type of such tessellations,called centroidal Voronoi tessellations or CVTs,are in big demand due to their optimality properties important for many applications.The availability of fast and reliable algorithms for their construction is crucial for their successful use in practical settings.This paper introduces a new multigrid algorithm for constructing CVTs that is based on the MG/Opt algorithm that was originally designed to solve large nonlinear optimization problems.Uniform convergence of the new method and its speedup comparing to existing techniques are demonstrated for linear and nonlinear densities for several 1d and 2d problems,and O(k)complexity estimation is provided for a problem with k generators.

An Improved Cooperative Team Spraying Control of a Diffusion Process With a Moving or Static Pollution Source

This paper is concerned with a control problem of a diffusion process with the help of static mesh sensor networks in a certain region of interest and a team of networked mobile actuators carrying chemical neutralizers.The major contribution of this paper can be divided into three parts:the first is the construction of a cyber-physical system framework based on centroidal Voronoi tessellations(CVTs),the second is the convergence analysis of the actuators location,and the last is a novel proportional integral(PI)control method for actuator motion planning and neutralizing control(e.g.,spraying)of a diffusion process with a moving or static pollution source,which is more effective than a proportional(P)control method.An optimal spraying control cost function is constructed.Then,the minimization problem of the spraying amount is addressed.Moreover,a new CVT algorithm based on the novel PI control method,henceforth called PI-CVT algorithm,is introduced together with the convergence analysis of the actuators location via a PI control law.Finally,a modified simulation platform called diffusion-mobile-actuators-sensors-2-dimension-proportional integral derivative(Diff-MAS2D-PID)is illustrated.In addition,a numerical simulation example for the diffusion process is presented to verify the effectiveness of our proposed controllers.

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.

A unified framework for isotropic meshing based on narrowband Euclidean distance transformation

In this paper, we propose a simpleyet-effective method for isotropic meshing relying on Euclidean distance transformation based centroidal Voronoi tessellation(CVT). Our approach improves the performance and robustness of computing CVT on curved domains while simultaneously providing highquality output meshes. While conventional extrinsic methods compute CVTs in the entire volume bounded by the input model, we restrict the computation to a 3D shell of user-controlled thickness. Taking voxels which contain surface samples as sites, we compute the exact Euclidean distance transform on the GPU. Our algorithm is parallel and memory-efficient,and can construct the shell space for resolutions up to 20483 at interactive speed. The 3D centroidal Voronoi tessellation and restricted Voronoi diagrams are also computed efficiently on the GPU. Since the shell space can bridge holes and gaps smaller than a certain tolerance, and tolerate non-manifold edges and degenerate triangles, our algorithm can handle models with such defects, which typically cause conventional remeshing methods to fail. Our method can process implicit surfaces, polyhedral surfaces, and point clouds in a unified framework. Computational results show that our GPU-based isotropic meshing algorithm produces results comparable to state-ofthe-art techniques, but is significantly faster than conventional CPU-based implementations.