A global interpolating meshless shape function based on the generalized moving least-square （GMLS） is formulated by the transformation technique. Both the shape function and its derivatives meet the Kronecker delta ...A global interpolating meshless shape function based on the generalized moving least-square （GMLS） is formulated by the transformation technique. Both the shape function and its derivatives meet the Kronecker delta function property. With the interpolating GMLS （IGMLS） shape function, an improved element-free Galerkin （EFG） method is proposed for the structural dynamic analysis. Compared with the conven- tional EFG method, the obvious advantage of the proposed method is that the essential boundary conditions including both displacements and derivatives can be imposed by the straightforward way. Meanwhile, it can greatly improve the ill-condition feature of the standard GMLS approximation, and provide good accuracy at low cost. The dynamic analyses of the Euler beam and Kirchhoff plate are performed to demonstrate the feasi- bility and effectiveness of the improved method. The comparison between the numerical results of the conventional method and the improved method shows that the proposed method has better stability, higher accuracy, and less time consumption.展开更多
In this paper, an iterative method is constructed to find the least-squares solutions of generalized Sylvester equation , where is real matrices group, and satisfies different linear constraint. By this iterative meth...In this paper, an iterative method is constructed to find the least-squares solutions of generalized Sylvester equation , where is real matrices group, and satisfies different linear constraint. By this iterative method, for any initial matrix group within a special constrained matrix set, a least squares solution group with satisfying different linear constraint can be obtained within finite iteration steps in the absence of round off errors, and the unique least norm least-squares solution can be obtained by choosing a special kind of initial matrix group. In addition, a minimization property of this iterative method is characterized. Finally, numerical experiments are reported to show the efficiency of the proposed method.展开更多
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.展开更多
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.展开更多
The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features becau...The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.展开更多
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.展开更多
The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of t...The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of the same type: In a breakdown situation the Galerkin method is unable to calculate an approximate solution, while the least-squares method, although does not really break down, is unsucessful in reducing the norm of its residual. In this paper we first establish a unified theorem which gives a relationship between breakdowns in the two methods. We further illustrate theoretically and experimentally that if the coefficient matrix of a lienar system is of high defectiveness with the associated eigenvalues less than 1, then the restarted Galerkin and least-squares methods will be in great risks of complete breakdowns. It appears that our findings may help to understand phenomena observed practically and to derive treatments for breakdowns of this type.展开更多
A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercriti...A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.展开更多
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.展开更多
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.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
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.展开更多
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.展开更多
The study of tidal circulation has a long history . The numerical simulation of tidal flow has been developed greatly with the development of computer techniques in the past two decades. The generalized wave equation ...The study of tidal circulation has a long history . The numerical simulation of tidal flow has been developed greatly with the development of computer techniques in the past two decades. The generalized wave equation finite-element method is a relatively new numerical model for studying shallow water flow . This method was used to simulate tidal waves of the Gulf of St. Lawrence in Canada . The very good agreement of the numerical results with the field data indicated that the model is an effective and promising numerical method for solving two-dimensional tidal wave problems .展开更多
基金Project supported by the National Natural Science Foundation of China(No.11176035)
文摘A global interpolating meshless shape function based on the generalized moving least-square （GMLS） is formulated by the transformation technique. Both the shape function and its derivatives meet the Kronecker delta function property. With the interpolating GMLS （IGMLS） shape function, an improved element-free Galerkin （EFG） method is proposed for the structural dynamic analysis. Compared with the conven- tional EFG method, the obvious advantage of the proposed method is that the essential boundary conditions including both displacements and derivatives can be imposed by the straightforward way. Meanwhile, it can greatly improve the ill-condition feature of the standard GMLS approximation, and provide good accuracy at low cost. The dynamic analyses of the Euler beam and Kirchhoff plate are performed to demonstrate the feasi- bility and effectiveness of the improved method. The comparison between the numerical results of the conventional method and the improved method shows that the proposed method has better stability, higher accuracy, and less time consumption.
文摘In this paper, an iterative method is constructed to find the least-squares solutions of generalized Sylvester equation , where is real matrices group, and satisfies different linear constraint. By this iterative method, for any initial matrix group within a special constrained matrix set, a least squares solution group with satisfying different linear constraint can be obtained within finite iteration steps in the absence of round off errors, and the unique least norm least-squares solution can be obtained by choosing a special kind of initial matrix group. In addition, a minimization property of this iterative method is characterized. Finally, numerical experiments are reported to show the efficiency of the proposed method.
文摘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.
基金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.
基金This project is supported by Research Foundation for Doctoral Program of Higher Education, China (No.98033532)
文摘The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.
基金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.
文摘The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of the same type: In a breakdown situation the Galerkin method is unable to calculate an approximate solution, while the least-squares method, although does not really break down, is unsucessful in reducing the norm of its residual. In this paper we first establish a unified theorem which gives a relationship between breakdowns in the two methods. We further illustrate theoretically and experimentally that if the coefficient matrix of a lienar system is of high defectiveness with the associated eigenvalues less than 1, then the restarted Galerkin and least-squares methods will be in great risks of complete breakdowns. It appears that our findings may help to understand phenomena observed practically and to derive treatments for breakdowns of this type.
基金the National Science Council ot Taiwan,China for funding this research(Project no.:NSC 94-2218-E-035-011)
文摘A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.
基金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.
基金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.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
文摘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.
文摘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.
文摘The study of tidal circulation has a long history . The numerical simulation of tidal flow has been developed greatly with the development of computer techniques in the past two decades. The generalized wave equation finite-element method is a relatively new numerical model for studying shallow water flow . This method was used to simulate tidal waves of the Gulf of St. Lawrence in Canada . The very good agreement of the numerical results with the field data indicated that the model is an effective and promising numerical method for solving two-dimensional tidal wave problems .
