A reduced finite element formulation based on POD method for two-dimensional solute transport problems

A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for two-dimensional solute transport problems with real practical applied background such that it is reduced into a reduced FE formulation with lower dimensions and high enough accuracy. The error estimates between the reduced POD FE solutions and the usual FE solutions are provided. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD FE method.     

A reduced finite volume element formulation and numerical simulations based on POD for parabolic problems

A proper orthogonal decomposition (POD) method is applied to a usual finite volume element (FVE) formulation for parabolic equations such that it is reduced to a POD FVE formulation with lower dimensions and high enough accuracy. The error estimates between the reduced POD FVE solution and the usual FVE solution are analyzed. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is also shown that the reduced POD FVE formulation based on POD method is both feasible and highly efficient.     

A reduced FE formulation based on POD method for hyperbolic equations

A proper orthogonal decomposition (POD) method was successfully used in the reduced-order modeling of complex systems. In this paper, we extend the applications of POD method, namely, apply POD method to a classical finite element (FE) formulation for second-order hyperbolic equations with real practical applied background, establish a reduced FE formulation with lower dimensions and high enough accuracy, and provide the error estimates between the reduced FE solutions and the classical FE solutions and the implementation of algorithm for solving reduced FE formulation so as to provide scientific theoretic basis for service applications. Some numerical examples illustrate the fact that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced FE formulation based on POD method is feasible and efficient for solving FE formulation for second-order hyperbolic equations.     

抛物型方程基于POD方法的时间二阶精度CN有限元降维格式

将特征正交分解(proper orthogonal decomposition, 简记为POD) 方法应用于抛物型方程通常的时间二阶精度Crank-Nicolson (简记为CN) 有限元格式, 简化其为一个自由度极少的时间二阶精度CN 有限元降维格式, 并给出简化的时间二阶精度CN 有限元解的误差分析. 数值例子表明在简化的时间二阶精度CN 有限元解和通常的时间二阶精度CN 有限元解之间的误差足够小的情况下, 简化的时间二阶精度CN 有限元格式能大大地节省自由度, 而且时间步长可以比时间一阶精度的格式取大10 倍, 以至能更快计算到所要时刻数值解, 减少计算机计算过程的截断误差, 提高计算速度和计算精度,从而验证降维时间二阶精度CN 有限元格式用于解类似于抛物型方程的时间依赖方程是很有效的.     

A reduced MFE formulation based on POD for the non-stationary conduction-convection problems

In this article, a reduced mixed finite element (MFE) formulation based on proper orthogonal decomposition (POD) for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.     

Convergence of one-leg methods for nonlinear neutral delay integro-differential equations

Some convergence results of one-leg methods for nonlinear neutral delay integro-differential equations (NDIDEs) are obtained. It is proved that a one-leg method is E (or EB) -convergent of order p for nonlinear NDIDEs if and only if it is A-stable and consistent of order p in classical sense for ODEs, where p = 1, 2. A numerical example that confirms the theoretical results is given in the end of this paper. This work was supported by National Natural Science Foundation of China (Grant No. 10871164), the Natural Science Foundation of Hunan Province (Grant No. 08JJ6002), and the Scientific Research Fund of Changsha University of Science and Technology (Grant No. 1004259)     

A level set method to reconstruct the discontinuity of the conductivity in EIT

In this paper, one level set method is applied to finding the interface of discontinuity of the conductivity in EIT(electrical impedance tomography) problem. By choosing one suitable velocity function, a level set reconstruction algorithm is proposed. The theoretical results for EIT problem and regularization are given. Finally the numerical examples demonstrate that the reconstruction algorithm is efficient and stable. The work was supported by the National Natural Science Foundation of China (Grant Nos. 10431030, 10771138), the Shanghai Natural Science Foundation (Grant No. 07JC14001), the National Basic Research Program (Grant No. 2005CB321701) and Ministry of Education of China and State Administration of Foreign Experts Affairs of China under an 111 project (Grant No. B08018).     

Consistency and asymptotic normality of profilekernel and backfitting estimators in semiparametric reproductive dispersion nonlinear models

