Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
Ion optics are crucial components of ion thrusters and the study of the different ion extraction solutions used since the beginning of the electric propulsion era is essential to understand the evolution of ion engine...Ion optics are crucial components of ion thrusters and the study of the different ion extraction solutions used since the beginning of the electric propulsion era is essential to understand the evolution of ion engines. This work describes ion engine grids＇ main functions, parameters and issues related to thermal expansion and sputter erosion, and then introduces a review of ion optics used for significant launched and tested ion thrusters since 1970. Configurations, geometries, materials and fabrication methods are analyzed to understand when typical ion thrusters use two or three grids, what are the thicknesses and aperture sizes of the screen, accelerator and decelerator grids, why molybdenum and carbon-based materials such as pyrolytic graphite and carbon–carbon are the best available options for ion optics and what is the manufacturing method for each material.展开更多
With the global economy integration and progress in energy transformation,it has become a general trend to surpass national boundaries to achieve wider and optimal energy resource allocations.Consequently,there is a c...With the global economy integration and progress in energy transformation,it has become a general trend to surpass national boundaries to achieve wider and optimal energy resource allocations.Consequently,there is a critical n eed to adopt scie ntific approaches in assessi ng cross-border power grid interconnection projects.First,con sidering the promotion of large-scale renewable energy resources and improvements in system adequacy,a comprehensive assessment index system,including costs,socio-economic benefits,environmental benefits,and technical benefits,is established in this study.Second,a synthetic assessment framework is proposed for cross-border power grid interconnection projects based on the index system comprising cost-benefit analysis,with market and network simulations,iterative methods for indicator weight evaluation,and technique for order preferenee by similarity to an ideal solution(TOPSIS)method for the project rankings.Fin ally,by assessi ng and comparing three cross-border projects betwee n Europe and Asia,the proposed index system and assessment framework have been proved to be effective and feasible;the results of this system can thus support investment decision-making related to such projects in the future.展开更多
In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of fini...In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.展开更多
With the accelerating urbanization process,the load demand of urban power grids is constantly increasing,giving rise to a batch of ultra-large urban power grids featuring large electricity demand,dense load distributi...With the accelerating urbanization process,the load demand of urban power grids is constantly increasing,giving rise to a batch of ultra-large urban power grids featuring large electricity demand,dense load distribution,and tight construction land constraints.This paper establishes a network planning method for urban power grids based on series reactors and MMC-MTEDC,focusing on four aspects:short-circuit current suppression,accommodation of external power supply,flexible inter-regional power support,and voltage stability enhancement in load centers.It proposes key indicators including node short-circuit current margin,line thermal stability margin,maximum fault-induced regional power loss,and voltage recovery time,thereby constructing an evaluation system for MMT-MTEDC network planning in urban power grids.Based on the Shenzhen power grid planning data,simulations using DSP software reveal that series reactors reduce short-circuit current by up to 5.0%,while the MMC-MTEDC system enhances node short-circuit margins by 4.212.9%and shortens voltage recovery time by 19.8%.Additionally,the MMC-MTEDC system maintains 3.34-6.76 percentage points higher thermal stability margins than conventional AC systems and enables complete avoidance of external power curtailment during N-2 faults via power reallocation between terminals.Compared with traditional AC or point-to-point HVDC schemes,the proposed hybrid planning method better adapts to the spatial and reliability demands of ultra-large receiving-end grids.This methodology provides practical insights into coordinated AC/DC development under high load density and strong external power reliance.Future work will extend the approach to include electromagnetic transient constraints and lightweight MMC station designs for urban applications.展开更多
A baroclinic typhoon model with a moving multi--nested grid is applied in marine environmental forecasts. This paper describes the numerical methods of the model including governing equations, finite differencing, spl...A baroclinic typhoon model with a moving multi--nested grid is applied in marine environmental forecasts. This paper describes the numerical methods of the model including governing equations, finite differencing, split scheme and time integration.展开更多
A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils.The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leave...A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils.The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leaves only small region adjacent to the aerodynamic bodies to be filled with the cloud of points used in the gridless methods,which results in a better combination of the computational efficiency of the Cartesian grid and the flexibility of the gridless method in handling complex geometries.The clouds of points in the local gridless region are implemented in an anisotropic way according to the features of the thin boundary layer of the viscous flows over the airfoils,and the clouds of points at the vicinity of the interface between the grid and the gridless regions are also controlled by using an adaptive refinement technique during the generation of the unstructured Cartesian grid.An implementation of the resulting hybrid method is presented for solving two-dimensional compressible Navier-Stokes(NS)equations.The simulations of the viscous flows over a RAE2822airfoil or a two-element airfoil are successfully carried out,and the obtained results agree well with the available experimental data.展开更多
Combined the purpose and requirements of security and stability economic operations of the smart grid, a more comprehensive, risk indicator system of smart grid is established from five aspects of the smart grid strat...Combined the purpose and requirements of security and stability economic operations of the smart grid, a more comprehensive, risk indicator system of smart grid is established from five aspects of the smart grid strategic risk, external risk, financial risk, compliance risk, operational risks. On this basis, it is conducted smart grid risk assessments by integrated of the use of Borda sequence value method of the original risk matrix and AHP (hierarchical analysis method), in order to assess risk facing the smart grid and enterprise development more comprehensively, objectively and systematically and conveniently, and prompt companies continue to reduce risk, improve economic efficiency, continuously improve the targeted improvement measures and continue to improve the level of grid development.展开更多
We present high order accurate numerical methods for the wave equation that combines efficient Hermite methods with geometrically flexible discontinuous Galerkin methods by using overset grids.Near boundaries we use t...We present high order accurate numerical methods for the wave equation that combines efficient Hermite methods with geometrically flexible discontinuous Galerkin methods by using overset grids.Near boundaries we use thin boundary fitted curvilinear grids and in the volume we use Cartesian grids so that the computational complexity of the solvers approaches a structured Cartesian Hermite method.Unlike many other overset methods we do not need to add artificial dissipation but we find that the built-in dissipation of the Hermite and discontinuous Galerkin methods is sufficient to maintain the stability.By numerical experiments we demonstrate the stability,accuracy,efficiency,and the applicability of the methods to forward and inverse problems.展开更多
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.展开更多
The multi-grid method has been known as an efficient iterative method for the linear systems and nonlinear systems that arise from finite difference approximations for partial differential equations. In this paper, th...The multi-grid method has been known as an efficient iterative method for the linear systems and nonlinear systems that arise from finite difference approximations for partial differential equations. In this paper, the multigrid method is extended to the application of solving integral equations which appear in electromagnetic scattering problems. The diakoptic theory is used for this purpose. Compared with other methods, the numerical results show that the multigrid method is powerful to solve electromagnetic scattering problems and can be used to compute electromagnetic scattering problems with electrically large bodies and complex structures.展开更多
In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the m...In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.展开更多
For a complex flow about multi-element airfoils a mixed grid method is set up. C-type grids are produced on each element′s body and in their wakes at first, O-type grids are given in the outmost area, and H-type grid...For a complex flow about multi-element airfoils a mixed grid method is set up. C-type grids are produced on each element′s body and in their wakes at first, O-type grids are given in the outmost area, and H-type grids are used in middle additional areas. An algebra method is used to produce the initial grids in each area. And the girds are optimized by elliptical differential equation method. Then C-O-H zonal patched grids around multi-element airfoils are produced automatically and efficiently. A time accurate finite-volume integration method is used to solve the compressible laminar and turbulent Navier-Stokes (N-S) equations on the grids. Computational results prove the method to be effective.展开更多
Rasterization is a conversion process accompanied with information loss, which includes the loss of features＇ shape, structure, position, attribute and so on. Two chief factors that affect estimating attribute accura...Rasterization is a conversion process accompanied with information loss, which includes the loss of features＇ shape, structure, position, attribute and so on. Two chief factors that affect estimating attribute accuracy loss in rasterization are grid cell size and evaluating method. That is, attribute accuracy loss in rasterization has a close relationship with grid cell size; besides, it is also influenced by evaluating methods. Therefore, it is significant to analyze these two influencing factors comprehensively. Taking land cover data of Sichuan at the scale of 1：250,000 in 2005 as a case, in view of data volume and its processing time of the study region, this study selects 16 spatial scales from 600 m to 30 km, uses rasterizing method based on the Rule of Maximum Area （RMA） in ArcGIS and two evaluating methods of attribute accuracy loss, which are Normal Analysis Method （NAM） and a new Method Based on Grid Cell （MBGC）, respectively, and analyzes the scale effect of attribute （it is area here） accuracy loss at 16 different scales by these two evaluating methods comparatively. The results show that： （1） At the same scale, average area accuracy loss of the entire study region evaluated by MBGC is significantly larger than the one estimated using NAM. Moreover, this discrepancy between the two is obvious in the range of 1 km to 10 km. When the grid cell is larger than 10 km, average area accuracy losses calculated by the two evaluating methods are stable, even tended to parallel. （2） MBGC can not only estimate RMA rasterization attribute accuracy loss accurately, but can express the spatial distribution of the loss objectively. （3） The suitable scale domain for RMA rasterization of land cover data of Sichuan at the scale of 1：250,000 in 2005 is better equal to or less than 800 m, in which the data volume is favorable and the processina time is not too Iona. as well as the area accuracv loss is less than 2.5%.展开更多
A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method o...A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.展开更多
The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Ro...The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.展开更多
On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the...On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software ＂Auto mapping system of submarine topography and geomorphology MBChart＂. This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China＇ s next Five-year Plan.展开更多
For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low vo...For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low voltage problems were highlighted in the rural power grid due to the characteristics of rural power grid. Using the distribution network flow calculation method, we evaluated the low voltage problems of the rural power grid which belongs to Chongqing Jiangbei Power Company. In addition, we collected the data of distribution transformers in electricity consumption peak period. Some practical management strategies were proposed by the analysis and evaluation of potential and appeared low voltage problems.展开更多
The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then sin...The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.展开更多
This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance ...This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.展开更多
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
文摘Ion optics are crucial components of ion thrusters and the study of the different ion extraction solutions used since the beginning of the electric propulsion era is essential to understand the evolution of ion engines. This work describes ion engine grids＇ main functions, parameters and issues related to thermal expansion and sputter erosion, and then introduces a review of ion optics used for significant launched and tested ion thrusters since 1970. Configurations, geometries, materials and fabrication methods are analyzed to understand when typical ion thrusters use two or three grids, what are the thicknesses and aperture sizes of the screen, accelerator and decelerator grids, why molybdenum and carbon-based materials such as pyrolytic graphite and carbon–carbon are the best available options for ion optics and what is the manufacturing method for each material.
基金the Science and Technology Project of Global Energy Interconnection Group Co.,Ltd.(No.524500180014).
文摘With the global economy integration and progress in energy transformation,it has become a general trend to surpass national boundaries to achieve wider and optimal energy resource allocations.Consequently,there is a critical n eed to adopt scie ntific approaches in assessi ng cross-border power grid interconnection projects.First,con sidering the promotion of large-scale renewable energy resources and improvements in system adequacy,a comprehensive assessment index system,including costs,socio-economic benefits,environmental benefits,and technical benefits,is established in this study.Second,a synthetic assessment framework is proposed for cross-border power grid interconnection projects based on the index system comprising cost-benefit analysis,with market and network simulations,iterative methods for indicator weight evaluation,and technique for order preferenee by similarity to an ideal solution(TOPSIS)method for the project rankings.Fin ally,by assessi ng and comparing three cross-border projects betwee n Europe and Asia,the proposed index system and assessment framework have been proved to be effective and feasible;the results of this system can thus support investment decision-making related to such projects in the future.
基金Project supported by the National Natural Science Foundation of China(Nos.11671157 and11826212)
文摘In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.
基金Shenzhen Power SupplyCo.,Ltd.Grant number 090000KC24040028.
文摘With the accelerating urbanization process,the load demand of urban power grids is constantly increasing,giving rise to a batch of ultra-large urban power grids featuring large electricity demand,dense load distribution,and tight construction land constraints.This paper establishes a network planning method for urban power grids based on series reactors and MMC-MTEDC,focusing on four aspects:short-circuit current suppression,accommodation of external power supply,flexible inter-regional power support,and voltage stability enhancement in load centers.It proposes key indicators including node short-circuit current margin,line thermal stability margin,maximum fault-induced regional power loss,and voltage recovery time,thereby constructing an evaluation system for MMT-MTEDC network planning in urban power grids.Based on the Shenzhen power grid planning data,simulations using DSP software reveal that series reactors reduce short-circuit current by up to 5.0%,while the MMC-MTEDC system enhances node short-circuit margins by 4.212.9%and shortens voltage recovery time by 19.8%.Additionally,the MMC-MTEDC system maintains 3.34-6.76 percentage points higher thermal stability margins than conventional AC systems and enables complete avoidance of external power curtailment during N-2 faults via power reallocation between terminals.Compared with traditional AC or point-to-point HVDC schemes,the proposed hybrid planning method better adapts to the spatial and reliability demands of ultra-large receiving-end grids.This methodology provides practical insights into coordinated AC/DC development under high load density and strong external power reliance.Future work will extend the approach to include electromagnetic transient constraints and lightweight MMC station designs for urban applications.
基金supported by the Chinese National Research Program of Science and Technology under Project! 85-903-03-04.
文摘A baroclinic typhoon model with a moving multi--nested grid is applied in marine environmental forecasts. This paper describes the numerical methods of the model including governing equations, finite differencing, split scheme and time integration.
基金Supported by the National Natural Science Foundation of China(11172134)the Funding of Jiangsu Innovation Program for Graduate Education(CXZZ110192)the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils.The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leaves only small region adjacent to the aerodynamic bodies to be filled with the cloud of points used in the gridless methods,which results in a better combination of the computational efficiency of the Cartesian grid and the flexibility of the gridless method in handling complex geometries.The clouds of points in the local gridless region are implemented in an anisotropic way according to the features of the thin boundary layer of the viscous flows over the airfoils,and the clouds of points at the vicinity of the interface between the grid and the gridless regions are also controlled by using an adaptive refinement technique during the generation of the unstructured Cartesian grid.An implementation of the resulting hybrid method is presented for solving two-dimensional compressible Navier-Stokes(NS)equations.The simulations of the viscous flows over a RAE2822airfoil or a two-element airfoil are successfully carried out,and the obtained results agree well with the available experimental data.
文摘Combined the purpose and requirements of security and stability economic operations of the smart grid, a more comprehensive, risk indicator system of smart grid is established from five aspects of the smart grid strategic risk, external risk, financial risk, compliance risk, operational risks. On this basis, it is conducted smart grid risk assessments by integrated of the use of Borda sequence value method of the original risk matrix and AHP (hierarchical analysis method), in order to assess risk facing the smart grid and enterprise development more comprehensively, objectively and systematically and conveniently, and prompt companies continue to reduce risk, improve economic efficiency, continuously improve the targeted improvement measures and continue to improve the level of grid development.
基金This work was supported in part by the National Science Foundation under Grant NSF-1913076.Any opinions,findings,and conclusions or recommendations expressed in this material are those of the author(s)and do not necessarily reflect the views of the National Science Foundation.
文摘We present high order accurate numerical methods for the wave equation that combines efficient Hermite methods with geometrically flexible discontinuous Galerkin methods by using overset grids.Near boundaries we use thin boundary fitted curvilinear grids and in the volume we use Cartesian grids so that the computational complexity of the solvers approaches a structured Cartesian Hermite method.Unlike many other overset methods we do not need to add artificial dissipation but we find that the built-in dissipation of the Hermite and discontinuous Galerkin methods is sufficient to maintain the stability.By numerical experiments we demonstrate the stability,accuracy,efficiency,and the applicability of the methods to forward and inverse problems.
基金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.
文摘The multi-grid method has been known as an efficient iterative method for the linear systems and nonlinear systems that arise from finite difference approximations for partial differential equations. In this paper, the multigrid method is extended to the application of solving integral equations which appear in electromagnetic scattering problems. The diakoptic theory is used for this purpose. Compared with other methods, the numerical results show that the multigrid method is powerful to solve electromagnetic scattering problems and can be used to compute electromagnetic scattering problems with electrically large bodies and complex structures.
文摘In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.
文摘For a complex flow about multi-element airfoils a mixed grid method is set up. C-type grids are produced on each element′s body and in their wakes at first, O-type grids are given in the outmost area, and H-type grids are used in middle additional areas. An algebra method is used to produce the initial grids in each area. And the girds are optimized by elliptical differential equation method. Then C-O-H zonal patched grids around multi-element airfoils are produced automatically and efficiently. A time accurate finite-volume integration method is used to solve the compressible laminar and turbulent Navier-Stokes (N-S) equations on the grids. Computational results prove the method to be effective.
基金The Independent Research of the State Key Laboratory of Resource and Environmental Information System,No.O88RA100SAThe Third Innovative and Cutting-edge Projects of Institute of Geographic Sciences andNatural Resources Research, CAS, No.O66U0309SZ
文摘Rasterization is a conversion process accompanied with information loss, which includes the loss of features＇ shape, structure, position, attribute and so on. Two chief factors that affect estimating attribute accuracy loss in rasterization are grid cell size and evaluating method. That is, attribute accuracy loss in rasterization has a close relationship with grid cell size; besides, it is also influenced by evaluating methods. Therefore, it is significant to analyze these two influencing factors comprehensively. Taking land cover data of Sichuan at the scale of 1：250,000 in 2005 as a case, in view of data volume and its processing time of the study region, this study selects 16 spatial scales from 600 m to 30 km, uses rasterizing method based on the Rule of Maximum Area （RMA） in ArcGIS and two evaluating methods of attribute accuracy loss, which are Normal Analysis Method （NAM） and a new Method Based on Grid Cell （MBGC）, respectively, and analyzes the scale effect of attribute （it is area here） accuracy loss at 16 different scales by these two evaluating methods comparatively. The results show that： （1） At the same scale, average area accuracy loss of the entire study region evaluated by MBGC is significantly larger than the one estimated using NAM. Moreover, this discrepancy between the two is obvious in the range of 1 km to 10 km. When the grid cell is larger than 10 km, average area accuracy losses calculated by the two evaluating methods are stable, even tended to parallel. （2） MBGC can not only estimate RMA rasterization attribute accuracy loss accurately, but can express the spatial distribution of the loss objectively. （3） The suitable scale domain for RMA rasterization of land cover data of Sichuan at the scale of 1：250,000 in 2005 is better equal to or less than 800 m, in which the data volume is favorable and the processina time is not too Iona. as well as the area accuracv loss is less than 2.5%.
基金supported by the National Basic Research Program of China (2009CB724104)the National Natural Science Foundation of China (90716010)
文摘A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.
基金This paper was supported bythe Natural Science Foundation of Shandong Province (Grant No.y2004f13)
文摘The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.
基金This study was supported by Projects under contract Nos DY105 China's 0cean-03-01-01 and DY105-03-01-07the National Natural Science Foundation of China under contract No.40506017the Youth Foundation of Marine High-tech Project of China under contract No.2002AA616010.
文摘On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software ＂Auto mapping system of submarine topography and geomorphology MBChart＂. This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China＇ s next Five-year Plan.
文摘For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low voltage problems were highlighted in the rural power grid due to the characteristics of rural power grid. Using the distribution network flow calculation method, we evaluated the low voltage problems of the rural power grid which belongs to Chongqing Jiangbei Power Company. In addition, we collected the data of distribution transformers in electricity consumption peak period. Some practical management strategies were proposed by the analysis and evaluation of potential and appeared low voltage problems.
基金supported by the National Natural Science Foundation of China under grant Nos.40474012,40874020 and 40821062
文摘The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.
文摘This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.
