在乙烯生产流程中,以脱丙烷塔与脱丁烷塔的复杂联塔系统作为研究对象,根据机理建立双塔联立的动态数学模型。该联塔数学模型是一个大规模的微分代数混合系统(DAEs),包括了精馏塔、控制器、阀门、泵等多种操作单元和设备。动态模拟实质...在乙烯生产流程中,以脱丙烷塔与脱丁烷塔的复杂联塔系统作为研究对象,根据机理建立双塔联立的动态数学模型。该联塔数学模型是一个大规模的微分代数混合系统(DAEs),包括了精馏塔、控制器、阀门、泵等多种操作单元和设备。动态模拟实质即为求解这个大规模的DAE系统。文章以动态模拟软件gPROMS为建模工具和模拟环境,通过调用大规模DAE求解算法,模拟联塔的动态操作过程,并与Aspen P lus稳态模拟结果进行了比较和分析。展开更多
We investigate an important relationship that exists between the Hopf bifurcation in the singularly perturbed nonlinear power systems and the singularity induced bifurcations (SIBs) in the corresponding different- tia...We investigate an important relationship that exists between the Hopf bifurcation in the singularly perturbed nonlinear power systems and the singularity induced bifurcations (SIBs) in the corresponding different- tial-algebraic equations (DAEs). In a generic case, the SIB phenomenon in a system of DAEs signals Hopf bifurcation in the singularly perturbed systems of ODEs. The analysis is based on the linear matrix pencil theory and polynomials with parameter dependent coefficients. A few numerical examples are included.展开更多
In this paper,we study the dynamics of an idealized benchmark bicycle moving on a surface of revolution.We employ symbolic manipulations to derive the contact constraint equations from an ordered process,and apply the...In this paper,we study the dynamics of an idealized benchmark bicycle moving on a surface of revolution.We employ symbolic manipulations to derive the contact constraint equations from an ordered process,and apply the Lagrangian equations of the first type to establish the nonlinear differential algebraic equations(DAEs),leaving nine coupled differential equations,six contact equations,two holonomic constraint equations and four nonholonomic constraint equations.We then present a complete description of hands-free circular motions,in which the time-dependent variables are eliminated through a rotation transformation.We find that the circular motions,similar to those of the bicycle moving on a horizontal surface,nominally fall into four solution families,characterized by four curves varying with the angular speed of the front wheel.Then,we numerically investigate how the topological profiles of these curves change with the parameter of the revolution surface.Furthermore,we directly linearize the nonlinear DAEs,from which a reduced linearized system is obtained by removing the dependent coordinates and counting the symmetries arising from cyclic coordinates.The stability of the circular motion is then analyzed according to the eigenvalues of the Jacobian matrix of the reduced linearized system around the equilibrium position.We find that a stable circular motion exists only if the curvature of the revolution surface is very small and it is limited in small sections of solution families.Finally,based on the numerical simulation of the original nonlinear DAEs system,we show that the stable circular motion is not asymptotically stable.展开更多
In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular inde...In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.展开更多
Studies the different types of multistep discretization of index 3 differential-algebraic equations in Hessenberg form. Existense, uniqueness and influence of perturbations; Local convergence of multistep discretizati...Studies the different types of multistep discretization of index 3 differential-algebraic equations in Hessenberg form. Existense, uniqueness and influence of perturbations; Local convergence of multistep discretization; Details on the numerical tests.展开更多
The simulation of a high-temperature gas-cooled reactor pebble-bed module(HTR-PM) plant is discussed.This lumped parameter model has the form of a set differential algebraic equations(DAEs) that include stiff equation...The simulation of a high-temperature gas-cooled reactor pebble-bed module(HTR-PM) plant is discussed.This lumped parameter model has the form of a set differential algebraic equations(DAEs) that include stiff equations to model point neutron kinetics.The nested approach is the most common method to solve DAE,but this approach is very expensive and time-consuming due to inner iterations.This paper deals with an alternative approach in which a simultaneous solution method is used.The DAEs are discretized over a time horizon using collocation on finite elements,and Radau collocation points are applied.The resulting nonlinear algebraic equations can be solved by existing solvers.The discrete algorithm is discussed in detail;both accuracy and stability issues are considered.Finally,the simulation results are presented to validate the efficiency and accuracy of the simultaneous approach that takes much less time than the nested one.展开更多
文摘在乙烯生产流程中,以脱丙烷塔与脱丁烷塔的复杂联塔系统作为研究对象,根据机理建立双塔联立的动态数学模型。该联塔数学模型是一个大规模的微分代数混合系统(DAEs),包括了精馏塔、控制器、阀门、泵等多种操作单元和设备。动态模拟实质即为求解这个大规模的DAE系统。文章以动态模拟软件gPROMS为建模工具和模拟环境,通过调用大规模DAE求解算法,模拟联塔的动态操作过程,并与Aspen P lus稳态模拟结果进行了比较和分析。
文摘We investigate an important relationship that exists between the Hopf bifurcation in the singularly perturbed nonlinear power systems and the singularity induced bifurcations (SIBs) in the corresponding different- tial-algebraic equations (DAEs). In a generic case, the SIB phenomenon in a system of DAEs signals Hopf bifurcation in the singularly perturbed systems of ODEs. The analysis is based on the linear matrix pencil theory and polynomials with parameter dependent coefficients. A few numerical examples are included.
基金National Natural Science Foundation of China(Grants 11932001 and 11702002).
文摘In this paper,we study the dynamics of an idealized benchmark bicycle moving on a surface of revolution.We employ symbolic manipulations to derive the contact constraint equations from an ordered process,and apply the Lagrangian equations of the first type to establish the nonlinear differential algebraic equations(DAEs),leaving nine coupled differential equations,six contact equations,two holonomic constraint equations and four nonholonomic constraint equations.We then present a complete description of hands-free circular motions,in which the time-dependent variables are eliminated through a rotation transformation.We find that the circular motions,similar to those of the bicycle moving on a horizontal surface,nominally fall into four solution families,characterized by four curves varying with the angular speed of the front wheel.Then,we numerically investigate how the topological profiles of these curves change with the parameter of the revolution surface.Furthermore,we directly linearize the nonlinear DAEs,from which a reduced linearized system is obtained by removing the dependent coordinates and counting the symmetries arising from cyclic coordinates.The stability of the circular motion is then analyzed according to the eigenvalues of the Jacobian matrix of the reduced linearized system around the equilibrium position.We find that a stable circular motion exists only if the curvature of the revolution surface is very small and it is limited in small sections of solution families.Finally,based on the numerical simulation of the original nonlinear DAEs system,we show that the stable circular motion is not asymptotically stable.
基金Project supported by the Foundation for the Authors of the National Excellent Doctoral Thesis Award of China (200720)
文摘In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.
文摘Studies the different types of multistep discretization of index 3 differential-algebraic equations in Hessenberg form. Existense, uniqueness and influence of perturbations; Local convergence of multistep discretization; Details on the numerical tests.
基金Project supported by the National Basic Research Program of China (No. 2009CB320603)the National Natural Science Foundation of China (Nos. 60974007 and 60934007)
文摘The simulation of a high-temperature gas-cooled reactor pebble-bed module(HTR-PM) plant is discussed.This lumped parameter model has the form of a set differential algebraic equations(DAEs) that include stiff equations to model point neutron kinetics.The nested approach is the most common method to solve DAE,but this approach is very expensive and time-consuming due to inner iterations.This paper deals with an alternative approach in which a simultaneous solution method is used.The DAEs are discretized over a time horizon using collocation on finite elements,and Radau collocation points are applied.The resulting nonlinear algebraic equations can be solved by existing solvers.The discrete algorithm is discussed in detail;both accuracy and stability issues are considered.Finally,the simulation results are presented to validate the efficiency and accuracy of the simultaneous approach that takes much less time than the nested one.
