Multiple Space Debris Collecting Mission—Debris Selection and Trajectory Optimization

This paper investigates the cost requirement for a space debris collecting mission aimed at removing heavy debris from low Earth orbits. The problem mixes combinatorial optimization to select the debris among a list of candidates and functional optimization to define the orbital manoeuvres. The solving methodology proceeds in two steps: Firstly, a specific transfer strategy with impulsive manoeuvres is defined so that the problem becomes of finite dimension; secondly the problem is linearized around an initial reference solution. A Branch and Bound algorithm is then applied iteratively to optimize simultaneously the debris selection and the orbital manoeuvres, yielding a new reference solution. The optimal solutions found are close to the initial guess despite a very complicated design space. The method is exemplified on a representative application case.     

泥石流固液分相流速计算方法研究

泥石流固液分相流速是泥石流对岸坡、防治结构冲击、磨损机理的核心问题．将泥石流体简化为具有相同粒径的固相和具有相同力学性质的液相,基于泥石流体为沿流动方向的一维两相流体,运用两相流理论建立了泥石流固液分相流速控制方程．构建了泥石流平均压力、彻体力及平均表面力的计算方法,尤其通过浆体的Binhanm体流变方程、Bagnold颗粒相互作用试验成果建立了控制体平均表面力计算方法;建立了固液两相流速比例系数,以及理论固相流速与实际流速的比例系数．据此求解控制方程得到了固液分相流速计算方法,该方法既可同时适用于粘性泥石流和稀性泥石流,也可在泥石流爆发以后通过现场采集沉积物分析反求泥石流爆发期间的分相流速．工程实例分析显示,该方法计算结果与实测结果吻合较好．     

用LSD法再探泥石流地貌要素之隐式关系

藉助MATLAB工具,对泥石流沟的六个地貌要素之流域面积、主沟长、切割度、沟口高程、流域高差和沟床比降等,一起用最小平方距离法(LSD)进行了初步处理分析之探讨:建立了线性相关的回归方程,并与用最小二乘法(LS)得到的线性回归方程进行比较,一些参数的估计值得到了改善。     

Modelling and remodelling of biological tissue in the framework of continuum biomechanics

