Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
The present work proposed a new method for the modeling by the finite element method of the acoustic propagation problems in infinite axisymmetric cylindrical guides lined with locally reacting absorbent materials wit...The present work proposed a new method for the modeling by the finite element method of the acoustic propagation problems in infinite axisymmetric cylindrical guides lined with locally reacting absorbent materials without flow. The method deals with the development of an efficient transparent boundary condition based on DtN operators. The method developed in this study is successfully applied to a straight axisymmetric lined guide by imposing a mode on one of the artificial boundaries of the truncated guide. The results are in good agreement with analytical solutions. Applying the method for a non-uniform axisymmetric lined guide which is a complex case, proved its effectiveness and the results compared to those of PML layers are in very good agreement.展开更多
Supply chain management is an essential part of an organisation's sustainable programme.Understanding the concentration of natural environment,public,and economic influence and feasibility of your suppliers and pu...Supply chain management is an essential part of an organisation's sustainable programme.Understanding the concentration of natural environment,public,and economic influence and feasibility of your suppliers and purchasers is becoming progressively familiar as all industries are moving towards a massive sustainable potential.To handle such sort of developments in supply chain management the involvement of fuzzy settings and their generalisations is playing an important role.Keeping in mind this role,the aim of this study is to analyse the role and involvement of complex q-rung orthopair normal fuzzy(CQRONF)information in supply chain management.The major impact of this theory is to analyse the notion of confidence CQRONF weighted averaging,confidence CQRONF ordered weighted averaging,confidence CQRONF hybrid averaging,confidence CQRONF weighted geometric,confidence CQRONF ordered weighted geometric,confidence CQRONF hybrid geometric operators and try to diagnose various properties and results.Furthermore,with the help of the CRITIC and VIKOR models,we diagnosed the novel theory of the CQRONF-CRITIC-VIKOR model to check the sensitivity analysis of the initiated method.Moreover,in the availability of diagnosed operators,we constructed a multi-attribute decision-making tool for finding a beneficial sustainable supplier to handle complex dilemmas.Finally,the initiated operator's efficiency is proved by comparative analysis.展开更多
We propose a symplectic partitioned Runge-Kutta (SPRK) method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this te...We propose a symplectic partitioned Runge-Kutta (SPRK) method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete (NAD) operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction (LWC) method and the eighth-order staggered-grid (SG)method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.展开更多
A novel nonlocal operator theory based on the variational principle is proposed for the solution of partial differential equations.Common differential operators as well as the variational forms are defined within the ...A novel nonlocal operator theory based on the variational principle is proposed for the solution of partial differential equations.Common differential operators as well as the variational forms are defined within the context of nonlocal operators.The present nonlocal formulation allows the assembling of the tangent stiffness matrix with ease and simplicity,which is necessary for the eigenvalue analysis such as the waveguide problem.The present formulation is applied to solve the differential electromagnetic vector wave equations based on electric fields.The governing equations are converted into nonlocal integral form.An hourglass energy functional is introduced for the elimination of zeroenergy modes.Finally,the proposed method is validated by testing three classical benchmark problems.展开更多
In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is...In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is obtained based on an a priori choice of the regularization parameter. Our analysis is not based on the sequential continuity of the normalized duality mapping.展开更多
The convergence of several Galerkin-Petrov methods, including polynomial collocation and analytic element collocation methods of Toeplitz operators on Dirichlet space, is established. In particular, it is shown that s...The convergence of several Galerkin-Petrov methods, including polynomial collocation and analytic element collocation methods of Toeplitz operators on Dirichlet space, is established. In particular, it is shown that such methods converge if the basis and test function own certain circular symmetry.展开更多
In this paper we develop two multilevel iteration methods for solving linear systems resulting from the Galerkin method and Tikhonov regularization for linear ill-posed problems. The two algorithms and their convergen...In this paper we develop two multilevel iteration methods for solving linear systems resulting from the Galerkin method and Tikhonov regularization for linear ill-posed problems. The two algorithms and their convergence analyses are presented in an abstract framework.展开更多
Unbounded operators can transform arbitrarily small vectors into arbitrarily large vectors—a phenomenon known as instability. Stabilization methods strive to approximate a value of an unbounded operator by applying a...Unbounded operators can transform arbitrarily small vectors into arbitrarily large vectors—a phenomenon known as instability. Stabilization methods strive to approximate a value of an unbounded operator by applying a family of bounded operators to rough approximate data that do not necessarily lie within the domain of unbounded operator. In this paper we shall be concerned with the stable method of computing values of unbounded operators having perturbations and the stability is established for this method.展开更多
In this article, we define the arithmetic operations of generalized trapezoidal picture fuzzy numbers by vertex method which is assembled on a combination of the (α, γ, β)-cut concept and standard interval analysis...In this article, we define the arithmetic operations of generalized trapezoidal picture fuzzy numbers by vertex method which is assembled on a combination of the (α, γ, β)-cut concept and standard interval analysis. Various related properties are explored. Finally, some computations of picture fuzzy functions over generalized picture fuzzy variables are illustrated by using our proposed technique.展开更多
We discuss the incomplete semi-iterative method (ISIM) for an approximate solution of a linear fixed point equations x=Tx+c with a bounded linear operator T acting on a complex Banach space X such that its resolvent h...We discuss the incomplete semi-iterative method (ISIM) for an approximate solution of a linear fixed point equations x=Tx+c with a bounded linear operator T acting on a complex Banach space X such that its resolvent has a pole of order k at the point 1. Sufficient conditions for the convergence of ISIM to a solution of x=Tx+c, where c belongs to the range space of R(I-T) k, are established. We show that the ISIM has an attractive feature that it is usually convergent even when the spectral radius of the operator T is greater than 1 and Ind 1T≥1. Applications in finite Markov chain is considered and illustrative examples are reported, showing the convergence rate of the ISIM is very high.展开更多
<Abstract>We show that the recently proposed invariant eigen-operator (lEO) method can be successfully applied to solving energy levels for SSH Hamiltonian describing Peierls phase transition.The electronic ener...<Abstract>We show that the recently proposed invariant eigen-operator (lEO) method can be successfully applied to solving energy levels for SSH Hamiltonian describing Peierls phase transition.The electronic energy band of compound lattice is also studied by IEO method.展开更多
The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ...The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ability of the discrete forms expressing to the element functions was talked about. In discrete operator difference method, the displacements of the elements can be reproduced exactly in the discrete forms whether the displacements are conforming or not. According to this point, discrete operator difference method is a method with good performance.展开更多
This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integ...This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.展开更多
We introduce multilevel augmentation methods for solving operator equations based on direct sum decompositions of the range space of the operator and the solution space of the operator equation and a matrix splitting ...We introduce multilevel augmentation methods for solving operator equations based on direct sum decompositions of the range space of the operator and the solution space of the operator equation and a matrix splitting scheme. We establish a general setting for the analysis of these methods, showing that the methods yield approximate solutions of the same convergence order as the best approximation from the subspace. These augmentation methods allow us to develop fast, accurate and stable nonconventional numerical algorithms for solving operator equations. In particular, for second kind equations, special splitting techniques are proposed to develop such algorithms. These algorithms are then applied to solve the linear systems resulting from matrix compression schemes using wavelet-like functions for solving Fredholm integral equations of the second kind. For this special case, a complete analysis for computational complexity and convergence order is presented. Numerical examples are included to demonstrate the efficiency and accuracy of the methods. In these examples we use the proposed augmentation method to solve large scale linear systems resulting from the recently developed wavelet Galerkin methods and fast collocation methods applied to integral equations of the secondkind. Our numerical results confirm that this augmentation method is particularly efficient for solving large scale linear systems induced from wavelet compression schemes.展开更多
Normal form theory is a very effective method when we study degenerate bifurcations of nonlinear dynamical systems. In this paper by using adjoint operator method, normal forms of order 3 and 4 for nonlinear dynamical...Normal form theory is a very effective method when we study degenerate bifurcations of nonlinear dynamical systems. In this paper by using adjoint operator method, normal forms of order 3 and 4 for nonlinear dynamical system with nilpotent linear part and Z(2)-asymmetry are computed. According to normal forms obtained, universal unfoldings for some degenerate bifurcation cases of codimension 3 and simple global characterizations, are studied.展开更多
For classical Hamiltonian with general form H = 1/2∑ijMijpipj+1/2∑ijLijqiqj we find a new convenient way to obtain its normal coordinates, namely, let H be quantised and then employ the invariant eigen-operator (...For classical Hamiltonian with general form H = 1/2∑ijMijpipj+1/2∑ijLijqiqj we find a new convenient way to obtain its normal coordinates, namely, let H be quantised and then employ the invariant eigen-operator (IEO) method (Fan et al. 2004 Phys. Lett. A 321 75) to derive them. The general matrix equation, which relies on M and L, for obtaining the normal coordinates of H is derived.展开更多
In this paper, we study some semi-closed 1-set-contractive operators A and investigate the boundary conditions under which the topological degrees of 1-set contractive fields, deg (I-A, Ω, p) are equal to 1. Correspo...In this paper, we study some semi-closed 1-set-contractive operators A and investigate the boundary conditions under which the topological degrees of 1-set contractive fields, deg (I-A, Ω, p) are equal to 1. Correspondingly, we can obtain some new fixed point theorems for 1-set-contractive operators which extend and improve many famous theorems such as the Leray-Schauder theorem, and operator equation, etc. Lemma 2.1 generalizes the famous theorem. The calculation of topological degrees and index are important things, which combine the existence of solution of for integration and differential equation and or approximation by iteration technique. So, we apply the effective modification of He’s variation iteration method to solve some nonlinear and linear equations are proceed to examine some a class of integral-differential equations, to illustrate the effectiveness and convenience of this method.展开更多
In this paper,we consider the inverse problem for identifying the source term of the time-fractional equation with a hyper-Bessel operator.First,we prove that this inverse problem is ill-posed,and give the conditional...In this paper,we consider the inverse problem for identifying the source term of the time-fractional equation with a hyper-Bessel operator.First,we prove that this inverse problem is ill-posed,and give the conditional stability.Then,we give the optimal error bound for this inverse problem.Next,we use the fractional Tikhonov regularization method and the fractional Landweber iterative regularization method to restore the stability of the ill-posed problem,and give corresponding error estimates under different regularization parameter selection rules.Finally,we verify the effectiveness of the method through numerical examples.展开更多
This paper studies the eigenfunction expansion method to solve the two dimensional (2D) elasticity problems based on the stress formulation. The fundamental system of partial differential equations of the 2D problem...This paper studies the eigenfunction expansion method to solve the two dimensional (2D) elasticity problems based on the stress formulation. The fundamental system of partial differential equations of the 2D problems is rewritten as an upper tri angular differential system based on the known results, and then the associated upper triangular operator matrix matrix is obtained. By further research, the two simpler com plete orthogonal systems of eigenfunctions in some space are obtained, which belong to the two block operators arising in the operator matrix. Then, a more simple and conve nient general solution to the 2D problem is given by the eigenfunction expansion method. Furthermore, the boundary conditions for the 2D problem, which can be solved by this method, are indicated. Finally, the validity of the obtained results is verified by a specific example.展开更多
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
文摘The present work proposed a new method for the modeling by the finite element method of the acoustic propagation problems in infinite axisymmetric cylindrical guides lined with locally reacting absorbent materials without flow. The method deals with the development of an efficient transparent boundary condition based on DtN operators. The method developed in this study is successfully applied to a straight axisymmetric lined guide by imposing a mode on one of the artificial boundaries of the truncated guide. The results are in good agreement with analytical solutions. Applying the method for a non-uniform axisymmetric lined guide which is a complex case, proved its effectiveness and the results compared to those of PML layers are in very good agreement.
文摘Supply chain management is an essential part of an organisation's sustainable programme.Understanding the concentration of natural environment,public,and economic influence and feasibility of your suppliers and purchasers is becoming progressively familiar as all industries are moving towards a massive sustainable potential.To handle such sort of developments in supply chain management the involvement of fuzzy settings and their generalisations is playing an important role.Keeping in mind this role,the aim of this study is to analyse the role and involvement of complex q-rung orthopair normal fuzzy(CQRONF)information in supply chain management.The major impact of this theory is to analyse the notion of confidence CQRONF weighted averaging,confidence CQRONF ordered weighted averaging,confidence CQRONF hybrid averaging,confidence CQRONF weighted geometric,confidence CQRONF ordered weighted geometric,confidence CQRONF hybrid geometric operators and try to diagnose various properties and results.Furthermore,with the help of the CRITIC and VIKOR models,we diagnosed the novel theory of the CQRONF-CRITIC-VIKOR model to check the sensitivity analysis of the initiated method.Moreover,in the availability of diagnosed operators,we constructed a multi-attribute decision-making tool for finding a beneficial sustainable supplier to handle complex dilemmas.Finally,the initiated operator's efficiency is proved by comparative analysis.
基金This research was supported by the National Natural Science Foundation of China (Nos. 41230210 and 41204074), the Science Foundation of the Education Department of Yunnan Province (No. 2013Z152), and Statoil Company (Contract No. 4502502663).
文摘We propose a symplectic partitioned Runge-Kutta (SPRK) method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete (NAD) operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction (LWC) method and the eighth-order staggered-grid (SG)method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.
文摘A novel nonlocal operator theory based on the variational principle is proposed for the solution of partial differential equations.Common differential operators as well as the variational forms are defined within the context of nonlocal operators.The present nonlocal formulation allows the assembling of the tangent stiffness matrix with ease and simplicity,which is necessary for the eigenvalue analysis such as the waveguide problem.The present formulation is applied to solve the differential electromagnetic vector wave equations based on electric fields.The governing equations are converted into nonlocal integral form.An hourglass energy functional is introduced for the elimination of zeroenergy modes.Finally,the proposed method is validated by testing three classical benchmark problems.
文摘In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is obtained based on an a priori choice of the regularization parameter. Our analysis is not based on the sequential continuity of the normalized duality mapping.
基金Supported by the National Natural Science Foundation of China (10371082)Chinese National Natural Science Foundation Committee Tianyuan Foundation (10526040)Guangzhou University Doctor Foundation (WXF-1001)
文摘The convergence of several Galerkin-Petrov methods, including polynomial collocation and analytic element collocation methods of Toeplitz operators on Dirichlet space, is established. In particular, it is shown that such methods converge if the basis and test function own certain circular symmetry.
基金The NSF(0611005)of Jiangxi Province and the SF(2007293)of Jiangxi Provincial Education Department.
文摘In this paper we develop two multilevel iteration methods for solving linear systems resulting from the Galerkin method and Tikhonov regularization for linear ill-posed problems. The two algorithms and their convergence analyses are presented in an abstract framework.
文摘Unbounded operators can transform arbitrarily small vectors into arbitrarily large vectors—a phenomenon known as instability. Stabilization methods strive to approximate a value of an unbounded operator by applying a family of bounded operators to rough approximate data that do not necessarily lie within the domain of unbounded operator. In this paper we shall be concerned with the stable method of computing values of unbounded operators having perturbations and the stability is established for this method.
文摘In this article, we define the arithmetic operations of generalized trapezoidal picture fuzzy numbers by vertex method which is assembled on a combination of the (α, γ, β)-cut concept and standard interval analysis. Various related properties are explored. Finally, some computations of picture fuzzy functions over generalized picture fuzzy variables are illustrated by using our proposed technique.
基金Project1 990 1 0 0 6 supported by National Natural Science Foundation of China,Doctoral Foundation of China,Chi-na Scholarship council and Laboratory of Computational Physics in Beijing of Chinathe second author is also supportedby the State Major Key
文摘We discuss the incomplete semi-iterative method (ISIM) for an approximate solution of a linear fixed point equations x=Tx+c with a bounded linear operator T acting on a complex Banach space X such that its resolvent has a pole of order k at the point 1. Sufficient conditions for the convergence of ISIM to a solution of x=Tx+c, where c belongs to the range space of R(I-T) k, are established. We show that the ISIM has an attractive feature that it is usually convergent even when the spectral radius of the operator T is greater than 1 and Ind 1T≥1. Applications in finite Markov chain is considered and illustrative examples are reported, showing the convergence rate of the ISIM is very high.
基金supported by the President Foundation of the Chinese Academy of Sciences and National Natural Science Foundation of China under Grant No.10475657
文摘<Abstract>We show that the recently proposed invariant eigen-operator (lEO) method can be successfully applied to solving energy levels for SSH Hamiltonian describing Peierls phase transition.The electronic energy band of compound lattice is also studied by IEO method.
文摘The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ability of the discrete forms expressing to the element functions was talked about. In discrete operator difference method, the displacements of the elements can be reproduced exactly in the discrete forms whether the displacements are conforming or not. According to this point, discrete operator difference method is a method with good performance.
基金The project supported by National Natural Science Foundation of China(9713008)Zhejiang Natural Science Foundation Special Funds No. RC.9601
文摘This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.
基金Supported in part by the Natural Science Foundation of China under grants 10371137and 10201034Foundation of Doctoral Program of National Higher Education of China under under grant 20030558008Guangdong Provincial Natural Science Foundation of China u
文摘We introduce multilevel augmentation methods for solving operator equations based on direct sum decompositions of the range space of the operator and the solution space of the operator equation and a matrix splitting scheme. We establish a general setting for the analysis of these methods, showing that the methods yield approximate solutions of the same convergence order as the best approximation from the subspace. These augmentation methods allow us to develop fast, accurate and stable nonconventional numerical algorithms for solving operator equations. In particular, for second kind equations, special splitting techniques are proposed to develop such algorithms. These algorithms are then applied to solve the linear systems resulting from matrix compression schemes using wavelet-like functions for solving Fredholm integral equations of the second kind. For this special case, a complete analysis for computational complexity and convergence order is presented. Numerical examples are included to demonstrate the efficiency and accuracy of the methods. In these examples we use the proposed augmentation method to solve large scale linear systems resulting from the recently developed wavelet Galerkin methods and fast collocation methods applied to integral equations of the secondkind. Our numerical results confirm that this augmentation method is particularly efficient for solving large scale linear systems induced from wavelet compression schemes.
文摘Normal form theory is a very effective method when we study degenerate bifurcations of nonlinear dynamical systems. In this paper by using adjoint operator method, normal forms of order 3 and 4 for nonlinear dynamical system with nilpotent linear part and Z(2)-asymmetry are computed. According to normal forms obtained, universal unfoldings for some degenerate bifurcation cases of codimension 3 and simple global characterizations, are studied.
基金supported by the National Natural Science Foundation of China (Grant No.10874174)the Specialized Research Fund for the Doctoral Program of Higher Education (Grant No.20070358009)
文摘For classical Hamiltonian with general form H = 1/2∑ijMijpipj+1/2∑ijLijqiqj we find a new convenient way to obtain its normal coordinates, namely, let H be quantised and then employ the invariant eigen-operator (IEO) method (Fan et al. 2004 Phys. Lett. A 321 75) to derive them. The general matrix equation, which relies on M and L, for obtaining the normal coordinates of H is derived.
文摘In this paper, we study some semi-closed 1-set-contractive operators A and investigate the boundary conditions under which the topological degrees of 1-set contractive fields, deg (I-A, Ω, p) are equal to 1. Correspondingly, we can obtain some new fixed point theorems for 1-set-contractive operators which extend and improve many famous theorems such as the Leray-Schauder theorem, and operator equation, etc. Lemma 2.1 generalizes the famous theorem. The calculation of topological degrees and index are important things, which combine the existence of solution of for integration and differential equation and or approximation by iteration technique. So, we apply the effective modification of He’s variation iteration method to solve some nonlinear and linear equations are proceed to examine some a class of integral-differential equations, to illustrate the effectiveness and convenience of this method.
基金supported by the National Natural Science Foundation of China(11961044)the Doctor Fund of Lan Zhou University of Technologythe Natural Science Foundation of Gansu Provice(21JR7RA214)。
文摘In this paper,we consider the inverse problem for identifying the source term of the time-fractional equation with a hyper-Bessel operator.First,we prove that this inverse problem is ill-posed,and give the conditional stability.Then,we give the optimal error bound for this inverse problem.Next,we use the fractional Tikhonov regularization method and the fractional Landweber iterative regularization method to restore the stability of the ill-posed problem,and give corresponding error estimates under different regularization parameter selection rules.Finally,we verify the effectiveness of the method through numerical examples.
基金supported by the Specialized Research Fund for the Doctoral Program of Higher Education of China (No. 20070126002)the National Natural Science Foundation of China (No. 10962004)
文摘This paper studies the eigenfunction expansion method to solve the two dimensional (2D) elasticity problems based on the stress formulation. The fundamental system of partial differential equations of the 2D problems is rewritten as an upper tri angular differential system based on the known results, and then the associated upper triangular operator matrix matrix is obtained. By further research, the two simpler com plete orthogonal systems of eigenfunctions in some space are obtained, which belong to the two block operators arising in the operator matrix. Then, a more simple and conve nient general solution to the 2D problem is given by the eigenfunction expansion method. Furthermore, the boundary conditions for the 2D problem, which can be solved by this method, are indicated. Finally, the validity of the obtained results is verified by a specific example.