In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the II...In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.展开更多
In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squar...In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.展开更多
In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-f...In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-free Galerkin （ICVEFG） method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square （CVMLS） approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.展开更多
In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved c...In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin（ICVEFG） method is presented for two-dimensional（2D） elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.展开更多
We first give a stabilized improved moving least squares （IMLS） approximation, which has better computational stability and precision than the IMLS approximation. Then, analysis of the improved element-free Galerkin...We first give a stabilized improved moving least squares （IMLS） approximation, which has better computational stability and precision than the IMLS approximation. Then, analysis of the improved element-free Galerkin method is provided theoretically for both linear and nonlinear elliptic boundary value problems. Finally, numerical examples are given to verify the theoretical analysis.展开更多
Based on the complex variable moving least-square （CVMLS） approximation, the complex variable element-free Galerkin （CVEFG） method for two-dimensional viscoelasticity problems under the creep condition is presente...Based on the complex variable moving least-square （CVMLS） approximation, the complex variable element-free Galerkin （CVEFG） method for two-dimensional viscoelasticity problems under the creep condition is presented in this paper. The Galerkin weak form is employed to obtain the equation system, and the penalty method is used to apply the essential boundary conditions, then the corresponding formulae of the CVEFG method for two-dimensional viscoelasticity problems under the creep condition are obtained. Compared with the element-free Galerkin （EFG） method, with the same node distribution, the CVEFG method has higher precision, and to obtain the similar precision, the CVEFG method has greater computational efficiency. Some numerical examples are given to demonstrate the validity and the efficiency of the method.展开更多
In this paper, an improved element-free Galerkin （IEFG） method is proposed to solve the generalized fifth-order Korteweg-de Vries （gfKdV） equation. When the traditional element-free Galerkin （EFG） method is used...In this paper, an improved element-free Galerkin （IEFG） method is proposed to solve the generalized fifth-order Korteweg-de Vries （gfKdV） equation. When the traditional element-free Galerkin （EFG） method is used to solve such an equation, unstable or even wrong numerical solutions may be obtained due to the violation of the consistency conditions of the moving least-squares （MLS） shape functions. To solve this problem, the EFG method is improved by employing the improved moving least-squares （IMLS） approximation based on the shifted polynomial basis functions. The effectiveness of the IEFG method for the gfKdV equation is investigated by using some numerical examples. Meanwhile, the motion of single solitary wave and the interaction of two solitons are simulated using the IEFG method.展开更多
The element-free Galerkin （EFG） method is used in this paper to find the numerical solution to a regularized long-wave （RLW） equation. The Galerkin weak form is adopted to obtain the discrete equations, and the es...The element-free Galerkin （EFG） method is used in this paper to find the numerical solution to a regularized long-wave （RLW） equation. The Galerkin weak form is adopted to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. The effectiveness of the EFG method of solving the RLW equation is investigated by two numerical examples in this paper.展开更多
The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the element-free Calerkin （EFG） method which is based on the moving least-squares approximation. A variational meth...The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the element-free Calerkin （EFG） method which is based on the moving least-squares approximation. A variational method is used to obtain discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for KdV equations needs only scattered nodes instead of meshing the domain of the problem. It does not require any element connectivity and does not suffer much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for the KdV equation is investigated by two numerical examples in this paper.展开更多
Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity proble...Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares （IMLS） method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin （IEFG） method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.展开更多
The improved element-free Galerkin （IEFG） method of elasticity is used to solve the topology optimization problems. In this method, the improved moving least-squares approximation is used to form the shape function....The improved element-free Galerkin （IEFG） method of elasticity is used to solve the topology optimization problems. In this method, the improved moving least-squares approximation is used to form the shape function. In a topology opti- mization process, the entire structure volume is considered as the constraint. From the solid isotropic microstructures with penalization, we select relative node density as a design variable. Then we choose the minimization of compliance to be an objective function, and compute its sensitivity with the adjoint method. The IEFG method in this paper can overcome the disadvantages of the singular matrices that sometimes appear in conventional element-free Galerkin （EFG） method. The central processing unit （CPU） time of each example is given to show that the IEFG method is more efficient than the EFG method under the same precision, and the advantage that the IEFG method does not form singular matrices is also shown.展开更多
The present paper deals with the numerical solution of a two-dimensional linear hyperbolic equation by using the element-free Galerkin （EFG） method which is based on the moving least-square approximation for the tes...The present paper deals with the numerical solution of a two-dimensional linear hyperbolic equation by using the element-free Galerkin （EFG） method which is based on the moving least-square approximation for the test and trial functions. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for hyperbolic problems needs only the scattered nodes instead of meshing the domain of the problem. It neither requires any element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for two-dimensional hyperbolic problems is investigated by two numerical examples in this paper.展开更多
A generalized Fisher equation （GFE） relates the time derivative of the average of the intrinsic rate of growth to its variance. The exact mathematical result of the GFE has been widely used in population dynamics an...A generalized Fisher equation （GFE） relates the time derivative of the average of the intrinsic rate of growth to its variance. The exact mathematical result of the GFE has been widely used in population dynamics and genetics, where it originated. Many researchers have studied the numerical solutions of the GFE, up to now. In this paper, we introduce an element-free Galerkin （EFG） method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics. Compared with other numerical methods, the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. In comparison with the traditional method, numerical solutions show that the new method has higher accuracy and better convergence. Several numerical examples are presented to demonstrate the effectiveness of the method.展开更多
The present paper deals with the numerical solution of the coupled Schrodinger-KdV equations using the elementfree Galerkin （EFG） method which is based on the moving least-square approximation. Instead of traditiona...The present paper deals with the numerical solution of the coupled Schrodinger-KdV equations using the elementfree Galerkin （EFG） method which is based on the moving least-square approximation. Instead of traditional mesh oriented methods such as the finite difference method （FDM） and the finite element method （FEM）, this method needs only scattered nodes in the domain. For this scheme, a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method. In numerical experiments, the results are presented and compared with the findings of the finite element method, the radial basis functions method, and an analytical solution to confirm the good accuracy of the presented scheme.展开更多
The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin （EFG） method, which is based on the moving least-square approximation. Compared w...The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin （EFG） method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α（0 〈 α≤ 1）. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.展开更多
The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （B...The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.展开更多
In this article,we propose and research a first-order,linearized discontinuous Galerkin method for the approximation of the hydrodynamic and sediment transport model.The method is decoupled and fully discrete,and is s...In this article,we propose and research a first-order,linearized discontinuous Galerkin method for the approximation of the hydrodynamic and sediment transport model.The method is decoupled and fully discrete,and is shown to be unconditionally stable.Furthermore,error estimates are proved.Finally,the theoretical analysis is confirmed by numerical examples.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations....A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.展开更多
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.
基金supported by the National Natural Science Foundation of China (Grant No.11026223)the Shanghai Leading Academic Discipline Project,China (Grant No.S30106)the Innovation Fund Project for Graduate Student of Shanghai University,China (Grant No.SHUCX112359)
文摘In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)the Innovation Fund Project for Graduate Student of Shanghai University,China (Grant No. SHUCX112359)
文摘In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-free Galerkin （ICVEFG） method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square （CVMLS） approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11171208 and U1433104)
文摘In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin（ICVEFG） method is presented for two-dimensional（2D） elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.
基金Project supported by the National Natural Science Foundation of China(Grant No.11471063)the Chongqing Research Program of Basic Research and Frontier Technology,China(Grant No.cstc2015jcyj BX0083)the Educational Commission Foundation of Chongqing City,China(Grant No.KJ1600330)
文摘We first give a stabilized improved moving least squares （IMLS） approximation, which has better computational stability and precision than the IMLS approximation. Then, analysis of the improved element-free Galerkin method is provided theoretically for both linear and nonlinear elliptic boundary value problems. Finally, numerical examples are given to verify the theoretical analysis.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘Based on the complex variable moving least-square （CVMLS） approximation, the complex variable element-free Galerkin （CVEFG） method for two-dimensional viscoelasticity problems under the creep condition is presented in this paper. The Galerkin weak form is employed to obtain the equation system, and the penalty method is used to apply the essential boundary conditions, then the corresponding formulae of the CVEFG method for two-dimensional viscoelasticity problems under the creep condition are obtained. Compared with the element-free Galerkin （EFG） method, with the same node distribution, the CVEFG method has higher precision, and to obtain the similar precision, the CVEFG method has greater computational efficiency. Some numerical examples are given to demonstrate the validity and the efficiency of the method.
基金the National Basic Research Program of China(Grant No.2012CB025903)
文摘In this paper, an improved element-free Galerkin （IEFG） method is proposed to solve the generalized fifth-order Korteweg-de Vries （gfKdV） equation. When the traditional element-free Galerkin （EFG） method is used to solve such an equation, unstable or even wrong numerical solutions may be obtained due to the violation of the consistency conditions of the moving least-squares （MLS） shape functions. To solve this problem, the EFG method is improved by employing the improved moving least-squares （IMLS） approximation based on the shifted polynomial basis functions. The effectiveness of the IEFG method for the gfKdV equation is investigated by using some numerical examples. Meanwhile, the motion of single solitary wave and the interaction of two solitons are simulated using the IEFG method.
基金supported by the Natural Science Foundation of Zhejiang Province of China (Grant No. Y6110007)
文摘The element-free Galerkin （EFG） method is used in this paper to find the numerical solution to a regularized long-wave （RLW） equation. The Galerkin weak form is adopted to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. The effectiveness of the EFG method of solving the RLW equation is investigated by two numerical examples in this paper.
基金Project supported by the Natural Science Foundation of Ningbo City (Grant No.2009A610014)the Natural Science Foundation of Zhejiang Province (Grant No.Y6090131)the Research Foundation of Ningbo University of Technology (Grant No.2008004)
文摘The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the element-free Calerkin （EFG） method which is based on the moving least-squares approximation. A variational method is used to obtain discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for KdV equations needs only scattered nodes instead of meshing the domain of the problem. It does not require any element connectivity and does not suffer much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for the KdV equation is investigated by two numerical examples in this paper.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)
文摘Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares （IMLS） method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin （IEFG） method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.
基金supported by the National Natural Science Foundation of China(Grant Nos.11571223 and U1433104)
文摘The improved element-free Galerkin （IEFG） method of elasticity is used to solve the topology optimization problems. In this method, the improved moving least-squares approximation is used to form the shape function. In a topology opti- mization process, the entire structure volume is considered as the constraint. From the solid isotropic microstructures with penalization, we select relative node density as a design variable. Then we choose the minimization of compliance to be an objective function, and compute its sensitivity with the adjoint method. The IEFG method in this paper can overcome the disadvantages of the singular matrices that sometimes appear in conventional element-free Galerkin （EFG） method. The central processing unit （CPU） time of each example is given to show that the IEFG method is more efficient than the EFG method under the same precision, and the advantage that the IEFG method does not form singular matrices is also shown.
基金Project supported by the Natural Science Foundation of Ningbo, China (Grant Nos 2009A610014, 2009A610154, 2008A610020 and 2007A610050)
文摘The present paper deals with the numerical solution of a two-dimensional linear hyperbolic equation by using the element-free Galerkin （EFG） method which is based on the moving least-square approximation for the test and trial functions. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for hyperbolic problems needs only the scattered nodes instead of meshing the domain of the problem. It neither requires any element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for two-dimensional hyperbolic problems is investigated by two numerical examples in this paper.
基金supported by the National Natural Science Foundation of China (Grant No. 11072117)the Natural Science Foundation of Ningbo City (Grant Nos. 2012A610038 and 2012A610152)+1 种基金the Scientific Research Fund of Education Department of Zhejiang Province,China (Grant No. Z201119278)K.C. Wong Magna Fund in Ningbo University
文摘A generalized Fisher equation （GFE） relates the time derivative of the average of the intrinsic rate of growth to its variance. The exact mathematical result of the GFE has been widely used in population dynamics and genetics, where it originated. Many researchers have studied the numerical solutions of the GFE, up to now. In this paper, we introduce an element-free Galerkin （EFG） method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics. Compared with other numerical methods, the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. In comparison with the traditional method, numerical solutions show that the new method has higher accuracy and better convergence. Several numerical examples are presented to demonstrate the effectiveness of the method.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11072117 and 61074142)the Natural Science Foundation of Zhejiang Province,China(Grant No.Y6110007)+3 种基金Scientific Research Fund of Zhejiang Provincial Education Department,China(Grant No.Z201119278)the Natural Science Foundation of Ningbo City(Grant Nos.2012A610152 and 2012A610038)the Disciplinary Project of Ningbo City,China(Grant No.SZXL1067)K.C.Wong Magna Fund in Ningbo University
文摘The present paper deals with the numerical solution of the coupled Schrodinger-KdV equations using the elementfree Galerkin （EFG） method which is based on the moving least-square approximation. Instead of traditional mesh oriented methods such as the finite difference method （FDM） and the finite element method （FEM）, this method needs only scattered nodes in the domain. For this scheme, a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method. In numerical experiments, the results are presented and compared with the findings of the finite element method, the radial basis functions method, and an analytical solution to confirm the good accuracy of the presented scheme.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Natural Science Foundationof Zhejiang Province,China(Grant Nos.Y6110007and Y6110502)the K.C.Wong Magna Fund in Ningbo University,China
文摘The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin （EFG） method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α（0 〈 α≤ 1）. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.
基金Project supported by the National Natural Science Foundation of China (Grant No 10871124)Innovation Program of Shanghai Municipal Education Commission (Grant No 09ZZ99)Shanghai Leading Academic Discipline Project (Grant No J50103)
文摘The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.
基金Supported by the National Natural Science Foundation of China(Grant No.41930643)the Natural Science Foundation of Henan Province(Grant No.232300420109).
文摘In this article,we propose and research a first-order,linearized discontinuous Galerkin method for the approximation of the hydrodynamic and sediment transport model.The method is decoupled and fully discrete,and is shown to be unconditionally stable.Furthermore,error estimates are proved.Finally,the theoretical analysis is confirmed by numerical examples.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
基金Supported by the National Natural Science Foundation of China(50976072,51106099,10902070)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)the Science Foundation for the Excellent Youth Scholar of Higher Education of Shanghai(slg09003)~~
文摘A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.
