Conservative numerical methods for model kinetic equations

A new conservative discrete ordinate method for nonlinear model kinetic equations is proposed. The conservation property with respect to the collision integral is achieved by satisfying at the discrete level approximation conditions used in deriving the model collision integrals. Additionally to the conservation property, the method ensures the correct approximation of the heat fluxes. Numerical examples of flows with large gradients are provided for the Shakhov and Rykov model kinetic equations.     

Gas kinetic algorithm using Boltzmann model equation

Based on the Boltzmann-BGK model equation, the unified simplified velocity distribution function equation adapted to various flow regimes can be presented by the aid of the basic characteristic on molecular movement and collision approaching to equilibrium. The optimum Golden Section principle is extended and applied to the discrete velocity ordinate method in order to discretize the corresponding velocity components, and then the molecular velocity distribution function equation will be cast into hyperbolic conservation laws form with non-linear source terms. In view of the unsteady characteristic of molecular convective movement and colliding relaxation, the time-splitting method is applied to decompose the velocity distribution function equations into the colliding relaxation equations with non-linear source terms and the convective motion equations. Based on the second-order Runge-Kutta method and the non-oscillatory, containing no free parameters, and dissipative (NND) finite difference method, the gas kinetic finite difference second-order scheme is constructed to directly solve the discrete velocity distribution functions. Four types of discrete velocity quadrature rules, such as the modified Gauss-Hermite formula and the Golden Section number-theoretic integral method, are developed and applied to evaluate the macroscopic flow moments of the distribution functions over the velocity space. As a result, a unified gas kinetic algorithm is established for the flows from rarefied transition to continuum regime. To test the reliability of the present method, the one-dimensional shock wave structures, the flows past two-dimensional circular cylinder and the three-dimensional flows over sphere with various Knudsen numbers are simulated. The computational results are found in high resolution of the flow fields and good agreement with the theoretical, DSMC, and experimental results.     

Gas-kinetic numerical study of complex flow problems covering various flow regimes

The Boltzmann simplified velocity distribution function equation, as adapted to various flow regimes, is described on the basis of the Boltzmann–Shakhov model from the kinetic theory of gases in this study. The discrete velocity ordinate method of gas-kinetic theory is studied and applied to simulate complex multi-scale flows. On the basis of using the uncoupling technique on molecular movements and collisions in the DSMC method, the gas-kinetic finite difference scheme is constructed by extending and applying the unsteady time-splitting method from computational fluid dynamics, which directly solves the discrete velocity distribution functions. The Gauss-type discrete velocity numerical quadrature technique for flows with different Mach numbers is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established for studying the three-dimensional complex flows with high Mach numbers from rarefied transition to continuum regimes. On the basis of the parallel characteristics of the respective independent discrete velocity points in the discretized velocity space, a parallel strategy suitable for the gas-kinetic numerical method is investigated and, then, the HPF (High Performance Fortran) parallel programming software is developed for simulating gas dynamical problems covering the full spectrum of flow regimes. To illustrate the feasibility of the present gas-kinetic numerical method and simulate gas transport phenomena covering various flow regimes, the gas flows around three-dimensional spheres and spacecraft-like shapes with different Knudsen numbers and Mach numbers are investigated to validate the accuracy of the numerical methods through HPF parallel computing. The computational results determine the flow fields in high resolution and agree well with the theoretical and experimental data. This computing, in practice, has confirmed that the present gas-kinetic algorithm probably provides a promising approach for resolving hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view for solving the mesoscopic Boltzmann model equation.     

Parallel three-dimensional direct simulation Monte Carlo method and its applications

The development of a parallel three-dimensional direct simulation Monte Carlo (DSMC) method using unstructured cells is reported. Variable hard sphere molecular model and no time counter method are used for the molecular collision kinetics, while the cell-by-cell ray-tracing technique is implemented for particle movement. Developed serial code has been verified by comparing the results of a supersonic corner flow with those of Bird’s three-dimensional structured DSMC code. In addition, a benchmark test is performed for an orifice expanding flow to verify the parallel implementation of DSMC method by comparing with available experimental data. Static physical domain decomposition is used to distribute the workload among multiple processors by considering the estimated particle weighting distribution. Two-step multi-level graph partitioning technique is used to perform the required domain decomposition. Completed code is then applied to compute a hypersonic flow over a sphere (external flow) and the flow field of a spiral drag pump (internal flow), respectively. Results of the former are in good agreement with previous numerical results using axisymmetric DSMC method and experimental data. Results of the latter also agree well with previous numerical results.     

