Genetic algorithms applied to the solution of hybrid optimal control problems in astrodynamics

Many space mission planning problems may be formulated as hybrid optimal control problems, i.e. problems that include both continuous-valued variables and categorical (binary) variables. There may be thousands to millions of possible solutions; a current practice is to pre-prune the categorical state space to limit the number of possible missions to a number that may be evaluated via total enumeration. Of course this risks pruning away the optimal solution. The method developed here avoids the need for pre-pruning by incorporating a new solution approach using nested genetic algorithms; an outer-loop genetic algorithm that optimizes the categorical variable sequence and an inner-loop genetic algorithm that can use either a shape-based approximation or a Lambert problem solver to quickly locate near-optimal solutions and return the cost to the outer-loop genetic algorithm. This solution technique is tested on three asteroid tour missions of increasing complexity and is shown to yield near-optimal, and possibly optimal, missions in many fewer evaluations than total enumeration would require.     

Software news and updates. Carma: a molecular dynamics analysis program

A computer program has been developed to aid the analysis of molecular dynamics trajectories. The program is tuned for macromolecular large-scale problems and supports features such as removal of global translations-rotations of the solute, calculation of average distance maps and their corresponding standard deviations, calculation of the variance-covariance and cross-correlation matrices, and principal component analysis of trajectories with the added ability to create artificial trajectories based on selected eigenvectors. Limited graphics (trajectory viewing) capabilities are also available.     

用准经典轨迹的理论方法研究H+ +H2反应

在新的势能面上,用准经典轨迹的方法对H H2反应体系进行了动力学研究．理论计算的结果发现,这个反应体系的前向和后向的角分布基本是对称的．同时还给出了在不同碰撞能下这个反应体系的转动态的分布情况．在碰撞能的从0．124到1．424 eV时,反应H H2的积分截面是随着碰撞能的增加而逐渐降低的．而且理论计算结果与实验结果也符合的非常好．     

Analysis of Cr atom focusing deposition properties in the double half Gaussian standing wave field

The use of the dipole force on atoms is a new technology that is used to build nanostructures. In this way, a high quality standard nano-grating can be obtained. Based on the semi-classical model, the motion equation is investigated and the trajectories of atoms in double half Gaussian standing wave field are simulated. Compared with the Gaussian standing wave field, the double half Gaussian standing wave can well focus the Cr atoms. In order to obtain this kind of beam, a prism is designed and the experimental result shows that the beam is well generated.     

THE EVOLUTION CHARACTERIZATION OF LIBRATIN POINT ORBITS AND APPLICATION RESEARCH

三体问题中, 轨道的受力和运动规律非常复杂. 对于特定的任务, 如何选择轨道的初始解是一大难题.针对平面三体问题, 利 用近拱点庞加莱映射, 对平动点顺行轨道和逆行轨道的长期和短期演化性质进行分析.根据轨道的初始状态将其分为逃逸轨道和捕获轨道.对于逃逸轨道, 给出了同宿轨道和异宿轨道的设计方法, 并利用两级微分修正法消除了拼接点处的位置不连续问题.对于捕获轨道, 得到了几类典型的周期和准周期轨道.对逆行轨道的演化性质进行分析时发现, 逆行轨道通常为准周期轨道, 比顺行轨道更加稳定.利用近拱点庞加莱映射可以快速确定不同类型轨道对应的初始状态, 为特定任务需求下的轨道设计提供了一种快速而有效的选择方案.     

Effect of reagent vibrational excitation and isotope substitution on the stereo-dynamics of the Ba+HF→BaF+H reaction

Based on an extended London-Eyring-Polanyi-Sato (LEPS) potential energy surface (PES), the Ba + HF reaction has been studied by the quasi-classical trajectory (QCT) method. The reaction integral cross section as a function of collision energy for the Ba + HF → BaF + H reaction is presented and the influence of isotope substitution on the differential cross sections (DCSs) and alignments of the product's rotational angular momentum have also been studied. The results suggest that the integral cross sections increase with increasing collision energy, and the vibrational excitation of the reagent has great influence on the DCS. In addition, the product's rotational polarization is very strong as a result of heavy-heavy-light (HHL) mass combination, and the distinct effect of isotope substitution on the stereodynamics is also revealed.     