Semiparametric reproductive dispersion nonlinear model (SRDNM) is an extension of nonlinear reproductive dispersion models and semiparametric nonlinear regression models, and includes semiparametric nonlinear model and semiparametric generalized linear model as its special cases. Based on the local kernel estimate of nonparametric component, profile-kernel and backfitting estimators of parameters of interest are proposed in SRDNM, and theoretical comparison of both estimators is also investigated in this paper. Under some regularity conditions, strong consistency and asymptotic normality of two estimators are proved. It is shown that the backfitting method produces a larger asymptotic variance than that for the profile-kernel method. A simulation study and a real example are used to illustrate the proposed methodologies. This work was supported by National Natural Science Foundation of China (Grant Nos. 10561008, 10761011), Natural Science Foundation of Department of Education of Zhejiang Province (Grant No. Y200805073), PhD Special Scientific Research Foundation of Chinese University (Grant No. 20060673002) and Program for New Century Excellent Talents in University (Grant No. NCET-07-0737)     

Numerical simulation based on POD for two-dimensional solute transport problems

A proper orthogonal decomposition (POD) method is applied to a usual finite element scheme for two-dimensional solute transport problems such that it is reduced into a reduced finite element formulation with lower dimensions and high enough accuracy. Numerical examples show that the results of numerical computations are consistent with accurate solutions. Moreover, this validates the feasibility and efficiency of POD method.     

Homoclinic and heteroclinic flows in global attractor for Davey-stewartson II equation

By the variable transformation and generalized Hirota method, exact homoclinic and heteroclinic solutions for Davey-Stewartson II （DSII） equation are obtained. For perturbed DSII equation, the existence of a global attractor is proved. The persistence of homoclinic and heteroclinic flows is investigated, and the special homoclinic and heteroclinic structure in attractors is shown.     

A reduced finite difference scheme based on singular value decomposition and proper orthogonal decomposition for Burgers equation

In this article, a reduced optimizing finite difference scheme (FDS) based on singular value decomposition (SVD) and proper orthogonal decomposition (POD) for Burgers equation is presented. Also the error estimates between the usual finite difference solution and the POD solution of reduced optimizing FDS are analyzed. It is shown by considering the results obtained for numerical simulations of cavity flows that the error between the POD solution of reduced optimizing FDS and the solution of the usual FDS is consistent with theoretical results. Moreover, it is also shown that the reduced optimizing FDS is feasible and efficient.     

Arrangement transformation approach to state-to-state quantum reactive scattering H + DH→DH + H, HH + D

Arrangement transformation approach (ATA) for doing state-to-state quantum reactive scattering calculations of atom-diatom systems is given. The state-to-state results of H + DH on ISTH potential energy surface are presented, and it can be deduced from the results that the ATA method can be used for the system with more than three atoms. Project supported by the National Natural Science Foundation of China (Grant No. 19774038) and Shandong Science Foundation (Grant No. Y96A08012).     

热传导对流方程基于POD的差分格式

本文用奇值分解和特征投影分解(proper orthogonal decomposition,简记为POD)研究热传导对流方程,导出其基于POD的一种简化的差分格式,并分析通常的差分格式的解和基于POD的简化的差分格式的解之间的误差估计.最后用方腔流数值例子验证本文的理论的正确性,从而验证了用基于POD的简化的差分格式解热传导对流方程的有效性.     

The truncatedSVD as a method for regularization

The truncated singular value decomposition (SVD) is considered as a method for regularization of ill-posed linear least squares problems. In particular, the truncated SVD solution is compared with the usual regularized solution. Necessary conditions are defined in which the two methods will yield similar results. This investigation suggests the truncated SVD as a favorable alternative to standard-form regularization in cases of ill-conditioned matrices with well-determined numerical rank.This work was carried out while the author visited the Dept. of Computer Science, Stanford University, California, U.S.A., and was supported in part by National Science Foundation Grant Number DCR 8412314, by a Fulbright Supplementary Grant, and by the Danish Space Board.     

A reduced FVE formulation based on POD method and error analysis for two-dimensional viscoelastic problem