光催化降解废气过程中紫外辐射照度模拟和优化

为在光催化设计中合理布置紫外灯,减少光在传输过程中的损失,概括总结径向辐射模型、多点源叠加近似辐射模型和离散坐标辐射模型等3种经典的紫外辐射物理模型;基于离散坐标辐射模型模拟某公司生产的光催化氧化废气处理装置的紫外辐射照度;通过调整灯管布置距离改变光强分布,并拟合光催化单元和光催化网表面的平均紫外辐射照度随灯管布置距离变化曲线.结果表明,距离参数在55～65 mm时,该光催化反应装置的光能利用率较好;紫外灯的合理布置可有效提高光能利用率.     

A Numerical Method for Simulation of Microflows by Solving Directly Kinetic Equations with WENO Schemes

A numerical method for simulation of transitional-regime gas flows in microdevices is presented. The method is based on solving relaxation-type kinetic equations using high-order shock capturing weighted essentially non-oscillatory (WENO) schemes in the coordinate space and the discrete ordinate techniques in the velocity space. In contrast to the direct simulation Monte Carlo (DSMC) method, this approach is not subject to statistical scattering and is equally efficient when simulating both steady and unsteady flows. The presented numerical method is used to simulate some classical problems of rarefied gas dynamics as well as some microflows of practical interest, namely shock wave propagation in a microchannel and steady and unsteady flows in a supersonic micronozzle. Computational results are compared with Navier–Stokes and DSMC solutions.     

Implementation of a transient adaptive sub-cell module for the parallel-DSMC code using unstructured grids

A method of transient adaptive sub-cells (TAS) suitable for unstructured grids that is modified from the existing one for the structured grids of DSMC is introduced. The TAS algorithm is implemented within the framework of a parallelized DSMC code (PDSC). Benchmarking tests are conducted for steady driven cavity flow, steady hypersonic flow over a two-dimensional cylinder, steady hypersonic flow over a cylinder/flare and the unsteady vortex shedding behind a two-dimensional cylinder. The use of TAS enables a reduction in the computational expense of the simulation since larger sampling cells and less simulation particles can be employed. Furthermore, the collision quality of the simulation is maintained or improved and the preservation of property gradients and vorticity at the scale of the sub-cells enables correct unsteady vortex shedding frequencies to be predicted. The use of TAS in a parallel-DSMC code allows simulations of unsteady processes at a level to be carried out efficiently, accurately and with acceptable computational time.     

A Hierarchical Grid Based Framework for Fast Collision Detection

We present a novel hierarchical grid based method for fast collision detection (CD) for deformable models on GPU architecture. A two‐level grid is employed to accommodate the non‐uniform distribution of practical scene geometry. A bottom‐to‐top method is implemented to assign the triangles into the hierarchical grid without any iteration while a deferred scheme is introduced to efficiently update the data structure. To address the issue of load balancing, which greatly influences the performance in SIMD parallelism, a propagation scheme which utilizes a parallel scan and a segmented scan is presented, distributing workloads evenly across all concurrent threads. The proposed method supports both discrete collision detection (DCD) and continuous collision detection (CCD) with self‐collision. Some typical benchmarks are tested to verify the effectiveness of our method. The results highlight our speedups over prior algorithms on different commodity GPUs.     

Interactive continuous collision detection for non-convex polyhedra

We present a highly interactive, continuous collision detection algorithm for rigid, general polyhedra. Given initial and final configurations of a moving polyhedral model, our algorithm creates a continuous motion with constant translational and angular velocities, thereby interpolating the initial and final configurations of the model. Then, our algorithm reports whether the model under the interpolated motion collides with other rigid polyhedral models in environments, and if it does, the algorithm reports its first time of contact (TOC) with the environment as well as its associated contact features at TOC. Our algorithm is a generalization of conservative advancement [19] to general polyhedra. In this approach, we calculate the motion bound of a moving polyhedral model and estimate the TOC based on this bound, and advance the model by the current TOC estimate. We iterate this process until the inter-distance between the moving model and the other objects in the environments is below a user-defined distance threshold. We pose the problem of calculating the motion bound as a linear programming problem and provide an efficient, novel solution based on the simplex method. Moreover, we also provide a hierarchical advancement technique based on a bounding volume traversal tree to generalize the conservative advancement for non-convex models. Our algorithm is relatively simple to implement and has very small computational overhead of merely performing discrete collision detection multiple times. We extensively benchmarked our algorithm in different scenarios, and in comparison to other known continuous collision detection algorithms, the performance improvement ranges by a factor of 1.4～45.5 depending on benchmarking scenarios. Moreover, our algorithm can perform CCD at 120～15460 frames per second on a 3.6 GHz Pentium 4 PC for complex models consisting of 10K～70K triangles.     

