With the integration of renewable power and electric vehicle,the power system stability is of increasing concern because the active power generated by the renewable energy and absorbed by the electric vehicle vary ran...With the integration of renewable power and electric vehicle,the power system stability is of increasing concern because the active power generated by the renewable energy and absorbed by the electric vehicle vary randomly.Based on the deterministic differential equation model,the nonlinear and linear stochastic differential equation models of power system under Gauss type random excitation are proposed in this paper.The angle curves under different random excitations were simulated using Euler-Maruyama(EM) numerical method.The numerical stability of EM method was proved.The mean stability and mean square stability of the power system under Gauss type of random small excitation were verified theoretically and illustrated with simulation sample.展开更多
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.展开更多
This paper concerns the development of high-order multidimensional gas kinetic schemes for the Navier-Stokes solutions.In the current approach,the state-of-the-art WENO-type initial reconstruction and the gas-kinetic ...This paper concerns the development of high-order multidimensional gas kinetic schemes for the Navier-Stokes solutions.In the current approach,the state-of-the-art WENO-type initial reconstruction and the gas-kinetic evolution model are used in the construction of the scheme.In order to distinguish the physical and numerical requirements to recover a physical solution in a discretized space,two particle collision times will be used in the current high-order gas-kinetic scheme(GKS).Different from the low order gas dynamic model of the Riemann solution in the Godunov type schemes,the current method is based on a high-order multidimensional gas evolution model,where the space and time variation of a gas distribution function along a cell interface from an initial piecewise discontinuous polynomial is fully used in the flux evaluation.The high-order flux function becomes a unification of the upwind and central difference schemes.The current study demonstrates that both the high-order initial reconstruction and high-order gas evolution model are important in the design of a high-order numerical scheme.Especially,for a compact method,the use of a high-order local evolution solution in both space and time may become even more important,because a short stencil and local low order dynamic evolution model,i.e.,the Riemann solution,are contradictory,where valid mechanism for the update of additional degrees of freedom becomes limited.展开更多
In this paper,we consider the initial-boundary value problem(IBVP)for the micropolar Naviers-Stokes equations(MNSE)and analyze a first order fully discrete mixed finite element scheme.We first establish some regularit...In this paper,we consider the initial-boundary value problem(IBVP)for the micropolar Naviers-Stokes equations(MNSE)and analyze a first order fully discrete mixed finite element scheme.We first establish some regularity results for the solution of MNSE,which seem to be not available in the literature.Next,we study a semi-implicit time-discrete scheme for the MNSE and prove L2-H1 error estimates for the time discrete solution.Furthermore,certain regularity results for the time discrete solution are establishes rigorously.Based on these regularity results,we prove the unconditional L2-H1 error estimates for the finite element solution of MNSE.Finally,some numerical examples are carried out to demonstrate both accuracy and efficiency of the fully discrete finite element scheme.展开更多
Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous...Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous, non-isotropic matter without using (or in the absence of) the mathematical models of the BVPs and the IVPs. These methods are also used for deriving mathematical models for BVPs and IVPs associated with isotropic, homogeneous as well as non-homogeneous, non-isotropic continuous matter. In energy methods when applied to IVPs, one constructs energy functional (<i>I</i>) consisting of kinetic energy, strain energy and the potential energy of loads. The first variation of this energy functional (<em>δI</em>) set to zero is a necessary condition for an extremum of <i>I</i>. In this approach one could use <i>δI</i> = 0 directly in constructing computational processes such as the finite element method or could derive Euler’s equations (differential or partial differential equations) from <i>δI</i> = 0, which is also satisfied by a solution obtained from <i>δI</i> = 0. The Euler’s equations obtained from <i>δI</i> = 0 indeed are the mathematical model associated with the energy functional <i>I</i>. In case of BVPs we follow the same approach except in this case, the energy functional <i>I</i> consists of strain energy and the potential energy of loads. In using the principle of virtual work for BVPs and the IVPs, we can also accomplish the same as described above using energy methods. In this paper we investigate consistency and validity of the mathematical models for isotropic, homogeneous and non-isotropic, non-homogeneous continuous matter for BVPs that are derived using energy functional consisting of strain energy and the potential energy of loads. Similar investigation is also presented for IVPs using energy functional consisting of kinetic energy, strain energy and the potential energy of loads. The computational approaches for BVPs and the IVPs designed using energy functional and principl展开更多
This paper is to study extension of high resolution kinetic flux-vector splitting (KFVS) methods. In this new method, two Maxwellians are first introduced to recover the Euler equations with an additional conservative...This paper is to study extension of high resolution kinetic flux-vector splitting (KFVS) methods. In this new method, two Maxwellians are first introduced to recover the Euler equations with an additional conservative equation. Next, based on the well-known connection between the Euler equations and Boltzmann equations, a class of high resolution KFVS methods are presented to solve numerically multicomponent flows. Our method does not solve any Riemann problems, and add any nonconservative corrections. The numerical results are also presented to show the accuracy and robustness of our methods. These include one-dimensional shock tube problem, and two-dimensional interface motion in compressible flows. The computed solutions are oscillation-free near material fronts, and produce correct shock speeds.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 51137002,51190102)the Fundamental Research Funds for the Central Universities (Grant No. BZX/09B101-32)
文摘With the integration of renewable power and electric vehicle,the power system stability is of increasing concern because the active power generated by the renewable energy and absorbed by the electric vehicle vary randomly.Based on the deterministic differential equation model,the nonlinear and linear stochastic differential equation models of power system under Gauss type random excitation are proposed in this paper.The angle curves under different random excitations were simulated using Euler-Maruyama(EM) numerical method.The numerical stability of EM method was proved.The mean stability and mean square stability of the power system under Gauss type of random small excitation were verified theoretically and illustrated with simulation sample.
基金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.
基金supported by Hong Kong Research Grant Council(Grant No.621011)HKUST research fund(Grant No.SRFI11SC05)
文摘This paper concerns the development of high-order multidimensional gas kinetic schemes for the Navier-Stokes solutions.In the current approach,the state-of-the-art WENO-type initial reconstruction and the gas-kinetic evolution model are used in the construction of the scheme.In order to distinguish the physical and numerical requirements to recover a physical solution in a discretized space,two particle collision times will be used in the current high-order gas-kinetic scheme(GKS).Different from the low order gas dynamic model of the Riemann solution in the Godunov type schemes,the current method is based on a high-order multidimensional gas evolution model,where the space and time variation of a gas distribution function along a cell interface from an initial piecewise discontinuous polynomial is fully used in the flux evaluation.The high-order flux function becomes a unification of the upwind and central difference schemes.The current study demonstrates that both the high-order initial reconstruction and high-order gas evolution model are important in the design of a high-order numerical scheme.Especially,for a compact method,the use of a high-order local evolution solution in both space and time may become even more important,because a short stencil and local low order dynamic evolution model,i.e.,the Riemann solution,are contradictory,where valid mechanism for the update of additional degrees of freedom becomes limited.
基金supported by the National Natural Science Foundation of China(Grant Nos.11871467,11471329).
文摘In this paper,we consider the initial-boundary value problem(IBVP)for the micropolar Naviers-Stokes equations(MNSE)and analyze a first order fully discrete mixed finite element scheme.We first establish some regularity results for the solution of MNSE,which seem to be not available in the literature.Next,we study a semi-implicit time-discrete scheme for the MNSE and prove L2-H1 error estimates for the time discrete solution.Furthermore,certain regularity results for the time discrete solution are establishes rigorously.Based on these regularity results,we prove the unconditional L2-H1 error estimates for the finite element solution of MNSE.Finally,some numerical examples are carried out to demonstrate both accuracy and efficiency of the fully discrete finite element scheme.
文摘Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous, non-isotropic matter without using (or in the absence of) the mathematical models of the BVPs and the IVPs. These methods are also used for deriving mathematical models for BVPs and IVPs associated with isotropic, homogeneous as well as non-homogeneous, non-isotropic continuous matter. In energy methods when applied to IVPs, one constructs energy functional (<i>I</i>) consisting of kinetic energy, strain energy and the potential energy of loads. The first variation of this energy functional (<em>δI</em>) set to zero is a necessary condition for an extremum of <i>I</i>. In this approach one could use <i>δI</i> = 0 directly in constructing computational processes such as the finite element method or could derive Euler’s equations (differential or partial differential equations) from <i>δI</i> = 0, which is also satisfied by a solution obtained from <i>δI</i> = 0. The Euler’s equations obtained from <i>δI</i> = 0 indeed are the mathematical model associated with the energy functional <i>I</i>. In case of BVPs we follow the same approach except in this case, the energy functional <i>I</i> consists of strain energy and the potential energy of loads. In using the principle of virtual work for BVPs and the IVPs, we can also accomplish the same as described above using energy methods. In this paper we investigate consistency and validity of the mathematical models for isotropic, homogeneous and non-isotropic, non-homogeneous continuous matter for BVPs that are derived using energy functional consisting of strain energy and the potential energy of loads. Similar investigation is also presented for IVPs using energy functional consisting of kinetic energy, strain energy and the potential energy of loads. The computational approaches for BVPs and the IVPs designed using energy functional and principl
文摘This paper is to study extension of high resolution kinetic flux-vector splitting (KFVS) methods. In this new method, two Maxwellians are first introduced to recover the Euler equations with an additional conservative equation. Next, based on the well-known connection between the Euler equations and Boltzmann equations, a class of high resolution KFVS methods are presented to solve numerically multicomponent flows. Our method does not solve any Riemann problems, and add any nonconservative corrections. The numerical results are also presented to show the accuracy and robustness of our methods. These include one-dimensional shock tube problem, and two-dimensional interface motion in compressible flows. The computed solutions are oscillation-free near material fronts, and produce correct shock speeds.
