In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the met...In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.展开更多
This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructe...This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.展开更多
Outgoing waves arising from high-velocity impacts between soil and structure can be reflected by the conventional truncated boundaries.Absorbing boundary conditions(ABCs),to attenuate the energy of the outward waves,a...Outgoing waves arising from high-velocity impacts between soil and structure can be reflected by the conventional truncated boundaries.Absorbing boundary conditions(ABCs),to attenuate the energy of the outward waves,are necessary to ensure the proper representation of the kinematic field and the accurate quantification of impact forces.In this paper,damping layer and dashpot ABCs are implemented in the material point method(MPM)with slight adjustments.Benchmark scenarios of different dynamic problems are modelled with the ABCs configured.Feasibility of the ABCs is assessed through the velocity fluctuations at specific observation points and the impact force fluctuations on the structures.The impact forces predicted by the MPM with ABCs are verified by comparison with those estimated using a computational fluid dynamics approach.展开更多
In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit...This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit scheme are obtained by discretizing the equivalent integral form of the equation.For the explicit scheme with nodal values,the algebraic system can be solved by the Thomas method.For the implicit scheme with both nodal values and their derivatives,the system can be implemented by a prediction-correction procedure,where in the correction stage,an implicit formula for recovering the nodal derivatives is introduced.Taking two point boundary value problem as an example,we prove that both the explicit and implicit schemes are convergent with fourth order accuracy with respect to some standard discrete norms using the energy method.Two numerical examples demonstrate the correctness and effectiveness of the schemes,as well as the indispensability of using non-uniform meshes.展开更多
Combining the radial point interpolation method(RPIM),the dual reciprocity method(DRM)and the hybrid boundary node method(HBNM),a dual reciprocity hybrid radial boundary node method(DHRBNM)is proposed for linear elast...Combining the radial point interpolation method(RPIM),the dual reciprocity method(DRM)and the hybrid boundary node method(HBNM),a dual reciprocity hybrid radial boundary node method(DHRBNM)is proposed for linear elasticity.Compared to DHBNM,RPIM is exploited to replace the moving least square(MLS)in DHRBNM,and it gets rid of the deficiency of MLS approximation,in which shape functions lack the delta function property,the boundary condition can not be applied easily and directly and it's computational expense is high.Besides,different approximate functions are discussed in DRM to get the interpolation property,in which the accuracy and efficiency for different basis functions are compared.Then RPIM is also applied in DRM to replace the conical function interpolation,which can greatly improve the accuracy of the present method.To demonstrate the effectiveness of the present method,DHBNM is applied for comparison,and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.展开更多
Volume parameter is the basic content of a spatial body object morphology analysis.However,the challenge lies in the volume calculation of irregular objects.The point cloud slicing method proposed in this study effect...Volume parameter is the basic content of a spatial body object morphology analysis.However,the challenge lies in the volume calculation of irregular objects.The point cloud slicing method proposed in this study effectively works in calculating the volume of the point cloud of the spatial object obtained through three-dimensional laser scanning(3DLS).In this method,a uniformly spaced sequent slicing process is first conducted in a specific direction on the point cloud of the spatial object obtained through 3DLS.A series of discrete point cloud slices corresponding to the point cloud bodies are then obtained.Subsequently,the outline boundary polygon of the point cloud slicing is searched one by one in accordance with the slicing sequence and areas of the polygon.The point cloud slice is also calculated.Finally,the individual point cloud section volume is calculated through the slicing areas and the adjacent slicing gap.Thus,the total volume of the scanned spatial object can be calculated by summing up the individual volumes.According to the results and analysis of the calculated examples,the slice-based volume-calculating method for the point cloud of irregular objects obtained through 3DLS is correct,concise in process,reliable in results,efficient in calculation methods,and controllable on accuracy.This method comes as a good solution to the volume calculation of irregular objects.展开更多
The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″（t） ＋ q（t）f（t,u（t）,u′（t）） = 0, 0 〈 t 〈 1, u（0） = u″（0） = u′（1） = 0 is studied....The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″（t） ＋ q（t）f（t,u（t）,u′（t）） = 0, 0 〈 t 〈 1, u（0） = u″（0） = u′（1） = 0 is studied. The iterative schemes for approximating the solutions are obtained by applying a monotone iterative method.展开更多
This paper deals with the existence of solutions to the p（t）-Laplacian equation with four-point boundary conditions. It is shown, by Leray-Schauder fixed point theorem and degree method, that under suitable conditio...This paper deals with the existence of solutions to the p（t）-Laplacian equation with four-point boundary conditions. It is shown, by Leray-Schauder fixed point theorem and degree method, that under suitable conditions, solutions of the problem exist. The interesting point is that p（t） is a general function.展开更多
Based on the analysis of the subjectivity of wetland boundary criteria and their causes at present, this paper suggested that, under the condition that the mechanism of wetland formation process has not been understoo...Based on the analysis of the subjectivity of wetland boundary criteria and their causes at present, this paper suggested that, under the condition that the mechanism of wetland formation process has not been understood, "black box" method of System Theory can be used to delineate wetland boundaries scientifically. After analyzing the difference of system construction among aquatic habitats, wetlands and uplands, the lower limit of rooted plants was chosen as the lower boundary criterion of wetlands. Because soil diagnostic horizon is the result of the long-term interaction among all environments, and it is less responsive than vegetation to short-term change, soil diagnostic horizon was chosen as the indicator to delineate wetland upper boundary, which lies at the thinning-out point of soil diagnostic horizon. Case study indicated that it was feasible using the lower limit of rooted plants and the thinning-out point of soil diagnostic horizon as criteria to delineate the lower and upper boundaries of wetland. In the study area, the thinning-out line of albic horizon was coincident with the 55.74m contour line, the maximum horizon error was less than 1m, and the maximum vertical error less than 0.04m. The problem on wetland definition always arises on the boundaries. Having delineated wetland boundaries, wetlands can be defined as follows: wetlands are the transitional zones between uplands and deepwater habitats, they are a kind of azonal complex that are inundated or saturated by surface or ground water, with the lower boundary lying at the lower limit of rooted plants, and the upper boundary at the thinning-out line of upland soil diagnostic horizon.展开更多
In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, w...In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.展开更多
In this paper, we discuss the existence of solution of a nonlinear two-point boundary value problem with a positive parameter Q arising in the study of surfacetension-induced flows of a liquid metal or semiconductor. ...In this paper, we discuss the existence of solution of a nonlinear two-point boundary value problem with a positive parameter Q arising in the study of surfacetension-induced flows of a liquid metal or semiconductor. By applying the Schauder's fixed-point theorem, we prove that the problem admits a solution for 0 ≤ Q ≤ 14.306.It improves the result of 0 ≤ Q ＜ 1 in [2] and 0 ≤ Q ≤ 13.213 in [3].展开更多
Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with...Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.展开更多
Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Br...Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Bratu’s equation, Troesch’s problems) occurs engineering and science, including the modeling of chemical reactions diffusion processes and heat transfer. An analytical expression pertaining to the concentration of substrate is obtained using Homotopy perturbation method for all values of parameters. These approximate analytical results were found to be in good agreement with the simulation results.展开更多
The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and...The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.展开更多
Two-dimensional tidal open boundary conditions of the M2 constituent in the Bohai and Yellow Seas(BYS) have been estimated by assimilating T/P altimeter data.During inversion,independent point(IP) strategy was used,in...Two-dimensional tidal open boundary conditions of the M2 constituent in the Bohai and Yellow Seas(BYS) have been estimated by assimilating T/P altimeter data.During inversion,independent point(IP) strategy was used,in which several IPs on the open boundary is assumed,values at these IPs can be optimized with an adjoint method,and those at other grid points are determined by linearly interpolating the values at IPs.The reasonability and feasibility of the model are tested by ideal twin experiments.In the practical experiment(PE) after assimilation,the cost function may reach 1% or less of its initial value.Mean absolute errors in amplitude and phase can be less than 5 cm and 5°,respectively,and the obtained co-chart can show the character of the M2 constituent in the BYS.The results of the PE indicate that using only two IPs on the open boundary can yield better simulated results.展开更多
基金supported by the National Natural Science Foundation of China (11132004 and 51078145)the Natural Science Foundation of Guangdong Province (9251064101000016)
文摘In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.
基金The authors thank the financial support of National Natural Science Foundation of China(NSFC)under Grant(No.11702238).
文摘This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.
基金the Key Science and Technology Plan of Power China Huadong Engineering Corporation(No.KY2018-ZD-01)China and the National Natural Science Foundations of China(No.51909248)。
文摘Outgoing waves arising from high-velocity impacts between soil and structure can be reflected by the conventional truncated boundaries.Absorbing boundary conditions(ABCs),to attenuate the energy of the outward waves,are necessary to ensure the proper representation of the kinematic field and the accurate quantification of impact forces.In this paper,damping layer and dashpot ABCs are implemented in the material point method(MPM)with slight adjustments.Benchmark scenarios of different dynamic problems are modelled with the ABCs configured.Feasibility of the ABCs is assessed through the velocity fluctuations at specific observation points and the impact force fluctuations on the structures.The impact forces predicted by the MPM with ABCs are verified by comparison with those estimated using a computational fluid dynamics approach.
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
文摘This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit scheme are obtained by discretizing the equivalent integral form of the equation.For the explicit scheme with nodal values,the algebraic system can be solved by the Thomas method.For the implicit scheme with both nodal values and their derivatives,the system can be implemented by a prediction-correction procedure,where in the correction stage,an implicit formula for recovering the nodal derivatives is introduced.Taking two point boundary value problem as an example,we prove that both the explicit and implicit schemes are convergent with fourth order accuracy with respect to some standard discrete norms using the energy method.Two numerical examples demonstrate the correctness and effectiveness of the schemes,as well as the indispensability of using non-uniform meshes.
基金Project supported by the National Basic Research Program of China(No.2010CB732006)the CAS/SAFEA International Partnership Program for Creative Research Teams(No.KZCX2-YW-T12)the National Natural Science Foundation of China(No.11002154)
文摘Combining the radial point interpolation method(RPIM),the dual reciprocity method(DRM)and the hybrid boundary node method(HBNM),a dual reciprocity hybrid radial boundary node method(DHRBNM)is proposed for linear elasticity.Compared to DHBNM,RPIM is exploited to replace the moving least square(MLS)in DHRBNM,and it gets rid of the deficiency of MLS approximation,in which shape functions lack the delta function property,the boundary condition can not be applied easily and directly and it's computational expense is high.Besides,different approximate functions are discussed in DRM to get the interpolation property,in which the accuracy and efficiency for different basis functions are compared.Then RPIM is also applied in DRM to replace the conical function interpolation,which can greatly improve the accuracy of the present method.To demonstrate the effectiveness of the present method,DHBNM is applied for comparison,and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.
文摘Volume parameter is the basic content of a spatial body object morphology analysis.However,the challenge lies in the volume calculation of irregular objects.The point cloud slicing method proposed in this study effectively works in calculating the volume of the point cloud of the spatial object obtained through three-dimensional laser scanning(3DLS).In this method,a uniformly spaced sequent slicing process is first conducted in a specific direction on the point cloud of the spatial object obtained through 3DLS.A series of discrete point cloud slices corresponding to the point cloud bodies are then obtained.Subsequently,the outline boundary polygon of the point cloud slicing is searched one by one in accordance with the slicing sequence and areas of the polygon.The point cloud slice is also calculated.Finally,the individual point cloud section volume is calculated through the slicing areas and the adjacent slicing gap.Thus,the total volume of the scanned spatial object can be calculated by summing up the individual volumes.According to the results and analysis of the calculated examples,the slice-based volume-calculating method for the point cloud of irregular objects obtained through 3DLS is correct,concise in process,reliable in results,efficient in calculation methods,and controllable on accuracy.This method comes as a good solution to the volume calculation of irregular objects.
基金Supported by the Natural Science Foundation of Zhejiang Province (Y605144)the XNF of Zhejiang University of Media and Communications (XN080012008034)
文摘The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″（t） ＋ q（t）f（t,u（t）,u′（t）） = 0, 0 〈 t 〈 1, u（0） = u″（0） = u′（1） = 0 is studied. The iterative schemes for approximating the solutions are obtained by applying a monotone iterative method.
基金The NSF(11271154)of Chinathe Key Lab of Symbolic Computation and Knowledge Engineering of Ministry of Education+1 种基金the 985 program of Jilin Universitythe DR Fund(150152)of Henan University of Technology
文摘This paper deals with the existence of solutions to the p（t）-Laplacian equation with four-point boundary conditions. It is shown, by Leray-Schauder fixed point theorem and degree method, that under suitable conditions, solutions of the problem exist. The interesting point is that p（t） is a general function.
基金Under the auspices of the Knowledge Innovation Program of Chinese Academy of Sciences(No.KZCX3-SW-NA-01)
文摘Based on the analysis of the subjectivity of wetland boundary criteria and their causes at present, this paper suggested that, under the condition that the mechanism of wetland formation process has not been understood, "black box" method of System Theory can be used to delineate wetland boundaries scientifically. After analyzing the difference of system construction among aquatic habitats, wetlands and uplands, the lower limit of rooted plants was chosen as the lower boundary criterion of wetlands. Because soil diagnostic horizon is the result of the long-term interaction among all environments, and it is less responsive than vegetation to short-term change, soil diagnostic horizon was chosen as the indicator to delineate wetland upper boundary, which lies at the thinning-out point of soil diagnostic horizon. Case study indicated that it was feasible using the lower limit of rooted plants and the thinning-out point of soil diagnostic horizon as criteria to delineate the lower and upper boundaries of wetland. In the study area, the thinning-out line of albic horizon was coincident with the 55.74m contour line, the maximum horizon error was less than 1m, and the maximum vertical error less than 0.04m. The problem on wetland definition always arises on the boundaries. Having delineated wetland boundaries, wetlands can be defined as follows: wetlands are the transitional zones between uplands and deepwater habitats, they are a kind of azonal complex that are inundated or saturated by surface or ground water, with the lower boundary lying at the lower limit of rooted plants, and the upper boundary at the thinning-out line of upland soil diagnostic horizon.
基金The Project was supported by National Natural Science Foundation of China
文摘In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.
基金The work was supported by National Natural Science Foundation(Grant No. 10471129) of China
文摘In this paper, we discuss the existence of solution of a nonlinear two-point boundary value problem with a positive parameter Q arising in the study of surfacetension-induced flows of a liquid metal or semiconductor. By applying the Schauder's fixed-point theorem, we prove that the problem admits a solution for 0 ≤ Q ≤ 14.306.It improves the result of 0 ≤ Q ＜ 1 in [2] and 0 ≤ Q ≤ 13.213 in [3].
文摘Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.
文摘Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Bratu’s equation, Troesch’s problems) occurs engineering and science, including the modeling of chemical reactions diffusion processes and heat transfer. An analytical expression pertaining to the concentration of substrate is obtained using Homotopy perturbation method for all values of parameters. These approximate analytical results were found to be in good agreement with the simulation results.
文摘The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.
基金Supported by the State Ministry of Science and Technology of China (Nos.2007AA09Z118,2008AA09A402)the National Natural Science Foundation of China(No.41076006)the Ministry of Education's 111 Project(No.B07036)
文摘Two-dimensional tidal open boundary conditions of the M2 constituent in the Bohai and Yellow Seas(BYS) have been estimated by assimilating T/P altimeter data.During inversion,independent point(IP) strategy was used,in which several IPs on the open boundary is assumed,values at these IPs can be optimized with an adjoint method,and those at other grid points are determined by linearly interpolating the values at IPs.The reasonability and feasibility of the model are tested by ideal twin experiments.In the practical experiment(PE) after assimilation,the cost function may reach 1% or less of its initial value.Mean absolute errors in amplitude and phase can be less than 5 cm and 5°,respectively,and the obtained co-chart can show the character of the M2 constituent in the BYS.The results of the PE indicate that using only two IPs on the open boundary can yield better simulated results.
