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 numericM integration scheme of nonholonomic Hamiltonian systems. The signal-stage numerical, multi-stage and parallel composition numerical integration sch...This paper focuses on studying a new energy-work relationship numericM integration scheme of nonholonomic Hamiltonian systems. The signal-stage numerical, multi-stage and parallel composition numerical integration schemes are presented. The high-order energy-work relation scheme of the system is constructed by a parallel connection of 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 nonholonomie Hamiltonian systems. This paper also gives that there is smaller error of the scheme when taking a large number of stages than a less one. 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.展开更多
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.展开更多
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.展开更多
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.展开更多
Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics(MD)simulations.Despite the massive use of the conventional leapfrog(LF)integrator in modern computational tools wit...Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics(MD)simulations.Despite the massive use of the conventional leapfrog(LF)integrator in modern computational tools within the framework of MD propagation,further development for better performance is still possible.The alternative version of LF in the middle thermostat scheme(LFmiddle)achieves a higher order of accuracy and efficiency and maintains stable dynamics even with the integration time stepsize extended by several folds.In this work,we perform a benchmark test of the two integrators(LF and LF-middle)in extensive conventional and enhanced sampling simulations,aiming at quantifying the time-stepsizeinduced variations of global properties(e.g.,detailed potential energy terms)as well as of local observables(e.g.,free energy changes or bondlengths)in practical simulations of complex systems.The test set is composed of six chemically and biologically relevant systems,including the conformational change of dihedral flipping in the N-methylacetamide and an AT(AdenineThymine)tract,the intra-molecular proton transfer inside malonaldehyde,the binding free energy calculations of benzene and phenol targeting T4 lysozyme L99A,the hydroxyl bond variations in ethaline deep eutectic solvent,and the potential energy of the blue-light using flavin photoreceptor.It is observed that the time-step-induced error is smaller for the LFmiddle scheme.The outperformance of LF-middle over the conventional LF integrator is much more significant for global properties than local observables.Overall,the current work demonstrates that the LF-middle scheme should be preferably applied to obtain accurate thermodynamics in the simulation of practical chemical and biological systems.展开更多
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%.展开更多
The solution of time-dependent hyperbolic conservation laws on cut cell meshes causes the small cell problem:standard schemes are not stable on the arbitrarily small cut cells if an explicit time stepping scheme is us...The solution of time-dependent hyperbolic conservation laws on cut cell meshes causes the small cell problem:standard schemes are not stable on the arbitrarily small cut cells if an explicit time stepping scheme is used and the time step size is chosen based on the size of the background cells.In May and Berger(J Sci Comput 71:919–943,2017),the mixed explicit-implicit approach in general and MUSCL-Trap(monotonic upwind scheme for conservation laws and trapezoidal scheme)in particular have been introduced to solve this problem by using implicit time stepping on the cut cells.Theoretical and numerical results have indicated that this might lead to a loss in accuracy when switching between the explicit and implicit time stepping.In this contribution,we examine this in more detail and will prove in one dimension that the specific combination MUSCL-Trap of an explicit second-order and an implicit second-order scheme results in a fully second-order mixed scheme.As this result is unlikely to hold in two dimensions,we also introduce two new versions of mixed explicit-implicit schemes based on exchanging the explicit scheme.We present numerical tests in two dimensions where we compare the new versions with the original MUSCL-Trap scheme.展开更多
We investigate the multisymplectic Euler box scheme for the Korteweg-de Vries （KdV） equation. A new completely explicit six-point scheme is derived. Numerical experiments of the new scheme with comparisons to the Za...We investigate the multisymplectic Euler box scheme for the Korteweg-de Vries （KdV） equation. A new completely explicit six-point scheme is derived. Numerical experiments of the new scheme with comparisons to the Zabusky-Kruskal scheme, the multisymplectic 12-point scheme, the narrow box scheme and the spectral method are made to show nice numerical stability and ability to preserve the integral invariant for long-time integration.展开更多
Using the Picard iteration method and treating the involved integration by numerical quadrature formulas, we propose a numerical scheme for the second kind nonlinear Volterra integral equations. For enlarging the conv...Using the Picard iteration method and treating the involved integration by numerical quadrature formulas, we propose a numerical scheme for the second kind nonlinear Volterra integral equations. For enlarging the convergence region of the Picard iteration method, multistage algorithm is devised. We also introduce an algorithm for problems with some singularities at the limits of integration including fractional integral equations. Numerical tests verify the validity of the proposed schemes.展开更多
The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain ...The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.展开更多
Considering the variable cross section thickness of longitudinal profiled plate and the dynamic reductions of straightening rolls,an analytical model combining curvature integral method with linear decreasing straight...Considering the variable cross section thickness of longitudinal profiled plate and the dynamic reductions of straightening rolls,an analytical model combining curvature integral method with linear decreasing straightening scheme was proposed to investigate the longitudinal profiled plate straightening process.Moreover,the calculation flow and solution algorithm of longitudinal profiled plate straightening process were presented.To verify the proposed model,calculated straightening forces were compared with the measured values,and very good agreements were achieved.Then,the reduction,contact angle,reverse bending curvature,residual curvature,straightening force and straightening moment of longitudinal profiled plate in the straightening process were calculated and analyzed,and the calculated results show that the curvature integral method can be used to reveal the mechanism of longitudinal profiled plate straightening.展开更多
We derive a new multisymplectic integrator for the Kawahara-type equation which is a fully explicit scheme and thus needs less computation cost. Multisympecticity of such scheme guarantees the long-time numerical beha...We derive a new multisymplectic integrator for the Kawahara-type equation which is a fully explicit scheme and thus needs less computation cost. Multisympecticity of such scheme guarantees the long-time numerical behaviors. Nu- merical experiments are presented to verify the accuracy of this scheme as well as the excellent performance on invariant preservation for three kinds of Kawahara-type equations.展开更多
This paper studies the limit distributions for discretization error of irregular sam- pling approximations of stochastic integral. The irregular sampling approximation was first presented in Hayashi et al.[3], which w...This paper studies the limit distributions for discretization error of irregular sam- pling approximations of stochastic integral. The irregular sampling approximation was first presented in Hayashi et al.[3], which was more general than the sampling approximation in Lindberg and Rootzen [10]. As applications, we derive the asymptotic distribution of hedging error and the Euler scheme of stochastic differential equation respectively.展开更多
This paper proposes an advanced substation integrated protection communication system based on the latest development in power system.The integrated protection communication system within a substation comprises a numb...This paper proposes an advanced substation integrated protection communication system based on the latest development in power system.The integrated protection communication system within a substation comprises a number of major components such as the merging unit, communication network and integrated protection unit.The design and capability evaluation of the communication network based on IEC61850 are focused on.Results show that huge information from 4 merging units can be transmitted real-time via adopting 100 Mbit/s optical fiber in a substation,and the proposed communication network is able to achieve required performance.展开更多
In the present paper,we derive the asymptotic expansion formula for the trapezoidal approximation of the fractional integral.We use the expansion formula to obtain approximations for the fractional integral of orders...In the present paper,we derive the asymptotic expansion formula for the trapezoidal approximation of the fractional integral.We use the expansion formula to obtain approximations for the fractional integral of ordersα,1+α,2+α,3+αand 4+α.The approximations are applied for computation of the numerical solutions of the ordinary fractional relaxation and the fractional oscillation equations expressed as fractional integral equations.展开更多
基金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 numericM integration scheme of nonholonomic Hamiltonian systems. The signal-stage numerical, multi-stage and parallel composition numerical integration schemes are presented. The high-order energy-work relation scheme of the system is constructed by a parallel connection of 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 nonholonomie Hamiltonian systems. This paper also gives that there is smaller error of the scheme when taking a large number of stages than a less one. 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.
基金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.
基金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.
文摘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.
基金supported by the National Natural Science Foundation of China(No.21961142017)the Ministry of Science and Technology of China(No.2017YFA0204901)。
文摘Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics(MD)simulations.Despite the massive use of the conventional leapfrog(LF)integrator in modern computational tools within the framework of MD propagation,further development for better performance is still possible.The alternative version of LF in the middle thermostat scheme(LFmiddle)achieves a higher order of accuracy and efficiency and maintains stable dynamics even with the integration time stepsize extended by several folds.In this work,we perform a benchmark test of the two integrators(LF and LF-middle)in extensive conventional and enhanced sampling simulations,aiming at quantifying the time-stepsizeinduced variations of global properties(e.g.,detailed potential energy terms)as well as of local observables(e.g.,free energy changes or bondlengths)in practical simulations of complex systems.The test set is composed of six chemically and biologically relevant systems,including the conformational change of dihedral flipping in the N-methylacetamide and an AT(AdenineThymine)tract,the intra-molecular proton transfer inside malonaldehyde,the binding free energy calculations of benzene and phenol targeting T4 lysozyme L99A,the hydroxyl bond variations in ethaline deep eutectic solvent,and the potential energy of the blue-light using flavin photoreceptor.It is observed that the time-step-induced error is smaller for the LFmiddle scheme.The outperformance of LF-middle over the conventional LF integrator is much more significant for global properties than local observables.Overall,the current work demonstrates that the LF-middle scheme should be preferably applied to obtain accurate thermodynamics in the simulation of practical chemical and biological systems.
文摘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%.
文摘The solution of time-dependent hyperbolic conservation laws on cut cell meshes causes the small cell problem:standard schemes are not stable on the arbitrarily small cut cells if an explicit time stepping scheme is used and the time step size is chosen based on the size of the background cells.In May and Berger(J Sci Comput 71:919–943,2017),the mixed explicit-implicit approach in general and MUSCL-Trap(monotonic upwind scheme for conservation laws and trapezoidal scheme)in particular have been introduced to solve this problem by using implicit time stepping on the cut cells.Theoretical and numerical results have indicated that this might lead to a loss in accuracy when switching between the explicit and implicit time stepping.In this contribution,we examine this in more detail and will prove in one dimension that the specific combination MUSCL-Trap of an explicit second-order and an implicit second-order scheme results in a fully second-order mixed scheme.As this result is unlikely to hold in two dimensions,we also introduce two new versions of mixed explicit-implicit schemes based on exchanging the explicit scheme.We present numerical tests in two dimensions where we compare the new versions with the original MUSCL-Trap scheme.
基金Supported by the National Baslc Research Programme under Grant No 2005CB321703, and the National Natural Science Foundation of China under Grant Nos 40221503, 10471067 and 40405019.
文摘We investigate the multisymplectic Euler box scheme for the Korteweg-de Vries （KdV） equation. A new completely explicit six-point scheme is derived. Numerical experiments of the new scheme with comparisons to the Zabusky-Kruskal scheme, the multisymplectic 12-point scheme, the narrow box scheme and the spectral method are made to show nice numerical stability and ability to preserve the integral invariant for long-time integration.
文摘Using the Picard iteration method and treating the involved integration by numerical quadrature formulas, we propose a numerical scheme for the second kind nonlinear Volterra integral equations. For enlarging the convergence region of the Picard iteration method, multistage algorithm is devised. We also introduce an algorithm for problems with some singularities at the limits of integration including fractional integral equations. Numerical tests verify the validity of the proposed schemes.
基金the National Natural Science Foundation of China(No.10571017)supported in part by the National Natural Science Foundation of China(No.60533020)supported in part by NSF DMS 0712744
文摘The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.
基金The authors are grateful for the supports of the National Key Research and Development Program of China(No.2017YFB0306404)Key Project of Hebei Education Department(No.ZD2018203).
文摘Considering the variable cross section thickness of longitudinal profiled plate and the dynamic reductions of straightening rolls,an analytical model combining curvature integral method with linear decreasing straightening scheme was proposed to investigate the longitudinal profiled plate straightening process.Moreover,the calculation flow and solution algorithm of longitudinal profiled plate straightening process were presented.To verify the proposed model,calculated straightening forces were compared with the measured values,and very good agreements were achieved.Then,the reduction,contact angle,reverse bending curvature,residual curvature,straightening force and straightening moment of longitudinal profiled plate in the straightening process were calculated and analyzed,and the calculated results show that the curvature integral method can be used to reveal the mechanism of longitudinal profiled plate straightening.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11271195 and 11271196)the Project of Graduate Education Innovation of Jiangsu Province,China(Grant No.CXZZ12-0385)
文摘We derive a new multisymplectic integrator for the Kawahara-type equation which is a fully explicit scheme and thus needs less computation cost. Multisympecticity of such scheme guarantees the long-time numerical behaviors. Nu- merical experiments are presented to verify the accuracy of this scheme as well as the excellent performance on invariant preservation for three kinds of Kawahara-type equations.
基金Supported by the National Natural Science Foundation of China(11371317)
文摘This paper studies the limit distributions for discretization error of irregular sam- pling approximations of stochastic integral. The irregular sampling approximation was first presented in Hayashi et al.[3], which was more general than the sampling approximation in Lindberg and Rootzen [10]. As applications, we derive the asymptotic distribution of hedging error and the Euler scheme of stochastic differential equation respectively.
基金National Natural Science Foundation of China(No.50677003).
文摘This paper proposes an advanced substation integrated protection communication system based on the latest development in power system.The integrated protection communication system within a substation comprises a number of major components such as the merging unit, communication network and integrated protection unit.The design and capability evaluation of the communication network based on IEC61850 are focused on.Results show that huge information from 4 merging units can be transmitted real-time via adopting 100 Mbit/s optical fiber in a substation,and the proposed communication network is able to achieve required performance.
基金This work was supported by the Bulgarian National Science Fund under Project KP-06-M32/2"Advanced Stochastic and Deterministic Approaches for Large-Scale Problems of Computational Mathematics"。
文摘In the present paper,we derive the asymptotic expansion formula for the trapezoidal approximation of the fractional integral.We use the expansion formula to obtain approximations for the fractional integral of ordersα,1+α,2+α,3+αand 4+α.The approximations are applied for computation of the numerical solutions of the ordinary fractional relaxation and the fractional oscillation equations expressed as fractional integral equations.
