Precise integration method for a class of singular two-point boundary value problems

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.

Machine Learning Enhanced Boundary Element Method:Prediction of Gaussian Quadrature Points

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.

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

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.

A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

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.

Three Dimensional Laser Point Cloud Slicing Method for Calculating Irregular Volume

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.

Existence and iteration of monotone positive solutions for a third-order two-point boundary value problem

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.

Existence of Solutions for a Four-point Boundary Value Problem with a p(t)-Laplacian

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.

THEORY AND METHOD FOR WETLAND BOUNDARY DELINEATION

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.

ASYMPTOTIC ERROR EXPANSIONS OF QUADRATIC SPLINE COLLOCATION SOLUTIONS FOR TWO-POINT BOUNDARY VALUE PROBLEMS

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.

ON THE EXISTENCE OF SOLUTION OF A NONLINEAR TWO-POINT BOUNDARY VALUE PROBLEM ARISING FROM A LIQUID METAL FLOW

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].

Analytical Solutions of Some Two-Point Non-Linear Elliptic Boundary Value Problems

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.

Efficient Decomposition Shooting Method for Solving Third-Order Boundary Value Problems

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.

基于边界表示模型的点核积分程序的开发与应用

点核积分方法是屏蔽设计中最常用的方法之一。为了克服传统点核积分方法因无法直接使用辐射防护现场基于边界表示(Brep,boundary representation)几何需重新建模而出现模型描述复杂且易出错的问题,本次研究开发了直接基于Brep几何的点核积分方法,可直接基于现场CAD模型进行辐射场计算,提升了点核积分方法的几何适应能力。使用某乏燃料运输容器基准问题对程序的准确性进行了验证,验证结果表明程序的计算结果与QAD程序吻合良好,使用华龙一号的主回路设备间屏蔽问题对程序的复杂几何适应能力进行验证,验证结果表明程序具备复杂场景的三维辐射场计算能力。

Inversion of two-dimensional tidal open boundary conditions of M_2 constituent in the Bohai and Yellow Seas

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.

Effects of slip condition on MHD stagnation-point flow over a power-law stretching sheet

The steady two-dimensional magnetohydrodynamic stagnation flow towards a nonlinear stretching surface is studied. The no-slip condition on the solid boundary is replaced with a partial slip condition. A scaling group transformation is used to get the invariants. Using the invariants, a third-order ordinary differential equation corresponding to the momentum is obtained. An analytical solution is obtained in a series form using a homotopy analysis method. Reliability and efficiency of series solutions are shown by the good agreement with numerical results presented in the literature. The effects of the slip parameter, the magnetic field parameter, the velocity ratio parameter, the suction velocity parameter, and the power law exponent on the flow are investigated. The results show that the velocity and shear stress profiles are greatly influenced by these parameters.

Landslide data mosaicking based on an airborne laser point cloud and multi-beam sonar images

Landslides are one of the most disastrous geological hazards in southwestern China.Once a landslide becomes unstable,it threatens the lives and safety of local residents.However,empirical studies on landslides have predominantly focused on landslides that occur on land.To this end,we aim to investigate ashore and underwater landslide data synchronously.This study proposes an optimized mosaicking method for ashore and underwater landslide data.This method fuses an airborne laser point cloud with multi-beam depth sounder images.Owing to their relatively high efficiency and large coverage area,airborne laser measurement systems are suitable for emergency investigations of landslides.Based on the airborne laser point cloud,the traversal of the point with the lowest elevation value in the point set can be used to perform rapid extraction of the crude channel boundaries.Further meticulous extraction of the channel boundaries is then implemented using the probability mean value optimization method.In addition,synthesis of the integrated ashore and underwater landslide data angle is realized using the spatial guide line between the channel boundaries and the underwater multibeam sonar images.A landslide located on the right bank of the middle reaches of the Yalong River is selected as a case study to demonstrate that the proposed method has higher precision thantraditional methods.The experimental results show that the mosaicking method in this study can meet the basic needs of landslide modeling and provide a basis for qualitative and quantitative analysis and stability prediction of landslides.

Analytical solution to stagnation-point flow and heat transfer over a stretching sheet based on homotopy analysis

This paper is concerned with two-dimensional stagnation-point steady flow of an incompressible viscous fluid towards a stretching sheet whose velocity is proportional to the distance from the slit. The governing system of partial differential equations is first transformed into a system of dimensionless ordinary differential equations. Analytical solutions of the velocity distribution and dimensionless temperature profiles are obtained for different ratios of free stream velocity and stretching velocity, Prandtl number, Eckert number and dimensionality index in series forms using homotopy analysis method（HAM）. It is shown that a boundary layer is formed when the free stream velocity exceeds the stretching velocity, and an inverted boundary layer is formed when the free stream velocity is less than the stretching velocity. Graphs are presented to show the effects of different parameters.

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.

Precise integration method for solving singular perturbation problems

This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.

Simple Singular Perturbation Problems with Turning Points

The paper considers the asymptotic solution of two-point boundary value problems εy” + A(x)y’ = 0, 0 ≤ x ≤ 1, when 0 1, A(x) is smooth with isolated zeros, y(0) = 0 and y(1) = 1. By using perturbation method, the limit asymptotic solutions of various cases are obtained. We provide a reliable and direct method for solving similar problems. The limiting solutions are constants in this paper, except in narrow boundary and interior layers of nonuniform convergence. These provide simple examples of boundary layer resonance.