This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the erro...This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.展开更多
The nonoverlapping domain decomposition method for parabolic partialdifferential equation on general domain is considered. A kind of domain decompositionthat uses the finite element procedure is given. The problem ove...The nonoverlapping domain decomposition method for parabolic partialdifferential equation on general domain is considered. A kind of domain decompositionthat uses the finite element procedure is given. The problem over the domains can beimplemented on parallel computer. Convergence analysis is also presented.展开更多
Abstract Two phase,incompressible miscible flow in porous media is governed by a system of nonlinear partial differential equations.The pressure equation,which is elliptic in appearance,is discretized by a standard fi...Abstract Two phase,incompressible miscible flow in porous media is governed by a system of nonlinear partial differential equations.The pressure equation,which is elliptic in appearance,is discretized by a standard five points difference method.The concentration equation is treated by an implicit finite difference method that applies a form of the method of characteristics to the transport terms.A class of biquadratic interpolation is introduced for the method of characteristics.Convergence rate is proved to be O(Δt+h 2).展开更多
A interior point scaling projected reduced Hessian method with combination of nonmonotonic backtracking technique and trust region strategy for nonlinear equality constrained optimization with nonegative constraint on...A interior point scaling projected reduced Hessian method with combination of nonmonotonic backtracking technique and trust region strategy for nonlinear equality constrained optimization with nonegative constraint on variables is proposed. In order to deal with large problems, a pair of trust region subproblems in horizontal and vertical subspaces is used to replace the general full trust region subproblem. The horizontal trust region subproblem in the algorithm is only a general trust region subproblem while the vertical trust region subproblem is defined by a parameter size of the vertical direction subject only to an ellipsoidal constraint. Both trust region strategy and line search technique at each iteration switch to obtaining a backtracking step generated by the two trust region subproblems. By adopting the l 1 penalty function as the merit function, the global convergence and fast local convergence rate of the proposed algorithm are established under some reasonable conditions. A nonmonotonic criterion and the second order correction step are used to overcome Maratos effect and speed up the convergence progress in some ill-conditioned cases. MR Subject Classification 90C30 - 65K05 - 49M40 Keywords trust region method - backtracking step - reduced Hessian - nonmonotonic technique - interior point Supported partially by the National Natural Science Foundation of China (10071050), Science Foundation (02ZA14070) of Shanghai Technical Sciences Committee and Science Foundation (02DK06) of Shanghai Education Committee.展开更多
Based on an asymptotic expansion of finite element,a new extrapolation formula and extrapolation cascadic multigrid method(EXCMG)are proposed,in which the new extrapolation and quadratic interpolation are used to prov...Based on an asymptotic expansion of finite element,a new extrapolation formula and extrapolation cascadic multigrid method(EXCMG)are proposed,in which the new extrapolation and quadratic interpolation are used to provide a better initial value on refined grid.In the case of triple grids,the error of the new initial value is analyzed in detail.A larger scale computation is completed in PC.展开更多
This paper discusses the order-preserving convergence for spectral approximation of the self-adjoint completely continuous operator T.Under the condition that the approximate operator Th converges to T in norm,it is p...This paper discusses the order-preserving convergence for spectral approximation of the self-adjoint completely continuous operator T.Under the condition that the approximate operator Th converges to T in norm,it is proven that the k-th eigenvalue of Th converges to the k-th eigenvalue of T.(We sorted the positive eigenvalues in decreasing order and negative eigenvalues in increasing order.) Then we apply this result to conforming elements,nonconforming elements and mixed elements of self-adjoint elliptic differential operators eigenvalue problems,and prove that the k-th approximate eigenvalue obtained by these methods converges to the k-th exact eigenvalue.展开更多
A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accura...A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.展开更多
In this paper,we consider the cascadic multigrid method for a parabolic type equation.Backward Euler approximation in time and linear finite element approximation in space are employed.A stability result is establishe...In this paper,we consider the cascadic multigrid method for a parabolic type equation.Backward Euler approximation in time and linear finite element approximation in space are employed.A stability result is established under some conditions on the smoother.Using new and sharper estimates for the smoothers that reflect the precise dependence on the time step and the spatial mesh parameter,these conditions are verified for a number of popular smoothers.Optimal error bound sare derived for both smooth and non-smooth data.Iteration strategies guaranteeing both the optimal accuracy and the optimal complexity are presented.展开更多
The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from t...The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.展开更多
Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which ...Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.展开更多
Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolatio...Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolation error of anisotropic elements based on Newton's formula of polynomial interpolation as well as its applications.展开更多
For numerical computations of multiple solutions of the nonlinear elliptic problem Δu + f(u = 0 in Ω, u = 0 on Γ, a search-extension method (SEM) was proposed and systematically studied by the authors. This paper s...For numerical computations of multiple solutions of the nonlinear elliptic problem Δu + f(u = 0 in Ω, u = 0 on Γ, a search-extension method (SEM) was proposed and systematically studied by the authors. This paper shall complete its theoretical analysis. It is assumed that the nonlinearity is non-convex and its solution is isolated, under some conditions the corresponding linearized problem has a unique solution. By use of the compactness of the solution family and the contradiction argument, in general conditions, the high order regularity of the solution u ∈ H 1+α, α > 0 is proved. Assume that some initial value searched by suitably many eigenbases is already fallen into the neighborhood of the isolated solution, then the optimal error estimates of its nonlinear finite element approximation are shown by the duality argument and continuation method.展开更多
This study discusses generalized Rayleigh quotient and high efficiency finite element discretization schemes. Some results are as follows: 1) Rayleigh quotient accelerate technique is extended to nonselfadjoint proble...This study discusses generalized Rayleigh quotient and high efficiency finite element discretization schemes. Some results are as follows: 1) Rayleigh quotient accelerate technique is extended to nonselfadjoint problems. Generalized Rayleigh quotients of operator form and weak form are defined and the basic relationship between approximate eigenfunction and its generalized Rayleigh quotient is established. 2) New error estimates are obtained by replacing the ascent of exact eigenvalue with the ascent of finite element approximate eigenvalue. 3) Based on the work of Xu Jinchao and Zhou Aihui, finite element two-grid discretization schemes are established to solve nonselfadjoint elliptic differential operator eigenvalue problems and these schemes are used in both conforming finite element and non-conforming finite element. Besides, the efficiency of the schemes is proved by both theoretical analysis and numerical experiments. 4) Iterated Galerkin method, interpolated correction method and gradient recovery for selfadjoint elliptic differential operator eigenvalue problems are extended to nonselfadjoint elliptic differential operator eigenvalue problems.展开更多
Consider an AMG for the linear system Au=f. Up to now, only the uniform convergence of two-level AMG is proved for symmetric and positive definite L-matrices with weak diagonal dominance. Using the new form (1), we ex...Consider an AMG for the linear system Au=f. Up to now, only the uniform convergence of two-level AMG is proved for symmetric and positive definite L-matrices with weak diagonal dominance. Using the new form (1), we extend the results in [1] to the case that A is a general symmetric and positive definite matrix with weak diagonal dominance. In the following, we shall use the same notations as in [1].展开更多
In this paper, we first introduce interpolation operator of projection type in three dimen- sions, from which we derive weak estimates for tricubic block finite elements. Then using the estimate for the W 2, 1-seminor...In this paper, we first introduce interpolation operator of projection type in three dimen- sions, from which we derive weak estimates for tricubic block finite elements. Then using the estimate for the W 2, 1-seminorm of the discrete derivative Green's function and the weak estimates, we show that the tricubic block finite element solution uh and the tricubic interpolant of projection type Πh3u have superclose gradient in the pointwise sense of the L∞-norm. Finally, this supercloseness is applied to superconvergence analysis, and the global superconvergence of the finite element approximation is derived.展开更多
For rectangular finite element, we give a superconvergence method by SPR technique based on the generalization of a new ultraconvergence record and the sharp Green function estimates, by which we prove that the deriva...For rectangular finite element, we give a superconvergence method by SPR technique based on the generalization of a new ultraconvergence record and the sharp Green function estimates, by which we prove that the derivative has ultra-convergence of order O(hk+3) (k 3 being odd) and displacement has order of O(hk+4) (k 4 being even) at the locally symmetry points.展开更多
This paper generalizes two nonconforming rectangular elements of the Reissner-Mindlin plate to the quadrilateral mesh. The first quadrilateral element uses the usual conforming bilinear element to approximate both com...This paper generalizes two nonconforming rectangular elements of the Reissner-Mindlin plate to the quadrilateral mesh. The first quadrilateral element uses the usual conforming bilinear element to approximate both components of the rotation, and the modified nonconforming rotated Q1 element enriched with the intersected term on each element to approximate the displacement, whereas the second one uses the enriched modified nonconforming rotated Q1 element to approximate both the rotation and the displacement. Both elements employ a more complicated shear force space to overcome the shear force locking, which will be described in detail in the introduction. We prove that both methods converge at optimal rates uniformly in the plate thickness t and the mesh distortion parameter in both the H1-and the L2-norms, and consequently they are locking free.展开更多
The main aim of this paper is to study the convergence properties of a low order mixed finite element for the Stokes problem under anisotropic meshes. We discuss the anisotropic convergence and superconvergence indepe...The main aim of this paper is to study the convergence properties of a low order mixed finite element for the Stokes problem under anisotropic meshes. We discuss the anisotropic convergence and superconvergence independent of the aspect ratio. Without the shape regularity assumption and inverse assumption on the meshes, the optimal error estimates and natural superconvergence at central points are obtained. The global superconvergence for the gradient of the velocity and the pressure is derived with the aid of a suitable postprocessing method. Furthermore, we develop a simple method to obtain the superclose properties which improves the results of the previous works .展开更多
文摘This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.
文摘The nonoverlapping domain decomposition method for parabolic partialdifferential equation on general domain is considered. A kind of domain decompositionthat uses the finite element procedure is given. The problem over the domains can beimplemented on parallel computer. Convergence analysis is also presented.
文摘Abstract Two phase,incompressible miscible flow in porous media is governed by a system of nonlinear partial differential equations.The pressure equation,which is elliptic in appearance,is discretized by a standard five points difference method.The concentration equation is treated by an implicit finite difference method that applies a form of the method of characteristics to the transport terms.A class of biquadratic interpolation is introduced for the method of characteristics.Convergence rate is proved to be O(Δt+h 2).
基金Supported partially by the National Natural Science Foundation of China( 1 0 0 71 0 5 0 ) ScienceFoundation ( 0 2 ZA1 4 0 70 ) of Shanghai Technical Sciences Committee and Science Foundation ( 0 2 DK0 6) ofShanghai Education Committee.
文摘A interior point scaling projected reduced Hessian method with combination of nonmonotonic backtracking technique and trust region strategy for nonlinear equality constrained optimization with nonegative constraint on variables is proposed. In order to deal with large problems, a pair of trust region subproblems in horizontal and vertical subspaces is used to replace the general full trust region subproblem. The horizontal trust region subproblem in the algorithm is only a general trust region subproblem while the vertical trust region subproblem is defined by a parameter size of the vertical direction subject only to an ellipsoidal constraint. Both trust region strategy and line search technique at each iteration switch to obtaining a backtracking step generated by the two trust region subproblems. By adopting the l 1 penalty function as the merit function, the global convergence and fast local convergence rate of the proposed algorithm are established under some reasonable conditions. A nonmonotonic criterion and the second order correction step are used to overcome Maratos effect and speed up the convergence progress in some ill-conditioned cases. MR Subject Classification 90C30 - 65K05 - 49M40 Keywords trust region method - backtracking step - reduced Hessian - nonmonotonic technique - interior point Supported partially by the National Natural Science Foundation of China (10071050), Science Foundation (02ZA14070) of Shanghai Technical Sciences Committee and Science Foundation (02DK06) of Shanghai Education Committee.
基金the National Natural Science Foundation of China(Grant Nos.10771063,10571053)Doctoral Programme of National Education Ministry of China(Grant No.20050542006)Programme for New Century Excellent Talents in University(Grant No.NCET-060712)
文摘Based on an asymptotic expansion of finite element,a new extrapolation formula and extrapolation cascadic multigrid method(EXCMG)are proposed,in which the new extrapolation and quadratic interpolation are used to provide a better initial value on refined grid.In the case of triple grids,the error of the new initial value is analyzed in detail.A larger scale computation is completed in PC.
基金supported by the National Natural Science Foundation of China (Grant No. 10761003)Guizhou Province Scientific Research for Senior Personnels
文摘This paper discusses the order-preserving convergence for spectral approximation of the self-adjoint completely continuous operator T.Under the condition that the approximate operator Th converges to T in norm,it is proven that the k-th eigenvalue of Th converges to the k-th eigenvalue of T.(We sorted the positive eigenvalues in decreasing order and negative eigenvalues in increasing order.) Then we apply this result to conforming elements,nonconforming elements and mixed elements of self-adjoint elliptic differential operators eigenvalue problems,and prove that the k-th approximate eigenvalue obtained by these methods converges to the k-th exact eigenvalue.
基金supported by National Natural Science Foundation of China (Grant Nos. 10871022, 10771065,and 60573158)Natural Science Foundation of Hebei Province (Grant No. A2007001027)
文摘A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.
基金the National Science Foundation(Grant Nos.DMS0409297,DMR0205232,CCF-0430349)US National Institute of Health-National Cancer Institute(Grant No.1R01CA125707-01A1)+2 种基金the National Natural Science Foundation of China(Grant No.10571172)the National Basic Research Program(Grant No.2005CB321704)the Youth's Innovative Program of Chinese Academy of Sciences(Grant Nos.K7290312G9,K7502712F9)
文摘In this paper,we consider the cascadic multigrid method for a parabolic type equation.Backward Euler approximation in time and linear finite element approximation in space are employed.A stability result is established under some conditions on the smoother.Using new and sharper estimates for the smoothers that reflect the precise dependence on the time step and the spatial mesh parameter,these conditions are verified for a number of popular smoothers.Optimal error bound sare derived for both smooth and non-smooth data.Iteration strategies guaranteeing both the optimal accuracy and the optimal complexity are presented.
基金the National Natural Science Foundation of China(Grant Nos.10471100,40437017,and 60573158)Beijing Jiaotong University Science and Technology Foundation
文摘The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.
基金the National Natural Science Foundation of China(Grant No.10771150)the National Basic Research Program of China(Grant No.2005CB321701)the Program for New Century Excellent Talents in University(Grant No.NCET-07-0584)
文摘Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.
基金the National Nutural Science Foundation of China(Grant Nos.10771198,10590353)
文摘Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolation error of anisotropic elements based on Newton's formula of polynomial interpolation as well as its applications.
基金This work was supported by the National Major Basic Research Projects (Grant No. G1999032804)the National Natural Science Foundation of China (Grant No.10471038, 10571053)+1 种基金the Research Fonds for Doctor Programme (Grant No. 20050542006)Programme for New Century Excellent Talent in University (GrantNo. NCET-06-0717)
文摘For numerical computations of multiple solutions of the nonlinear elliptic problem Δu + f(u = 0 in Ω, u = 0 on Γ, a search-extension method (SEM) was proposed and systematically studied by the authors. This paper shall complete its theoretical analysis. It is assumed that the nonlinearity is non-convex and its solution is isolated, under some conditions the corresponding linearized problem has a unique solution. By use of the compactness of the solution family and the contradiction argument, in general conditions, the high order regularity of the solution u ∈ H 1+α, α > 0 is proved. Assume that some initial value searched by suitably many eigenbases is already fallen into the neighborhood of the isolated solution, then the optimal error estimates of its nonlinear finite element approximation are shown by the duality argument and continuation method.
基金supported by National Natural Science Foundation of China (Grant No.10761003) the Governor's Special Foundation of Guizhou Province for Outstanding Scientific Education Personnel (Grant No.[2005]155)
文摘This study discusses generalized Rayleigh quotient and high efficiency finite element discretization schemes. Some results are as follows: 1) Rayleigh quotient accelerate technique is extended to nonselfadjoint problems. Generalized Rayleigh quotients of operator form and weak form are defined and the basic relationship between approximate eigenfunction and its generalized Rayleigh quotient is established. 2) New error estimates are obtained by replacing the ascent of exact eigenvalue with the ascent of finite element approximate eigenvalue. 3) Based on the work of Xu Jinchao and Zhou Aihui, finite element two-grid discretization schemes are established to solve nonselfadjoint elliptic differential operator eigenvalue problems and these schemes are used in both conforming finite element and non-conforming finite element. Besides, the efficiency of the schemes is proved by both theoretical analysis and numerical experiments. 4) Iterated Galerkin method, interpolated correction method and gradient recovery for selfadjoint elliptic differential operator eigenvalue problems are extended to nonselfadjoint elliptic differential operator eigenvalue problems.
基金Project supported by the National Natural Science Foundation of China.
文摘Consider an AMG for the linear system Au=f. Up to now, only the uniform convergence of two-level AMG is proved for symmetric and positive definite L-matrices with weak diagonal dominance. Using the new form (1), we extend the results in [1] to the case that A is a general symmetric and positive definite matrix with weak diagonal dominance. In the following, we shall use the same notations as in [1].
基金supported by Natural Science Foundation of Ningbo City (Grant No. 2008A610020)National Natural Science Foundation of China (Grant No. 10671065)the Scientific Research Fund of Hunan Provincial Education Department (Grant Nos. 07C576, 03C212)
文摘In this paper, we first introduce interpolation operator of projection type in three dimen- sions, from which we derive weak estimates for tricubic block finite elements. Then using the estimate for the W 2, 1-seminorm of the discrete derivative Green's function and the weak estimates, we show that the tricubic block finite element solution uh and the tricubic interpolant of projection type Πh3u have superclose gradient in the pointwise sense of the L∞-norm. Finally, this supercloseness is applied to superconvergence analysis, and the global superconvergence of the finite element approximation is derived.
文摘For rectangular finite element, we give a superconvergence method by SPR technique based on the generalization of a new ultraconvergence record and the sharp Green function estimates, by which we prove that the derivative has ultra-convergence of order O(hk+3) (k 3 being odd) and displacement has order of O(hk+4) (k 4 being even) at the locally symmetry points.
基金the National Natural Science Foundation of China (Grant No. 10601003)National Excellent Doctoral Dissertation of China (Grant No. 200718)
文摘This paper generalizes two nonconforming rectangular elements of the Reissner-Mindlin plate to the quadrilateral mesh. The first quadrilateral element uses the usual conforming bilinear element to approximate both components of the rotation, and the modified nonconforming rotated Q1 element enriched with the intersected term on each element to approximate the displacement, whereas the second one uses the enriched modified nonconforming rotated Q1 element to approximate both the rotation and the displacement. Both elements employ a more complicated shear force space to overcome the shear force locking, which will be described in detail in the introduction. We prove that both methods converge at optimal rates uniformly in the plate thickness t and the mesh distortion parameter in both the H1-and the L2-norms, and consequently they are locking free.
基金the National Natural Science Foundation of China under the grant 10771198
文摘The main aim of this paper is to study the convergence properties of a low order mixed finite element for the Stokes problem under anisotropic meshes. We discuss the anisotropic convergence and superconvergence independent of the aspect ratio. Without the shape regularity assumption and inverse assumption on the meshes, the optimal error estimates and natural superconvergence at central points are obtained. The global superconvergence for the gradient of the velocity and the pressure is derived with the aid of a suitable postprocessing method. Furthermore, we develop a simple method to obtain the superclose properties which improves the results of the previous works .
