The multi-resolution adaptive grids method is proposed to solve the problems of inefficiency in the previous grid-based methods,and it can be used in clouds simulation as well as the interactive simulation between obj...The multi-resolution adaptive grids method is proposed to solve the problems of inefficiency in the previous grid-based methods,and it can be used in clouds simulation as well as the interactive simulation between objects and clouds.Oriented bounding box(OBB)hierarchical trees of objects are established,and the resolutions of global and local grids can be selected automatically.The motion equations of fluid dynamics are simplified.Upwind difference is applied to ensure the stability of the simulation process during the discrete process of partial differential equations.To solve the speed problem of existed phase functions,the improved phase function is applied to the illumination calculation of clouds.Experimental results show that the proposed methods can promote the simulation efficiency and meet the need for the simulation of large-scale clouds scene.Real-time rendering of clouds and the interaction between clouds and objects have been realized without preprocessing stage.展开更多
High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of th...High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.展开更多
In this paper the algebraic multi-grid principle is applied to the multilevel moment method, which makes the new multilevel method easier to implement and more adaptive to structure. Moreover, the error spectrum is an...In this paper the algebraic multi-grid principle is applied to the multilevel moment method, which makes the new multilevel method easier to implement and more adaptive to structure. Moreover, the error spectrum is analyzed, and the reason why conjugate gradient iteration is not a good relaxation scheme for multi-grid algorithm is explored. The numerical results show that our algebraic block Gauss Seidel multi-grid algorithm is very effective.展开更多
In order to break through the limitationof thelatitude/longitudegrid and hexagon grid, a new subdivision unit, Half-honeycomb Trapezoid, is proposed. Based on the summarization of the geometric properties and subdivis...In order to break through the limitationof thelatitude/longitudegrid and hexagon grid, a new subdivision unit, Half-honeycomb Trapezoid, is proposed. Based on the summarization of the geometric properties and subdivision performance of Half-honeycomb Trapezoid, a new discrete global topographic grid system is established, and its compatibility with hexagonal grid is analyzed. At last, the visualization of multi-resolution global grid is achieved.展开更多
A model of a hypertorus communication grid has been constructed in the form of an infinite Petri net. A grid cell represents either a packet switching device or a bioplast cell. A parametric expression is obtained to ...A model of a hypertorus communication grid has been constructed in the form of an infinite Petri net. A grid cell represents either a packet switching device or a bioplast cell. A parametric expression is obtained to allow a finite specification of an infinite Petri net. To prove properties of an ideal communication protocol, we derive an infinite Diophantine system of equations from it, which is subsequently solved. Then we present the programs htgen and ht-mcrl2-gen, developed in the C language, which generate Petri net and process algebra models of a hypertorus with a given number of dimensions and grid size. These are the inputs for the respective modeling tools Tina and mCRL2, which provide model visualization, step simulation, state space generation and reduction, and structural analysis techniques. Benchmarks to compare the two approaches are obtained. An ad-hoc induction-like technique on invariants,obtained for a series of generated models, allows the calculation of a solution of the Diophantine system in a parametric form.It is proven that the basic solutions of the infinite system have been found and that the infinite Petri net is bounded and conservative. Some remarks regarding liveness and liveness enforcing techniques are also presented.展开更多
Quality and robustness of grid deformation is of the most importance in the field of aircraft design, and grid in high quality is essential for improving the precision of numerical simulation. In order to maintain the...Quality and robustness of grid deformation is of the most importance in the field of aircraft design, and grid in high quality is essential for improving the precision of numerical simulation. In order to maintain the orthogonality of deformed grid, the displacement of grid points is divided into rotational and translational parts in this paper, and inverse distance weighted interpolation is used to transfer the changing location from boundary grid to the spatial grid. Moreover, the deformation of rotational part is implemented in combination with the exponential space mapping that improves the certainty and stability of quaternion interpolation. Furthermore, the new grid deformation technique named ‘‘layering blend deformation＇＇ is built based on the basic quaternion technique, which combines the layering arithmetic with transfinite interpolation（TFI） technique. Then the proposed technique is applied in the movement of airfoil, parametric modeling, and the deformation of complex configuration, in which the robustness of grid quality is tested. The results show that the new method has the capacity to deal with the problems with large deformation, and the ‘‘layering blend deformation＇＇ improves the efficiency and quality of the basic quaternion deformation method significantly.展开更多
Discrete Global Grid Systems(DGGSs) are spatial references that use a hierarchical tessellation of cells to partition and address the entire globe. They provide an organizational structure that permits fast integratio...Discrete Global Grid Systems(DGGSs) are spatial references that use a hierarchical tessellation of cells to partition and address the entire globe. They provide an organizational structure that permits fast integration between multiple sources of large and variable geospatial data sufficient for visualization and analysis. Despite a significant body of research supporting hexagonal DGGSs as the superior choice, the application thereof has been hindered owing in part to the lack of a rational hierarchy with an efficient addressing system. This paper presents an algebraic model of encoding scheme for the Aperture 3 Hexagonal(A3H) DGGS. Firstly, the definition of a grid cell, which is composed of vertices, edges, and a center, is introduced to describe fundamental elements of grids. Secondly, by identifying the grid cell with its center, this paper proves that cell centers at different levels can be represented exactly using a mixed positional number system in the complex plane through the recursive geometric relationship between two successive levels, which reveals that grid cells are essentially special complex radix numbers. Thirdly, it is shown that through the recursive geometric relationship of successive odd or even levels, the mixed positional number system can also be applied to uniquely represent cell centers at different levels under specific constraint conditions, according to which the encoding scheme is designed. Finally, it is shown that by extending the scheme to 20 triangular faces of the regular icosahedron,multi-resolution grids on closed surfaces of the icosahedron are addressed perfectly. Contrast experiments show that the proposed encoding scheme has the advantages of theoretical rigor and high programming efficiency and that the efficiency of cross-face adjacent cell searching is 242.9 times that of a similar scheme. Moreover, the proposed complex radix number representation is an ideal formalized description tool for grid systems. The research ideas introduced herein can be used to create a universal theoretical framework for DGGSs.展开更多
Analyzing network robustness under various circumstances is generally regarded as a challenging problem.Robustness against failure is one of the essential properties of large-scale dynamic network systems such as powe...Analyzing network robustness under various circumstances is generally regarded as a challenging problem.Robustness against failure is one of the essential properties of large-scale dynamic network systems such as power grids,transportation systems,communication systems,and computer networks.Due to the network diversity and complexity,many topological features have been proposed to capture specific system properties.For power grids,a popular process for improving a network’s structural robustness is via the topology design.However,most of existing methods focus on localized network metrics,such as node connectivity and edge connectivity,which do not encompass a global perspective of cascading propagation in a power grid.In this paper,we use an informative global metric algebraic connectivity because it is sensitive to the connectedness in a broader spectrum of graphs.Our process involves decreasing the average propagation in a power grid by minimizing the increase in its algebraic connectivity.We propose a topology-based greedy strategy to optimize the robustness of the power grid.To evaluate the network robustness,we calculate the average propagation using MATCASC to simulate cascading line outages in power grids.Experimental results illustrate that our proposed method outperforms existing techniques.展开更多
基金supported by the National Natural Science Foundation of China(No.61102167)
文摘The multi-resolution adaptive grids method is proposed to solve the problems of inefficiency in the previous grid-based methods,and it can be used in clouds simulation as well as the interactive simulation between objects and clouds.Oriented bounding box(OBB)hierarchical trees of objects are established,and the resolutions of global and local grids can be selected automatically.The motion equations of fluid dynamics are simplified.Upwind difference is applied to ensure the stability of the simulation process during the discrete process of partial differential equations.To solve the speed problem of existed phase functions,the improved phase function is applied to the illumination calculation of clouds.Experimental results show that the proposed methods can promote the simulation efficiency and meet the need for the simulation of large-scale clouds scene.Real-time rendering of clouds and the interaction between clouds and objects have been realized without preprocessing stage.
文摘High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.
基金Supported by the Natlonal Natural Science Foundation of China
文摘In this paper the algebraic multi-grid principle is applied to the multilevel moment method, which makes the new multilevel method easier to implement and more adaptive to structure. Moreover, the error spectrum is analyzed, and the reason why conjugate gradient iteration is not a good relaxation scheme for multi-grid algorithm is explored. The numerical results show that our algebraic block Gauss Seidel multi-grid algorithm is very effective.
基金Supported by Key Scientific and Technological Project of Anhui Province(No.1401b042009)Provincal Natural Science Foundation of the Higher Education Institutions of Anhui(No.KJ2014ZD27)
文摘In order to break through the limitationof thelatitude/longitudegrid and hexagon grid, a new subdivision unit, Half-honeycomb Trapezoid, is proposed. Based on the summarization of the geometric properties and subdivision performance of Half-honeycomb Trapezoid, a new discrete global topographic grid system is established, and its compatibility with hexagonal grid is analyzed. At last, the visualization of multi-resolution global grid is achieved.
文摘A model of a hypertorus communication grid has been constructed in the form of an infinite Petri net. A grid cell represents either a packet switching device or a bioplast cell. A parametric expression is obtained to allow a finite specification of an infinite Petri net. To prove properties of an ideal communication protocol, we derive an infinite Diophantine system of equations from it, which is subsequently solved. Then we present the programs htgen and ht-mcrl2-gen, developed in the C language, which generate Petri net and process algebra models of a hypertorus with a given number of dimensions and grid size. These are the inputs for the respective modeling tools Tina and mCRL2, which provide model visualization, step simulation, state space generation and reduction, and structural analysis techniques. Benchmarks to compare the two approaches are obtained. An ad-hoc induction-like technique on invariants,obtained for a series of generated models, allows the calculation of a solution of the Diophantine system in a parametric form.It is proven that the basic solutions of the infinite system have been found and that the infinite Petri net is bounded and conservative. Some remarks regarding liveness and liveness enforcing techniques are also presented.
文摘Quality and robustness of grid deformation is of the most importance in the field of aircraft design, and grid in high quality is essential for improving the precision of numerical simulation. In order to maintain the orthogonality of deformed grid, the displacement of grid points is divided into rotational and translational parts in this paper, and inverse distance weighted interpolation is used to transfer the changing location from boundary grid to the spatial grid. Moreover, the deformation of rotational part is implemented in combination with the exponential space mapping that improves the certainty and stability of quaternion interpolation. Furthermore, the new grid deformation technique named ‘‘layering blend deformation＇＇ is built based on the basic quaternion technique, which combines the layering arithmetic with transfinite interpolation（TFI） technique. Then the proposed technique is applied in the movement of airfoil, parametric modeling, and the deformation of complex configuration, in which the robustness of grid quality is tested. The results show that the new method has the capacity to deal with the problems with large deformation, and the ‘‘layering blend deformation＇＇ improves the efficiency and quality of the basic quaternion deformation method significantly.
基金supported by the National Natural Science Foundation of China (Grant No. 41671410)the Postdoctoral Science Foundation of China (Grant No. 2013T60161)the Excellent Young Scholar Foundation of Information Engineering University (Grant No. 2016610802)
文摘Discrete Global Grid Systems(DGGSs) are spatial references that use a hierarchical tessellation of cells to partition and address the entire globe. They provide an organizational structure that permits fast integration between multiple sources of large and variable geospatial data sufficient for visualization and analysis. Despite a significant body of research supporting hexagonal DGGSs as the superior choice, the application thereof has been hindered owing in part to the lack of a rational hierarchy with an efficient addressing system. This paper presents an algebraic model of encoding scheme for the Aperture 3 Hexagonal(A3H) DGGS. Firstly, the definition of a grid cell, which is composed of vertices, edges, and a center, is introduced to describe fundamental elements of grids. Secondly, by identifying the grid cell with its center, this paper proves that cell centers at different levels can be represented exactly using a mixed positional number system in the complex plane through the recursive geometric relationship between two successive levels, which reveals that grid cells are essentially special complex radix numbers. Thirdly, it is shown that through the recursive geometric relationship of successive odd or even levels, the mixed positional number system can also be applied to uniquely represent cell centers at different levels under specific constraint conditions, according to which the encoding scheme is designed. Finally, it is shown that by extending the scheme to 20 triangular faces of the regular icosahedron,multi-resolution grids on closed surfaces of the icosahedron are addressed perfectly. Contrast experiments show that the proposed encoding scheme has the advantages of theoretical rigor and high programming efficiency and that the efficiency of cross-face adjacent cell searching is 242.9 times that of a similar scheme. Moreover, the proposed complex radix number representation is an ideal formalized description tool for grid systems. The research ideas introduced herein can be used to create a universal theoretical framework for DGGSs.
基金supported by the National Natural Science Foundation of China(No.U1866602)the National Key R&D Program of China(Nos.2019YFB1600700 and 2018AAA0101505)。
文摘Analyzing network robustness under various circumstances is generally regarded as a challenging problem.Robustness against failure is one of the essential properties of large-scale dynamic network systems such as power grids,transportation systems,communication systems,and computer networks.Due to the network diversity and complexity,many topological features have been proposed to capture specific system properties.For power grids,a popular process for improving a network’s structural robustness is via the topology design.However,most of existing methods focus on localized network metrics,such as node connectivity and edge connectivity,which do not encompass a global perspective of cascading propagation in a power grid.In this paper,we use an informative global metric algebraic connectivity because it is sensitive to the connectedness in a broader spectrum of graphs.Our process involves decreasing the average propagation in a power grid by minimizing the increase in its algebraic connectivity.We propose a topology-based greedy strategy to optimize the robustness of the power grid.To evaluate the network robustness,we calculate the average propagation using MATCASC to simulate cascading line outages in power grids.Experimental results illustrate that our proposed method outperforms existing techniques.
