This paper considers the multi-symplectic formulations of the generalized fifth-order KdV equation in Hamiltonian space. Recurring to the midpoint rule, it presents an implicit multi-symplectic scheme with discrete mu...This paper considers the multi-symplectic formulations of the generalized fifth-order KdV equation in Hamiltonian space. Recurring to the midpoint rule, it presents an implicit multi-symplectic scheme with discrete multi-symplectic conservation law to solve the partial differential equations which are derived from the generalized fifth-order KdV equation numerically. The results of the numerical experiments show that this multi-symplectic algorithm is good in accuracy and its long-time numerical behaviour is also perfect.展开更多
The generalized Boussinesq equation that represents a group of important nonlinear equations possesses many interesting properties. Multi-symplectic formulations of the generalized Boussinesq equation in the Hamilton ...The generalized Boussinesq equation that represents a group of important nonlinear equations possesses many interesting properties. Multi-symplectic formulations of the generalized Boussinesq equation in the Hamilton space are introduced in this paper. And then an implicit multi-symplectic scheme equivalent to the multi-symplectic Box scheme is constructed to solve the partial differential equations （PDEs） derived from the generalized Boussinesq equation. Finally, the numerical experiments on the soliton solutions of the generalized Boussinesq equation are reported. The results show that the multi-symplectic method is an efficient algorithm with excellent long-time numerical behaviors for nonlinear partial differential equations.展开更多
In the present paper, a general solution involv- ing three arbitrary functions for the generalized （2＋1）- dimensional KdV-mKdV equation, which is derived from the generalized （1＋1）-dimensional KdV-mKdV equa- tio...In the present paper, a general solution involv- ing three arbitrary functions for the generalized （2＋1）- dimensional KdV-mKdV equation, which is derived from the generalized （1＋1）-dimensional KdV-mKdV equa- tion, is first introduced by means of the Wiess, Tabor, Carnevale （WTC） truncation method. And then multi- symplectic formulations with several conservation laws taken into account are presented for the generalized （2＋1）- dimensional KdV-mKdV equation based on the multi- symplectic theory of Bridges. Subsequently, in order to simulate the periodic wave solutions in terms of rational functions of the Jacobi elliptic functions derived from thegeneral solution, a semi-implicit multi-symplectic scheme is constructed that is equivalent 1：o the Preissmann scheme. From the results of the numerical experiments, we can con- clude that the multi-symplectic schemes can accurately sim- ulate the periodic wave solutions of the generalized （2＋1）- dimensional KdV-mKdV equation while preserve approxi- mately the conservation laws.展开更多
The parallel multisection method for solving algebraic eigenproblem has been presented in recent years with the development of the parallel computers, but all the research work is limited in standard eigenproblems of ...The parallel multisection method for solving algebraic eigenproblem has been presented in recent years with the development of the parallel computers, but all the research work is limited in standard eigenproblems of symmetric tridiagonal matrix. The multisection method for solving the generalized eigenproblem applied significantly in many science and engineering domains has not been studied. The parallel region preserving multisection method (PRM for short) for solving generalized eigenproblems of large sparse and real symmetric matrix is presented in this paper. This method not only retains the advantages of the conventional determinant search method (DS for short), but also overcomes its disadvantages such as leaking roots and disconvergence. We have tested the method on the YH 1 vector computer, and compared it with the parallel region preserving determinant search method the parallel region preserving bisection method (PRB for short). The numerical results show that PRM has a higher speed up, for instance, it attains the speed up of 7.7 when the scale of the problem is 2 114 and the eigenpair found is 3, and PRM is superior to PRB when the scale of the problem is large.展开更多
Applying the generalized method, which is a direct and unified algebraic method for constructing multipletravelling wave solutions of nonlinear partial differential equations (PDEs), and implementing in a computer alg...Applying the generalized method, which is a direct and unified algebraic method for constructing multipletravelling wave solutions of nonlinear partial differential equations (PDEs), and implementing in a computer algebraicsystem, we consider the generalized Zakharov-Kuzentsov equation with nonlinear terms of any order. As a result, wecan not only successfully recover the previously known travelling wave solutions found by existing various tanh methodsand other sophisticated methods, but also obtain some new formal solutions. The solutions obtained include kink-shapedsolitons, bell-shaped solitons, singular solitons, and periodic solutions.展开更多
We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler ...We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler method for temporal discretization. It is verified that the proposed method has corresponding discrete multi-symplectic conservation laws. Numerical simulations indicate that the proposed scheme is characterized by excellent conservation.展开更多
In this paper, we consider the convergence of the generalized alternating direction method of multipliers（GADMM） for solving linearly constrained nonconvex minimization model whose objective contains coupled functio...In this paper, we consider the convergence of the generalized alternating direction method of multipliers（GADMM） for solving linearly constrained nonconvex minimization model whose objective contains coupled functions. Under the assumption that the augmented Lagrangian function satisfies the Kurdyka-Lojasiewicz inequality, we prove that the sequence generated by the GADMM converges to a critical point of the augmented Lagrangian function when the penalty parameter in the augmented Lagrangian function is sufficiently large. Moreover, we also present some sufficient conditions guaranteeing the sublinear and linear rate of convergence of the algorithm.展开更多
The modified AOR method for solving linear complementarity problem(LCP(M,p))was proposed in literature,with some convergence results.In this paper,we considered the MAOR method for generalized-order linear complementa...The modified AOR method for solving linear complementarity problem(LCP(M,p))was proposed in literature,with some convergence results.In this paper,we considered the MAOR method for generalized-order linear complementarity problem(ELCP(M,N,p,q)),where M,N are nonsingular matrices of the following form:M=[D11H1K1D2],N=[D12H2K2D22],D11,D12,D21 and D22 are square nonsingular diagonal matrices.展开更多
In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseu...In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseudo-contractive mappings, and the set of solutions of the variational inclusion problem with multi-valued maximal monotone mappings and inverse-strongly monotone mappings in Hilbert space. Under suitable conditions, some strong convergence theorems are proved. Our results extends the recent results in G.L.Acedo and H.K.Xu [2], Zhang, Lee and Chan [8], Wakahashi and Toyoda [9], Takahashi and Takahashi [I0] and S. S. Chang, H. W. Joseph Lee and C. K. Chan [II], S.Takahashi and W.Takahashi [12]. Moreover, the method of proof adopted in this article is different from those of [4] and [12].展开更多
Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that fo...Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that for finite element method.展开更多
In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-...In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws-a multi-symplectic conservation law （CLS）, a local energy conservation law （ECL） as well as a local momentum conservation law （MCL） --is constructed to discrete the PDEs that are derived from the membrane free vibration equation. The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior,展开更多
As the main unconventional natural gas reservoirs,shale gas reservoirs and coalbed methane(CBM)reservoirs belong to adsorptive gas reservoirs,i.e.,gas reservoirs containing adsorbed gas.Shale gas and CBM reservoirs us...As the main unconventional natural gas reservoirs,shale gas reservoirs and coalbed methane(CBM)reservoirs belong to adsorptive gas reservoirs,i.e.,gas reservoirs containing adsorbed gas.Shale gas and CBM reservoirs usually have the characteristics of rich adsorbed gas and obvious dynamic changes of porosity and permeability.A generalized material balance equation and the corresponding reserve evaluation method considering all the mechanisms for both shale gas reservoirs and CBM reservoirs are necessary.In this work,a generalized material balance equation(GMBE)considering the effects of critical desorption pressure,stress sensitivity,matrix shrinkage,water production,water influx,and solubility of natural gas in water is established.Then,by converting the GMBE to a linear relationship between two parameter groups related with known formation/fluid properties and dynamic performance data,the straight-line reserve evaluation method is proposed.By using the slope and the y-intercept of this straight line,the original adsorbed gas in place(OAGIP),original free gas in place(OFGIP),original dissolved gas in place(ODGIP),and the original gas in place(OGIP)can be quickly calculated.Third,two validation cases for shale gas reservoir and CBM reservoir are conducted using commercial reservoir simulator and the coalbed methane dynamic performance analysis software,respectively.Finally,two field studies in the Fuling shale gas field and the Baode CBM field are presented.Results show that the GMBE and the corresponding straight-line reserve evaluation method are rational,accurate,and effective for both shale gas reservoirs and CBM reservoirs.More detailed information about reserves of shale gas and CBM reservoirs can be clarified,and only the straight-line fitting approach is used to determine all kinds of reserves without iteration,proving that the proposed method has great advantages compared with other current methods.展开更多
In this paper, the generalized extended tanh-function method is used for constructing the traveling wave solutions of nonlinear evolution equations. We choose Fisher＇s equation, the nonlinear schr＆#168;odinger equat...In this paper, the generalized extended tanh-function method is used for constructing the traveling wave solutions of nonlinear evolution equations. We choose Fisher＇s equation, the nonlinear schr＆#168;odinger equation to illustrate the validity and ad-vantages of the method. Many new and more general traveling wave solutions are obtained. Furthermore, this method can also be applied to other nonlinear equations in physics.展开更多
In this paper, the trial function method is extended to study the generalized nonlinear Schrodinger equation with time- dependent coefficients. On the basis of a generalized traveling wave transformation and a trial f...In this paper, the trial function method is extended to study the generalized nonlinear Schrodinger equation with time- dependent coefficients. On the basis of a generalized traveling wave transformation and a trial function, we investigate the exact envelope traveling wave solutions of the generalized nonlinear Schrodinger equation with time-dependent coefficients. Taking advantage of solutions to trial function, we successfully obtain exact solutions for the generalized nonlinear Schrodinger equation with time-dependent coefficients under constraint conditions.展开更多
Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic ...Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations （PDEs） derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.展开更多
A generalized finite spectral method is proposed. The method is of highorder accuracy. To attain high accuracy in time discretization, the fourth-order AdamsBashforth-Moulton predictor and corrector scheme was used. T...A generalized finite spectral method is proposed. The method is of highorder accuracy. To attain high accuracy in time discretization, the fourth-order AdamsBashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Hermite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation （nonlinear convection- diffusion problem） and KdV equation（single solitary and 2-solitary wave problems）, where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.展开更多
Dynamic response of multispan viscoelastic thin beams subjected to a moving mass is studied by an efficient numerical method in some detail. To this end, the unknown parameters of the problem are discretized in spatia...Dynamic response of multispan viscoelastic thin beams subjected to a moving mass is studied by an efficient numerical method in some detail. To this end, the unknown parameters of the problem are discretized in spatial domain using generalized moving least square method （GMLSM） and then, discrete equations of motion based on Lagrange＇s equation are obtained. Maximum deflection and bending moments are considered as the important design parameters. The design parameter spectra in terms of mass weight and velocity of the moving mass are presented for multispan viscoelastic beams as well as various values of relaxation rate and beam span number. A reasonable good agreement is achieved between the results of the proposed solution and those obtained by other researchers. The results indicate that, although the load inertia effects in beams with higher span number would be intensified for higher levels of moving mass velocity, the maximum values of design parameters would increase either. Moreover, the possibility of mass separation is shown to be more critical as the span number of the beam increases. This fact also violates the linear relation between the mass weight of the moving load and the associated design parameters, especially for high moving mass velocities. However, as the relaxation rate of the beam material increases, the load inertia effects as well as the possibility of moving mass separation reduces.展开更多
Semi-inverse method, which is an integration and an extension of Hu's try-and-error method, Chien's veighted residual method and Liu's systematic method, is proposed to establish generalized variational pr...Semi-inverse method, which is an integration and an extension of Hu's try-and-error method, Chien's veighted residual method and Liu's systematic method, is proposed to establish generalized variational principles with multi-variables without arty variational crisis phenomenon. The method is to construct an energy trial-functional with an unknown function F, which can be readily identified by making the trial-functional stationary and using known constraint equations. As a result generalized variational principles with two kinds of independent variables (such as well-known Hellinger-Reissner variational principle and Hu-Washizu principle) and generalized variational principles with three kinds of independent variables (such as Chien's generalized variational principles) in elasticity have been deduced without using Lagrange multiplier method. By semi-inverse method, the author has also proved that Hu-Washizu principle is actually a variational principle with only two kinds of independent variables, stress-strain relations are still its constraints.展开更多
Using the idea of splitting numerical methods and the multi-symplectic methods, we propose a multisymplectic splitting （MSS） method to solve the two-dimensional nonlinear Schrodinger equation （2D-NLSE） in this pap...Using the idea of splitting numerical methods and the multi-symplectic methods, we propose a multisymplectic splitting （MSS） method to solve the two-dimensional nonlinear Schrodinger equation （2D-NLSE） in this paper. It is further shown that the method constructed in this way preserve the global symplectieity exactly. Numerical experiments for the plane wave solution and singular solution of the 2D-NLSE show the accuracy and effectiveness of the proposed method.展开更多
The generalized differential quadrature method （GDQM） is employed to con- sider the free vibration and critical speed of moderately thick rotating laminated compos- ite conical shells with different boundary conditi...The generalized differential quadrature method （GDQM） is employed to con- sider the free vibration and critical speed of moderately thick rotating laminated compos- ite conical shells with different boundary conditions developed from the first-order shear deformation theory （FSDT）. The equations of motion are obtained applying Hamilton＇s concept, which contain the influence of the centrifugal force, the Coriolis acceleration, and the preliminary hoop stress. In addition, the axial load is applied to the conical shell as a ratio of the global critical buckling load. The governing partial differential equations are given in the expressions of five components of displacement related to the points ly- ing on the reference surface of the shell. Afterward, the governing differential equations are converted into a group of algebraic equations by using the GDQM. The outcomes are achieved considering the effects of stacking sequences, thickness of the shell, rotating velocities, half-vertex cone angle, and boundary conditions. Furthermore, the outcomes indicate that the rate of the convergence of frequencies is swift, and the numerical tech- nique is superior stable. Three comparisons between the selected outcomes and those of other research are accomplished, and excellent agreement is achieved.展开更多
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10572119, 10772147 and 10632030)the Doctoral Program Foundation of Education Ministry of China (Grant No 20070699028)+1 种基金the National Natural Science Foundation of Shaanxi Province of China (Grant No 2006A07)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment
文摘This paper considers the multi-symplectic formulations of the generalized fifth-order KdV equation in Hamiltonian space. Recurring to the midpoint rule, it presents an implicit multi-symplectic scheme with discrete multi-symplectic conservation law to solve the partial differential equations which are derived from the generalized fifth-order KdV equation numerically. The results of the numerical experiments show that this multi-symplectic algorithm is good in accuracy and its long-time numerical behaviour is also perfect.
基金Project supported by the National Natural Science Foundation of China (Nos.10572119,10772147,10632030)the Doctoral Program Foundation of Education Ministry of China (No.20070699028)+1 种基金the Natural Science Foundation of Shaanxi Province of China (No.2006A07)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment,Dalian University of Technology.
文摘The generalized Boussinesq equation that represents a group of important nonlinear equations possesses many interesting properties. Multi-symplectic formulations of the generalized Boussinesq equation in the Hamilton space are introduced in this paper. And then an implicit multi-symplectic scheme equivalent to the multi-symplectic Box scheme is constructed to solve the partial differential equations （PDEs） derived from the generalized Boussinesq equation. Finally, the numerical experiments on the soliton solutions of the generalized Boussinesq equation are reported. The results show that the multi-symplectic method is an efficient algorithm with excellent long-time numerical behaviors for nonlinear partial differential equations.
基金supported by the National Natural Science Foundation of China (11002115,10972182,11172239)the Science Foundation of Aviation of China (2010ZB53021)+5 种基金the China Postdoctoral Science Special Foundation (201003682)111 project(B07050) to the Northwestern Polytechnical Universitythe NPU Foundation for Fundamental Research (JC200938,JC20110259)the Doctoral Program Foundation of Education Ministry of China(20106102110019)the Open Foundation of State Key Laboratory of Mechanical System & Vibration (MSV-2011-21)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment (GZ0802)
文摘In the present paper, a general solution involv- ing three arbitrary functions for the generalized （2＋1）- dimensional KdV-mKdV equation, which is derived from the generalized （1＋1）-dimensional KdV-mKdV equa- tion, is first introduced by means of the Wiess, Tabor, Carnevale （WTC） truncation method. And then multi- symplectic formulations with several conservation laws taken into account are presented for the generalized （2＋1）- dimensional KdV-mKdV equation based on the multi- symplectic theory of Bridges. Subsequently, in order to simulate the periodic wave solutions in terms of rational functions of the Jacobi elliptic functions derived from thegeneral solution, a semi-implicit multi-symplectic scheme is constructed that is equivalent 1：o the Preissmann scheme. From the results of the numerical experiments, we can con- clude that the multi-symplectic schemes can accurately sim- ulate the periodic wave solutions of the generalized （2＋1）- dimensional KdV-mKdV equation while preserve approxi- mately the conservation laws.
文摘The parallel multisection method for solving algebraic eigenproblem has been presented in recent years with the development of the parallel computers, but all the research work is limited in standard eigenproblems of symmetric tridiagonal matrix. The multisection method for solving the generalized eigenproblem applied significantly in many science and engineering domains has not been studied. The parallel region preserving multisection method (PRM for short) for solving generalized eigenproblems of large sparse and real symmetric matrix is presented in this paper. This method not only retains the advantages of the conventional determinant search method (DS for short), but also overcomes its disadvantages such as leaking roots and disconvergence. We have tested the method on the YH 1 vector computer, and compared it with the parallel region preserving determinant search method the parallel region preserving bisection method (PRB for short). The numerical results show that PRM has a higher speed up, for instance, it attains the speed up of 7.7 when the scale of the problem is 2 114 and the eigenpair found is 3, and PRM is superior to PRB when the scale of the problem is large.
基金The project supported by National Natural Science Foundation of China under Grant No.10072013the National Key Basic Research Development Program under Grant No.G1998030600
文摘Applying the generalized method, which is a direct and unified algebraic method for constructing multipletravelling wave solutions of nonlinear partial differential equations (PDEs), and implementing in a computer algebraicsystem, we consider the generalized Zakharov-Kuzentsov equation with nonlinear terms of any order. As a result, wecan not only successfully recover the previously known travelling wave solutions found by existing various tanh methodsand other sophisticated methods, but also obtain some new formal solutions. The solutions obtained include kink-shapedsolitons, bell-shaped solitons, singular solitons, and periodic solutions.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 10971226, 91130013, and 11001270)the National Basic Research Program of China (Grant No. 2009CB723802)
文摘We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler method for temporal discretization. It is verified that the proposed method has corresponding discrete multi-symplectic conservation laws. Numerical simulations indicate that the proposed scheme is characterized by excellent conservation.
基金Supported by the National Natural Science Foundation of China(Grant Nos.1157117811801455)the Fundamental Research Funds of China West Normal University(Grant No.17E084)
文摘In this paper, we consider the convergence of the generalized alternating direction method of multipliers（GADMM） for solving linearly constrained nonconvex minimization model whose objective contains coupled functions. Under the assumption that the augmented Lagrangian function satisfies the Kurdyka-Lojasiewicz inequality, we prove that the sequence generated by the GADMM converges to a critical point of the augmented Lagrangian function when the penalty parameter in the augmented Lagrangian function is sufficiently large. Moreover, we also present some sufficient conditions guaranteeing the sublinear and linear rate of convergence of the algorithm.
文摘The modified AOR method for solving linear complementarity problem(LCP(M,p))was proposed in literature,with some convergence results.In this paper,we considered the MAOR method for generalized-order linear complementarity problem(ELCP(M,N,p,q)),where M,N are nonsingular matrices of the following form:M=[D11H1K1D2],N=[D12H2K2D22],D11,D12,D21 and D22 are square nonsingular diagonal matrices.
基金supported by Scientific Research Fund of Sichuan Provincial Education Department (09ZB102)Scientific Research Fund of Science and Technology Deportment of Sichuan Provincial (2011JYZ011)
文摘In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseudo-contractive mappings, and the set of solutions of the variational inclusion problem with multi-valued maximal monotone mappings and inverse-strongly monotone mappings in Hilbert space. Under suitable conditions, some strong convergence theorems are proved. Our results extends the recent results in G.L.Acedo and H.K.Xu [2], Zhang, Lee and Chan [8], Wakahashi and Toyoda [9], Takahashi and Takahashi [I0] and S. S. Chang, H. W. Joseph Lee and C. K. Chan [II], S.Takahashi and W.Takahashi [12]. Moreover, the method of proof adopted in this article is different from those of [4] and [12].
基金This work is supported by the Foundatiorl of Zhongshan University Advanced Research Centre
文摘Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that for finite element method.
基金the National Natural Science Foundation of China(Nos.10632030 and 10572119)Program for New Century Excellent Talents of Ministry of Education of China(No.NCET-04-0958)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment
文摘In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws-a multi-symplectic conservation law （CLS）, a local energy conservation law （ECL） as well as a local momentum conservation law （MCL） --is constructed to discrete the PDEs that are derived from the membrane free vibration equation. The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior,
基金supported by Science and Technology Major Project of Shanxi Province,China(No.20201101002)Science and Technology Major Project of China,China(No.2016ZX05043002)+1 种基金National Natural Science Foundation Project of China,China(No.51874319)Science Foundation of China University of Petroleum(Beijing),China(No.2462020QNXZ003)to support part of this work
文摘As the main unconventional natural gas reservoirs,shale gas reservoirs and coalbed methane(CBM)reservoirs belong to adsorptive gas reservoirs,i.e.,gas reservoirs containing adsorbed gas.Shale gas and CBM reservoirs usually have the characteristics of rich adsorbed gas and obvious dynamic changes of porosity and permeability.A generalized material balance equation and the corresponding reserve evaluation method considering all the mechanisms for both shale gas reservoirs and CBM reservoirs are necessary.In this work,a generalized material balance equation(GMBE)considering the effects of critical desorption pressure,stress sensitivity,matrix shrinkage,water production,water influx,and solubility of natural gas in water is established.Then,by converting the GMBE to a linear relationship between two parameter groups related with known formation/fluid properties and dynamic performance data,the straight-line reserve evaluation method is proposed.By using the slope and the y-intercept of this straight line,the original adsorbed gas in place(OAGIP),original free gas in place(OFGIP),original dissolved gas in place(ODGIP),and the original gas in place(OGIP)can be quickly calculated.Third,two validation cases for shale gas reservoir and CBM reservoir are conducted using commercial reservoir simulator and the coalbed methane dynamic performance analysis software,respectively.Finally,two field studies in the Fuling shale gas field and the Baode CBM field are presented.Results show that the GMBE and the corresponding straight-line reserve evaluation method are rational,accurate,and effective for both shale gas reservoirs and CBM reservoirs.More detailed information about reserves of shale gas and CBM reservoirs can be clarified,and only the straight-line fitting approach is used to determine all kinds of reserves without iteration,proving that the proposed method has great advantages compared with other current methods.
基金The NSF(11001042) of ChinaSRFDP(20100043120001)FRFCU(09QNJJ002)
文摘In this paper, the generalized extended tanh-function method is used for constructing the traveling wave solutions of nonlinear evolution equations. We choose Fisher＇s equation, the nonlinear schr＆#168;odinger equation to illustrate the validity and ad-vantages of the method. Many new and more general traveling wave solutions are obtained. Furthermore, this method can also be applied to other nonlinear equations in physics.
基金Project supported in part by the National Natural Science Foundation of China(Grant No.11071177)
文摘In this paper, the trial function method is extended to study the generalized nonlinear Schrodinger equation with time- dependent coefficients. On the basis of a generalized traveling wave transformation and a trial function, we investigate the exact envelope traveling wave solutions of the generalized nonlinear Schrodinger equation with time-dependent coefficients. Taking advantage of solutions to trial function, we successfully obtain exact solutions for the generalized nonlinear Schrodinger equation with time-dependent coefficients under constraint conditions.
基金supported by the National Natural Science Foundation of China (Nos. 10772147 and10632030)the Ph. D. Program Foundation of Ministry of Education of China (No. 20070699028)+2 种基金the Natural Science Foundation of Shaanxi Province of China (No. 2006A07)the Open Foundationof State Key Laboratory of Structural Analysis of Industrial Equipment (No. GZ0802)the Foundation for Fundamental Research of Northwestern Polytechnical University
文摘Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations （PDEs） derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.
基金Project supported by the National Natural Science Foundation of China (No.10272118) the Hong Kong Polytechnic University Research Grant (No.A-PE28) the Research Fund for the Doctoral Program of Ministry of Education of China (No.20020558013)
文摘A generalized finite spectral method is proposed. The method is of highorder accuracy. To attain high accuracy in time discretization, the fourth-order AdamsBashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Hermite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation （nonlinear convection- diffusion problem） and KdV equation（single solitary and 2-solitary wave problems）, where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.
文摘Dynamic response of multispan viscoelastic thin beams subjected to a moving mass is studied by an efficient numerical method in some detail. To this end, the unknown parameters of the problem are discretized in spatial domain using generalized moving least square method （GMLSM） and then, discrete equations of motion based on Lagrange＇s equation are obtained. Maximum deflection and bending moments are considered as the important design parameters. The design parameter spectra in terms of mass weight and velocity of the moving mass are presented for multispan viscoelastic beams as well as various values of relaxation rate and beam span number. A reasonable good agreement is achieved between the results of the proposed solution and those obtained by other researchers. The results indicate that, although the load inertia effects in beams with higher span number would be intensified for higher levels of moving mass velocity, the maximum values of design parameters would increase either. Moreover, the possibility of mass separation is shown to be more critical as the span number of the beam increases. This fact also violates the linear relation between the mass weight of the moving load and the associated design parameters, especially for high moving mass velocities. However, as the relaxation rate of the beam material increases, the load inertia effects as well as the possibility of moving mass separation reduces.
文摘Semi-inverse method, which is an integration and an extension of Hu's try-and-error method, Chien's veighted residual method and Liu's systematic method, is proposed to establish generalized variational principles with multi-variables without arty variational crisis phenomenon. The method is to construct an energy trial-functional with an unknown function F, which can be readily identified by making the trial-functional stationary and using known constraint equations. As a result generalized variational principles with two kinds of independent variables (such as well-known Hellinger-Reissner variational principle and Hu-Washizu principle) and generalized variational principles with three kinds of independent variables (such as Chien's generalized variational principles) in elasticity have been deduced without using Lagrange multiplier method. By semi-inverse method, the author has also proved that Hu-Washizu principle is actually a variational principle with only two kinds of independent variables, stress-strain relations are still its constraints.
基金Supported by the Natural Science Foundation of China under Grant No.0971226the 973 Project of China under Grant No.2009CB723802+1 种基金the Research Innovation Fund of Hunan Province under Grant No.CX2011B011the Innovation Fund of NUDT under Grant No.B110205
文摘Using the idea of splitting numerical methods and the multi-symplectic methods, we propose a multisymplectic splitting （MSS） method to solve the two-dimensional nonlinear Schrodinger equation （2D-NLSE） in this paper. It is further shown that the method constructed in this way preserve the global symplectieity exactly. Numerical experiments for the plane wave solution and singular solution of the 2D-NLSE show the accuracy and effectiveness of the proposed method.
文摘The generalized differential quadrature method （GDQM） is employed to con- sider the free vibration and critical speed of moderately thick rotating laminated compos- ite conical shells with different boundary conditions developed from the first-order shear deformation theory （FSDT）. The equations of motion are obtained applying Hamilton＇s concept, which contain the influence of the centrifugal force, the Coriolis acceleration, and the preliminary hoop stress. In addition, the axial load is applied to the conical shell as a ratio of the global critical buckling load. The governing partial differential equations are given in the expressions of five components of displacement related to the points ly- ing on the reference surface of the shell. Afterward, the governing differential equations are converted into a group of algebraic equations by using the GDQM. The outcomes are achieved considering the effects of stacking sequences, thickness of the shell, rotating velocities, half-vertex cone angle, and boundary conditions. Furthermore, the outcomes indicate that the rate of the convergence of frequencies is swift, and the numerical tech- nique is superior stable. Three comparisons between the selected outcomes and those of other research are accomplished, and excellent agreement is achieved.