碰撞C20-C20富勒烯分子的聚合与分裂

针对碳富勒烯制备过程中普遍存在的碳原子团簇与碳原子团簇的相撞、聚合与生长现象,采用分子动力学与量子力学相结合的方法,对C20富勒烯分子之间碰撞后的聚合与分裂过程进行了理论模拟。模拟结果表明,(1)当碰撞速度较小时,分子间发生弹性碰撞;当碰撞速度较高时,分子间发生聚合;当碰撞速度非常高时,分子间碰撞后发生分裂。(2)C20分子碰撞后的聚合与分裂行为,不仅与分子之间的初始碰撞速度有关,而且还与分子的碰撞形式有关。     

土壤碰撞问题的刚-散耦合动力学分析

本文以探测器着陆行星土壤为背景,对土壤碰撞问题进行刚-散耦合动力学建模与仿真分析研究.结合离散元方法和多体动力学方法,对半球壳装置土壤跌落问题进行耦合动力学仿真.通过与实验结果及有限元仿真结果对比,验证所采用离散元方法的有效性.分析了颗粒场中颗粒尺寸、恢复系数、静摩擦系数等参数,对碰撞中物体和颗粒场的碰撞加速度、碰撞持续时间、振动波形等动力学响应的影响.本研究将拓展对刚-散耦合动力学问题的理论认识,为探测器着陆系统的设计提供技术支持.     

基于冲突分类模型的冲突解析算法

IEEE 802.11 网络中的传统退避算法利用竞争窗口机制进行冲突解析,冲突节点被重新分布在一个更大的窗口范围内以避免冲突.然而,只要这些分布窗口之间存在着交集,就仍有可能引发冲突.为了解决一个问题,提出了一种冲突分类模型,将网络中的分组冲突分成交叉冲突和同级冲突,并且提出针对这两种不同类型的冲突应该采取不同的策略进行解析.利用顺序离散窗口分布机制(sequential discrete window distribution,简称SDWD)解析交叉冲突,通过节点分布窗口的离散化避免交叉冲突;同时,通过     

高超声速飞行器基于Back-stepping的离散控制器设计

根据高超声速飞行器纵向模型的特点,提出了基于Back-stepping的离散控制器设计方法.首先通过合理的简化,将飞行器的模型转化为连续非线性系统的严格反馈形式;然后采用欧拉法得到其近似的离散模型,根据近似离散模型并结合Back-stepping和反馈线性化方法,设计了高超声速飞行器的离散控制器.利用高超声速飞行器的纵向模型对算法进行仿真验证,得到了较为满意的控制效果.     

A model for retrieval of surface spectral reflectance from satellite radiance measurements using realistic atmospheric aerosol profiles

Abstract

An effort has been made to retrieve surface spectral reflectance from satellite radiance measurements. For this purpose, a simple cubic equation with surface reflectance as one of its roots has been derived. Besides satellite measured radiance, the other required quantities in the equation are calculated by using the actual solution of radiative transfer equation by a discrete ordinate method, taking into account the inhomogeneous aerosol distribution. The method has been tested for internal consistency by numerical simulations using realistic aerosol profiles. Finally, it has been tested with a Landsat image, giving good agreement with ground observations.     

交互式乐器演奏的六自由度力觉渲染方法

