With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the...With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the influence on elastic rock properties. We start with a discussion and an analysis about carbonate rock pore structure utilizing rock slices. Then, given appropriate assumptions, we introduce a new approach to modeling carbonate rocks and construct a pore structure algorithm to identify pore structure mutation with a basis on the Gassmann equation and the Eshelby-Walsh ellipsoid inclusion crack theory. Finally, we compute a single well＇s porosity using this new approach with full wave log data and make a comparison with the predicted result of traditional method and simultaneously invert for reservoir parameters. The study results reveal that the rock pore structure can significantly influence the rocks＇ elastic properties and the predicted porosity error of the new modeling approach is merely 0.74%. Therefore, the approach we introduce can effectively decrease the predicted error of reservoir parameters.展开更多
To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat s...To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.展开更多
Many applications require the solution of large nonsymmetric linear systems with multiple right hand sides. Instead of applying an iterative method to each of these systems individually, it is often more efficient to...Many applications require the solution of large nonsymmetric linear systems with multiple right hand sides. Instead of applying an iterative method to each of these systems individually, it is often more efficient to use a block version of the method that generates iterates for all the systems simultaneously. In this paper, we propose a block version of generalized minimum backward (GMBACK) for solving large multiple nonsymmetric linear systems. The new method employs the block Arnoldi process to construct a basis for the Krylov subspace K m(A, R 0) and seeks X m∈X 0+K m(A, R 0) to minimize the norm of the perturbation to the data given in A.展开更多
Ebulliometric method based on the error analysis equation is presented for systems with the large phase equilibrium constant. Application is given for the determination and calculation of binary vapor-liquid equilibri...Ebulliometric method based on the error analysis equation is presented for systems with the large phase equilibrium constant. Application is given for the determination and calculation of binary vapor-liquid equilibrium data for the ethanol+n-hexane system. It is also given the comparison results between the ebulliometric method based on the error analysis equation and the quasi-static method.展开更多
With the increasing researches on geotechnical properties of the diesel contaminated soil( DCS),the water content measured is indispensable part during the early period. In this study,the relative error of water conte...With the increasing researches on geotechnical properties of the diesel contaminated soil( DCS),the water content measured is indispensable part during the early period. In this study,the relative error of water content measurement using the traditional method is as high as 20. 78%,which is no longer suitable for contaminated soil. Through a series of tests to measure the loss coefficient of diesel in the drying time,the authors finally proposed a modified calculation formula for test samples. The results show that the maximum relative error calculated by using the modified formula is 0. 96%,far lower than that of traditional formula,which can provide accurate data for further study of diesel contaminated soil.展开更多
Feature initialization is an important issue in the monocular simultaneous locahzation ana mapping （SLAM） literature as the feature depth can not be obtained at one observation. In this paper, we present a new featu...Feature initialization is an important issue in the monocular simultaneous locahzation ana mapping （SLAM） literature as the feature depth can not be obtained at one observation. In this paper, we present a new feature initialization method named modified homogeneous parameterization （MHP）, which allows undelayed initialization with scale invariant representation of point features located at various depths. The linearization error of the measurement equation is quantified using a depth estimation model and the feature initialization process is described. In order to verify the performance of the proposed method, the simulation is carried out. Results show that with the proposed method, the SLAM algorithm can achieve better consistency as compared with the existing inverse depth parameterization （IDP） method.展开更多
The artificial boundary method is one of the most popular and effective numerical methods tor solving partial differential equations on unbounded domains, with more than thirty years development. The artificiM boundar...The artificial boundary method is one of the most popular and effective numerical methods tor solving partial differential equations on unbounded domains, with more than thirty years development. The artificiM boundary method has reached maturity in recent years. It has been applied to various problems in scientific and engineering computations, and the theoretical issues such as the convergence and error estimates of the artificial boundary method have been solved gradually. This paper reviews the development and discusses different forms of the artificial boundary method.展开更多
Due to the difficulty in obtaining the a priori estimate,it is very hard to establish the optimal point-wise error bound of a finite difference scheme for solving a nonlinear partial differential equation in high dime...Due to the difficulty in obtaining the a priori estimate,it is very hard to establish the optimal point-wise error bound of a finite difference scheme for solving a nonlinear partial differential equation in high dimensions(2D or 3D).We here propose and analyze finite difference methods for solving the coupled GrossPitaevskii equations in two dimensions,which models the two-component Bose-Einstein condensates with an internal atomic Josephson junction.The methods which we considered include two conservative type schemes and two non-conservative type schemes.Discrete conservation laws and solvability of the schemes are analyzed.For the four proposed finite difference methods,we establish the optimal convergence rates for the error at the order of O(h^2+τ~2)in the l~∞-norm(i.e.,the point-wise error estimates)with the time stepτand the mesh size h.Besides the standard techniques of the energy method,the key techniques in the analysis is to use the cut-off function technique,transformation between the time and space direction and the method of order reduction.All the methods and results here are also valid and can be easily extended to the three-dimensional case.Finally,numerical results are reported to confirm our theoretical error estimates for the numerical methods.展开更多
We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation （NLDE） in the nonrelativistic limit regime, involving a small dimensionless parameter 0...We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation （NLDE） in the nonrelativistic limit regime, involving a small dimensionless parameter 0 〈 ε〈〈1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O（ ε^2） and O（1） in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference （CNFD） method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ- as well as the small parameter 0 〈 ε≤1 Based on the error bound, in order to obtain ＇correct＇ numerical solutions in the nonrelativistic limit regime, i.e., 0 〈 ε≤1 , the CNFD method requests the ε-scalability： τ- = O（ε3） and h = O（√ε）. Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time- splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O（ε2） and h = O（1） when 0 〈 ε 〈〈 1. Extensive numerical results are reported to confirm our error estimates.展开更多
A Crank-Nicolson scheme based on nonconforming finite element with moving grids is dis- cussed for a class of parabolic integro-differential equations under anisotropic meshes. The corresponding convergence analysis i...A Crank-Nicolson scheme based on nonconforming finite element with moving grids is dis- cussed for a class of parabolic integro-differential equations under anisotropic meshes. The corresponding convergence analysis is presented and the error estimates are obtained by using the interpolation operator instead of the conventional elliptic projection which is an indispensable tool in the convergence analysis of traditional finite element methods in previous literature.展开更多
A nonconforming finite element method for the nonlinear parabolic equations is studied inthis paper.The convergence analysis is presented and the optimal error estimate in L^2(‖·‖_h)norm isobtained through Ritz...A nonconforming finite element method for the nonlinear parabolic equations is studied inthis paper.The convergence analysis is presented and the optimal error estimate in L^2(‖·‖_h)norm isobtained through Ritz projection technique,where ‖·‖_h is a norm over the finite element space.展开更多
基金sponsored by the National Nature Science Foundation of China (Grant No.40904034 and 40839905)
文摘With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the influence on elastic rock properties. We start with a discussion and an analysis about carbonate rock pore structure utilizing rock slices. Then, given appropriate assumptions, we introduce a new approach to modeling carbonate rocks and construct a pore structure algorithm to identify pore structure mutation with a basis on the Gassmann equation and the Eshelby-Walsh ellipsoid inclusion crack theory. Finally, we compute a single well＇s porosity using this new approach with full wave log data and make a comparison with the predicted result of traditional method and simultaneously invert for reservoir parameters. The study results reveal that the rock pore structure can significantly influence the rocks＇ elastic properties and the predicted porosity error of the new modeling approach is merely 0.74%. Therefore, the approach we introduce can effectively decrease the predicted error of reservoir parameters.
基金Supported by China Postdoctoral Science Foundation under Grant No.20090460102
文摘To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.
文摘Many applications require the solution of large nonsymmetric linear systems with multiple right hand sides. Instead of applying an iterative method to each of these systems individually, it is often more efficient to use a block version of the method that generates iterates for all the systems simultaneously. In this paper, we propose a block version of generalized minimum backward (GMBACK) for solving large multiple nonsymmetric linear systems. The new method employs the block Arnoldi process to construct a basis for the Krylov subspace K m(A, R 0) and seeks X m∈X 0+K m(A, R 0) to minimize the norm of the perturbation to the data given in A.
基金Supported by the National Natural Science Foundation of China(No. 29976035) and the Natural Science Foundation of Zhejiang Province.
文摘Ebulliometric method based on the error analysis equation is presented for systems with the large phase equilibrium constant. Application is given for the determination and calculation of binary vapor-liquid equilibrium data for the ethanol+n-hexane system. It is also given the comparison results between the ebulliometric method based on the error analysis equation and the quasi-static method.
文摘With the increasing researches on geotechnical properties of the diesel contaminated soil( DCS),the water content measured is indispensable part during the early period. In this study,the relative error of water content measurement using the traditional method is as high as 20. 78%,which is no longer suitable for contaminated soil. Through a series of tests to measure the loss coefficient of diesel in the drying time,the authors finally proposed a modified calculation formula for test samples. The results show that the maximum relative error calculated by using the modified formula is 0. 96%,far lower than that of traditional formula,which can provide accurate data for further study of diesel contaminated soil.
文摘Feature initialization is an important issue in the monocular simultaneous locahzation ana mapping （SLAM） literature as the feature depth can not be obtained at one observation. In this paper, we present a new feature initialization method named modified homogeneous parameterization （MHP）, which allows undelayed initialization with scale invariant representation of point features located at various depths. The linearization error of the measurement equation is quantified using a depth estimation model and the feature initialization process is described. In order to verify the performance of the proposed method, the simulation is carried out. Results show that with the proposed method, the SLAM algorithm can achieve better consistency as compared with the existing inverse depth parameterization （IDP） method.
基金supported by National National Science Foundation of China(Grant No.10971116)FRG of Hong Kong Baptist University(Grant No.FRG1/11-12/051)
文摘The artificial boundary method is one of the most popular and effective numerical methods tor solving partial differential equations on unbounded domains, with more than thirty years development. The artificiM boundary method has reached maturity in recent years. It has been applied to various problems in scientific and engineering computations, and the theoretical issues such as the convergence and error estimates of the artificial boundary method have been solved gradually. This paper reviews the development and discusses different forms of the artificial boundary method.
基金supported by National Natural Science Foundation of China(Grant No.11201239)the Singapore A*STAR SERC PSF(Grant No.1321202067)
文摘Due to the difficulty in obtaining the a priori estimate,it is very hard to establish the optimal point-wise error bound of a finite difference scheme for solving a nonlinear partial differential equation in high dimensions(2D or 3D).We here propose and analyze finite difference methods for solving the coupled GrossPitaevskii equations in two dimensions,which models the two-component Bose-Einstein condensates with an internal atomic Josephson junction.The methods which we considered include two conservative type schemes and two non-conservative type schemes.Discrete conservation laws and solvability of the schemes are analyzed.For the four proposed finite difference methods,we establish the optimal convergence rates for the error at the order of O(h^2+τ~2)in the l~∞-norm(i.e.,the point-wise error estimates)with the time stepτand the mesh size h.Besides the standard techniques of the energy method,the key techniques in the analysis is to use the cut-off function technique,transformation between the time and space direction and the method of order reduction.All the methods and results here are also valid and can be easily extended to the three-dimensional case.Finally,numerical results are reported to confirm our theoretical error estimates for the numerical methods.
基金supported by the Ministry of Education of Singapore(Grant No.R146-000-196-112)National Natural Science Foundation of China(Grant No.91430103)
文摘We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation （NLDE） in the nonrelativistic limit regime, involving a small dimensionless parameter 0 〈 ε〈〈1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O（ ε^2） and O（1） in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference （CNFD） method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ- as well as the small parameter 0 〈 ε≤1 Based on the error bound, in order to obtain ＇correct＇ numerical solutions in the nonrelativistic limit regime, i.e., 0 〈 ε≤1 , the CNFD method requests the ε-scalability： τ- = O（ε3） and h = O（√ε）. Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time- splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O（ε2） and h = O（1） when 0 〈 ε 〈〈 1. Extensive numerical results are reported to confirm our error estimates.
基金This research is supported by the National Natural Science Foundation of China under Grant Nos. 10671184 and 10971203.
文摘A Crank-Nicolson scheme based on nonconforming finite element with moving grids is dis- cussed for a class of parabolic integro-differential equations under anisotropic meshes. The corresponding convergence analysis is presented and the error estimates are obtained by using the interpolation operator instead of the conventional elliptic projection which is an indispensable tool in the convergence analysis of traditional finite element methods in previous literature.
基金supported by the Natural Science Foundation of China under Grant Nos.10671184 and 10971203
文摘A nonconforming finite element method for the nonlinear parabolic equations is studied inthis paper.The convergence analysis is presented and the optimal error estimate in L^2(‖·‖_h)norm isobtained through Ritz projection technique,where ‖·‖_h is a norm over the finite element space.
