The aim of this paper is to obtain the numerical solutions of generalized space-fractional Burgers' equations with initial-boundary conditions by the Jacobi spectral collocation method using the shifted Jacobi-Gau...The aim of this paper is to obtain the numerical solutions of generalized space-fractional Burgers' equations with initial-boundary conditions by the Jacobi spectral collocation method using the shifted Jacobi-Gauss-Lobatto collocation points. By means of the simplifed Jacobi operational matrix, we produce the diferentiation matrix and transfer the space-fractional Burgers' equation into a system of ordinary diferential equations that can be solved by the fourth-order Runge-Kutta method. The numerical simulations indicate that the Jacobi spectral collocation method is highly accurate and fast convergent for the generalized space-fractional Burgers' equation.展开更多
Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of...Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.展开更多
In this paper, Radial point collocation method (RPCM), a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local suppo...In this paper, Radial point collocation method (RPCM), a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local supported domains based on radial basis functions. As a result, this method is local and hence the system matrix is banded which is very attractive for practical engineering problems. In the numerical examination, RPCM is applied to solve non-linear convection-diffusion 2D Burgers equations. The results obtained by RPCM demonstrate the accuracy and efficiency of the proposed method for solving transient fluid dynamic problems. A fictitious point scheme is adopted to improve the solution accuracy while Neumann boundary conditions exist. The meshfree feature of the nresent method is verv attractive in solving comnutational fluid nroblems.展开更多
A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a...A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM (finite element method) and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.展开更多
This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredepend...This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.展开更多
In this article,we approximate the solution of high order linear Fredholm integro-differential equations with a variable coefficient under the initial-boundary conditions by Bell polynomials.Using collocation points a...In this article,we approximate the solution of high order linear Fredholm integro-differential equations with a variable coefficient under the initial-boundary conditions by Bell polynomials.Using collocation points and treating the solution as a linear combination of Bell polynomials,the problem is reduced to linear system of equations whose unknown variables are Bell coefficients.The solution to this algebraic system determines the approximate solution.Error estimation of approximate solution is done.Some examples are provided to illustrate the performance of the method.The numerical results are compared with the collocation method based on Legendre polynomials and the other two methods based on Taylor polynomials.It is observed that the method is better than Legendre collocation method and as accurate as the methods involving Taylor polynomials.展开更多
With the terms of the exact series solution taken as trial functions, the method of point collocation was used to calculate the large deflection of a circular plate. The axisymmetrical bending formulae were developed ...With the terms of the exact series solution taken as trial functions, the method of point collocation was used to calculate the large deflection of a circular plate. The axisymmetrical bending formulae were developed for the calculation of a circular plate subjected to polynomial distributed loads, a concentrated load at the center, uniform radial forces and moments along the edge or their combinations. The support may be elastic. The buckling load was calculated. Under action of uniformly distributed load, central load or their compound load, solutions were compared with those obtained by other methods. Buckling beyond critical thrust was compared with that calculated by the power series method. The method presented in this paper has advantages of wide convergent range, high precision and short computing time. Moreover, the computing time is nearly independent of the complexity of the loads.展开更多
In this paper, we propose to replace the Chebyshev series used in pseudospectral methods with the equivalent Chebyshev economized power series that can be evaluated more rapidly. We keep the rest of the implementation...In this paper, we propose to replace the Chebyshev series used in pseudospectral methods with the equivalent Chebyshev economized power series that can be evaluated more rapidly. We keep the rest of the implementation the same as the spectral method so that there is no new mathematical principle involved. We show by numerical examples that the new approach works well and there is indeed no significant loss of solution accuracy. The advantages of using power series also include simplicity in its formulation and implementation such that it could be used for complex systems. We investigate the important issue of collocation point selection. Our numerical results indicate that there is a clear accuracy advantage of using collocation points corresponding to roots of the Chebyshev polynomial.展开更多
In this paper, the nonlinear equations of motion for shallow spherical shells with axisymmetric deformation including transverse shear are derived. The nonlinear static and dynamic response and dynamic buckling of sha...In this paper, the nonlinear equations of motion for shallow spherical shells with axisymmetric deformation including transverse shear are derived. The nonlinear static and dynamic response and dynamic buckling of shallow spherical shells with circular hole on elastically restrained edge are investigated. By using the orthogonal point collocation method for space and Newmarh-β scheme for time, the displacement functions are separated and the nonlinear differential equations are replaced by linear algebraic equations to seek solutions. The numerical results are presented for different cases and compared with available data.展开更多
In this paper, an exponential method is presented for the approximate solutions of the HIV infection model of CD4+T. The method is based on exponential polynomi- als and collocation points. This model problem corresp...In this paper, an exponential method is presented for the approximate solutions of the HIV infection model of CD4+T. The method is based on exponential polynomi- als and collocation points. This model problem corresponds to a system of nonlinear ordinary differential equations. Matrix relations are constructed for the exponential functions. By aid of these matrix relations and the collocation points, the proposed technique transforms the model problem into a system of nonlinear algebraic equations. By solving the system of the algebraic equations, the unknown coefficients are com- puted and thus the approximate solutions are obtained. The applications of the method for the considered problem are given and the comparisons are made with the other methods.展开更多
基金This work is supported by the National Natural Science Foundation of China(Grant Nos.11701358,11774218)。
文摘The aim of this paper is to obtain the numerical solutions of generalized space-fractional Burgers' equations with initial-boundary conditions by the Jacobi spectral collocation method using the shifted Jacobi-Gauss-Lobatto collocation points. By means of the simplifed Jacobi operational matrix, we produce the diferentiation matrix and transfer the space-fractional Burgers' equation into a system of ordinary diferential equations that can be solved by the fourth-order Runge-Kutta method. The numerical simulations indicate that the Jacobi spectral collocation method is highly accurate and fast convergent for the generalized space-fractional Burgers' equation.
文摘Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.
基金Project (No. 10572128) supported by the National Natural ScienceFoundation of China
文摘In this paper, Radial point collocation method (RPCM), a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local supported domains based on radial basis functions. As a result, this method is local and hence the system matrix is banded which is very attractive for practical engineering problems. In the numerical examination, RPCM is applied to solve non-linear convection-diffusion 2D Burgers equations. The results obtained by RPCM demonstrate the accuracy and efficiency of the proposed method for solving transient fluid dynamic problems. A fictitious point scheme is adopted to improve the solution accuracy while Neumann boundary conditions exist. The meshfree feature of the nresent method is verv attractive in solving comnutational fluid nroblems.
基金the National Natural Science Foundation of China (No. 50275059).
文摘A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM (finite element method) and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.
基金supported by a grant from the National Science and Technology Council of the Republic of China(Grant Number:MOST 112-2221-E-006-048-MY2).
文摘This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.
文摘In this article,we approximate the solution of high order linear Fredholm integro-differential equations with a variable coefficient under the initial-boundary conditions by Bell polynomials.Using collocation points and treating the solution as a linear combination of Bell polynomials,the problem is reduced to linear system of equations whose unknown variables are Bell coefficients.The solution to this algebraic system determines the approximate solution.Error estimation of approximate solution is done.Some examples are provided to illustrate the performance of the method.The numerical results are compared with the collocation method based on Legendre polynomials and the other two methods based on Taylor polynomials.It is observed that the method is better than Legendre collocation method and as accurate as the methods involving Taylor polynomials.
文摘With the terms of the exact series solution taken as trial functions, the method of point collocation was used to calculate the large deflection of a circular plate. The axisymmetrical bending formulae were developed for the calculation of a circular plate subjected to polynomial distributed loads, a concentrated load at the center, uniform radial forces and moments along the edge or their combinations. The support may be elastic. The buckling load was calculated. Under action of uniformly distributed load, central load or their compound load, solutions were compared with those obtained by other methods. Buckling beyond critical thrust was compared with that calculated by the power series method. The method presented in this paper has advantages of wide convergent range, high precision and short computing time. Moreover, the computing time is nearly independent of the complexity of the loads.
文摘In this paper, we propose to replace the Chebyshev series used in pseudospectral methods with the equivalent Chebyshev economized power series that can be evaluated more rapidly. We keep the rest of the implementation the same as the spectral method so that there is no new mathematical principle involved. We show by numerical examples that the new approach works well and there is indeed no significant loss of solution accuracy. The advantages of using power series also include simplicity in its formulation and implementation such that it could be used for complex systems. We investigate the important issue of collocation point selection. Our numerical results indicate that there is a clear accuracy advantage of using collocation points corresponding to roots of the Chebyshev polynomial.
文摘In this paper, the nonlinear equations of motion for shallow spherical shells with axisymmetric deformation including transverse shear are derived. The nonlinear static and dynamic response and dynamic buckling of shallow spherical shells with circular hole on elastically restrained edge are investigated. By using the orthogonal point collocation method for space and Newmarh-β scheme for time, the displacement functions are separated and the nonlinear differential equations are replaced by linear algebraic equations to seek solutions. The numerical results are presented for different cases and compared with available data.
文摘In this paper, an exponential method is presented for the approximate solutions of the HIV infection model of CD4+T. The method is based on exponential polynomi- als and collocation points. This model problem corresponds to a system of nonlinear ordinary differential equations. Matrix relations are constructed for the exponential functions. By aid of these matrix relations and the collocation points, the proposed technique transforms the model problem into a system of nonlinear algebraic equations. By solving the system of the algebraic equations, the unknown coefficients are com- puted and thus the approximate solutions are obtained. The applications of the method for the considered problem are given and the comparisons are made with the other methods.