AN EFFICIENT METHOD FOR MULTIOBJECTIVE OPTIMAL CONTROL AND OPTIMAL CONTROL SUBJECT TO INTEGRAL CONSTRAINTS

We introduce a new and efficient numerical method for multicriterion optimal control and single criterion optimal control under integral constraints. The approach is based on extending the state space to include information on a ＂budget＂ remaining to satisfy each constraint; the augmented Hamilton-Jacobi-Bellman PDE is then solved numerically. The efficiency of our approach hinges on the causality in that PDE, i.e., the monotonicity of characteristic curves in one of the newly added dimensions. A semi-Lagrangian ＂marching＂ method is used to approximate the discontinuous viscosity solution efficiently. We compare this to a recently introduced ＂weighted sum＂ based algorithm for the same problem [25]. We illustrate our method using examples from flight path planning and robotic navigation in the presence of friendly and adversarial observers.     

燃气射流起始冲击波场形态的实验研究

介绍利用激光摩尔偏折仪和压阻式压力传感器对发动机喷流起始冲击波的基本形态进行的研究． 研究结果验证了在燃气流团前存在着起始冲击波及发展初期具有非球面的特性,并提出了喷口起始冲击波波阵面发展的新结构,对前人有关起始冲击波发展结构的研究结果进行了修正．     

用铁塔资料研究地面冷锋的结构

利用铁塔的梯度观测资料研究夜间干冷锋过境进时的边界层结构和部分地面锋面和特征。资料表明近地面，锋面与地表面近于垂直，同附近垂直运动很强，只比水平风速小一个量级，表征湍流变换强度的总体拖是曳系数ＣＤ和热量输送系数ＣＨ在锋面过境前后有明显的响应。并存在着明显的日变化，陆面上白天的ＣＤ和ＣＨ值较海面上大一个量级。     

螺旋桨运输机进气道阻力特性分析

以某型螺旋桨运输机带动力装置飞行试验台为研究对象，建立了三维流场模型，采用CFD进行了网格划分及数值计算。以前期飞行试验数据和发动机热力循环模型计算结果为边界条件，开展了不同飞行马赫数、高度及桨叶角的数值计算。选用桨叶角表征发动机状态，分析了不同飞行高度、马赫数、螺旋滑流对进气道前罩作用力和附加前体力影响的变化规律。计算结果表明：不同高度、桨叶角及马赫数均会改变进气道阻力，在马赫数0.5，高度4 km状态下，桨叶角由40°增大至50°，前罩作用力减小1 700 N，附加前体力减小3 600 N,而附加前体力系数及前罩作用力系数随飞行高度基本不变。进气道阻力特性研究为后续试飞中确定涡桨发动机外部阻力特性分析提供了数据支持。     

纳米多晶铜中冲击波阵面的分子动力学研究

冲击波阵面反映材料在冲击压缩下的弹塑性变形行为以及屈服强度、应变率条件等宏观量, 还与冲击压缩后的强度变化联系. 本文使用分子动力学方法, 模拟研究了冲击压缩下纳米多晶铜中的动态塑性变形过程, 考察了冲击波阵面和弹塑性机理对晶界存在的依赖, 并与纳米多晶铝的冲击压缩进行了比较. 研究发现: 相比晶界对纳米多晶铝的贡献而言, 纳米多晶铜中晶界对冲击波阵面宽度的影响较小; 并且其塑性变形机理主要以不全位错的发射和传播为主, 很少观察到全位错和形变孪晶的出现. 模拟还发现纳米多晶铜的冲击波阵面宽度随着冲击应力的增加而减小, 并得到了冲击波阵面宽度与冲击应力之间的定量反比关系, 该定量关系与他人纳米多晶铜模拟结果相近, 而与粗晶铜的冲击压缩实验结果相差较大.     

Estimation of the degree of asphericity of a glass sphere using a vectorial shearing interferometer

The degree of asphericity is estimated by determining the average radius of curvature in different sections, at various points on the surface of a sphere, and the deviation from it. We employ the vectorial shearing interferometer (VSI) as the instrument to determine the radius of curvature from two subapertures of the transparent glass sphere. We incorporate the sphere as a thick lens into the interferometric setup, illuminating it with an expanded beam. The spherical aberration, introduced by the sphere in the wave front, depends on the local sphere radius, on the refraction index of the glass, and on the cone angle of the source. The wave front aberrated by the sphere impinges on the VSI. Here, the wave front is divided in two in amplitude, it is sheared vectorially, and it is superimposed with itself. The fringe pattern is formed in the intersection of the wave fronts. The shape of the resulting fringe pattern is directly related to spherical aberration. We estimate qualitatively the degree of asphericity, comparing the phase gradients in different sections of the sphere. Here, we report on the experimental setup to test the asphericity, the results with different vectorial shearing (magnitude and direction). Finally, we perform a comparison with the theoretical predictions.     

基于元胞蚁群的主动队列管理算法研究

针对网络拥塞现象,基于弃头方式提出了一种新的主动队列管理算法.该算法首先利用元胞蚁群建立了实际网络队长最大值的计算方法,同时通过判断网络队长与阈值的关系,采取从队列头部丢弃N个数据包的方法.最后,以长相关数据进行仿真实验,对比分析了DFCA与RED、DROP-TAIL之间的优劣,结果表明该算法具有较好的适应性.     

一种特殊平面装药爆轰波阵面的摄影与结果分析

对由导爆索构成的平面装药的爆轰波阵面进行了高速摄影,并将其摄影照片同理论作图做了对比,结果表明此理论作图的方法是正确的。     

分岔管道内瓦斯爆炸火焰传播规律实验及数值模拟

为了研究煤矿井下瓦斯爆炸火焰在分岔巷道内的传播规律,自制45°分岔管道实验装置开展甲烷体积分数为9.5%的瓦斯爆炸火焰传播实验,用Fluent 16.0软件模拟分岔管道内瓦斯爆炸火焰传播过程。对比分析实验数据与模拟结果,得到分岔管道瓦斯爆炸火焰传播的变化规律。研究结果表明:1)分岔管道内瓦斯爆炸火焰在分岔处产生漩涡,加速管道内爆炸火焰湍流化,火焰冲击反射现象明显;2)分岔支管截面处爆炸火焰温度、传播速度、冲击波超压与离子电流峰值最大;3)瓦斯爆炸火焰传播的模拟结果与实验数据在数值上存在一定差异,但各参量总体变化趋势相同。研究结果为深入认识井下瓦斯爆炸传播机制和在巷道分岔处采取瓦斯爆炸火焰传播抑制措施提供一定参考。     

一类含时滞反应扩散方程波前解的存在性

利用J.Wu和X.Zou(J.Dynam.Diff.Eqns.,2001,13(3):651～687.)建立的解的存在性理论,研究 2u1(x,t) u1(x,t) t=D1b1+a1u2(x,t-τ2)], x2+r1u1(x,t)[1-u1(x,t-τ1) u2(x,t) 2u2(x,t) t=D2b2+a2u1(x,t-τ4)], x2+r2u2(x,t)[1-u2(x,t-τ3)的行波解,其中x∈R,t∈R,ui(x,t)∈R,Di>0,ri>0,ai>0,bi>0,i=1,2,a1a2<1,τj>0,j=1,2,3,4,得到了这个系统波前解存在的充分条件.