Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement...
详细信息
Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp://lsec. cc. ac. cn/phg/), a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.
This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding diffe...
详细信息
This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential *** particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical *** an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.
For large sparse non-Hermitian positive definite system of linear equations, we present several variants of the Hermitian and skew-Hermitian splitting (HSS) about the coefficient matrix and establish correspondingly s...
详细信息
For large sparse non-Hermitian positive definite system of linear equations, we present several variants of the Hermitian and skew-Hermitian splitting (HSS) about the coefficient matrix and establish correspondingly several HSS-based iterative schemes. Theoretical analyses show that these methods are convergent unconditionally to the exact solution of the referred system of linear equations, and they may show advantages on problems that the HSS method is ineffective.
A discrete total variation calculus with variable time steps is presented for mechanico-electrical systems where there exist non-potential and dissipative forces. By using this discrete variation calculus, the symplec...
详细信息
A discrete total variation calculus with variable time steps is presented for mechanico-electrical systems where there exist non-potential and dissipative forces. By using this discrete variation calculus, the symplectic-energy-first integrators for mechanico-electrical systems are derived. To do this, the time step adaptation is employed. The discrete variational principle and the Euler-Lagrange equation are derived for the systems. By using this discrete algorithm it is shown that mechanico-electrical systems are not symplectic and their energies are not conserved unless they are Lagrange mechanico-electrical systems. A practical example is presented to illustrate these results.
Starting from the variable separation approach, the algebraic soliton solution and the solution describing the interaction between line soliton and algebraic soliton are obtained by selecting appropriate seed solution...
详细信息
Starting from the variable separation approach, the algebraic soliton solution and the solution describing the interaction between line soliton and algebraic soliton are obtained by selecting appropriate seed solution for (2+1)-dimensional ANNV equation. The behaviors of interactions are discussed in detail both analytically and graphically. It is shown that there are two kinds of singular interactions between line soliton and algebraic soliton: 1) the resonant interaction where the algebraic soliton propagates together with the line soliton and persists infinitely; 2) the extremely repulsive interaction where the algebraic soliton affects the motion of the line soliton infinitely apart.
In this paper,we study a posteriori error estimates of the edge stabilization Galerkin method for the constrained optimal control problem governed by convection-dominated diffusion *** residual-type a posteriori error...
详细信息
In this paper,we study a posteriori error estimates of the edge stabilization Galerkin method for the constrained optimal control problem governed by convection-dominated diffusion *** residual-type a posteriori error estimators yield both upper and lower bounds for control u measured in L2-norm and for state y and costate p measured in energy *** numerical examples are presented to illustrate the effectiveness of the error estimators provided in this paper.
Stochastic approximation problem is to find some root or extremum of a nonlinear function for which only noisy measurements of the function are available. The classical algorithm for stochastic approximation problem i...
详细信息
Stochastic approximation problem is to find some root or extremum of a nonlinear function for which only noisy measurements of the function are available. The classical algorithm for stochastic approximation problem is the Robbins-Monro (RM) algorithm, which uses the noisy evaluation of the negative gradient direction as the iterative direction. In order to accelerate the RM algorithm, this paper gives a flame algorithm using adaptive iterative directions. At each iteration, the new algorithm goes towards either the noisy evaluation of the negative gradient direction or some other directions under some switch criterions. Two feasible choices of the criterions are proposed and two corresponding frame algorithms are formed. Different choices of the directions under the same given switch criterion in the frame can also form different algorithms. We also proposed the simultanous perturbation difference forms for the two frame algorithms. The almost surely convergence of the new algorithms are all established. The numerical experiments show that the new algorithms are promising.
In this paper the work on implementing two mesh partitioning algorithms, the refinement-tree based partitioning algorithm and the space-filling curve partitioning algorithm, in the parallel adaptive finite element too...
详细信息
In this paper the work on implementing two mesh partitioning algorithms, the refinement-tree based partitioning algorithm and the space-filling curve partitioning algorithm, in the parallel adaptive finite element toolbox PHG (parallel hierarchical grid) is presented. These algorithms are used for both initial mesh partitioning and mesh repartitioning for dynamical load balancing in adaptive finite element computations. In the implementations improved algorithms are designed. Partitioning time and quality of our code are compared with existing publicly available mesh or graph partitioners, including ParMETIS and Zoltan, through some numerical examples.
PHG (parallel hierarchical grid) is a scalable parallel adaptive finite element toolbox under active developmentat the statekeylaboratory of scientific and engineeringcomputing, Chinese Academy of Sciences. This pa...
详细信息
PHG (parallel hierarchical grid) is a scalable parallel adaptive finite element toolbox under active developmentat the statekeylaboratory of scientific and engineeringcomputing, Chinese Academy of Sciences. This paper demonstrates its application to adaptive finite element computations of electromagnetic problems. Two examples on solving the time harmonic Maxwell's equations are shown. Results of some large scale adaptive finite element simulations with up to 1 billion degrees of freedom and using up to 2048 CPUs are presented.
暂无评论