Proper orthogonal decomposition (POD) method has been successfully used in the reduced-order modeling of complex systems. In this paper, we extend the applications of POD method, i.e., combine a classical finite volume element (FVE) method with POD method to establish a reduced FVE formulation with lower dimensions and sufficiently high accuracy for two-dimensional viscoelastic problem with real practical applied background, and analyze the errors between the reduced POD FVE solution and the classical FVE solution so as to provide scientific theoretic basis for service applications. Some numerical examples illustrate the fact that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is also shown that the reduced FVE formulation based on POD method is feasible and efficient for solving two-dimensional viscoelastic problem.     

Numerical simulation of a modified KdV equation on the whole real axis

The numerical simulation of the solution to a modified KdV equation on the whole real axis is considered in this paper. Based on the work of Fokas (Comm Pure Appl Math 58(5):639–670, 2005), a kind of exact nonreflecting boundary conditions which are suitable for numerical purposes are presented with the inverse scattering theory. With these boundary conditions imposed on the artificially introduced boundary points, a reduced problem defined on a finite computational interval is formulated. The discretization of the nonreflecting boundary conditions is studied in detail, and a dual-Petrov–Galerkin spectral method is proposed for the numerical solution to the reduced problem. Some numerical tests are given, which validate the effectiveness, and suggest the stability of the proposed scheme.Supported by the National Natural Science Foundation of China under Grant No. 10401020, the Alexander von Humboldt Foundation, and the Key Project of China High Performance Scientific Computation Research.     

How to decrease the levels of modular integrals

A method to decrease the level of a modular integral with a given level is advanced. Project supported by the National Natural Science Foundation of China (Grant No. 19501009) and the Natural Science Foundation of Guangdong Province.     

A convergence analysis of the inexact Rayleigh quotient iteration and simplified Jacobi-Davidson method for the large Hermitian matrix eigenproblem

The inexact Rayleigh quotient iteration (RQI) is used for computing the smallest eigenpair of a large Hermitian matrix. Under certain condition, the method was proved to converge quadratically in literature. However, it is shown in this paper that under the original given condition the inexact RQI may not quadratically converge to the desired eigenpair and even may misconverge to some other undesired eigenpair. A new condition, called the uniform positiveness condition, is given that can fix misconvergence problem and ensure the quadratic convergence of the inexact RQI. An alternative to the inexact RQI is the Jacobi-Davidson (JD) method without subspace acceleration. A new proof of its linear convergence is presented and a sharper bound is established in the paper. All the results are verified and analyzed by numerical experiments. This work was supported by the National Natural Science Foundation of China (Grant Nos. 10471074, 10771116) and the Doctoral Program of the Ministry of Education of China (Grant No. 20060003003)     

On the Finite Termination of an Entropy Function Based Non-Interior Continuation Method for Vertical Linear Complementarity Problems

By using a smooth entropy function to approximate the non-smooth max-type function, a vertical linear complementarity problem (VLCP) can be treated as a family of parameterized smooth equations. A Newton-type method with a testing procedure is proposed to solve such a system. We show that under some milder than usual assumptions the proposed algorithm finds an exact solution of VLCP in a finite number of iterations. Some computational results are included to illustrate the potential of this approach. This author’s work was partially supported by the National Natural Science Foundation of China (Grant Nos. 10271002 and 10401038). This author’s work was partially supported by the Scientific Research Foundation of Tianjin University for the Returned Overseas Chinese Scholars and the Scientific Research Foundation of Liu Hui Center for Applied Mathematics, Nankai University-Tianjin University.     

A Trust Region Method for Optimization Problem with Singular Solutions

In this paper, we propose a trust region method for minimizing a function whose Hessian matrix at the solutions may be singular. The global convergence of the method is obtained under mild conditions. Moreover, we show that if the objective function is LC ² function, the method possesses local superlinear convergence under the local error bound condition without the requirement of isolated nonsingular solution. This is the first regularized Newton method with trust region technique which possesses local superlinear (quadratic) convergence without the assumption that the Hessian of the objective function at the solution is nonsingular. Preliminary numerical experiments show the efficiency of the method. This work is partly supported by the National Natural Science Foundation of China (Grant Nos. 70302003, 10571106, 60503004, 70671100) and Science Foundation of Beijing Jiaotong University (2007RC014).