In recent years,meshless methods have been increasingly applied to the simulation of various engineering problems due to their inherent advantages over traditional mesh-based approaches,including greater flexibility,i...In recent years,meshless methods have been increasingly applied to the simulation of various engineering problems due to their inherent advantages over traditional mesh-based approaches,including greater flexibility,independence from predefined meshing,simpler adaptive analysis,improved automation,and suitability for complex problems.Several meshless methods have been used for porous media simulation,and are broadly categorized into collocation,global weak form and local weak form methods.In this study,a comprehensive comparison of the applicability of these three categories of meshless methods for simulating coupled flow and transport problems in porous media is presented.The Radial Point Collocation Method(RPCM)(strong form),the Element Free Galerkin Method(EFGM)(global weak form)and the Meshless Local Petrov Galerkin(MLPG)method(local weak form)are implemented and systematically compared.These methods are applied to the analysis of flow in a synthetic regular domain aquifer,flow and non-reactive contaminant transport in a synthetic irregular boundary porous media problem and groundwater flow in a field aquifer located in India.The simulated groundwater heads are compared with analytical solution,observed field data and results obtained from widely used MODFLOW-MT3DMS models.The deviation of the solutions from the analytical solution is in the range of O.67%to O.l6%for the hypothetical case study.For the field-scale case study,mean absolute error of 0.183%,0.18l%and 0.188%are obtained for the RPCM,EFGM and MLPG models,respectively,outperforming MODFLOW,which exhibits a deviation of 0.254%from observed values.Overall,the present study reaffirms the practical applicability of these meshless methods for real-world groundwater problems and provides valuable insights into the utilization of each category of meshless method,with respect to problem type,computational efficiency and accuracy requirements.展开更多
After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir develo...After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.展开更多
In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equi...In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.展开更多
In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is pres...In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.展开更多
The objectives of this study are to employ the meshless local Petrov-Galerkin method (MLPGM) to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many fact...The objectives of this study are to employ the meshless local Petrov-Galerkin method (MLPGM) to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors' sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.展开更多
Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as t...Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as the basis of the moving least square method to construct the meshless interpolation function. Multi-resolution analysis is used to decompose the field variables into high and low scales and the high scale component can commonly represent the gradient of the solution according to inherent characteristics of wavelets. The high scale component in the present method can directly detect high gradient regions of the field variables. The developed adaptive refinement scheme has been applied to simulate actual examples, and the effectiveness of the present adaptive refinement scheme has been verified.展开更多
A meshiess local discontinuous Petrov-Galerkin (MLDPG) method based on the local symmetric weak form (LSWF) is presented with the application to blasting problems. The derivation is similar to that of mesh-based R...A meshiess local discontinuous Petrov-Galerkin (MLDPG) method based on the local symmetric weak form (LSWF) is presented with the application to blasting problems. The derivation is similar to that of mesh-based Runge-Kutta Discontinuous Galerkin (RKDG) method. The solutions are reproduced in a set of overlapped spherical sub-domains, and the test functions are employed from a partition of unity of the local basis functions. There is no need of any traditional nonoverlapping mesh either for local approximation purpose or for Galerkin integration purpose in the presented method. The resulting MLDPG method is a meshless, stable, high-order accurate and highly parallelizable scheme which inherits both the advantages of RKDG and meshless method (MM), and it can handle the problems with extremely complicated physics and geometries easily. Three numerical examples of the one-dimensional Sod shock-tube problem, the blast-wave problem and the Woodward-Colella interacting shock wave problem are given. All the numerical results are in good agreement with the closed solutions. The higher-order MLDPG schemes can reproduce more accurate solution than the lower-order schemes.展开更多
Based on the complex variable moving least-square(CVMLS) approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin(CVMLPG) method of solving two-dimensional potential proble...Based on the complex variable moving least-square(CVMLS) approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin(CVMLPG) method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square(MLS) approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin(MLPG) method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.展开更多
Traditional topology optimization methods often introduce weak artificial material to mimic voids to avoid the singularity of the global stiffness matrix and carry out topology optimization with a fixed finite element...Traditional topology optimization methods often introduce weak artificial material to mimic voids to avoid the singularity of the global stiffness matrix and carry out topology optimization with a fixed finite element(FE)mesh.This treatment,however,may not only increase the computational cost for structural analysis but also lead to unfavorable numerical instabilities,especially when large deformations and dynamic/buckling behaviors are involved.In the present work,a new meshless moving morphable component-based method(ML-MMC),which structural analysis is carried out only on the solid region occupied by components,is proposed.In this approach,the coupling of discrete components is achieved through the adaptively constructed influence domain of the meshless shape function.Therefore,the singularity problem of the stiffness matrix can be naturally avoided without introducing weak artificial material.Compared with traditional methods,the number of degrees of freedoms(DOFs)can be reduced substantially under this treatment.The effectiveness of the proposed approach is also illustrated by some representative examples.展开更多
This paper presents a meshless method for the nonlinear generalized regularized long wave (GRLW) equation based on the moving least-squares approximation. The nonlinear discrete scheme of the GRLW equation is obtain...This paper presents a meshless method for the nonlinear generalized regularized long wave (GRLW) equation based on the moving least-squares approximation. The nonlinear discrete scheme of the GRLW equation is obtained and is solved using the iteration method. A theorem on the convergence of the iterative process is presented and proved using theorems of the infinity norm. Compared with numerical methods based on mesh, the meshless method for the GRLW equation only requires the.scattered nodes instead of meshing the domain of the problem. Some examples, such as the propagation of single soliton and the interaction of two solitary waves, are given to show the effectiveness of the meshless method.展开更多
The element-free Galerkin (EFG) method for numerically solving the compound Korteweg-de Vries-Burgers (KdVB) equation is discussed in this paper. The Galerkin weak form is used to obtain the discrete equation and ...The element-free Galerkin (EFG) method for numerically solving the compound Korteweg-de Vries-Burgers (KdVB) equation is discussed in this paper. The Galerkin weak form is used to obtain the discrete equation and the essential boundary conditions are enforced by the penalty method. The effectiveness of the EFG method of solving the compound Korteweg-de Vries-Burgers (KdVB) equation is illustrated by three numerical examples.展开更多
We propose a novel machine learning approach to reconstruct meshless surface wind speed fields,i.e.,to reconstruct the surface wind speed at any location,based on meteorological background fields and geographical info...We propose a novel machine learning approach to reconstruct meshless surface wind speed fields,i.e.,to reconstruct the surface wind speed at any location,based on meteorological background fields and geographical information.The random forest method is selected to develop the machine learning data reconstruction model(MLDRM-RF)for wind speeds over Beijing from 2015-19.We use temporal,geospatial attribute and meteorological background field features as inputs.The wind speed field can be reconstructed at any station in the region not used in the training process to cross-validate model performance.The evaluation considers the spatial distribution of and seasonal variations in the root mean squared error(RMSE)of the reconstructed wind speed field across Beijing.The average RMSE is 1.09 m s−1,considerably smaller than the result(1.29 m s−1)obtained with inverse distance weighting(IDW)interpolation.Finally,we extract the important feature permutations by the method of mean decrease in impurity(MDI)and discuss the reasonableness of the model prediction results.MLDRM-RF is a reasonable approach with excellent potential for the improved reconstruction of historical surface wind speed fields with arbitrary grid resolutions.Such a model is needed in many wind applications,such as wind energy and aviation safety assessments.展开更多
Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attr...Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson's equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.展开更多
Based on our previous study,the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs)are used,so that it could re...Based on our previous study,the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs)are used,so that it could result in significant error in solving partial differential equations with Neumann boundary conditions.To overcome this drawback,the Consistent Compactly Supported Radial Basis Functions(CCSRBFs)are developed,which satisfy the predetermined consistency con- ditions.Meshless method based on point collocation with CCSRBFs is developed for solving partial differential equations.Numerical studies show that the proposed method improves the accuracy of approximation significantly.展开更多
A meshless local radial point interpolation method (LRPIM) for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using ...A meshless local radial point interpolation method (LRPIM) for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function (RBF) coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.展开更多
Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature (PUQ),for Galerkin meshless methods was presented. The technique is b...Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature (PUQ),for Galerkin meshless methods was presented. The technique is based on finite covering and partition of unity. There is no need to decompose the physical domain into small cell. It possesses remarkable integration accuracy. Using Element-free Galerkin methods as example, Galerkin meshless methods based on PUQ were studied in detail. Meshing is always not required in the procedure of constitution of approximate function or numerical quadrature, so Galerkin meshless methods based on PUQ are “truly” meshless methods.展开更多
Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the movi...Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.展开更多
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.展开更多
The meshless weighted least-square (MWLS) method was developed based on the weighted least-square method. The method possesses several advantages, such as high accuracy, high stability and high e?ciency. Moreover, the...The meshless weighted least-square (MWLS) method was developed based on the weighted least-square method. The method possesses several advantages, such as high accuracy, high stability and high e?ciency. Moreover, the coe?cient matrix obtained is symmetric and semi- positive de?nite. In this paper, the method is further examined critically. The e?ects of several parameters on the results of MWLS are investigated systematically by using a cantilever beam and an in?nite plate with a central circular hole. The numerical results are compared with those obtained by using the collocation-based meshless method (CBMM) and Galerkin-based meshless method (GBMM). The investigated parameters include the type of approximations, the type of weight functions, the number of neighbors of an evaluation point, as well as the manner in which the neighbors of an evaluation point are determined. This study shows that the displacement accuracy and convergence rate obtained by MWLS is comparable to that of the GBMM while the stress accuracy and convergence rate yielded by MWLS is even higher than that of GBMM. Furthermore, MWLS is much more e?cient than GBMM. This study also shows that the instability of CBMM is mainly due to the neglect of the equi- librium residuals at boundary nodes. In MWLS, the residuals of all the governing equations are minimized in a weighted least-square sense.展开更多
A simple direct space-time meshless scheme,based on the radial or non-radial basis function,is proposed for the onedimensional Klein-Gordon equations.Since these equations are time-dependent,it is worthwhile to presen...A simple direct space-time meshless scheme,based on the radial or non-radial basis function,is proposed for the onedimensional Klein-Gordon equations.Since these equations are time-dependent,it is worthwhile to present two schemes for the basis functions from radial and non-radial aspects.The first scheme is fulfilled by considering time variable as normal space variable,to construct an"isotropic"space-time radial basis function.The other scheme considered a realistic relationship between space variable and time variable which is not radial.The timedependent variable is treated regularly during the whole solution process and the Klein-Gordon equations can be solved in a direct way.Numerical results show that the proposed meshless schemes are simple,accurate,stable,easy-to-program and efficient for the Klein-Gordon equations.展开更多
摘要In recent years,meshless methods have been increasingly applied to the simulation of various engineering problems due to their inherent advantages over traditional mesh-based approaches,including greater flexibility,independence from predefined meshing,simpler adaptive analysis,improved automation,and suitability for complex problems.Several meshless methods have been used for porous media simulation,and are broadly categorized into collocation,global weak form and local weak form methods.In this study,a comprehensive comparison of the applicability of these three categories of meshless methods for simulating coupled flow and transport problems in porous media is presented.The Radial Point Collocation Method(RPCM)(strong form),the Element Free Galerkin Method(EFGM)(global weak form)and the Meshless Local Petrov Galerkin(MLPG)method(local weak form)are implemented and systematically compared.These methods are applied to the analysis of flow in a synthetic regular domain aquifer,flow and non-reactive contaminant transport in a synthetic irregular boundary porous media problem and groundwater flow in a field aquifer located in India.The simulated groundwater heads are compared with analytical solution,observed field data and results obtained from widely used MODFLOW-MT3DMS models.The deviation of the solutions from the analytical solution is in the range of O.67%to O.l6%for the hypothetical case study.For the field-scale case study,mean absolute error of 0.183%,0.18l%and 0.188%are obtained for the RPCM,EFGM and MLPG models,respectively,outperforming MODFLOW,which exhibits a deviation of 0.254%from observed values.Overall,the present study reaffirms the practical applicability of these meshless methods for real-world groundwater problems and provides valuable insights into the utilization of each category of meshless method,with respect to problem type,computational efficiency and accuracy requirements.
基金funded by the 14th Five-Year Plan Major Science and Technology Project of CNOOC project number KJGG2021-0506.
摘要After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.
基金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 for Graduate Student of Shanghai University,China (Grant No. SHUCX120125)
摘要In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.
基金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 for Graduate Student of Shanghai University of China (Grant No.SHUCX120125)
摘要In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.
基金the Scientific Foundation of National Outstanding Youth of China(No.50225520)the Science Foundation of Shandong University of Technology of China(No.2006KJM33).
摘要The objectives of this study are to employ the meshless local Petrov-Galerkin method (MLPGM) to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors' sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.
基金supported by the Scientific Foundation of National Outstanding Youth of China(No.50225520)Science Foundation of Shandong University of Technology of China(No.2006KJM33).
摘要Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as the basis of the moving least square method to construct the meshless interpolation function. Multi-resolution analysis is used to decompose the field variables into high and low scales and the high scale component can commonly represent the gradient of the solution according to inherent characteristics of wavelets. The high scale component in the present method can directly detect high gradient regions of the field variables. The developed adaptive refinement scheme has been applied to simulate actual examples, and the effectiveness of the present adaptive refinement scheme has been verified.
基金Supported by New Century Excellent Talents in University in China(NCET),National"973" Program(No.61338)Innovative Research Project of Xi'an Hi-Tech Institute(EPXY0806)
摘要A meshiess local discontinuous Petrov-Galerkin (MLDPG) method based on the local symmetric weak form (LSWF) is presented with the application to blasting problems. The derivation is similar to that of mesh-based Runge-Kutta Discontinuous Galerkin (RKDG) method. The solutions are reproduced in a set of overlapped spherical sub-domains, and the test functions are employed from a partition of unity of the local basis functions. There is no need of any traditional nonoverlapping mesh either for local approximation purpose or for Galerkin integration purpose in the presented method. The resulting MLDPG method is a meshless, stable, high-order accurate and highly parallelizable scheme which inherits both the advantages of RKDG and meshless method (MM), and it can handle the problems with extremely complicated physics and geometries easily. Three numerical examples of the one-dimensional Sod shock-tube problem, the blast-wave problem and the Woodward-Colella interacting shock wave problem are given. All the numerical results are in good agreement with the closed solutions. The higher-order MLDPG schemes can reproduce more accurate solution than the lower-order schemes.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant No. 11102125)
摘要Based on the complex variable moving least-square(CVMLS) approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin(CVMLPG) method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square(MLS) approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin(MLPG) method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.
基金supported by the National Natural Science Foundation(Grant Nos.11821202,11732004,12002077 and 12002073)the National Key Research and Development Plan(Grant No.2020YFB1709401)+2 种基金the Fundamental Research Funds for the Central Universities(Grant Nos.DUT21-RC(3)076 and DUT20RC(3)020)the Doctoral Scientific Research Foundation of Liaoning Province(Grant No.2021-BS-063)111 Project(Grant No.B14013).
摘要Traditional topology optimization methods often introduce weak artificial material to mimic voids to avoid the singularity of the global stiffness matrix and carry out topology optimization with a fixed finite element(FE)mesh.This treatment,however,may not only increase the computational cost for structural analysis but also lead to unfavorable numerical instabilities,especially when large deformations and dynamic/buckling behaviors are involved.In the present work,a new meshless moving morphable component-based method(ML-MMC),which structural analysis is carried out only on the solid region occupied by components,is proposed.In this approach,the coupling of discrete components is achieved through the adaptively constructed influence domain of the meshless shape function.Therefore,the singularity problem of the stiffness matrix can be naturally avoided without introducing weak artificial material.Compared with traditional methods,the number of degrees of freedoms(DOFs)can be reduced substantially under this treatment.The effectiveness of the proposed approach is also illustrated by some representative examples.
基金supported by the National Natural Science Foundation of China (Grant No. 10871124)the Innovation Program of the Shanghai Municipal Education Commission,China (Grant No. 09ZZ99)
摘要This paper presents a meshless method for the nonlinear generalized regularized long wave (GRLW) equation based on the moving least-squares approximation. The nonlinear discrete scheme of the GRLW equation is obtained and is solved using the iteration method. A theorem on the convergence of the iterative process is presented and proved using theorems of the infinity norm. Compared with numerical methods based on mesh, the meshless method for the GRLW equation only requires the.scattered nodes instead of meshing the domain of the problem. Some examples, such as the propagation of single soliton and the interaction of two solitary waves, are given to show the effectiveness of the meshless method.
基金Project supported by the National Natural Science Foundation of China (Grant No.10871124)the Natural Science Foundation of Zhejiang Province of China (Grant No.Y6110007)
摘要The element-free Galerkin (EFG) method for numerically solving the compound Korteweg-de Vries-Burgers (KdVB) equation is discussed in this paper. The Galerkin weak form is used to obtain the discrete equation and the essential boundary conditions are enforced by the penalty method. The effectiveness of the EFG method of solving the compound Korteweg-de Vries-Burgers (KdVB) equation is illustrated by three numerical examples.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA19030402)the Key Special Projects for International Cooperation in Science and Technology Innovation between Governments(Grant No.2017YFE0133600the Beijing Municipal Natural Science Foundation Youth Project 8214066:Application Research of Beijing Road Visibility Prediction Based on Machine Learning Methods.
摘要We propose a novel machine learning approach to reconstruct meshless surface wind speed fields,i.e.,to reconstruct the surface wind speed at any location,based on meteorological background fields and geographical information.The random forest method is selected to develop the machine learning data reconstruction model(MLDRM-RF)for wind speeds over Beijing from 2015-19.We use temporal,geospatial attribute and meteorological background field features as inputs.The wind speed field can be reconstructed at any station in the region not used in the training process to cross-validate model performance.The evaluation considers the spatial distribution of and seasonal variations in the root mean squared error(RMSE)of the reconstructed wind speed field across Beijing.The average RMSE is 1.09 m s−1,considerably smaller than the result(1.29 m s−1)obtained with inverse distance weighting(IDW)interpolation.Finally,we extract the important feature permutations by the method of mean decrease in impurity(MDI)and discuss the reasonableness of the model prediction results.MLDRM-RF is a reasonable approach with excellent potential for the improved reconstruction of historical surface wind speed fields with arbitrary grid resolutions.Such a model is needed in many wind applications,such as wind energy and aviation safety assessments.
基金supported by the Natural Science Foundation of Ningbo,China (Grant Nos.2009A610014 and 2009A610154)the Natural Science Foundation of Zhejiang Province,China (Grant No.Y6090131)
摘要Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson's equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.
基金The project supported by the National Natural Science Foundation of China (10172052)
摘要Based on our previous study,the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs)are used,so that it could result in significant error in solving partial differential equations with Neumann boundary conditions.To overcome this drawback,the Consistent Compactly Supported Radial Basis Functions(CCSRBFs)are developed,which satisfy the predetermined consistency con- ditions.Meshless method based on point collocation with CCSRBFs is developed for solving partial differential equations.Numerical studies show that the proposed method improves the accuracy of approximation significantly.
基金supported by the National 973 Scientific and Technological Innovation Project (No. 2004CB719402)National Natural Science Foundation of China (No. 10672055)+3 种基金Key Project of NSFC (No. 60635020)Natural Science Foundation for Out standing Youth of China (No. 50625519)Hunan Provincial Natural Science Foundation of China (No. 07JJ6002)Scientific Research Fund of Hunan Provincial Education Department of China (No. 08C230)
摘要A meshless local radial point interpolation method (LRPIM) for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function (RBF) coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.
摘要Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature (PUQ),for Galerkin meshless methods was presented. The technique is based on finite covering and partition of unity. There is no need to decompose the physical domain into small cell. It possesses remarkable integration accuracy. Using Element-free Galerkin methods as example, Galerkin meshless methods based on PUQ were studied in detail. Meshing is always not required in the procedure of constitution of approximate function or numerical quadrature, so Galerkin meshless methods based on PUQ are “truly” meshless methods.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Natural Science Foundation of Ningbo City,China(GrantNo.2013A610103)+2 种基金the Natural Science Foundation of Zhejiang Province,China(Grant No.Y6090131)the Disciplinary Project of Ningbo City,China(GrantNo.SZXL1067)the K.C.Wong Magna Fund in Ningbo University,China
摘要Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.
基金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 (No.10172052).
摘要The meshless weighted least-square (MWLS) method was developed based on the weighted least-square method. The method possesses several advantages, such as high accuracy, high stability and high e?ciency. Moreover, the coe?cient matrix obtained is symmetric and semi- positive de?nite. In this paper, the method is further examined critically. The e?ects of several parameters on the results of MWLS are investigated systematically by using a cantilever beam and an in?nite plate with a central circular hole. The numerical results are compared with those obtained by using the collocation-based meshless method (CBMM) and Galerkin-based meshless method (GBMM). The investigated parameters include the type of approximations, the type of weight functions, the number of neighbors of an evaluation point, as well as the manner in which the neighbors of an evaluation point are determined. This study shows that the displacement accuracy and convergence rate obtained by MWLS is comparable to that of the GBMM while the stress accuracy and convergence rate yielded by MWLS is even higher than that of GBMM. Furthermore, MWLS is much more e?cient than GBMM. This study also shows that the instability of CBMM is mainly due to the neglect of the equi- librium residuals at boundary nodes. In MWLS, the residuals of all the governing equations are minimized in a weighted least-square sense.
基金Supported by Anhui Provincial Natural Science Foundation(1908085QA09)
摘要A simple direct space-time meshless scheme,based on the radial or non-radial basis function,is proposed for the onedimensional Klein-Gordon equations.Since these equations are time-dependent,it is worthwhile to present two schemes for the basis functions from radial and non-radial aspects.The first scheme is fulfilled by considering time variable as normal space variable,to construct an"isotropic"space-time radial basis function.The other scheme considered a realistic relationship between space variable and time variable which is not radial.The timedependent variable is treated regularly during the whole solution process and the Klein-Gordon equations can be solved in a direct way.Numerical results show that the proposed meshless schemes are simple,accurate,stable,easy-to-program and efficient for the Klein-Gordon equations.