Gradient Recovery Based Two-Grid Finite Element Method for Parabolic Integro-Differential Optimal Control Problems

In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and co-state variables, and piecewise constant function is used to approximate control variables. Generally, the optimal conditions for the problem are solved iteratively until the control variable reaches error tolerance. In order to calculate all the variables individually and parallelly, we introduce a gradient recovery based two-grid method. First, we solve the small scaled optimal control problem on coarse grids. Next, we use the gradient recovery technique to recover the gradients of state and co-state variables. Finally, using the recovered variables, we solve the large scaled optimal control problem for all variables independently. Moreover, we estimate priori error for the proposed scheme, and use an example to validate the theoretical results.

Hermite WENO-based limiters for high order discontinuous Galerkin method on unstructured grids

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.

A Finite Volume Method with Unstructured Triangular Grids for Numerical Modeling of Tidal Current

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.

Resources calculation of cobalt-rich crusts with the grid subdivision and integral method

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.

Analysis and management strategies for the low voltage problems of rural power grid based on the distribution network flow calculation method

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.

Unstructured Grid Immersed Boundary Method for Numerical Simulation of Fluid Structure Interaction

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.

A scheme to treat the singularity in global seismic wavefield simulation using pseudospectral method with staggered grids

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.

An explicit finite volume element method for solving characteristic level set equation on triangular grids

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.

New Immersed Boundary Method on the Adaptive Cartesian Grid Applied to the Local Discontinuous Galerkin Method

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.

A finite volume method for global electromagnetic induction forward modeling on collocated unstructured grids

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.

Index system and methods for comprehensive assessment of cross-border power grid interconnection projects

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.

Thermal Comfort Measurement and Landscape Reconstruction Strategy of Street Space Based on Decomposition Grid Method:A Case Study of Yangmeizhu Street

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.

Two-grid methods for semi-linear elliptic interface problems by immersed finite element methods

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.

A baroclinic typhoon model with a moving multi-nested grid and variational adjustment initialisationI. Numerical method

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.

Chebyshev finite spectral method with extended moving grids

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.

HYBRID CARTESIAN GRID/GRIDLESS METHOD FOR CALCULATING VISCOUS FLOWS OVER MULTI-ELEMENT AIRFOILS

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.

Preliminary Analysis of Smart Grid Risk Index System and Evaluation Methods

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 the Penalty Method of the Steady Navier-Stokes Equations

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.

A Staggered Grid Method for Solving Incompressible Flow on Unstructured Meshes

A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid technique.The computational procedure can handle cells of arbitrary shapes,although solutions presented in this paper were only involved with triangular and quadrilateral cells.The pressure or pressure-correction value was stored on the vertex of cells.The mass conservation equation was discretized on the dual cells surrounding the vertex of primary cells,while the velocity components and other scale variables were saved on the central of primary cells.Since the semi-staggered arrangement can’t guarantee a strong coupling relationship between pressure and velocity,thus a weak coupling relationship leads to the oscillations for pressure and velocity.In order to eliminate such an oscillation,a special interpolation scheme was used to construct the pressure-correction equation.Computational results of several viscous flow problems show good agreement with the analytical or numerical results in previous literature.This semi-staggered grid method can be applied to arbitrary shape elements,while it has the most efficiency for triangular cells.

3D numerical simulation in acoustic-elastic coupled media with staggered-grid finite-difference method

Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.