This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method (FEM) by recursive application of the one-dimensional (1D) element energy pro...This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method (FEM) by recursive application of the one-dimensional (1D) element energy projection (EEP) technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element (FE) model is reached. This conceptual dimension-by- dimension (D-by-D) discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional (2D) and three-dimensional (3D) problems can be obtained over the domain. Numerical examples of 3D Poisson's equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.展开更多
Singular initial value problems arise in solving one-dimensional steady transonic flow of dualmode scramjet. The existing solution method has the problems of large initial value errors in principles. This paper puts f...Singular initial value problems arise in solving one-dimensional steady transonic flow of dualmode scramjet. The existing solution method has the problems of large initial value errors in principles. This paper puts forward an improved algorithm based on variable transformation, and constructs a nonsingular one-dimensional steady transonic flow equation by defining a new variable. The improved algorithm can eliminate the singularity of the differential equation, and can solve the singular initial value problems of one-dimensional steady transonic flow of dual-mode scramjet.展开更多
By using Mawhin's continuation theorem, the existence of a solution for a class of m-point boundary value problem at resonance with one-dimensional p-Laplacian is obtained. An example is given to demonstrate the main...By using Mawhin's continuation theorem, the existence of a solution for a class of m-point boundary value problem at resonance with one-dimensional p-Laplacian is obtained. An example is given to demonstrate the main result of this paper.展开更多
In this paper, we present an explicit one-step method for solving periodic initial value problems of second order ordinary differential equations. The method is P-stable, and of first algebraic order and high phase-la...In this paper, we present an explicit one-step method for solving periodic initial value problems of second order ordinary differential equations. The method is P-stable, and of first algebraic order and high phase-lag order. To improve the algebraic order, we give a composition second order scheme with the proposed method and its adjoint. We report some numerical results to illustrate the efficiency of our methods.展开更多
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.展开更多
The Kuhn-Tucker theorem in nondifferential form is a well-known classical optimality criterion for a convex programming problems which is true for a convex problem in the case when a Kuhn-Tucker vector exists. It is n...The Kuhn-Tucker theorem in nondifferential form is a well-known classical optimality criterion for a convex programming problems which is true for a convex problem in the case when a Kuhn-Tucker vector exists. It is natural to extract two features connected with the classical theorem. The first of them consists in its possible “impracticability” (the Kuhn-Tucker vector does not exist). The second feature is connected with possible “instability” of the classical theorem with respect to the errors in the initial data. The article deals with the so-called regularized Kuhn-Tucker theorem in nondifferential sequential form which contains its classical analogue. A proof of the regularized theorem is based on the dual regularization method. This theorem is an assertion without regularity assumptions in terms of minimizing sequences about possibility of approximation of the solution of the convex programming problem by minimizers of its regular Lagrangian, that are constructively generated by means of the dual regularization method. The major distinctive property of the regularized Kuhn-Tucker theorem consists that it is free from two lacks of its classical analogue specified above. The last circumstance opens possibilities of its application for solving various ill-posed problems of optimization, optimal control, inverse problems.展开更多
Since the time step of the traditional finite-difference time-domain(FDTD) method is limited by the small grid size, it is inefficient when dealing with the electromagnetic problems of multi-scale structures.Therefore...Since the time step of the traditional finite-difference time-domain(FDTD) method is limited by the small grid size, it is inefficient when dealing with the electromagnetic problems of multi-scale structures.Therefore, the explicit and unconditionally stable FDTD(US-FDTD) approach has been developed to break through the limitation of Courant–Friedrich–Levy(CFL) condition.However, the eigenvalues and eigenvectors of the system matrix must be calculated before the time iteration in the explicit US-FDTD.Moreover, the eigenvalue decomposition is also time consuming, especially for complex electromagnetic problems in practical application.In addition, compared with the traditional FDTD method, the explicit US-FDTD method is more difficult to introduce the absorbing boundary and plane wave.To solve the drawbacks of the traditional FDTD and the explicit US-FDTD, a new hybrid FDTD algorithm is proposed in this paper.This combines the explicit US-FDTD with the traditional FDTD, which not only overcomes the limitation of CFL condition but also reduces the system matrix dimension, and introduces the plane wave and the perfectly matched layer(PML) absorption boundary conveniently.With the hybrid algorithm, the calculation of the eigenvalues is only required in the fine mesh region and adjacent coarse mesh region.Therefore, the calculation efficiency is greatly enhanced.Furthermore, the plane wave and the absorption boundary introduction of the traditional FDTD method can be directly utilized.Numerical results demonstrate the effectiveness, accuracy, stability, and convenience of this hybrid algorithm.展开更多
The stability problem of the Rindler spacetime is carefully studies by using the scalar wave perturbation. Using two different coordinate systems, the scalar wave equation is investigated. The results are different in...The stability problem of the Rindler spacetime is carefully studies by using the scalar wave perturbation. Using two different coordinate systems, the scalar wave equation is investigated. The results are different in the two cases. They are analysed and compared with each other in detail. The following conclusions are obtained: (a) the Rindler spacetime as a whole is not stable; (b) the Rindler spacetime can exist stably only as part of the Minkowski spacetime, and the Minkowski spacetime can be a real entity independently; (c) there are some defects for the scalar wave equation written by the Rindler coordinates, and it is unsuitable for the investigation of the stability properties of the Rindler spacetime. All these results may shed some light on the stability properties of the Schwarzschild black hole. It is natural and reasonable for one to infer that: (a) perhaps the Regge-Wheeler equation is not sufficient to determine the stable properties; (b) the Schwarzschild black hole as a whole might be really unstable; (c) the Kruskal spacetime is stable and can exist as a real physical entity; whereas the Schwarzschild black hole can occur only as part of the Kruskal spacetime.展开更多
This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing t...This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing the order of approximation in fixed number of subdomains,rather than by resorting to a further partitioning.The stability and the convergence of this method are proved.展开更多
In this paper, the authors consider the inverse piston problem for the system of one-dimensional isentropic flow and obtain that, under suitable conditions, the piston velocity can be uniquely determined by the initia...In this paper, the authors consider the inverse piston problem for the system of one-dimensional isentropic flow and obtain that, under suitable conditions, the piston velocity can be uniquely determined by the initial state of the gas on the right side of the piston and the position of the forward shock.展开更多
In this paper, we consider a multipoint boundary value problem for one-dimensional p-Laplacian. Using a fixed point theorem due to Bai and Ge, we study the existence of at least three positive solutions to the boundar...In this paper, we consider a multipoint boundary value problem for one-dimensional p-Laplacian. Using a fixed point theorem due to Bai and Ge, we study the existence of at least three positive solutions to the boundary value problem. In this problem, the nonlinear term explicitly involves a first-order derivative, which is different from some previous ones.展开更多
The stable marriage problem stands for a class of assignment problems and has been attracting high attention from operational research community. Various versions of the problem and corresponding algorithms have been ...The stable marriage problem stands for a class of assignment problems and has been attracting high attention from operational research community. Various versions of the problem and corresponding algorithms have been published. This paper introduces the characteristics of GPMS in object\|oriented simulation modelling through constructing a simulation model of the problem.展开更多
This paper disccusses the inverse scattering problem for one-dimensional Schrodinger operatorsrelated to the general Stark effect. We provide a ganeral framework which can be applied both to theStark-effect operator a...This paper disccusses the inverse scattering problem for one-dimensional Schrodinger operatorsrelated to the general Stark effect. We provide a ganeral framework which can be applied both to theStark-effect operator and the ordinary Schrodinger operator.展开更多
By using Mawhin's continuation theorem, the existence of solution for a class of m-point boundary value problem at resonance with one-dimensional p-Laplacian is obtained. An example is given to demonstrate the main r...By using Mawhin's continuation theorem, the existence of solution for a class of m-point boundary value problem at resonance with one-dimensional p-Laplacian is obtained. An example is given to demonstrate the main result of the paper.展开更多
基金supported by the National Natural Science Foundation of China(Nos.51378293 and 51078199)
文摘This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method (FEM) by recursive application of the one-dimensional (1D) element energy projection (EEP) technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element (FE) model is reached. This conceptual dimension-by- dimension (D-by-D) discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional (2D) and three-dimensional (3D) problems can be obtained over the domain. Numerical examples of 3D Poisson's equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.
基金Hi TechResearchandDevelopmentProgramofChina(2002AA723011),OutstandingYouthFoundationofHeilongjiang Province
文摘Singular initial value problems arise in solving one-dimensional steady transonic flow of dualmode scramjet. The existing solution method has the problems of large initial value errors in principles. This paper puts forward an improved algorithm based on variable transformation, and constructs a nonsingular one-dimensional steady transonic flow equation by defining a new variable. The improved algorithm can eliminate the singularity of the differential equation, and can solve the singular initial value problems of one-dimensional steady transonic flow of dual-mode scramjet.
基金The NSF (Kj2007b055) of Anhui Educational Departmentthe Youth Project Foundation (2007jqL101,2007jqL102) of Anhui Educational Department.
文摘By using Mawhin's continuation theorem, the existence of a solution for a class of m-point boundary value problem at resonance with one-dimensional p-Laplacian is obtained. An example is given to demonstrate the main result of this paper.
基金The project is supported by NSF of Anhui Province(No.2005jk218), China.
文摘In this paper, we present an explicit one-step method for solving periodic initial value problems of second order ordinary differential equations. The method is P-stable, and of first algebraic order and high phase-lag order. To improve the algebraic order, we give a composition second order scheme with the proposed method and its adjoint. We report some numerical results to illustrate the efficiency of our methods.
文摘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.
文摘The Kuhn-Tucker theorem in nondifferential form is a well-known classical optimality criterion for a convex programming problems which is true for a convex problem in the case when a Kuhn-Tucker vector exists. It is natural to extract two features connected with the classical theorem. The first of them consists in its possible “impracticability” (the Kuhn-Tucker vector does not exist). The second feature is connected with possible “instability” of the classical theorem with respect to the errors in the initial data. The article deals with the so-called regularized Kuhn-Tucker theorem in nondifferential sequential form which contains its classical analogue. A proof of the regularized theorem is based on the dual regularization method. This theorem is an assertion without regularity assumptions in terms of minimizing sequences about possibility of approximation of the solution of the convex programming problem by minimizers of its regular Lagrangian, that are constructively generated by means of the dual regularization method. The major distinctive property of the regularized Kuhn-Tucker theorem consists that it is free from two lacks of its classical analogue specified above. The last circumstance opens possibilities of its application for solving various ill-posed problems of optimization, optimal control, inverse problems.
基金Project supported by the National Natural Science Foundation of China(Grant No.61571348)the Equipment Pre-Research Foundation of China(Grant No.61405180202)
文摘Since the time step of the traditional finite-difference time-domain(FDTD) method is limited by the small grid size, it is inefficient when dealing with the electromagnetic problems of multi-scale structures.Therefore, the explicit and unconditionally stable FDTD(US-FDTD) approach has been developed to break through the limitation of Courant–Friedrich–Levy(CFL) condition.However, the eigenvalues and eigenvectors of the system matrix must be calculated before the time iteration in the explicit US-FDTD.Moreover, the eigenvalue decomposition is also time consuming, especially for complex electromagnetic problems in practical application.In addition, compared with the traditional FDTD method, the explicit US-FDTD method is more difficult to introduce the absorbing boundary and plane wave.To solve the drawbacks of the traditional FDTD and the explicit US-FDTD, a new hybrid FDTD algorithm is proposed in this paper.This combines the explicit US-FDTD with the traditional FDTD, which not only overcomes the limitation of CFL condition but also reduces the system matrix dimension, and introduces the plane wave and the perfectly matched layer(PML) absorption boundary conveniently.With the hybrid algorithm, the calculation of the eigenvalues is only required in the fine mesh region and adjacent coarse mesh region.Therefore, the calculation efficiency is greatly enhanced.Furthermore, the plane wave and the absorption boundary introduction of the traditional FDTD method can be directly utilized.Numerical results demonstrate the effectiveness, accuracy, stability, and convenience of this hybrid algorithm.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10475013, 10375087 and 10373003), the National Basic Research Program (Grant No 2004CB318000) and National Science Foundation for Post-Doctoral Scientists of China.
文摘The stability problem of the Rindler spacetime is carefully studies by using the scalar wave perturbation. Using two different coordinate systems, the scalar wave equation is investigated. The results are different in the two cases. They are analysed and compared with each other in detail. The following conclusions are obtained: (a) the Rindler spacetime as a whole is not stable; (b) the Rindler spacetime can exist stably only as part of the Minkowski spacetime, and the Minkowski spacetime can be a real entity independently; (c) there are some defects for the scalar wave equation written by the Rindler coordinates, and it is unsuitable for the investigation of the stability properties of the Rindler spacetime. All these results may shed some light on the stability properties of the Schwarzschild black hole. It is natural and reasonable for one to infer that: (a) perhaps the Regge-Wheeler equation is not sufficient to determine the stable properties; (b) the Schwarzschild black hole as a whole might be really unstable; (c) the Kruskal spacetime is stable and can exist as a real physical entity; whereas the Schwarzschild black hole can occur only as part of the Kruskal spacetime.
文摘This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing the order of approximation in fixed number of subdomains,rather than by resorting to a further partitioning.The stability and the convergence of this method are proved.
文摘In this paper, the authors consider the inverse piston problem for the system of one-dimensional isentropic flow and obtain that, under suitable conditions, the piston velocity can be uniquely determined by the initial state of the gas on the right side of the piston and the position of the forward shock.
文摘In this paper, we consider a multipoint boundary value problem for one-dimensional p-Laplacian. Using a fixed point theorem due to Bai and Ge, we study the existence of at least three positive solutions to the boundary value problem. In this problem, the nonlinear term explicitly involves a first-order derivative, which is different from some previous ones.
文摘The stable marriage problem stands for a class of assignment problems and has been attracting high attention from operational research community. Various versions of the problem and corresponding algorithms have been published. This paper introduces the characteristics of GPMS in object\|oriented simulation modelling through constructing a simulation model of the problem.
文摘This paper disccusses the inverse scattering problem for one-dimensional Schrodinger operatorsrelated to the general Stark effect. We provide a ganeral framework which can be applied both to theStark-effect operator and the ordinary Schrodinger operator.
基金The work is sponsored by Youth Project Foundation of Anhui Educational Dept.(2007jqL101,2007jqL102)Natural Science Foundation of Anhui Educational Dept.(kj2007b055).
文摘By using Mawhin's continuation theorem, the existence of solution for a class of m-point boundary value problem at resonance with one-dimensional p-Laplacian is obtained. An example is given to demonstrate the main result of the paper.