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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow pheno...Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.展开更多
Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and ...Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and this method employs the adaptive Cartesian grid to improve the adaptability to complex shapes and the immersed boundary to increase computational e ciency. The new immersed boundary method employs different boundary cells(the physical cell and ghost cell) to impose the boundary condition and the reconstruction algorithm of the ghost cell is the key for this method. The classical model elliptic equation is used to test the method. This method is tested and analyzed from the viewpoints of boundary cell type, error distribution and accuracy. The numerical result shows that the presented method has low error and a good rate of the convergence and works well in complex geometries. The method has good prospect for practical application research of the numerical calculation research.展开更多
Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellit...Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellites,inverting the Earth’s three-dimensional conductivity distribution on a global scale becomes attainable.A key requirement in the global conductivity inversion is to have a forward solver with high-accuracy and efficiency.In this study,a finite volume method for global electromagnetic induction forward modeling is developed based on unstructured grids.Arbitrary polyhedral grids are supported in our algorithms to obtain high geometric adaptability.We employ a cell-centered collocated variable arrangement which allows convenient discretization for complex geometries and straightforward implementation of multigrid technique.To validate the method,we test our code with two synthetic models and compare our finite volume results with an analytical solution and a finite element numerical solution.Good agreements are observed between our solution and other results,indicating acceptable accuracy of the proposed method.展开更多
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.展开更多
As an important public space in the city,street space is the second residence of people.Taking Yangmeizhu Street as the research object,this study established a street microclimate model in Beijing through microclimat...As an important public space in the city,street space is the second residence of people.Taking Yangmeizhu Street as the research object,this study established a street microclimate model in Beijing through microclimate measurement and simulation,and then calculated thermal comfort index of the street and evaluated thermal comfort of the street space and distribution model.The evaluation method used grid method to decompose and study and made several 10m×10m grids for inner space of the street,so as to compare deeply difference of microclimate environment in the street and analyze quantitatively relevance between street space elements(street pavement,street greening,building shadow coverage)and thermal comfort(physiological equivalent temperature PET)in the data.Finally,adaptive strategy of microclimate on street level was put forward to improve the environmental quality of public space in the old city by means of landscape design,aiming at creating a resilient and humanized street space and improving adaptability of microclimate and thermal comfort for living.展开更多
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.展开更多
In the context of the transformation method, we propose a general approach to construct numerically the mapping generated by imposing specific boundary conditions with a targeted function, and the necessary material a...In the context of the transformation method, we propose a general approach to construct numerically the mapping generated by imposing specific boundary conditions with a targeted function, and the necessary material and heat source spatial distributions are then derived with the help of transformation method. The construction of mapping by grid generation method through solving partial differential equations circumvents the limitation of device geometry, which paves the way for designing more complex heat flow control devices. Two numerical examples are also given to show how to design material properties and heat source in order to control temperature patterns.展开更多
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 Chebyshev finite spectral method on non-uniform meshes is proposed. An equidistribution scheme for two types of extended moving grids is used to generate grids. One type is designed to provide better resolution for ...A Chebyshev finite spectral method on non-uniform meshes is proposed. An equidistribution scheme for two types of extended moving grids is used to generate grids. One type is designed to provide better resolution for the wave surface, and the other type is for highly variable gradients. The method has high-order accuracy because of the use of the Chebyshev polynomial as the basis function. The polynomial is used to interpolate the values between the two non-uniform meshes from a previous time step to the current time step. To attain high accuracy in the time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme is used. To avoid numerical oscillations caused by the dispersion term in the Korteweg-de Vries （KdV） equation, a numerical technique on non-uniform meshes is introduced. The proposed numerical scheme is validated by the applications to the Burgers equation （nonlinear convectiondiffusion problems） and the KdV equation （single solitary and 2-solitary wave problems）, where analytical solutions are available for comparisons. Numerical results agree very well with the corresponding analytical solutions in all cases.展开更多
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.展开更多
A two grid technique for solving the steady incompressible Navier Stokes equations in a penalty method was presented and the convergence of numerical solutions was analyzed. If a coarse size H and a fine size ...A two grid technique for solving the steady incompressible Navier Stokes equations in a penalty method was presented and the convergence of numerical solutions was analyzed. If a coarse size H and a fine size h satisfy H=O(h 13-s )(s=0(n=2);s=12(n=3), where n is a space dimension), this method has the same convergence accuracy as the usual finite element method. But the two grid method can save a lot of computation time for its brief calculation. Moreover, a numerical test was couducted in order to verify the correctness of above theoretical analysis.展开更多
基金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.
文摘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.
基金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.
文摘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 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.
基金King Mongkut’s University of Technology North Bangkok (KMUTNB)the Office of the Higher Education Commission (OHEC)the National Metal and Materials Technology Center (MTEC) for supporting this research work
文摘Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.
基金Supported by National Natural Science Foundation of China(Grant No.51405375)National Key Basic Research and Development Program of China(973 Program,Grant No.2011CB706606)
文摘Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and this method employs the adaptive Cartesian grid to improve the adaptability to complex shapes and the immersed boundary to increase computational e ciency. The new immersed boundary method employs different boundary cells(the physical cell and ghost cell) to impose the boundary condition and the reconstruction algorithm of the ghost cell is the key for this method. The classical model elliptic equation is used to test the method. This method is tested and analyzed from the viewpoints of boundary cell type, error distribution and accuracy. The numerical result shows that the presented method has low error and a good rate of the convergence and works well in complex geometries. The method has good prospect for practical application research of the numerical calculation research.
基金supported by the National Natural Science Foundation of China(41922027,4214200052)by the Macao Foundation+1 种基金by the Pre-research Project on Civil Aerospace Technologies No.D020308/D020303 funded by China National Space Administrationby the Macao Science and Technology Development Fund,grant No.0001/2019/A1。
文摘Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellites,inverting the Earth’s three-dimensional conductivity distribution on a global scale becomes attainable.A key requirement in the global conductivity inversion is to have a forward solver with high-accuracy and efficiency.In this study,a finite volume method for global electromagnetic induction forward modeling is developed based on unstructured grids.Arbitrary polyhedral grids are supported in our algorithms to obtain high geometric adaptability.We employ a cell-centered collocated variable arrangement which allows convenient discretization for complex geometries and straightforward implementation of multigrid technique.To validate the method,we test our code with two synthetic models and compare our finite volume results with an analytical solution and a finite element numerical solution.Good agreements are observed between our solution and other results,indicating acceptable accuracy of the proposed method.
基金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.
基金Beijing Natural Science Foundation(8202017)2018 Beijing Municipal University Academic Human Resources Development-Youth Talent Support Program(PXM2018–014212–000043)National Natural Science Foundation of China(51708004).
文摘As an important public space in the city,street space is the second residence of people.Taking Yangmeizhu Street as the research object,this study established a street microclimate model in Beijing through microclimate measurement and simulation,and then calculated thermal comfort index of the street and evaluated thermal comfort of the street space and distribution model.The evaluation method used grid method to decompose and study and made several 10m×10m grids for inner space of the street,so as to compare deeply difference of microclimate environment in the street and analyze quantitatively relevance between street space elements(street pavement,street greening,building shadow coverage)and thermal comfort(physiological equivalent temperature PET)in the data.Finally,adaptive strategy of microclimate on street level was put forward to improve the environmental quality of public space in the old city by means of landscape design,aiming at creating a resilient and humanized street space and improving adaptability of microclimate and thermal comfort for living.
基金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.
基金supported by the National Natural Science Foundation of China(Nos.11372035,10832002 and 11221202)China Postdoctoral Science Foundation(No.2014M550054)
文摘In the context of the transformation method, we propose a general approach to construct numerically the mapping generated by imposing specific boundary conditions with a targeted function, and the necessary material and heat source spatial distributions are then derived with the help of transformation method. The construction of mapping by grid generation method through solving partial differential equations circumvents the limitation of device geometry, which paves the way for designing more complex heat flow control devices. Two numerical examples are also given to show how to design material properties and heat source in order to control temperature patterns.
基金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 Research Grants Council of Hong Kong (No. 522007)the National Marine Public Welfare Research Projects of China (No. 201005002)
文摘A Chebyshev finite spectral method on non-uniform meshes is proposed. An equidistribution scheme for two types of extended moving grids is used to generate grids. One type is designed to provide better resolution for the wave surface, and the other type is for highly variable gradients. The method has high-order accuracy because of the use of the Chebyshev polynomial as the basis function. The polynomial is used to interpolate the values between the two non-uniform meshes from a previous time step to the current time step. To attain high accuracy in the time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme is used. To avoid numerical oscillations caused by the dispersion term in the Korteweg-de Vries （KdV） equation, a numerical technique on non-uniform meshes is introduced. The proposed numerical scheme is validated by the applications to the Burgers equation （nonlinear convectiondiffusion problems） and the KdV equation （single solitary and 2-solitary wave problems）, where analytical solutions are available for comparisons. Numerical results agree very well with the corresponding analytical solutions in all cases.
基金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.
文摘A two grid technique for solving the steady incompressible Navier Stokes equations in a penalty method was presented and the convergence of numerical solutions was analyzed. If a coarse size H and a fine size h satisfy H=O(h 13-s )(s=0(n=2);s=12(n=3), where n is a space dimension), this method has the same convergence accuracy as the usual finite element method. But the two grid method can save a lot of computation time for its brief calculation. Moreover, a numerical test was couducted in order to verify the correctness of above theoretical analysis.