A biological tissue in general is formed by cells, extracellular matrix (ECM) and fluids. Consequently, its overall material behaviour results from its components and their interaction among each other. Furthermore, in case of living tissues, the material properties do not remain constant but naturally change due to adaptation processes or diseases. In the context of the Theory of Porous Media (TPM), a continuum-mechanical model is introduced to describe the complex fluid-structure interaction in biological tissue on a macroscopic scale. The tissue is treated as an aggregate of two immiscible constituents, where the cells and the ECM are summarised to a solid phase, whereas the fluid phase represents the extracellular and interstitial liquids as well as necrotic debris and cell or matrix precursors in solution. The growth and remodelling processes are described by a distinct mass exchange between the fluid and solid phase, which also results in a change of the constituent material behaviour. To furthermore guarantee the compliance with the entropy principle, the growth energy is introduced as an additional quantity. It measures the average of chemical energy available for cell metabolism, and thus, controls the growth and remodelling processes. To set an example, the presented model is applied for the simulation of the early stages of avascular tumour growth in the framework of the finite element method (FEM). (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

目标拦截模拟系统中的造型方法

对目标拦截模拟系统中的造型技术作了研究.包括飞行器造型、地形造型、火焰喷射模拟等.提出了蒙皮技术造型中存在的三个问题及相应改进措施.给出了一个基于蒙皮技术的地形生成新算法.利用嵌套圆台扰动的思想较好地实现了导弹飞行时的火焰喷射模拟.     

Application of metal laminates to aircraft structures: Prediction of penetration performance

A major thrust of the transportation industries in the US is the incorporation of advanced structural materials in airplanes and automobiles. These advanced materials include metal matrix composites, where particulate or whisker reinforcements stiffen and strengthen a ductile phase matrix. In aerospace structures, specific stiffness, improved strength, and weight reduction are key factors. Both in the certain parts of the airframe structure and in the engine fan containment area, increased emphasis is being given to the ability of the material to resist penetration from engine debris and other projectiles which might impact the aircraft structures. Experimental measurements of the ballistic limit velocity of a material versus material thickness gives a method to rank the relative penetration performance of aircraft structural materials. Dynamic finite element analysis aids in understanding the experimental results and in predicting the aircraft debris containment response. For certain aluminum alloys and metal laminates, the relationship between the ballistic limit velocity and plate thickness is linear, while for an aerospace titanium alloy, the ballistic penetration response is more complex.     

泥石流对速流结构的磨损作用研究

泥石流对砼防治结构的磨损是泥石流研究极其薄弱的环节,泥石流磨损大大缩短了防治结构的有效使用年限．作为大型特大型公路泥石流有效防治结构的速流结构,因泥石流的磨损使其有效使用寿命缩短了近三分之一．基于泥石流两相流流速计算理论,研究了泥石流对速流结构的磨损机理,包括均质浆体对速流槽的磨蚀机理和固相颗粒对速流槽的切削机理;基于牛顿运动原理建立了泥石流液相浆体及固相颗粒的磨损动力学方程;运用量纲分析法构建了速流槽磨损量计算式,计算结果与实际情况吻合．研究成果为速流结构的有效设计提供了理论依据．     

随机模拟化学反应系统的加速L-leap算法

提出了随机模拟化学反应系统的加速L-leap算法,该算法根据leap条件确定具有最大倾向函数的反应通道的反应次数,并利用二项分布随机数生成其它反应通道在当前leap时间区间内的反应次数．L-leap算法可更好地满足leap条件．数值模拟实验表明该算法能取得更好的模拟性能．     

A statistical process control approach to selecting a warm-up period for a discrete-event simulation

The selection of a warm-up period for a discrete-event simulation continues to be problematic. A variety of selection methods have been devised, and are briefly reviewed. It is apparent that no one method can be recommended above any other. A new approach, based upon the principles of statistical process control, is described (SPC method). Because simulation output data are often highly autocorrelated and potentially non-normal, the batch means method is employed in constructing the control chart. The SPC method is tested on seven data sets and encouraging results are obtained concerning its accuracy. The approach is also discussed with respect to its ease of implementation, simplicity, generality of use and requirements for parameter estimation.     

Numerical simulation of gas-assisted injection molding using CLSVOF method

It is a typical gas-liquid two phase flow phenomenon that gas penetrates the polymer melt in gas-assisted injection molding (GAIM) process. Numerical simulation is now playing an important role in GAIM, in which the accurate simulation of moving interface is of great importance. The level set (LS) method is a popular interface tracking method, but it does not ensure naturally mass-conservation. In order to improve the mass-conservation of LS method, a coupled level-set and volume-of-fluid (CLSVOF) method with mass-correction is presented for the numerical simulations of interfacial flows in GAIM. The performance of this CLSVOF method is demonstrated by two numerical tests including the three-dimensional deformation field test and the dam break problems. Finally the CLSVOF method is employed to simulate the 3D moving interfaces in GAIM, including gas-melt interface and the melt-front interface. The influences of melt temperature and gas delay time are also analyzed detailedly. As a case study, the processes that gas penetrates the polymer melt in complex cavities are also simulated using this method, and the simulation results are in agreement with those obtained by other researchers.     

上证50ETF期权定价有效性的研究:基于B-S-M模型和蒙特卡罗模拟

上证50ETF期权作为中国资本市场上股票期权的第一个试点产品,其定价问题尤为重要。本文分别运用B-S-M期权定价模型和蒙特卡罗模拟方法对其定价进行实证研究,分析结果表明:1)IGARCH模型比传统的GARCH模型更能较好地拟合上证50ETF的波动率;2)当模拟次数为1000时,蒙特卡罗方法的效率一致地高于B-S-M模型,并且除了对偶变量技术的拟蒙特卡罗其他模型的精确度也都高于B-S-M模型;3)B-S-M模型和蒙特卡罗模拟方法都可以较为准确地、有效地模拟出上证50ETF期权价格。这些研究将为今后期权定价模型的发展和完善提供必要的参考和指引。     

Digital Simulation in the Teaching of Management and Technology

The application of mathematical models to computer simulation laboratories covering management and technological situations is discussed. Laboratory simulation is illustrated in relation to a production‐control laboratory and other teaching applications are discussed. The desirability of ‘stretching’ students during a limited teaching period demands that time on the laboratory model should be speeded up by a factor of perhaps 20,000 over real‐time. Ways of achieving this and some other teaching objectives are described, including some points in the application of hill‐climbing techniques to simulation problems. Finally, the cost effectiveness of this method is briefly covered.     

Application-driven sequential designs for simulation experiments: Kriging metamodelling

This paper proposes a novel method to select an experimental design for interpolation in simulation. Although the paper focuses on Kriging in deterministic simulation, the method also applies to other types of metamodels (besides Kriging), and to stochastic simulation. The paper focuses on simulations that require much computer time, so it is important to select a design with a small number of observations. The proposed method is therefore sequential. The novelty of the method is that it accounts for the specific input/output function of the particular simulation model at hand; that is, the method is application-driven or customized. This customization is achieved through cross-validation and jackknifing. The new method is tested through two academic applications, which demonstrate that the method indeed gives better results than either sequential designs based on an approximate Kriging prediction variance formula or designs with prefixed sample sizes.     

“清除太空垃圾”是否存在商机呢?

对太空轨道碎片清理的商业运营模式和是否存在商机问题建模进行研究。首先,提出了太空轨道碎片清理公司的商业运营模式,即以清理公司、保险公司和卫星公司为主体的三位一体的系统运营模式;其次,给出私营清理公司商机的定义,通过对运营系统资金费用的分析,定量建立系统正常运营的保险公司最低保单价格和卫星公司最高保险价格求解的模型,并分析了模型的可计算性,给出了商机存在的条件以及商机大小的计算方法;最后,借助于清除效率,建立3种不同轨道碎片清除模式下清除成本的模型,给出了计算不同清除方式下保险公司最低保单价格和卫星公司最高保险价格的模型,验证了定义的商机是合理的并可计算的。     

