Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a...Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.展开更多
Dynamic exploration for oil and gas requires careful monitoring of reservoir contents for safety and efficiency of oil extraction. This paper proposes a multi-source and multi-azimuth walk-around vertical electromagne...Dynamic exploration for oil and gas requires careful monitoring of reservoir contents for safety and efficiency of oil extraction. This paper proposes a multi-source and multi-azimuth walk-around vertical electromagnetic profiling （MM-VEP） technique for surface-to-borehole electromagnetic surveying. Based on the difference in conductivities between reservoirs with different concentrations of oil and water, MM-VEP can be used to monitor reservoirs as they are injected with water. The MM-VEP response in five azimuth planes is modeled with three-dimensional （3D） integral equation calculations. The progress of waterflooding in four stages for enhanced oil recovery is shown to be indicated by field anomalies MM-VEP caused by variations in the reservoir resistivity. Numerical modeling demonstrates that MM-VEP measurements provides enough quantitative information from an underground reservoir to accurately detect oil deposits and monitor the progress of waterflooding.展开更多
We study the existence of solutions to the second order three-point boundary value problem：{x″（t）＋f（t,x（t）,x′（t））=0,t≠ti,△x（ti）=Ii（x（ti）,x′（ti））,i=1,2,…,k,△x′（ti）=Ji（x（ti）,x′（t））,i=1...We study the existence of solutions to the second order three-point boundary value problem：{x″（t）＋f（t,x（t）,x′（t））=0,t≠ti,△x（ti）=Ii（x（ti）,x′（ti））,i=1,2,…,k,△x′（ti）=Ji（x（ti）,x′（t））,i=1,2,…,k,x（0）=0=x（1）-αx（η）,where 0〈η〈1,α∈R,and f：[0,1]×R×R→R,Ii：R×R→R,Ji：R×R→R（i=1,2,…,k）are continuous. Our results is new and different from previous results. In particular, we obtain the Green function of the problem, which makes the problem simpler.展开更多
The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the ...The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the turbulent coherent structure in the channel flow. Th is method contains three dimensional coupling difference scheme with high accur acy and high resolution, and the high order time splitting methods. Compared with the general spectral method, the method can be used to research turbule nt coherent structure under more general boundary conditions and in flow domains . In this paper, the generation and evolution of the turbulent coherent structur es ind uced by wall pulse in the channel flow are simulated, and the basic characterist ics and rules of the turbulent coherent structure are shown. Computational r esults indicate that a wall negative pulse is more convenient than the resonant three wave model.展开更多
In this paper, we extend the mapping deformation method proposed by Lou. It is used to find new exacttravelling wave solutions of nonlinear partial differential equation or coupled nonlinear partial differential equat...In this paper, we extend the mapping deformation method proposed by Lou. It is used to find new exacttravelling wave solutions of nonlinear partial differential equation or coupled nonlinear partial differential equations(PDEs). Based on the idea of the homogeneous balance method, we construct the general mapping relation betweenthe solutions of the PDEs and those of the cubic nonlinear Klein-Gordon (NKG) equation. By using this relation andthe abundant solutions of the cubic NKG equation, many explicit and exact travelling wave solutions of three systemsof coupled PDEs, which contain solitary wave solutions, trigonometric function solutions, Jacobian elliptic functionsolutions, and rational solutions, are obtained.展开更多
In this paper, the solution of Chebyshev equation with its argument being greater than 1 is obtained. The initial value of the derivative of the solution is the expression of magnetization, which is valid for any spin...In this paper, the solution of Chebyshev equation with its argument being greater than 1 is obtained. The initial value of the derivative of the solution is the expression of magnetization, which is valid for any spin quantum number S. The Chebyshev equation is transformed from an ordinary differential equation obtained when we dealt with Heisenberg model, in order to calculate all three components of magnetization, by many-body Green＇s function under random phase approximation. The Chebyshev functions with argument being greater than 1 are discussed. This paper shows that the Chebyshev polynomials with their argument being greater than 1 have their physical application.展开更多
A new three-component Camassa-Holm equation is introduced. This system is endowed with a structuresimilar to the Camassa-Holm equation. It has peakon solitons and conserves H^1-norm conservation law.
Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the invest...Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.展开更多
In this paper we apply fractional calculus to solve the 3rd order ordinary differential equation of the following form: (z-a)(z-b)(z-c)φ 3+(βz 2+γz+D)φ 2+(α(2β-3α-3)z+αγ+α(α+1)(a+b+c))φ 1+α(α-...In this paper we apply fractional calculus to solve the 3rd order ordinary differential equation of the following form: (z-a)(z-b)(z-c)φ 3+(βz 2+γz+D)φ 2+(α(2β-3α-3)z+αγ+α(α+1)(a+b+c))φ 1+α(α-1)(β-2α-2)φ=f.展开更多
By using cone expansion-compression theorem in this paper, we study boundary value problems for a coupled system of nonlinear third-order differential equation. Some sufficient conditions are obtained which guarantee ...By using cone expansion-compression theorem in this paper, we study boundary value problems for a coupled system of nonlinear third-order differential equation. Some sufficient conditions are obtained which guarantee the boundary value problems for a coupled system of nonlinear third-order differential equation has at least one positive solution. Some examples are given to verify our results.展开更多
The symmetries, symmetry reductions, and exact solutions of a coupled nonlinear Schrodinger （CNLS） equation derived from the governing system for atmospheric gravity waves are researched by means of classical Lie gr...The symmetries, symmetry reductions, and exact solutions of a coupled nonlinear Schrodinger （CNLS） equation derived from the governing system for atmospheric gravity waves are researched by means of classical Lie group approach in this paper. Calculation shows the CNLS equation is invariant under some Galilean transformations, scaling transformations, phase shifts, and space-time translations. Some ordinary differential equations are derived from the CNLS equation. Several exact solutions including envelope cnoidal waves, solitary waves and trigonometric function solutions for the CNLS equation are also obtained by making use of symmetries.展开更多
The objective of this paper is to provide a provable solution of the ancient Greek problem of trisecting an arbitrary angle employing only compass and straightedge （ruler）. （Pierre Laurent Wantzel, 1837） obscurely...The objective of this paper is to provide a provable solution of the ancient Greek problem of trisecting an arbitrary angle employing only compass and straightedge （ruler）. （Pierre Laurent Wantzel, 1837） obscurely presented a proof based on ideas from Galois field showing that, the solution of angle trisection corresponds to solution of the cubic equation; x3 - 3x - 1 = 0, which is geometrically irreducible [1]. The focus of this work is to show the possibility to solve the trisection of an angle by correcting some flawed methods meant for general construction of angles, and exemplify why the stated trisection impossible proof is not geometrically valid. The revealed proof is based on a concept from the Archimedes proposition of straightedge construction [2, 3].展开更多
The dynamic behavior of two collinear cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied using the Schmidt method for the permeable crack surface conditions. By using the Fou...The dynamic behavior of two collinear cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied using the Schmidt method for the permeable crack surface conditions. By using the Fourier transform, the problem can be solved with a set of triple integral equations in which the unknown variable is the jump of the displacements across the crack surfaces. In solving the triple integral equations, the jump of the displacements across the crack surface is expanded in a series of Jacobi polynomials. It can be obtained that the stress field is independent of the electric field and the magnetic flux.展开更多
基金supported by the National Natural Science Foundation of China(No.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.
基金supported by the National Science and Technology Major Project(No.2011ZX05019-007)National Natural Science Foundation of China(No.41604097)+1 种基金China Postdoctoral Science Foundation(No.2016M592611)Project(Nos.002401003503 and 002401003514)from Guilin University of Technology
文摘Dynamic exploration for oil and gas requires careful monitoring of reservoir contents for safety and efficiency of oil extraction. This paper proposes a multi-source and multi-azimuth walk-around vertical electromagnetic profiling （MM-VEP） technique for surface-to-borehole electromagnetic surveying. Based on the difference in conductivities between reservoirs with different concentrations of oil and water, MM-VEP can be used to monitor reservoirs as they are injected with water. The MM-VEP response in five azimuth planes is modeled with three-dimensional （3D） integral equation calculations. The progress of waterflooding in four stages for enhanced oil recovery is shown to be indicated by field anomalies MM-VEP caused by variations in the reservoir resistivity. Numerical modeling demonstrates that MM-VEP measurements provides enough quantitative information from an underground reservoir to accurately detect oil deposits and monitor the progress of waterflooding.
基金Supported by the National Natural Science Foundation of China(10371006)
文摘We study the existence of solutions to the second order three-point boundary value problem：{x″（t）＋f（t,x（t）,x′（t））=0,t≠ti,△x（ti）=Ii（x（ti）,x′（ti））,i=1,2,…,k,△x′（ti）=Ji（x（ti）,x′（t））,i=1,2,…,k,x（0）=0=x（1）-αx（η）,where 0〈η〈1,α∈R,and f：[0,1]×R×R→R,Ii：R×R→R,Ji：R×R→R（i=1,2,…,k）are continuous. Our results is new and different from previous results. In particular, we obtain the Green function of the problem, which makes the problem simpler.
文摘The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the turbulent coherent structure in the channel flow. Th is method contains three dimensional coupling difference scheme with high accur acy and high resolution, and the high order time splitting methods. Compared with the general spectral method, the method can be used to research turbule nt coherent structure under more general boundary conditions and in flow domains . In this paper, the generation and evolution of the turbulent coherent structur es ind uced by wall pulse in the channel flow are simulated, and the basic characterist ics and rules of the turbulent coherent structure are shown. Computational r esults indicate that a wall negative pulse is more convenient than the resonant three wave model.
文摘In this paper, we extend the mapping deformation method proposed by Lou. It is used to find new exacttravelling wave solutions of nonlinear partial differential equation or coupled nonlinear partial differential equations(PDEs). Based on the idea of the homogeneous balance method, we construct the general mapping relation betweenthe solutions of the PDEs and those of the cubic nonlinear Klein-Gordon (NKG) equation. By using this relation andthe abundant solutions of the cubic NKG equation, many explicit and exact travelling wave solutions of three systemsof coupled PDEs, which contain solitary wave solutions, trigonometric function solutions, Jacobian elliptic functionsolutions, and rational solutions, are obtained.
基金The project supported by the State Key Project of Fundamental Research of China under Grant No. G2000067101
文摘In this paper, the solution of Chebyshev equation with its argument being greater than 1 is obtained. The initial value of the derivative of the solution is the expression of magnetization, which is valid for any spin quantum number S. The Chebyshev equation is transformed from an ordinary differential equation obtained when we dealt with Heisenberg model, in order to calculate all three components of magnetization, by many-body Green＇s function under random phase approximation. The Chebyshev functions with argument being greater than 1 are discussed. This paper shows that the Chebyshev polynomials with their argument being greater than 1 have their physical application.
基金Supported by National Natural Science Foundation of China under Grant Nos. 10671156 and 10671153
文摘A new three-component Camassa-Holm equation is introduced. This system is endowed with a structuresimilar to the Camassa-Holm equation. It has peakon solitons and conserves H^1-norm conservation law.
基金Project(2009ZX04004-031-04) supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China
文摘Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.
文摘In this paper we apply fractional calculus to solve the 3rd order ordinary differential equation of the following form: (z-a)(z-b)(z-c)φ 3+(βz 2+γz+D)φ 2+(α(2β-3α-3)z+αγ+α(α+1)(a+b+c))φ 1+α(α-1)(β-2α-2)φ=f.
基金Foundation item： Supported by the National Natural Science Foundation of China（10801001） Supported by the Natural Science Foundation of Anhui Province（1208085MA13）
文摘By using cone expansion-compression theorem in this paper, we study boundary value problems for a coupled system of nonlinear third-order differential equation. Some sufficient conditions are obtained which guarantee the boundary value problems for a coupled system of nonlinear third-order differential equation has at least one positive solution. Some examples are given to verify our results.
基金supported by the Scientific Research Foundation for the Doctors of University of Electronic Science and Technology of China Zhongshan Institutethe National Natural Science Foundation of China under Grant Nos. 10735030 and 90503006
文摘The symmetries, symmetry reductions, and exact solutions of a coupled nonlinear Schrodinger （CNLS） equation derived from the governing system for atmospheric gravity waves are researched by means of classical Lie group approach in this paper. Calculation shows the CNLS equation is invariant under some Galilean transformations, scaling transformations, phase shifts, and space-time translations. Some ordinary differential equations are derived from the CNLS equation. Several exact solutions including envelope cnoidal waves, solitary waves and trigonometric function solutions for the CNLS equation are also obtained by making use of symmetries.
文摘The objective of this paper is to provide a provable solution of the ancient Greek problem of trisecting an arbitrary angle employing only compass and straightedge （ruler）. （Pierre Laurent Wantzel, 1837） obscurely presented a proof based on ideas from Galois field showing that, the solution of angle trisection corresponds to solution of the cubic equation; x3 - 3x - 1 = 0, which is geometrically irreducible [1]. The focus of this work is to show the possibility to solve the trisection of an angle by correcting some flawed methods meant for general construction of angles, and exemplify why the stated trisection impossible proof is not geometrically valid. The revealed proof is based on a concept from the Archimedes proposition of straightedge construction [2, 3].
文摘The dynamic behavior of two collinear cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied using the Schmidt method for the permeable crack surface conditions. By using the Fourier transform, the problem can be solved with a set of triple integral equations in which the unknown variable is the jump of the displacements across the crack surfaces. In solving the triple integral equations, the jump of the displacements across the crack surface is expanded in a series of Jacobi polynomials. It can be obtained that the stress field is independent of the electric field and the magnetic flux.
