RETROSPECTIVE TIME INTEGRAL SCHEME AND ITS APPLICATIONS TO THE ADVECTION EQUATION

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.

A STEP-BY-STEP INTEGRAL METHOD FOR ACTIVELY CONTROLLED STRUCTURES

A new active control method was proposed, in which the analytical control law was deduced by using a step by step integral method to differential equation of motion under the condition of static error being zero. This control law is terse in mathematical expression and convenient for practical use. The simulation results demonstrate that the proposed method can provide much more remarkable peak response reduction of seismically excited structures than the classical LQR method.

OPTIMAL SCHEME FOR SEQUENTIAL COMPUTATIONS OF F_m(z) INTEGRALS IN AB INITIO CALCULATIONS----COMBINATORY USE OF UPWARD AND DOWNWARD RECURSIONS

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%.

Multistage Numerical Picard Iteration Methods for Nonlinear Volterra Integral Equations of the Second Kind

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.

Topology Discovery Sub-Layer for Integrated Terrestrial-Satellite Network Routing Schemes

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.

Discretization error of irregular sampling approximations of stochastic integrals

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.

DIRECT INTEGRATION METHODS WITH INTEGRAL MODEL FOR DYNAMIC SYSTEMS

A new approach which is a direct integration method with integral model (DIM-IM) to solve dynamic governing equations is presented. The governing equations ore integrated into the integral equations. An algorithm with explicit and predict-correct and self-starting and fourth-order accuracy to integrate the integral equations is given. Theoretical analysis and numerical examples shaw that DIM-IM described in this paper suitable for strong nonlinear and non-conservative system have higher accuracy than central difference, Houbolt, Newmark and Wilson-Theta methods.

COORDINATION OF THE MOTION PARAMETERS WITHIN STEP-BY-STEP INTEGRATION FOR DYNAMIC EQUATION

A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved using an energy equation and the calculating data of the original method. The method presented is better than the original method in terms of calculating postulations and is in better conformity with the system's movement. Take the Wilson-θ method as an example. By using the coordination process, the calculation precision has been greatly im proved (reducing the errors by approximately 90% ), and the greater part of overshooting of the calculation result has been eliminated. The study suggests that the mal-coordination of the motion parameters within the time-step is the major factor that contributes to the result errors of step-by-step integration for the dynamic equation.

Approximations of the Fractional Integral and Numerical Solutions of Fractional Integral Equations

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.

Simple Program for Step-by-Step Time Integration in Chemical Kinetics, Applied to Simple Model for Hydrogen Combustion

A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by accurate computational studies of [3]. 3) In the first step of the algorithm, the only nonzero species concentrations are those of reactants [H2], [O2]. So, the maximum reaction rate is that of the forward initiation reaction max R = Rif = kif[H2] [O2], where the rate constant kif is very slow. Thus, the first time steplength Δt1 = 10-6mol·L-1/max R results long in sec. After the first step, the sequences of all the following Δt’s are very short, in μsec. So, the first time steplength Δt1 can be considered as ignition delay time. 4) It is assumed that explosion corresponds to ignition delay time Δt1 t1 = 10 sec., can be considered as explosion limit curve. This curve is compared to the corresponding one obtained by the accurate computational studies of [2].

Incompressible Flow and Heat Transfer over a Plate: A Hybrid Integral Domain-Discretized Numerical Procedure

This work deals with incompressible two-dimensional viscous flow over a semi-infinite plate ac-cording to the approximations resulting from Prandtl boundary layer theory. The governing non-linear coupled partial differential equations describing laminar flow are converted to a self-simi- lar type third order ordinary differential equation known as the Falkner-Skan equation. For the purposes of a numerical solution, the Falkner-Skan equation is converted to a system of first order ordinary differential equations. These are numerically addressed by the conventional shooting and bisection methods coupled with the Runge-Kutta technique. However the accompanying energy equation lends itself to a hybrid numerical finite element-boundary integral application. An appropriate complementary differential equation as well as the Green second identity paves the way for the integral representation of the energy equation. This is followed by a finite element-type discretization of the problem domain. Based on the quality of the results obtained herein, a strong case is made for a hybrid numerical scheme as a useful approach for the numerical resolution of boundary layer flows and species transport. Thanks to the sparsity of the resulting coefficient matrix, the solution profiles not only agree with those of similar problems in literature but also are in consonance with the physics they represent.