差异驱动型评价方法的稳定性及差异凸显能力比较

针对综合评价中方法种类繁多且无统一比较标准的问题,选取了四种突出被评价对象之间差异的评价方法,采用随机模拟的方式分别从评价方法的稳定性及对差异的凸显能力两个方面进行了比较分析。得出了四种方法稳定性由高到低分别为均方差法、最大离差法、熵值法、拉开档次法,且评价方法的稳定性越高,则其对差异的凸显能力反而越差的结论。该研究不仅验证了差异驱动型评价方法的相关特性,为评价者关于评价方法的选取提供了参考意见,而且随机模拟方法的应用,可为类似的多评价方法的比较问题提供技术参考。     

Starting and stopping rules for simulations using a priori information

This paper proposes a method of interval estimation (based on the autoregressive method) that exploits a priori information about activity level (traffic intensity) and data skewness in a queueing simulation. The method relies on two rules: One determines when to start collecting data during a run and the other determines when to stop collection. The rules are designed to use the a priori information in a way that mitigates the effects of the initial conditions in the simulation and ensures representative congested behavior in the collected data.Experiments with a simulation of the M/M//c queue, with c = 1, 2, 4 and = 0.7, 0.8, 0.9, 0.95 produce favorable results. For = 0.7, 0.8, 0.9, the coverage rates for interval estimates are close to the specified theoretical coverage rates and are higher than those reported in the literature for other methods of interval estimation. The sample sizes to obtain the coverage rates are moderate and are insensitive to variation in the number of servers and the activity level. Experiments with a fixed truncation starting rule and a fixed sample-size stopping rule clearly demonstrate the effectiveness of the proposed method. It is anticipated that a priori information also exists in more complex simulations, and that such information can be used in an analogous way to achieve desired coverage levels.     

A hybrid prediction for wind buffeting noises of vehicle rear window based on LES-LAA method

Based on the theories of computational fluid dynamics (CFD) and aeroacoustics, a hybrid simulation technique, the so-called LES-LAA method is proposed in this paper, for predicting the wind buffeting noises generated by opening rear windows of a running vehicle. The LES-LAA is developed by combining the large eddy simulation (LES) and the Lighthill acoustic analogy (LAA) methods. Based on the established vehicle and wind tunnel models, the wind buffeting noises from rear windows are predicted by using the proposed LES-LAA method and the obtained results are compared with experimental data. The results show that the calculation error of sound pressure level (SPL) from the LES-LAA method is less than 2%, which suggests that the proposed method has good accuracy in predicting the wind noise of the rear window of a vehicle. The wind noise when both sides of the rear window are opened at the same time is much lower than the case when just one window is opened. In conclusion, the hybrid LES-LAA technique presented in this paper is effective and feasible for predicting the wind buffeting noise, which can be applied to other types of vehicle and is a promising approach for solving other aero-acoustical engineering problems.     

A random model approach for the LASSO

The least absolute selection and shrinkage operator (LASSO) is a method of estimation for linear models similar to ridge regression. It shrinks the effect estimates, potentially shrinking some to be identically zero. The amount of shrinkage is governed by a single parameter. Using a random model formulation of the LASSO, this parameter can be specified as the ratio of dispersion parameters. These parameters are estimated using an approximation to the marginal likelihood of the observed data. The observed score equations from the approximation are biased and hence are adjusted by subtracting an empirical estimate of the expected value. After estimation, the model effects can be tested (via simulation) as the distribution of the observed data given that all model effects are zero is known. Two related simulation studies are presented that show that dispersion parameter estimation results in effect estimates that are competitive with other estimation methods (including other LASSO methods).     

Adaptive Local Linear Quantile Regression

In this paper we propose a new method of local linear adaptive smoothing for nonparametric conditional quantile regression. Some theoretical properties of the procedure are investigated. Then we demonstrate the performance of the method on a simulated example and compare it with other methods. The simulation results demonstrate a reasonable performance of our method proposed especially in situations when the underlying image is piecewise linear or can be approximated by such images. Generally speaking, our ...     

Customized sequential designs for random simulation experiments: Kriging metamodeling and bootstrapping

This paper proposes a novel method to select an experimental design for interpolation in random simulation, especially discrete event simulation. (Though the paper focuses on Kriging, this design approach may also apply to other types of metamodels such as non-linear regression models and splines.) Assuming that simulation requires much computer time, it is important to select a design with a small number of observations (or simulation runs). The proposed method is therefore sequential. Its novelty is that it accounts for the specific input/output behavior (or response function) of the particular simulation at hand; i.e., the method is customized or application-driven. A tool for this customization is bootstrapping, which enables the estimation of the variances of predictions for inputs not yet simulated. The method is tested through two classic simulation models, namely the expected steady-state waiting time of the M/M/1 queuing model, and the mean costs of a terminating (s, S) inventory simulation. For these two simulation models the novel design indeed gives better results than a popular alternative design, namely Latin Hypercube Sampling (LHS) with a prefixed sample.