The energy approach is used to theoretically verify that the average acceleration method (AAM), which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical ...The energy approach is used to theoretically verify that the average acceleration method (AAM), which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.展开更多
An explicit unconditionally stable algorithm for hybrid tests,which is developed from the traditional HHT-α algorithm,is proposed.The unconditional stability is first proven by the spectral radius method for a linear...An explicit unconditionally stable algorithm for hybrid tests,which is developed from the traditional HHT-α algorithm,is proposed.The unconditional stability is first proven by the spectral radius method for a linear system.If the value of α is selected within [-0.5,0],then the algorithm is shown to be unconditionally stable.Next,the root locus method for a discrete dynamic system is applied to analyze the stability of a nonlinear system.The results show that the proposed method is conditionally stable for dynamic systems with stiffness hardening.To improve the stability of the proposed method,the structure stiffness is then identified and updated.Both numerical and pseudo-dynamic tests on a structure with the collision effect prove that the stiffness updating method can effectively improve stability.展开更多
By means of the frequency domain method and the inequality analysis, we discuss the unconditional stability problem for the hyperneutral type constant linear control system with delays, and obtain some precise suffici...By means of the frequency domain method and the inequality analysis, we discuss the unconditional stability problem for the hyperneutral type constant linear control system with delays, and obtain some precise sufficient, sufficient and necessary conditions.展开更多
In this paper, by means of the frequency domain method and the inequality analysis, unconditional stability problem for the hyperneutral type constant linear control system with delays are discussed, and some precise ...In this paper, by means of the frequency domain method and the inequality analysis, unconditional stability problem for the hyperneutral type constant linear control system with delays are discussed, and some precise sufficient, sufficient and necessary conditions are obtained.展开更多
Applying the frequency domain method and the inequality method, we discussed the unconditional stability problem of the multigroup multidelays neutral type linear constant continuous control system, and obtained some ...Applying the frequency domain method and the inequality method, we discussed the unconditional stability problem of the multigroup multidelays neutral type linear constant continuous control system, and obtained some sufficient conditions.展开更多
In this paper, the sufficient and necessary conditions of the unconditional stability, and the delay bound of the third-order neutral delay differential equation with real constant coefficients are given. The conditio...In this paper, the sufficient and necessary conditions of the unconditional stability, and the delay bound of the third-order neutral delay differential equation with real constant coefficients are given. The conditions are brief and practical algebraic criterions Furthermore, we get the delay bound.展开更多
The ABE-I (Alternating Block Explicit-lmplicit) method for diffusion problem is extended to solve the variable coefficient problem and the unconditional stability of the ABE-I method is proved by the energy method.
A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a...A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a structure does not need to be considered. This is because this subfamily is unconditionally stable for any instantaneous stiffness softening system, linear elastic system and instantaneous stiffness hardening system that might occur in the pseudodynamic testing of a real structure. In addition, it also offers good accuracy when compared to a general second-order accurate method for both linear elastic and nonlinear systems.展开更多
Two explicit integration algorithms with unconditional stability for linear elastic systems have been successfully developed for pseudodynamic testing. Their numerical properties in the solution of a linear elastic sy...Two explicit integration algorithms with unconditional stability for linear elastic systems have been successfully developed for pseudodynamic testing. Their numerical properties in the solution of a linear elastic system have been well explored and their applications to the pseudodynamic testing of a nonlinear system have been shown to be feasible. However, their numerical properties in the solution of a nonlinear system are not apparent. Therefore, the performance of both algorithms for use in the solution of a nonlinear system has been analytically evaluated after introducing an instantaneous degree of nonlinearity. The two algorithms have roughly the same accuracy for a small value of the product of the natural frequency and step size. Meanwhile, the first algorithm is unconditionally stable when the instantaneous degree of nonlinearity is less than or equal to 1, and it becomes conditionally stable when it is greater than 1. The second algorithm is conditionally stable as the instantaneous degree of nonlinearity is less than 1/9, and becomes unstable when it is greater than 1. It can have unconditional stability for the range between 1/9 and 1. Based on these evaluations, it was concluded that the first algorithm is superior to the second one. Also, both algorithms were found to require commensurate computational efforts, which are much less than needed for the Newmark explicit method in general structural dynamic problems.展开更多
Some new Saul'yev type asymmetric difference schemes for Burgers' equation is given, by the use of the schemes, a kind of alternating group four points method for solving nonlinear Burgers' equation is con...Some new Saul'yev type asymmetric difference schemes for Burgers' equation is given, by the use of the schemes, a kind of alternating group four points method for solving nonlinear Burgers' equation is constructed here. The basic idea of the method is that the grid points on the same time level is divided into a number of groups, the difference equations of each group can be solved independently, hence the method with intrinsic parallelism can be used directly on parallel computer. The method is unconditionally stable by analysis of linearization procedure. The numerical experiments show that the method has good stability and accuracy.展开更多
Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in trac...Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.展开更多
In this paper, a new alternating group explicit-implicit (nAGEI) scheme for dispersive equations with a periodic boundary condition is derived. This new unconditionally stable scheme has a fourth-order truncation er...In this paper, a new alternating group explicit-implicit (nAGEI) scheme for dispersive equations with a periodic boundary condition is derived. This new unconditionally stable scheme has a fourth-order truncation error in space and a convergence ratio faster than some known alternating methods such as ASEI and AGE. Comparison in accuracy with the AGEI and AGE methods is presented in the numerical experiment.展开更多
The main aim of this paper is to analyze the numerical method based upon the spectral element technique for the numerical solution of the fractional advection-diffusion equa-tion.The time variable has been discretized...The main aim of this paper is to analyze the numerical method based upon the spectral element technique for the numerical solution of the fractional advection-diffusion equa-tion.The time variable has been discretized by a second-order finite difference procedure.The stability and the convergence of the semi-discrete formula have been proven.Then,the spatial variable of the main PDEs is approximated by the spectral element method.The convergence order of the fully discrete scheme is studied.The basis functions of the spectral element method are based upon a class of Legendre polynomials.The numerical experiments confirm the theoretical results.展开更多
In this paper, we propose a fully decoupled and linear scheme for the magnetohydrodynamic (MHD) equation with the backward differential formulation (BDF) and finite element method (FEM). To solve the system, we adopt ...In this paper, we propose a fully decoupled and linear scheme for the magnetohydrodynamic (MHD) equation with the backward differential formulation (BDF) and finite element method (FEM). To solve the system, we adopt a technique based on the “zero-energy-contribution” contribution, which separates the magnetic and fluid fields from the coupled system. Additionally, making use of the pressure projection methods, the pressure variable appears explicitly in the velocity field equation, and would be computed in the form of a Poisson equation. Therefore, the total system is divided into several smaller sub-systems that could be simulated at a significantly low cost. We prove the unconditional energy stability, unique solvability and optimal error estimates for the proposed scheme, and present numerical results to verify the accuracy, efficiency and stability of the scheme.展开更多
For solving nonlinear parabolic equation on massive parallel computers, the construction of parallel difference schemes with simple design, high parallelism and unconditional stability and second order global accuracy...For solving nonlinear parabolic equation on massive parallel computers, the construction of parallel difference schemes with simple design, high parallelism and unconditional stability and second order global accuracy in space, has long been desired. In the present work, a new kind of general parallel difference schemes for the nonlinear parabolic system is proposed. The general parallel difference schemes include, among others, two new parallel schemes. In one of them, to obtain the interface values on the interface of sub-domains an explicit scheme of Jacobian type is employed, and then the fully implicit scheme is used in the sub-domains. Here, in the explicit scheme of Jacobian type, the values at the points being adjacent to the interface points are taken as the linear combination of values of previous two time layers at the adjoining points of the inner interface. For the construction of another new parallel difference scheme, the main procedure is as follows. Firstly the linear combination of values of previous two time layers at the interface points among the sub-domains is used as the (Dirichlet) boundary condition for solving the sub-domain problems. Then the values in the sub-domains are calculated by the fully implicit scheme. Finally the interface values are computed by the fully implicit scheme, and in fact these calculations of the last step are explicit since the values adjacent to the interface points have been obtained in the previous step. The existence, uniqueness, unconditional stability and the second order accuracy of the discrete vector solutions for the parallel difference schemes are proved. Numerical results are presented to examine the stability, accuracy and parallelism of the parallel schemes.展开更多
Proposes an explicit fully discrete three-level pseudo-spectral scheme with unconditional stability for the Cahn-Hilliard equation. Equations for pseudo-spectral scheme; Analysis of linear stability of critical points.
The paper is concerned with the unconditional stability and error estimates of fully discrete Galerkin-Galerkin FEMs for the equations of incompressible miscible flows in porous media.We prove that the optimal L 2 err...The paper is concerned with the unconditional stability and error estimates of fully discrete Galerkin-Galerkin FEMs for the equations of incompressible miscible flows in porous media.We prove that the optimal L 2 error estimates hold without any time-step(convergence)conditions,while all previous works require certain time-step restrictions.Theoretical analysis is based on a splitting of the error into two parts:the error from the time discretization of the PDEs and the error from the finite element discretization of the corresponding time-discrete PDEs,which was proposed in our previous work[26,27].Numerical results for both two and three-dimensional flow models are presented to confirm our theoretical analysis.展开更多
In this paper, we are interested in the multigroup multidelays interval coefficient constant and time varying linear continuous control systems, by means of the equivalence method and the differential inequality in t...In this paper, we are interested in the multigroup multidelays interval coefficient constant and time varying linear continuous control systems, by means of the equivalence method and the differential inequality in the time domain, we obtain some unconditional robust stability results for the multigroup multidelays constant and time-varying interval coefficient linear continuous control systems, respectively.展开更多
Thermal phase change problems are widespread in mathematics,nature,and science.They are particularly useful in simulating the phenomena of melting and solidification in materials science.In this paper we propose a nov...Thermal phase change problems are widespread in mathematics,nature,and science.They are particularly useful in simulating the phenomena of melting and solidification in materials science.In this paper we propose a novel class of arbitrarily high-order and unconditionally energy stable schemes for a thermal phase changemodel,which is the coupling of a heat transfer equation and a phase field equation.The unconditional energy stability and consistency error estimates are rigorously proved for the proposed schemes.A detailed implementation demonstrates that the proposed method requires only the solution of a system of linear elliptic equations at each time step,with an efficient scheme of sufficient accuracy to calculate the solution at the first step.It is observed from the comparison with the classical explicit Runge-Kutta method that the new schemes allow to use larger time steps.Adaptive time step size strategies can be applied to further benefit from this unconditional stability.Numerical experiments are presented to verify the theoretical claims and to illustrate the accuracy and effectiveness of our method.展开更多
基金National Natural Science Foundation of ChinaUnder Grant No. 50578047, 50338020 China Ministry ofEducation (Program for New Century Excellent Talents inUniversity) China Ministry of Science and Technology UnderGrant No.2003AA602150
文摘The energy approach is used to theoretically verify that the average acceleration method (AAM), which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.
基金Scientific Research Fund of the Institute of Engineering Mechanics,CEA under Grant Nos.2017A02,2016B09 and 2016A06the National Science-technology Support Plan Projects under Grant No.2015BAK17B02the National Natural Science Foundation of China under Grant Nos.51378478,51408565,51678538 and 51161120360
文摘An explicit unconditionally stable algorithm for hybrid tests,which is developed from the traditional HHT-α algorithm,is proposed.The unconditional stability is first proven by the spectral radius method for a linear system.If the value of α is selected within [-0.5,0],then the algorithm is shown to be unconditionally stable.Next,the root locus method for a discrete dynamic system is applied to analyze the stability of a nonlinear system.The results show that the proposed method is conditionally stable for dynamic systems with stiffness hardening.To improve the stability of the proposed method,the structure stiffness is then identified and updated.Both numerical and pseudo-dynamic tests on a structure with the collision effect prove that the stiffness updating method can effectively improve stability.
基金Supported by the Natural Science Foundation of Hubei Province ( 2 0 0 0 A490 0 5 )
文摘By means of the frequency domain method and the inequality analysis, we discuss the unconditional stability problem for the hyperneutral type constant linear control system with delays, and obtain some precise sufficient, sufficient and necessary conditions.
基金Supported by the natural science Foundation of Hubei Provincec Education Committee
文摘In this paper, by means of the frequency domain method and the inequality analysis, unconditional stability problem for the hyperneutral type constant linear control system with delays are discussed, and some precise sufficient, sufficient and necessary conditions are obtained.
文摘Applying the frequency domain method and the inequality method, we discussed the unconditional stability problem of the multigroup multidelays neutral type linear constant continuous control system, and obtained some sufficient conditions.
文摘In this paper, the sufficient and necessary conditions of the unconditional stability, and the delay bound of the third-order neutral delay differential equation with real constant coefficients are given. The conditions are brief and practical algebraic criterions Furthermore, we get the delay bound.
基金Project supported by Special Funds for Major State Basic Research Projects (G1999032802) in China and NNSF (10076006) of China.
文摘The ABE-I (Alternating Block Explicit-lmplicit) method for diffusion problem is extended to solve the variable coefficient problem and the unconditional stability of the ABE-I method is proved by the energy method.
基金Science Council, Chinese Taipei Under Grant No. NSC-95-2221-E-027-099
文摘A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a structure does not need to be considered. This is because this subfamily is unconditionally stable for any instantaneous stiffness softening system, linear elastic system and instantaneous stiffness hardening system that might occur in the pseudodynamic testing of a real structure. In addition, it also offers good accuracy when compared to a general second-order accurate method for both linear elastic and nonlinear systems.
基金Science Council,Chinese Taipei,Under Grant No. NSC-96-2211-E-027-030
文摘Two explicit integration algorithms with unconditional stability for linear elastic systems have been successfully developed for pseudodynamic testing. Their numerical properties in the solution of a linear elastic system have been well explored and their applications to the pseudodynamic testing of a nonlinear system have been shown to be feasible. However, their numerical properties in the solution of a nonlinear system are not apparent. Therefore, the performance of both algorithms for use in the solution of a nonlinear system has been analytically evaluated after introducing an instantaneous degree of nonlinearity. The two algorithms have roughly the same accuracy for a small value of the product of the natural frequency and step size. Meanwhile, the first algorithm is unconditionally stable when the instantaneous degree of nonlinearity is less than or equal to 1, and it becomes conditionally stable when it is greater than 1. The second algorithm is conditionally stable as the instantaneous degree of nonlinearity is less than 1/9, and becomes unstable when it is greater than 1. It can have unconditional stability for the range between 1/9 and 1. Based on these evaluations, it was concluded that the first algorithm is superior to the second one. Also, both algorithms were found to require commensurate computational efforts, which are much less than needed for the Newmark explicit method in general structural dynamic problems.
文摘Some new Saul'yev type asymmetric difference schemes for Burgers' equation is given, by the use of the schemes, a kind of alternating group four points method for solving nonlinear Burgers' equation is constructed here. The basic idea of the method is that the grid points on the same time level is divided into a number of groups, the difference equations of each group can be solved independently, hence the method with intrinsic parallelism can be used directly on parallel computer. The method is unconditionally stable by analysis of linearization procedure. The numerical experiments show that the method has good stability and accuracy.
基金supported by the National Natural Science Foundation of China(Grant Numbers 11872090,11672019,11472035).
文摘Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.
基金National Natural Science Foundation of China (No.10671113)
文摘In this paper, a new alternating group explicit-implicit (nAGEI) scheme for dispersive equations with a periodic boundary condition is derived. This new unconditionally stable scheme has a fourth-order truncation error in space and a convergence ratio faster than some known alternating methods such as ASEI and AGE. Comparison in accuracy with the AGEI and AGE methods is presented in the numerical experiment.
基金The authors are grateful to the two reviewers for carefully reading this paper and for their comments and suggestions which have highly improved the paper.
文摘The main aim of this paper is to analyze the numerical method based upon the spectral element technique for the numerical solution of the fractional advection-diffusion equa-tion.The time variable has been discretized by a second-order finite difference procedure.The stability and the convergence of the semi-discrete formula have been proven.Then,the spatial variable of the main PDEs is approximated by the spectral element method.The convergence order of the fully discrete scheme is studied.The basis functions of the spectral element method are based upon a class of Legendre polynomials.The numerical experiments confirm the theoretical results.
文摘In this paper, we propose a fully decoupled and linear scheme for the magnetohydrodynamic (MHD) equation with the backward differential formulation (BDF) and finite element method (FEM). To solve the system, we adopt a technique based on the “zero-energy-contribution” contribution, which separates the magnetic and fluid fields from the coupled system. Additionally, making use of the pressure projection methods, the pressure variable appears explicitly in the velocity field equation, and would be computed in the form of a Poisson equation. Therefore, the total system is divided into several smaller sub-systems that could be simulated at a significantly low cost. We prove the unconditional energy stability, unique solvability and optimal error estimates for the proposed scheme, and present numerical results to verify the accuracy, efficiency and stability of the scheme.
基金The project is supported by the Special Funds for Major State Basic Research Projects 2005CB321703, the National Nature Science Foundation of China (No. 10476002, 60533020).
文摘For solving nonlinear parabolic equation on massive parallel computers, the construction of parallel difference schemes with simple design, high parallelism and unconditional stability and second order global accuracy in space, has long been desired. In the present work, a new kind of general parallel difference schemes for the nonlinear parabolic system is proposed. The general parallel difference schemes include, among others, two new parallel schemes. In one of them, to obtain the interface values on the interface of sub-domains an explicit scheme of Jacobian type is employed, and then the fully implicit scheme is used in the sub-domains. Here, in the explicit scheme of Jacobian type, the values at the points being adjacent to the interface points are taken as the linear combination of values of previous two time layers at the adjoining points of the inner interface. For the construction of another new parallel difference scheme, the main procedure is as follows. Firstly the linear combination of values of previous two time layers at the interface points among the sub-domains is used as the (Dirichlet) boundary condition for solving the sub-domain problems. Then the values in the sub-domains are calculated by the fully implicit scheme. Finally the interface values are computed by the fully implicit scheme, and in fact these calculations of the last step are explicit since the values adjacent to the interface points have been obtained in the previous step. The existence, uniqueness, unconditional stability and the second order accuracy of the discrete vector solutions for the parallel difference schemes are proved. Numerical results are presented to examine the stability, accuracy and parallelism of the parallel schemes.
文摘Proposes an explicit fully discrete three-level pseudo-spectral scheme with unconditional stability for the Cahn-Hilliard equation. Equations for pseudo-spectral scheme; Analysis of linear stability of critical points.
基金supported in part by a grant from National Science Foundation(Project No.11301262)a grant from the Research Grants Council of the Hong Kong Special Administrative Region,China(Project No.CityU 102613)The work of J.Wang and W.Sun was supported in part by a grant from the Research Grants Council of the Hong Kong Special Administrative Region,China(Project No.CityU 102613).
文摘The paper is concerned with the unconditional stability and error estimates of fully discrete Galerkin-Galerkin FEMs for the equations of incompressible miscible flows in porous media.We prove that the optimal L 2 error estimates hold without any time-step(convergence)conditions,while all previous works require certain time-step restrictions.Theoretical analysis is based on a splitting of the error into two parts:the error from the time discretization of the PDEs and the error from the finite element discretization of the corresponding time-discrete PDEs,which was proposed in our previous work[26,27].Numerical results for both two and three-dimensional flow models are presented to confirm our theoretical analysis.
基金the Natural Science Foundation of Hubei Province Education Committee.
文摘In this paper, we are interested in the multigroup multidelays interval coefficient constant and time varying linear continuous control systems, by means of the equivalence method and the differential inequality in the time domain, we obtain some unconditional robust stability results for the multigroup multidelays constant and time-varying interval coefficient linear continuous control systems, respectively.
文摘Thermal phase change problems are widespread in mathematics,nature,and science.They are particularly useful in simulating the phenomena of melting and solidification in materials science.In this paper we propose a novel class of arbitrarily high-order and unconditionally energy stable schemes for a thermal phase changemodel,which is the coupling of a heat transfer equation and a phase field equation.The unconditional energy stability and consistency error estimates are rigorously proved for the proposed schemes.A detailed implementation demonstrates that the proposed method requires only the solution of a system of linear elliptic equations at each time step,with an efficient scheme of sufficient accuracy to calculate the solution at the first step.It is observed from the comparison with the classical explicit Runge-Kutta method that the new schemes allow to use larger time steps.Adaptive time step size strategies can be applied to further benefit from this unconditional stability.Numerical experiments are presented to verify the theoretical claims and to illustrate the accuracy and effectiveness of our method.