Influence of Isotope Effects on Product Polarizations of N(2D)+D2, N(2D)+H2 and N(2D)+HD Reactive Systems

To figure out the influence of isotope effect on product polarizations of the N(²D)+D₂ reactive system and its isotope variants, quasi-classical trajectory(QCT) calculation was performed on Ho’s potential energy surface(PES) of ²A″ state. Product polarizations such as product distributions of P(θ_r), P(φ_r) and P(θ_r,φ_r), as well as the generalized polarization-dependent differential cross sections(PDDCSs) were discussed and compared in detail among the four product channels of the title reactions. Both the intermolecular and intramolecular isotope effects were proved to be influential on product polarizations.     

Effects of reagent's rotational and vibrational excitations on reaction O(3P) + H2 H2(ν = 0, 3, j = 0, 3; 5; 7; 9; 12; 15) → OH + H

To investigate the effect of reagent's rotational and vibrational excitations on the stereo-dynamics of reaction product, the title reaction is theoretically simulated using the quasi-classical trajectory (QCT) method on the ³A" and ³A′ potential energy surfaces (PESs). The reaction cross section is considered as the only scalar property in this work at four different collision energies. Furthermore the vector properties including two polarization-dependent differential cross sections (PDDCSs), the angular distributions of product' rotational momentum are discussed at one fixed collision energy. Effects of reagents' rotational excitation on the reaction do exist regularly.     

Two-phase micro- and macro-time scales in particle-laden turbulent channel flows

The micro-and macro-time scales in two-phaseturbulent channel flows are investigated using the direct numerical simulation and the Lagrangian particle trajectorymethods for the fluid-and the particle-phases,respectively.Lagrangian and Eulerian time scales of both phases are calculated using velocity correlation functions.Due to flowanisotropy,micro-time scales are not the same with the theoretical estimations in large Reynolds number(isotropic) turbulence.Lagrangian macro-time scales of particle-phaseand of fluid-phase seen by particles are both dependent onparticle Stokes number.The fluid-phase Lagrangian integral time scales increase with distance from the wall,longerthan those time scales seen by particles.The Eulerian integral macro-time scales increase in near-wall regions but decrease in out-layer regions.The moving Eulerian time scalesare also investigated and compared with Lagrangian integraltime scales,and in good agreement with previous measurements and numerical predictions.For the fluid particles themicro Eulerian time scales are longer than the Lagrangianones in the near wall regions,while away from the walls themicro Lagrangian time scales are longer.The Lagrangianintegral time scales are longer than the Eulerian ones.Theresults are useful for further understanding two-phase flowphysics and especially for constructing accurate predictionmodels of inertial particle dispersion.     

复合材料层板的抗贯穿机理与模拟研究

为了研究树脂基纤维增强复合材料层板的抗侵彻贯穿机理和动态力学行为与抗侵彻毁伤的关系, 通过球形破片模拟弹贯穿实验表征了复合材料层板抗高速侵彻的吸能特性;通过高速摄影技术分析了层板 贯穿过程的瞬态变形失效特点;采用CT扫描成像及SEM 电镜分析等手段研究了复合材料层板的抗贯穿破 坏耗能模式。实验结果显示,高速冲击下层板抗贯穿吸能与入射速度成正比;高速侵彻过程是复合材料层板 高应变率变形的动态过程,高应变率动态力学行为对复合层板抗贯穿吸能特性影响显著;冲击波在层板中的 传播特性决定了不同破坏模式阶段的划分以及损伤区域的范围。基于复合层板高速贯穿下的动力学瞬态分 析,建立了复合层板抗高速侵彻吸能的两阶段动态破坏模型,模型计算值与实验值符合良好。研究结果表明, 应变率效应与惯性效应在复合材料层板抗侵彻性能分析中是不可忽视的2个关键因素。