Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, p...Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.展开更多
In this paper, supose Γ be a boundary of a Jordan domain D and Γ satisfied Альпер condition, the order that rational type interpolating operators at Fejer's points of f(z)∈C(Γ) converge to f(z) in the se...In this paper, supose Γ be a boundary of a Jordan domain D and Γ satisfied Альпер condition, the order that rational type interpolating operators at Fejer's points of f(z)∈C(Γ) converge to f(z) in the sense of uniformly convergence is obtained.展开更多
Suppose that the outer mapping function of domain D has its second continuous derivatives. In this paper, the order proximation by (0,1,…,q) Hermite-Fejer interpolating polynomials at nearly Fejer's points of fun...Suppose that the outer mapping function of domain D has its second continuous derivatives. In this paper, the order proximation by (0,1,…,q) Hermite-Fejer interpolating polynomials at nearly Fejer's points of function of class A(D) are presented. Moreover in general the order of approximation is sharp.展开更多
An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in th...An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.展开更多
The results of accurate order of uniform approximation and simultaneous approximation in the early work "Jackson Type Theorems on Complex Curves" are improved from Fejer points to disturbed Fejer points in this arti...The results of accurate order of uniform approximation and simultaneous approximation in the early work "Jackson Type Theorems on Complex Curves" are improved from Fejer points to disturbed Fejer points in this article. Furthermore, the stability of convergence of Tn,∈(f,z) with disturbed sample values f(z^*) + Sk are also proved in this article.展开更多
The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their ...The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.展开更多
In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral c...In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.展开更多
Parameter identification problem is one of essential problem in order to model effectively experimental data by fractal interpolation function.In this paper,we first present an example to explain a relationship betwee...Parameter identification problem is one of essential problem in order to model effectively experimental data by fractal interpolation function.In this paper,we first present an example to explain a relationship between iteration procedure and fractal function.Then we discuss conditions that vertical scaling factors must obey in one typical case.展开更多
Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materia...Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materials constituting the Gobi result in notable differences in saltation processes across various Gobi surfaces.It is challenging to describe these processes according to a uniform morphology.Therefore,it becomes imperative to articulate surface characteristics through parameters such as the three-dimensional(3D)size and shape of gravel.Collecting morphology information for Gobi gravels is essential for studying its genesis and sand saltation.To enhance the efficiency and information yield of gravel parameter measurements,this study conducted field experiments in the Gobi region across Dunhuang City,Guazhou County,and Yumen City(administrated by Jiuquan City),Gansu Province,China in March 2023.A research framework and methodology for measuring 3D parameters of gravel using point cloud were developed,alongside improved calculation formulas for 3D parameters including gravel grain size,volume,flatness,roundness,sphericity,and equivalent grain size.Leveraging multi-view geometry technology for 3D reconstruction allowed for establishing an optimal data acquisition scheme characterized by high point cloud reconstruction efficiency and clear quality.Additionally,the proposed methodology incorporated point cloud clustering,segmentation,and filtering techniques to isolate individual gravel point clouds.Advanced point cloud algorithms,including the Oriented Bounding Box(OBB),point cloud slicing method,and point cloud triangulation,were then deployed to calculate the 3D parameters of individual gravels.These systematic processes allow precise and detailed characterization of individual gravels.For gravel grain size and volume,the correlation coefficients between point cloud and manual measurements all exceeded 0.9000,confirming the feasibility of the proposed methodology for measuring 3D parameters of individual gravels.The proposed workflow yields accurate calculations of relevant parameters for Gobi gravels,providing essential data support for subsequent studies on Gobi environments.展开更多
In this paper, we consider dynamical system, in the presence of parameter uncertainties. We apply max-min principles to determine the saddle point solution for the class of differential game arising from the associate...In this paper, we consider dynamical system, in the presence of parameter uncertainties. We apply max-min principles to determine the saddle point solution for the class of differential game arising from the associated dynamical system. We also provide sufficient condition for the existence of this saddle point.展开更多
In this paper an error in[4]is pointed out and a method for constructingsurface interpolating scattered data points is presented.The main feature of the methodin this paper is that the surface so constructed is polyno...In this paper an error in[4]is pointed out and a method for constructingsurface interpolating scattered data points is presented.The main feature of the methodin this paper is that the surface so constructed is polynomial,which makes the construction simple and the calculation easy.展开更多
In computer aided geometric design(CAGD) ,it is often needed to produce a convexity-preserving interpolating curve according to the given planar data points. However,most existing pertinent methods cannot generate con...In computer aided geometric design(CAGD) ,it is often needed to produce a convexity-preserving interpolating curve according to the given planar data points. However,most existing pertinent methods cannot generate convexity-preserving in-terpolating transcendental curves;even constructing convexity-preserving interpolating polynomial curves,it is required to solve a system of equations or recur to a complicated iterative process. The method developed in this paper overcomes the above draw-backs. The basic idea is:first to construct a kind of trigonometric polynomial curves with a shape parameter,and interpolating trigonometric polynomial parametric curves with C2(or G1) continuity can be automatically generated without having to solve any system of equations or do any iterative computation. Then,the convexity of the constructed curves can be guaranteed by the appropriate value of the shape parameter. Performing the method is easy and fast,and the curvature distribution of the resulting interpolating curves is always well-proportioned. Several numerical examples are shown to substantiate that our algorithm is not only correct but also usable.展开更多
The simulating exactly compared with realty of ground surface to run a model is more and more highly required. In the real, terrain of the earth surface is always complicated by the natural and human made ground objec...The simulating exactly compared with realty of ground surface to run a model is more and more highly required. In the real, terrain of the earth surface is always complicated by the natural and human made ground objects. Because of limitation of collecting and storing technologies in the past time, data are usually not detailed so that the data can not be full for the simulation. Besides computing power and simulation increase more day by day, the increasing requirements more detailed of topography surface simulation is a demand. In simulated flooding phenomenon or phenomena related to energy and momentum of water flow, the linear objects of ground surface such as roads, dikes, dams, etc. need to have their vertical dimension along continuously. However, these datas have often no height information alternately, there are only discrete elevation points that are extracted from topographic maps. Consequently, the demand of a suitable method for linear objects height interpolation is necessary. This paper aims to provide a method and evaluate its accuracy to meet this requirement.展开更多
This paper considered the optimal control problem for distributed parameter systems with mixed phase-control constraints and end-point constraints. Pontryagin's maximum principle for optimal control are derived vi...This paper considered the optimal control problem for distributed parameter systems with mixed phase-control constraints and end-point constraints. Pontryagin's maximum principle for optimal control are derived via Duboviskij-Milujin theorem.展开更多
A method determining vertical scaling parameters of fractal interpolation is given in this paper. By computer experiments, it is clear that this method is very effective.
Light detection and ranging(LiDAR)has contributed immensely to forest mapping and 3D tree modelling.From the perspective of data acquisition,the integration of LiDAR data from different platforms would enrich forest i...Light detection and ranging(LiDAR)has contributed immensely to forest mapping and 3D tree modelling.From the perspective of data acquisition,the integration of LiDAR data from different platforms would enrich forest information at the tree and plot levels.This research develops a general framework to integrate ground-based and UAV-LiDAR(ULS)data to better estimate tree parameters based on quantitative structure modelling(QSM).This is accomplished in three sequential steps.First,the ground-based/ULS LiDAR data were co-registered based on the local density peaks of the clustered canopy.Next,redundancy and noise were removed for the ground-based/ULS LiDAR data fusion.Finally,tree modeling and biophysical parameter retrieval were based on QSM.Experiments were performed for Backpack/Handheld/UAV-based multi-platform mobile LiDAR data of a subtropical forest,including poplar and dawn redwood species.Generally,ground-based/ULS LiDAR data fusion outperforms ground-based LiDAR with respect to tree parameter estimation compared to field data.The fusion-derived tree height,tree volume,and crown volume significantly improved by up to 9.01%,5.28%,and 18.61%,respectively,in terms of rRMSE.By contrast,the diameter at breast height(DBH)is the parameter that has the least benefits from fusion,and rRMSE remains approximately the same,because stems are already well sampled from ground data.Additionally,particularly for dense forests,the fusion-derived tree parameters were improved compared to those derived from ground-based LiDAR.Ground-based LiDAR can potentially be used to estimate tree parameters in low-stand-density forests,whereby the improvement owing to fusion is not significant.展开更多
This paper presents a method for creating modificable quartic and quintic curves with shape parameters. The curves can achieve C 2 even C 3 continuity and unify both interpolation and approximation to the control poin...This paper presents a method for creating modificable quartic and quintic curves with shape parameters. The curves can achieve C 2 even C 3 continuity and unify both interpolation and approximation to the control points without solving a system of equations or inserting additional control points. They have the local properties like the cubic B spline. Besides, the quintic curve would be able globally to tend the control polygon.展开更多
Spatial interpolation has been frequently encountered in earth sciences and engineering.A reasonable appraisal of subsurface heterogeneity plays a significant role in planning,risk assessment and decision making for g...Spatial interpolation has been frequently encountered in earth sciences and engineering.A reasonable appraisal of subsurface heterogeneity plays a significant role in planning,risk assessment and decision making for geotechnical practice.Geostatistics is commonly used to interpolate spatially varying properties at un-sampled locations from scatter measurements.However,successful application of classic geostatistical models requires prior characterization of spatial auto-correlation structures,which poses a great challenge for unexperienced engineers,particularly when only limited measurements are available.Data-driven machine learning methods,such as radial basis function network(RBFN),require minimal human intervention and provide effective alternatives for spatial interpolation of non-stationary and non-Gaussian data,particularly when measurements are sparse.Conventional RBFN,however,is direction independent(i.e.isotropic)and cannot quantify prediction uncertainty in spatial interpolation.In this study,an ensemble RBFN method is proposed that not only allows geotechnical anisotropy to be properly incorporated,but also quantifies uncertainty in spatial interpolation.The proposed method is illustrated using numerical examples of cone penetration test(CPT)data,which involve interpolation of a 2D CPT cross-section from limited continuous 1D CPT soundings in the vertical direction.In addition,a comparative study is performed to benchmark the proposed ensemble RBFN with two other non-parametric data-driven approaches,namely,Multiple Point Statistics(MPS)and Bayesian Compressive Sensing(BCS).The results reveal that the proposed ensemble RBFN provides a better estimation of spatial patterns and associated prediction uncertainty at un-sampled locations when a reasonable amount of data is available as input.Moreover,the prediction accuracy of all the three methods improves as the number of measurements increases,and vice versa.It is also found that BCS prediction is less sensitive to the number of measurement data and outperforms RBFN and MPS when only limited point observations are available.展开更多
Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and the...Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.展开更多
The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the ...The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method (FEM) as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline (TPS) radial basis functions.展开更多
基金supported by the National Natural Science Foundation of China (Grant No.40874052)the Key Laboratory of Geo-detection (China University of Geosciences,Beijing),Ministry of Education
文摘Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.
文摘In this paper, supose Γ be a boundary of a Jordan domain D and Γ satisfied Альпер condition, the order that rational type interpolating operators at Fejer's points of f(z)∈C(Γ) converge to f(z) in the sense of uniformly convergence is obtained.
文摘Suppose that the outer mapping function of domain D has its second continuous derivatives. In this paper, the order proximation by (0,1,…,q) Hermite-Fejer interpolating polynomials at nearly Fejer's points of function of class A(D) are presented. Moreover in general the order of approximation is sharp.
基金supported by the National Natural Science Foundation of China(Grant No.11171208)the Natural Science Foundation of Shanxi Province,China(Grant No.2013011022-6)
文摘An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.
基金Supported by NSF of Henan Province of China(20001110001)
文摘The results of accurate order of uniform approximation and simultaneous approximation in the early work "Jackson Type Theorems on Complex Curves" are improved from Fejer points to disturbed Fejer points in this article. Furthermore, the stability of convergence of Tn,∈(f,z) with disturbed sample values f(z^*) + Sk are also proved in this article.
文摘The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.
基金The research for this paper was supported by(1)the National Natural Science Foundation of China(Grants Nos.51708429,51708428)the Open Projects Foundation(Grant No.2017-04-GF)of State Key Laboratory for Health and Safety of Bridge Structures+1 种基金Wuhan Institute of Technology Science Found(Grant No.K201734)the science and technology projects of Wuhan Urban and Rural Construction Bureau(Grants Nos.201831,201919).
文摘In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.
文摘Parameter identification problem is one of essential problem in order to model effectively experimental data by fractal interpolation function.In this paper,we first present an example to explain a relationship between iteration procedure and fractal function.Then we discuss conditions that vertical scaling factors must obey in one typical case.
基金funded by the National Natural Science Foundation of China(42071014).
文摘Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materials constituting the Gobi result in notable differences in saltation processes across various Gobi surfaces.It is challenging to describe these processes according to a uniform morphology.Therefore,it becomes imperative to articulate surface characteristics through parameters such as the three-dimensional(3D)size and shape of gravel.Collecting morphology information for Gobi gravels is essential for studying its genesis and sand saltation.To enhance the efficiency and information yield of gravel parameter measurements,this study conducted field experiments in the Gobi region across Dunhuang City,Guazhou County,and Yumen City(administrated by Jiuquan City),Gansu Province,China in March 2023.A research framework and methodology for measuring 3D parameters of gravel using point cloud were developed,alongside improved calculation formulas for 3D parameters including gravel grain size,volume,flatness,roundness,sphericity,and equivalent grain size.Leveraging multi-view geometry technology for 3D reconstruction allowed for establishing an optimal data acquisition scheme characterized by high point cloud reconstruction efficiency and clear quality.Additionally,the proposed methodology incorporated point cloud clustering,segmentation,and filtering techniques to isolate individual gravel point clouds.Advanced point cloud algorithms,including the Oriented Bounding Box(OBB),point cloud slicing method,and point cloud triangulation,were then deployed to calculate the 3D parameters of individual gravels.These systematic processes allow precise and detailed characterization of individual gravels.For gravel grain size and volume,the correlation coefficients between point cloud and manual measurements all exceeded 0.9000,confirming the feasibility of the proposed methodology for measuring 3D parameters of individual gravels.The proposed workflow yields accurate calculations of relevant parameters for Gobi gravels,providing essential data support for subsequent studies on Gobi environments.
文摘In this paper, we consider dynamical system, in the presence of parameter uncertainties. We apply max-min principles to determine the saddle point solution for the class of differential game arising from the associated dynamical system. We also provide sufficient condition for the existence of this saddle point.
文摘In this paper an error in[4]is pointed out and a method for constructingsurface interpolating scattered data points is presented.The main feature of the methodin this paper is that the surface so constructed is polynomial,which makes the construction simple and the calculation easy.
基金Project supported by the National Basic Research Program (973) of China (No. 2004CB719400)the National Natural Science Founda-tion of China (Nos. 60673031 and 60333010) the National Natural Science Foundation for Innovative Research Groups of China (No. 60021201)
文摘In computer aided geometric design(CAGD) ,it is often needed to produce a convexity-preserving interpolating curve according to the given planar data points. However,most existing pertinent methods cannot generate convexity-preserving in-terpolating transcendental curves;even constructing convexity-preserving interpolating polynomial curves,it is required to solve a system of equations or recur to a complicated iterative process. The method developed in this paper overcomes the above draw-backs. The basic idea is:first to construct a kind of trigonometric polynomial curves with a shape parameter,and interpolating trigonometric polynomial parametric curves with C2(or G1) continuity can be automatically generated without having to solve any system of equations or do any iterative computation. Then,the convexity of the constructed curves can be guaranteed by the appropriate value of the shape parameter. Performing the method is easy and fast,and the curvature distribution of the resulting interpolating curves is always well-proportioned. Several numerical examples are shown to substantiate that our algorithm is not only correct but also usable.
文摘The simulating exactly compared with realty of ground surface to run a model is more and more highly required. In the real, terrain of the earth surface is always complicated by the natural and human made ground objects. Because of limitation of collecting and storing technologies in the past time, data are usually not detailed so that the data can not be full for the simulation. Besides computing power and simulation increase more day by day, the increasing requirements more detailed of topography surface simulation is a demand. In simulated flooding phenomenon or phenomena related to energy and momentum of water flow, the linear objects of ground surface such as roads, dikes, dams, etc. need to have their vertical dimension along continuously. However, these datas have often no height information alternately, there are only discrete elevation points that are extracted from topographic maps. Consequently, the demand of a suitable method for linear objects height interpolation is necessary. This paper aims to provide a method and evaluate its accuracy to meet this requirement.
文摘This paper considered the optimal control problem for distributed parameter systems with mixed phase-control constraints and end-point constraints. Pontryagin's maximum principle for optimal control are derived via Duboviskij-Milujin theorem.
基金This paper is partly supported by the funds of Beijing Education Commission (00KG-125)and Xi'an University of Technology.
文摘A method determining vertical scaling parameters of fractal interpolation is given in this paper. By computer experiments, it is clear that this method is very effective.
基金supported by the National Natural Science Foundation of China(Project No.42171361)the Research Grants Council of the Hong Kong Special Administrative Region,China,under Project PolyU 25211819the Hong Kong Polytechnic University under Projects 1-ZE8E and 1-ZVN6.
文摘Light detection and ranging(LiDAR)has contributed immensely to forest mapping and 3D tree modelling.From the perspective of data acquisition,the integration of LiDAR data from different platforms would enrich forest information at the tree and plot levels.This research develops a general framework to integrate ground-based and UAV-LiDAR(ULS)data to better estimate tree parameters based on quantitative structure modelling(QSM).This is accomplished in three sequential steps.First,the ground-based/ULS LiDAR data were co-registered based on the local density peaks of the clustered canopy.Next,redundancy and noise were removed for the ground-based/ULS LiDAR data fusion.Finally,tree modeling and biophysical parameter retrieval were based on QSM.Experiments were performed for Backpack/Handheld/UAV-based multi-platform mobile LiDAR data of a subtropical forest,including poplar and dawn redwood species.Generally,ground-based/ULS LiDAR data fusion outperforms ground-based LiDAR with respect to tree parameter estimation compared to field data.The fusion-derived tree height,tree volume,and crown volume significantly improved by up to 9.01%,5.28%,and 18.61%,respectively,in terms of rRMSE.By contrast,the diameter at breast height(DBH)is the parameter that has the least benefits from fusion,and rRMSE remains approximately the same,because stems are already well sampled from ground data.Additionally,particularly for dense forests,the fusion-derived tree parameters were improved compared to those derived from ground-based LiDAR.Ground-based LiDAR can potentially be used to estimate tree parameters in low-stand-density forests,whereby the improvement owing to fusion is not significant.
文摘This paper presents a method for creating modificable quartic and quintic curves with shape parameters. The curves can achieve C 2 even C 3 continuity and unify both interpolation and approximation to the control points without solving a system of equations or inserting additional control points. They have the local properties like the cubic B spline. Besides, the quintic curve would be able globally to tend the control polygon.
基金supported by grants from the Research Grants Council of Hong Kong Special Administrative Region,China(Project No.City U 11213119 and T22-603/15N)The financial support is gratefully acknowledgedfinancial support from the Hong Kong Ph.D.Fellowship Scheme funded by the Research Grants Council of Hong Kong,China。
文摘Spatial interpolation has been frequently encountered in earth sciences and engineering.A reasonable appraisal of subsurface heterogeneity plays a significant role in planning,risk assessment and decision making for geotechnical practice.Geostatistics is commonly used to interpolate spatially varying properties at un-sampled locations from scatter measurements.However,successful application of classic geostatistical models requires prior characterization of spatial auto-correlation structures,which poses a great challenge for unexperienced engineers,particularly when only limited measurements are available.Data-driven machine learning methods,such as radial basis function network(RBFN),require minimal human intervention and provide effective alternatives for spatial interpolation of non-stationary and non-Gaussian data,particularly when measurements are sparse.Conventional RBFN,however,is direction independent(i.e.isotropic)and cannot quantify prediction uncertainty in spatial interpolation.In this study,an ensemble RBFN method is proposed that not only allows geotechnical anisotropy to be properly incorporated,but also quantifies uncertainty in spatial interpolation.The proposed method is illustrated using numerical examples of cone penetration test(CPT)data,which involve interpolation of a 2D CPT cross-section from limited continuous 1D CPT soundings in the vertical direction.In addition,a comparative study is performed to benchmark the proposed ensemble RBFN with two other non-parametric data-driven approaches,namely,Multiple Point Statistics(MPS)and Bayesian Compressive Sensing(BCS).The results reveal that the proposed ensemble RBFN provides a better estimation of spatial patterns and associated prediction uncertainty at un-sampled locations when a reasonable amount of data is available as input.Moreover,the prediction accuracy of all the three methods improves as the number of measurements increases,and vice versa.It is also found that BCS prediction is less sensitive to the number of measurement data and outperforms RBFN and MPS when only limited point observations are available.
基金Project(2010CB732103)supported by the National Basic Research Program of ChinaProject(51179092)supported by the National Natural Science Foundation of ChinaProject(2012-KY-02)supported by the State Key Laboratory of Hydroscience and Engineering,China
文摘Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.
基金supported by the Key Program of the National Natural Science Foundation of China (Grand No. 51138001)the China-German Cooperation Project (Grand No. GZ566)+1 种基金the Innovative Research Groups Funded by the National Natural Science Foundation of China (Grand No. 51121005)the Special Funds for the Basic Scientific Research Expenses for the Central University (Grant No. DUT13LK16)
文摘The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method (FEM) as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline (TPS) radial basis functions.