Many engineering applications need to analyse the system dynamics on the macro and micro level,which results in a larger computational effort.An explicit-implicit asynchronous step algorithm is introduced to solve the...Many engineering applications need to analyse the system dynamics on the macro and micro level,which results in a larger computational effort.An explicit-implicit asynchronous step algorithm is introduced to solve the structural dynamics in multi-scale both the space domain and time domain.The discrete FEA model is partitioned into explicit and implicit parts using the nodal partition method.Multiple boundary node method is adopted to handle the interface coupled problem.In coupled region,the implicit Newmark coupled with an explicit predictor corrector Newmark whose predictive wave propagates into the implicit mesh.During the explicit subcycling process,the variables of boundary nodes are solved directly by dynamics equilibrium equation.The dissipation energy is dynamically determined in accordance with the energy balance checking.A cantilever beam and a building two numerical examples are proposed to verify that the method can greatly reduce the computing time while maintaining a high accuracy.展开更多
This paper focuses on studying a new energy-work relationship numerical integration scheme of nonholo-nomie Hamiltonian systems.The signal-stage numerical,multi-stage and parallel composition numerical integrationsche...This paper focuses on studying a new energy-work relationship numerical integration scheme of nonholo-nomie Hamiltonian systems.The signal-stage numerical,multi-stage and parallel composition numerical integrationschemes are presented.The high-order energy-work relation scheme of the system is constructed by a parallel connectionof n multi-stage schemes of order 2,its order of accuracy is 2n.The connection,which is discrete analogue of usual case,between the change of energy and work of nonholonomic constraint forces is obtained for nonholonomic Hamiltoniansystems.This paper also gives that there is smaller error of the scheme when taking a large number of stages than a lessone.Finally,an applied example is discussed to illustrate these results.展开更多
To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of...To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.展开更多
By introducing the s-parameterized generalized Wigner operator into phase-space quantum mechanics we invent the technique of integration within s-ordered product of operators (which considers normally ordered, antino...By introducing the s-parameterized generalized Wigner operator into phase-space quantum mechanics we invent the technique of integration within s-ordered product of operators (which considers normally ordered, antinormally ordered and Weyl ordered product of operators as its special cases). The s-ordered operator expansion (denoted by s…s ) formula of density operators is derived, which isρ=2/1-s∫d^2β/π〈-β|ρ|β〉sexp{2/s-1(s|β|^2-β*α+βa-αα)}s The s-parameterized quantization scheme is thus completely established.展开更多
A group of asymmetric difference schemes to approach the Korteweg-de Vries (KdV) equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segme...A group of asymmetric difference schemes to approach the Korteweg-de Vries (KdV) equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.展开更多
With the booming development of terrestrial network, scaling terrestrial network over satellite network to build Integrated Terrestrial-Satellite Network(ITSN) and meanwhile to provide the global Internet access, has ...With the booming development of terrestrial network, scaling terrestrial network over satellite network to build Integrated Terrestrial-Satellite Network(ITSN) and meanwhile to provide the global Internet access, has become ever more attractive. Naturally, the widely and successfully used terrestrial routing protocols are the promising protocols to integrate the terrestrial and satellite networks. However, the terrestrial routing protocols, which rely on propagating routing messages to discover New Network Topology(NNT) in the terrestrial network with rare topology changes, will suffer from overly numerous routing messages in satellite network whose topology frequently changes as satellites move. In this paper, a Topology Discovery Sub-layer for ITSN Routing Schemes(TDS-IRS) is firstly proposed to avoid the propagation of numerous routing messages by taking advantage of the movement predictability of satellite and the requirements of routing schemes to discover NNT in advance of topology change. Secondly, a Weighted Perfect Matching based Topology Discovery(WPM-TD) model is designed to conduct the NNT discovery on the ground. Thirdly, this paper builds a testbed with real network devices and meanwhile interconnect that testbed with real Internet, to validate that RS-TDS can discover NNT immediately with the less on-board overhead compared with optimized routing schemes. Finally, different network scenarios are applied to validate the WPM-TD, i.e., the core module of TDS-IRS. Extensive experiments show WPM-TD can work efficiently, avoiding the invalid NNT discovery and decreasing 20% ~ 57% of potential topology changes, which can also improve up to 47% ~ 105% of network throughput.展开更多
A direct time integration scheme based on Gauss-Legendre quadrature is proposed to solve problems in linear structural dynamics.The proposed method is a oneparameter non-dissipative scheme.Improved stability,accuracy,...A direct time integration scheme based on Gauss-Legendre quadrature is proposed to solve problems in linear structural dynamics.The proposed method is a oneparameter non-dissipative scheme.Improved stability,accuracy,and dispersion characteristics are achieved using appropriate values of the parameter.The proposed scheme has second-order accuracy with and without physical damping.Moreover,its stability,accuracy,and dispersion are analyzed.In addition,its performance is demonstrated by the two-dimensional scalar wave problem,the single-degree-of-freedom problem,two degrees-of-freedom spring system,and beam with boundary constraints.The wave propagation problem is solved in the high frequency wave regime to demonstrate the advantage of the proposed scheme.When the proposed scheme is applied to solve the wave problem,more accurate solutions than those of other methods are obtained by using the appropriate value of the parameter.For the single-degree-offreedom system,two degrees-of-freedom system,and the time responses of beam,the proposed scheme can be used effectively owing to its high accuracy and lower computational cost.展开更多
In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubi...In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.展开更多
This paper presents a ZUC-256 stream cipher algorithm hardware system in order to prevent the advanced security threats for 5 G wireless network.The main innovation of the hardware system is that a six-stage pipeline ...This paper presents a ZUC-256 stream cipher algorithm hardware system in order to prevent the advanced security threats for 5 G wireless network.The main innovation of the hardware system is that a six-stage pipeline scheme comprised of initialization and work stage is employed to enhance the solving speed of the critical logical paths.Moreover,the pipeline scheme adopts a novel optimized hardware structure to fast complete the Mod(231-1)calculation.The function of the hardware system has been validated experimentally in detail.The hardware system shows great superiorities.Compared with the same type system in recent literatures,the logic delay reduces by 47%with an additional hardware resources of only 4 multiplexers,the throughput rate reaches 5.26 Gbps and yields at least 45%better performance,the throughput rate per unit area increases 14.8%.The hardware system provides a faster and safer encryption module for the 5G wireless network.展开更多
The quantitative rules of the transfer and variation of errors,when the Gaussian integral functions F.(z) are evaluated sequentially by recurring,have been expounded.The traditional viewpoint to negate the applicabili...The quantitative rules of the transfer and variation of errors,when the Gaussian integral functions F.(z) are evaluated sequentially by recurring,have been expounded.The traditional viewpoint to negate the applicability and reliability of upward recursive formula in principle is amended.An optimal scheme of upward-and downward-joint recursions has been developed for the sequential F(z) computations.No additional accuracy is needed with the fundamental term of recursion because the absolute error of Fn(z) always decreases with the recursive approach.The scheme can be employed in modifying any of existent subprograms for Fn<z> computations.In the case of p-d-f-and g-type Gaussians,combining this method with Schaad's formulas can reduce,at least,the additive operations by a factor 40%;the multiplicative and exponential operations by a factor 60%.展开更多
基金supported by the National Key Research and Development Program of China(2016YFB0201800)the National Natural Science Foundation of China(No.51475287 and No.11772192).
文摘Many engineering applications need to analyse the system dynamics on the macro and micro level,which results in a larger computational effort.An explicit-implicit asynchronous step algorithm is introduced to solve the structural dynamics in multi-scale both the space domain and time domain.The discrete FEA model is partitioned into explicit and implicit parts using the nodal partition method.Multiple boundary node method is adopted to handle the interface coupled problem.In coupled region,the implicit Newmark coupled with an explicit predictor corrector Newmark whose predictive wave propagates into the implicit mesh.During the explicit subcycling process,the variables of boundary nodes are solved directly by dynamics equilibrium equation.The dissipation energy is dynamically determined in accordance with the energy balance checking.A cantilever beam and a building two numerical examples are proposed to verify that the method can greatly reduce the computing time while maintaining a high accuracy.
基金supported by National Natural Science Foundation of China under Grant No.10672143the Natural Science Foundation of Henan Province under Grant No.0511022200
文摘This paper focuses on studying a new energy-work relationship numerical integration scheme of nonholo-nomie Hamiltonian systems.The signal-stage numerical,multi-stage and parallel composition numerical integrationschemes are presented.The high-order energy-work relation scheme of the system is constructed by a parallel connectionof n multi-stage schemes of order 2,its order of accuracy is 2n.The connection,which is discrete analogue of usual case,between the change of energy and work of nonholonomic constraint forces is obtained for nonholonomic Hamiltoniansystems.This paper also gives that there is smaller error of the scheme when taking a large number of stages than a lessone.Finally,an applied example is discussed to illustrate these results.
基金The project supported by the National Key Program for Developing Basic Sciences (G1999043408 and G1998040901-1)the National Natural Sciences Foundation of China (40175024 and 40035010)
文摘To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 10775097 and 10874174)
文摘By introducing the s-parameterized generalized Wigner operator into phase-space quantum mechanics we invent the technique of integration within s-ordered product of operators (which considers normally ordered, antinormally ordered and Weyl ordered product of operators as its special cases). The s-ordered operator expansion (denoted by s…s ) formula of density operators is derived, which isρ=2/1-s∫d^2β/π〈-β|ρ|β〉sexp{2/s-1(s|β|^2-β*α+βa-αα)}s The s-parameterized quantization scheme is thus completely established.
基金Project supported by the National Natural Science Foundation of China(No.10671113)the Natural Science Foundation of Shandong Province of China(No.Y2003A04)
文摘A group of asymmetric difference schemes to approach the Korteweg-de Vries (KdV) equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.
基金supported by State Key Program of National Natural Science of China (91738202)Science &Technology Program of Beijing (Z171100005217001)
文摘With the booming development of terrestrial network, scaling terrestrial network over satellite network to build Integrated Terrestrial-Satellite Network(ITSN) and meanwhile to provide the global Internet access, has become ever more attractive. Naturally, the widely and successfully used terrestrial routing protocols are the promising protocols to integrate the terrestrial and satellite networks. However, the terrestrial routing protocols, which rely on propagating routing messages to discover New Network Topology(NNT) in the terrestrial network with rare topology changes, will suffer from overly numerous routing messages in satellite network whose topology frequently changes as satellites move. In this paper, a Topology Discovery Sub-layer for ITSN Routing Schemes(TDS-IRS) is firstly proposed to avoid the propagation of numerous routing messages by taking advantage of the movement predictability of satellite and the requirements of routing schemes to discover NNT in advance of topology change. Secondly, a Weighted Perfect Matching based Topology Discovery(WPM-TD) model is designed to conduct the NNT discovery on the ground. Thirdly, this paper builds a testbed with real network devices and meanwhile interconnect that testbed with real Internet, to validate that RS-TDS can discover NNT immediately with the less on-board overhead compared with optimized routing schemes. Finally, different network scenarios are applied to validate the WPM-TD, i.e., the core module of TDS-IRS. Extensive experiments show WPM-TD can work efficiently, avoiding the invalid NNT discovery and decreasing 20% ~ 57% of potential topology changes, which can also improve up to 47% ~ 105% of network throughput.
文摘A direct time integration scheme based on Gauss-Legendre quadrature is proposed to solve problems in linear structural dynamics.The proposed method is a oneparameter non-dissipative scheme.Improved stability,accuracy,and dispersion characteristics are achieved using appropriate values of the parameter.The proposed scheme has second-order accuracy with and without physical damping.Moreover,its stability,accuracy,and dispersion are analyzed.In addition,its performance is demonstrated by the two-dimensional scalar wave problem,the single-degree-of-freedom problem,two degrees-of-freedom spring system,and beam with boundary constraints.The wave propagation problem is solved in the high frequency wave regime to demonstrate the advantage of the proposed scheme.When the proposed scheme is applied to solve the wave problem,more accurate solutions than those of other methods are obtained by using the appropriate value of the parameter.For the single-degree-offreedom system,two degrees-of-freedom system,and the time responses of beam,the proposed scheme can be used effectively owing to its high accuracy and lower computational cost.
文摘In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.
基金supported in part by the National R&D Program for Major Research Instruments of China(Grant No:62027814)the National Natural Science Foundation of China(Grant No:62104054)+2 种基金the Natural Science Foundation of Heilongjiang Province(Grant No:F2018010)the Postdoctoral Science Foundation of Heilongjiang Province,China(No:LBH-Z20133)the Fundamental Research Funds for The Central Universities,China(3072021CF0806)。
文摘This paper presents a ZUC-256 stream cipher algorithm hardware system in order to prevent the advanced security threats for 5 G wireless network.The main innovation of the hardware system is that a six-stage pipeline scheme comprised of initialization and work stage is employed to enhance the solving speed of the critical logical paths.Moreover,the pipeline scheme adopts a novel optimized hardware structure to fast complete the Mod(231-1)calculation.The function of the hardware system has been validated experimentally in detail.The hardware system shows great superiorities.Compared with the same type system in recent literatures,the logic delay reduces by 47%with an additional hardware resources of only 4 multiplexers,the throughput rate reaches 5.26 Gbps and yields at least 45%better performance,the throughput rate per unit area increases 14.8%.The hardware system provides a faster and safer encryption module for the 5G wireless network.
文摘The quantitative rules of the transfer and variation of errors,when the Gaussian integral functions F.(z) are evaluated sequentially by recurring,have been expounded.The traditional viewpoint to negate the applicability and reliability of upward recursive formula in principle is amended.An optimal scheme of upward-and downward-joint recursions has been developed for the sequential F(z) computations.No additional accuracy is needed with the fundamental term of recursion because the absolute error of Fn(z) always decreases with the recursive approach.The scheme can be employed in modifying any of existent subprograms for Fn<z> computations.In the case of p-d-f-and g-type Gaussians,combining this method with Schaad's formulas can reduce,at least,the additive operations by a factor 40%;the multiplicative and exponential operations by a factor 60%.