A Closed-Form Solution of a Kinetic Integral Equation for Rarefied Gas Flow in a Cylindrical Duct

A spectral method based on Hermite cubic splines expansions combined with a collocation scheme is used to develop a solution for the vector form integral S-model kinetic equation describing rarefied gas flows in cylindrical geometry. Some manipulations are made to facilitate the computational treatment of the singularities inherent to the kernel. Numerical results for the simulation of flows generated by pressure and thermal gradients, Poiseuille and thermal-creep problems, are presented.

Energy-Work Connection Integration Scheme for Nonholonomic Hamiltonian Systems

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.

制导控制系统半实物仿真总体方案设计及集成联试

半实物仿真作为一个大型分布式仿真系统,其总体方案的设计优劣和集成联试的是否顺利,直接影响了系统性能指标的表现和建设目标的实现。归纳总结了总体方案设计的工作内容,分析了实时性、兼容性、扩展性和安全性的性能要求;给出了包括功能层次、运行机制和结构组成在内的典型半实物仿真系统总体方案;概括了仿真系统的集成联试内容,总结了典型系统的集成联试步骤,给出了系统验收测试的内容和系统可信度评估方法。

高速铁路列车运行图一体化自动编制方法研究

高速铁路列车运行图编制不仅要在图面上生成列车运行线,还要考虑客运需求、车底交路运用、车站股道运用等问题,影响因素众多,且各影响因素相关性较强,因而我国高速铁路列车运行图一直以来都采用人工铺画方式,缺少计算机一体化自动编制方法。为进一步提高高速铁路运行图编制效率,在借鉴人工编图经验的基础上,提出基于车底交路方案、列车结构方案及调整策略的高速铁路列车运行图一体化自动编制流程,同时生成列车时刻表、车底交路方案、车站股道运用方案。针对不易处理的列车运行图结构优化问题,设计微观站场进路占用策略与运行图优化编制模型,并以京沪高速铁路为例进行验证,用时48 h重新铺画运行图,各项指标与真实运行图相近,大幅减少了现有运行图编制模式对人力物力的消耗,为高速铁路列车运行图编制提供一种自动化的理论与技术支撑。

偏好异质性、农户参与式方案创设与养殖废弃物资源化利用——来自养殖户种养结合的选择实验

生态补偿政策构成的农业环境方案对于废弃物资源化利用具有激励作用,但现行补偿政策方案可能存在推广瓶颈,突破瓶颈的关键在于了解农户的政策选择偏好。以养殖废弃物肥料化还田为例,通过选择实验法构建了由4个政策属性组合而成的备选政策情境,结合江苏省346户生猪养殖户实验数据,运用Mixed Logit模型检验了补偿政策的选择偏好及协同互补效应。结果表明:(1)养殖户对选择实验方案中的服务类、收入类和价格类补偿政策的偏好程度更高,而对于技术支持类补偿政策不具备明显偏好;(2)养殖户环保认知与风险特征是上述政策偏好异质性的重要来源,风险偏好型和环保认知水平高的养殖户对各类补偿政策方案的接受意愿更高;(3)技术支持、外包服务同收入类、价格类补偿政策之间呈现显著的互补性。根据“重视多元化补偿,强化正协同效应”原则,本文创设了针对养殖废弃物还田利用的推广方案,旨在优化和完善现行的生态补偿政策体系,打好相关政策措施的“组合拳”。

基于积分分析方法对聚变堆启动氚量与所需氚增殖比的评估