目的在进行虚拟乐器交互演奏时,需要模拟触力觉-视觉-听觉多通道同步反馈,其中触力觉反馈的难点在于模拟人手操作乐器的六自由度(6-Do F)力觉交互过程。方法提出一种基于混合模型和单边约束优化的六自由度力觉合成方法,实现了虚拟人手和琴弦的多点多区域接触力觉模拟。虚拟人手采用层次化球树模型表达,古琴采用混合模型表达,其中琴体和琴弦分别采用层次化球树模型和直线解析模型。提出了基于混合模型的离散碰撞检测算法,实时检测虚拟手和琴弦是否产生碰撞;基于发生碰撞的几何元素对建立单边不可穿透约束方程,通过Active Set方法求解约束优化后方程,获得6维位姿变量保证图形显示场景中的虚拟手不会和琴弦产生穿透。为模拟琴弦变形,提出变直径的圆柱体模型来模拟琴弦在不同振动幅度下的动力学响应;提出交互状态敏感的力计算模型以刻画人手在弹奏不同状态琴弦(静态、振动态)的力觉感受差异。结果基于力觉交互设备Phantom Premium 3.0建立了实验平台,实验结果表明,本文算法可以模拟单点、多点等不同接触状态,并能模拟6维力和力矩,操作者可以感受到琴弦振动时的细腻力感觉,力觉交互过程稳定,算法计算效率在1 k Hz以上。结论算法可模拟针对琴弦一类的超薄形状物体的多点接触力觉交互过程,算法计算效率高,包含碰撞检测、约束优化、琴弦变形仿真等计算回路的更新频率也能达到要求,该混合模型能为后续复杂形状物体的碰撞响应研究提供思路。     

Generalization of the collision cone approach for motion safety in 3-D environments

Avoidance of collision between moving objects in a 3-D environment is fundamental to the problem of planning safe trajectories in dynamic environments. This problem appears in several diverse fields including robotics, air vehicles, underwater vehicles and computer animation. Most of the existing literature on collision prediction assumes objects to be modelled as spheres. While the conservative spherical bounding box is valid in many cases, in many other cases, where objects operate in close proximity, a less conservative approach, that allows objects to be modelled using analytic surfaces that closely mimic the shape of the object, is more desirable. In this paper, a collision cone approach (previously developed only for objects moving on a plane) is used to determine collision between objects, moving in 3-D space, whose shapes can be modelled by general quadric surfaces. Exact collision conditions for such quadric surfaces are obtained and used to derive dynamic inversion based avoidance strategies.     

Multi-core-CPU and GPU-accelerated radiative transfer models based on the discrete ordinate method

The operational processing of remote sensing data usually requires high-performance radiative transfer model (RTM) simulations. To date, multi-core CPUs and also Graphical Processing Units (GPUs) have been used for highly intensive parallel computations. In this paper, we have compared multi-core and GPU implementations of an RTM based on the discrete ordinate solution method. To implement GPUs, the original CPU code has been redesigned using the C-oriented Compute Unified Device Architecture (CUDA) developed by NVIDIA.     

A second-order TVD implicit–explicit finite volume method for time-dependent convection-reaction equations

A class of conservative methods is developed in the more general framework of cell-centered upwind differences to approximate numerically the solution of one-dimensional non-linear conservation laws with (possibly) stiff reaction source terms. These methods are based on a non-oscillatory piecewise linear polynomial representation of the discrete solution within any mesh interval to compute pointwise solution values. The piecewise linear approximate solution is obtained by approximating the cell average of the analytical solution and the solution slope in every mesh cell. These two quantities are evolved in time by solving a set of discrete equations that are suitably designed to ensure formal second-order consistency. Several numerical tests which are taken from literature illustrate the performance of the method in solving non-stiff and stiff convection-reaction equations in conservative form.     

Stabilization of Discrete‐time Switched Systems with State Constraints Based on Mode‐Dependent Average Dwell Time

This paper considers the stabilization problem for a class of switched systems with state constraints based on mode‐dependent average dwell time (MDADT) in discrete‐time context. An improved average dwell time method is proposed, which is less conservative than the common average dwell time method. The sufficient conditions and stabilizing state feedback controllers for stabilization of discrete‐time switched systems with state constraints under MDADT switching are derived. Finally, the simulation results show that the approach designed by this paper is effective.     

Multi-layer structure and stability of discrete large-scale systems

The structure of a discrete large-scale system is analyzed. It is shown that a discrete large-scale system can be considered with a multi-layer structure, each layer of which can be considered as a minimal strongly connected system (MSCS), except the highest layer which is a simply connected system. In studying the stability, a direct method is firstly presented for MSCS. Analysis and examples show that the stability criteria given here for MSCS are less conservative than those in other literature. When the isolated subsystems in each layer (which are also MSCS) are stable, the large-scale system is stable.