托卡马克聚变堆的主要发展方式包括混合堆、纯聚变堆。关于托卡马克聚变堆氚自持的研究,国内外主要采用平均滞留时间方法进行研究,并且针对聚变功率较低的混合堆的氚自持研究较少。本工作采用更符合实际的积分分析方法分析了混合堆、纯聚变堆氚自持的启动氚量、氚增殖比(TBR)要求。研究结果表明:启动氚量、备用氚量与聚变功率具有线性关系,所需TBR与聚变功率呈反比例关系;混合堆聚变功率较低,所需TBR较高,工程实现所需TBR挑战较大,需要通过限制长期氚滞留量以降低所需TBR要求;纯聚变堆聚变功率高,所需TBR较低,工程实现所需TBR挑战较小,但备用氚需求达数十千克,应考虑氚系统的冗余设计或提高氚系统的可靠性、可维护性,以降低备用氚的使用规模;运行因子是聚变堆的一个重要设计指标,在此着重分析了运行因子对所需TBR的影响,并重新定义了一个聚变堆氚自持的关系式,以突出运行因子对氚自持的重要影响。

基于天窗基本单元的编组站天窗方案优化研究

编组站接发列车方向多、数量大,机车换挂和调车作业频繁,直接影响施工维修作业。既有固定天窗时长的天窗方案,难以兼顾维修施工作业和运输工作连续性需求,导致天窗兑现率低,增加安全隐患,研究提出基于天窗基本单元的编组站天窗方案优化方法。根据各维修养护工种作业特性、编组站工作人员数量及效率等因素,划分天窗基本单元;考虑编组站维修养护与运输组织需求,以效用、作业完成率和天窗阈值利用率最优为目标,以时序关系为约束,建立基于阈值的天窗方案多目标优化模型,并利用数学软件进行求解;以成都北编组站为例,验证模型的有效性,结果表明,通过合理优化天窗方案,天窗兑现率达100%,有助于缓解编组站维修与运输之间的冲突。

改进蜉蝣算法求解工艺规划与调度集成问题

针对工艺规划与调度集成问题(IPPS),基于AND/OR析取图和“组合”的概念以最小化最大完工时间为优化目标建立数学模型;设计一种整数编码方案,该方案无需预先为工件生成加工路径,可同时处理工艺规划与调度问题,由于初始解对算法的寻优能力的影响,在种群初始化时考虑机器的负荷以提高初始种群的质量;解码时采用插空的方法生成活跃调度,以缩短工件整体的加工时间;离散蜉蝣算法实现求解IPPS,采用Metropolis准则接收较差的解,引入自适应因子提高算法的收敛速度,设计局部搜索算法提高算法的精度。最后通过大规模基准算例进行实验并与其他算法对比分析,验证该算法有较好寻优性能。

Enhancement of Solar PV Panel Efficiency Using Double Integral Sliding Mode MPPT Control

The extraction of maximum power from the solar panels,using the sliding mode control scheme,becomes popular for partial weather atmospheric conditions due to its effective dynamic duty cycle ratio.However,the sliding mode control scheme was sophisticated with single integral and double integral sliding mode control scheme,which offer enhanced maximum power extraction and support enhanced solar panel efficiency in partial weather conditions.The operation of the sliding mode control scheme depends on the selection of a sliding surface selection based on the atmospheric weather condition,which enables the effective sliding duty cycle ratio operation for the DC/DC boost converter.The duty cycle ratio of the sliding mode control resembles the usual dynamic behavior to achieve enhanced efficiency compared to the various maximum power point tracking(MPPT)schemes.The major limitation of the sliding mode control scheme is to achieve the steady state voltage error of the solar panel in minimum settling time duration.The single integral sliding mode control scheme achieves the expected steady state voltage error limit but fails to achieve minimum settling time duration.Hence,the single integral sliding mode control is extended to a double integral sliding mode control scheme to achieve both steady state voltage error limits within the minimum settling time duration.This double integral sliding mode control scheme allows us to obtain the higher sliding surface duty cycle ratio which acts as the input signal to the boost converter.This activates the enhanced stable and reliable system operation,and nullifies the lacuna of maximum solar panel efficiency under partial weather conditions.Hence,this paper aims to present the design and performance operation of the double integral sliding mode(DISM)MPPT control scheme.To validate the performance analysis of the proposed DISM MPPT control scheme,the MATLAB/Simulink model is designed and verified.Also,the performance analysis of the proposed DISM MPPT control scheme is compared with the sliding mode controller(SMC)scheme and single integral sliding mode controller(SiSMC)scheme.The performance analysis of the proposed double integral sliding mode controller(DISMC)scheme attains 99.10%of efficiency and a very less setting time of 0.035s when compared to other existingmethods.