分区拼接网格算法数值模拟超声速复杂流场

以分区拼接网格数值模拟方法为研究对象,对二阶Godunov方法在拼接网格中流场的数值模拟进行了研究,发展了适用于Godunov格式的通量守恒型算法,结合二阶Godunov有限体积法离散非定常Euler方程,数值模拟了捆绑火箭及航天飞机的超声速流场,计算结果正确描述了流场中的激波相交、反射等干扰特性.     

含激波跨声速S_2流面反问题的椭圆型方程间断解

在充分的理论依据之上,通过找到叶轮机械S_2流面反问题运动方程的适当守恒形式,把激波关系嵌入主方程,提出了求子午面速度为亚声速的含激波跨声速S_2反问题的椭圆型方程间断解的数值方法.由于方程是统一的椭圆型的,算法简单可靠.编制了计算机程序并试算了例题.本方法可用于与S_1“激波拟合——分区计算”的跨声速三元流动迭代求解方法中.     

基于熵条件二阶差分格式的嵌套网格分区算法

以带熵条件的二阶TVD格式为基本格式,设计了三维曲线坐标下的Euler方程数值计算方法.对多个物体使用分区嵌套网格设计了相应的分区算法.所建立的算法适用于处理物体相互分离以及相互接触的定常与非定常流动问题.给出了绕两圆柱的超声速流动在两圆柱处于不同位置时的详细结果,并模拟了装配鸭舵火箭弹在超声速旋转飞行时的气动特性.     

基于非结构变形网格的间断装配法原理

在激波捕捉法计算得到的流场基础上采用辨识算法得到初始间断位置, 从ALE方程出发, 考虑离散几何守恒律, 采用变形网格和网格重构技术解决计算过程中间断运动和变形, 新旧网格之间流场采用高精度信息传递方法保持时间精度, 建立了基于非结构动网格技术的间断装配方法.通过激波管问题的二维模拟, 模拟了初始间断分解为激波和接触间断激波遇到固壁反射后与接触间断相交的非定常流动过程, 对这种新方法的基本原理进行了介绍.      

流体力学分区算法及相关理论问题

分区计算是计算流体力学中处理复杂几何区域、多块物体相互运动和进行并行计算的一类重要方法.对分区算法、相关理论问题和一些数值现象进行了综述.     

流动/运动松耦合与紧耦合计算方法及稳定性分析

为了更加精确地模拟流动/运动耦合问题, 建立了耦合动态混合网格生成非定常流场计算和六自由度运动方程求解的一体化计算方法, 并在统一框架内同时实现了松耦合与紧耦合方法.通过圆柱涡致自激振荡(vortex induced vibration, VIV)的模拟, 对不同时间精度的松耦合和紧耦合算法的优劣及适用范围进行了评估和分析; 通过引入附加质量的概念, 对耦合算法的稳定性进行了理论分析.研究表明:在流体的密度与物体的密度接近时, 松耦合方法是不稳定的, 必须采用紧耦合方法.最后利用耦合算法对二维鱼体的自主游动和钝锥三自由度自由飞过程进行了数值模拟, 证实了理论分析的结论.      

拉氏数值模拟中的人工粘性方法

较好的人工粘性需要满足较小的计算开销、不能去除真实具有的涡运动等条件.提出一种应用于拉氏数值模拟中基于Lew人工粘性,同时增加了限制器的人工粘性方法.可以有效减少数值模拟结果对网格的依赖;采用特征值限制器控制施加的人工粘性大小,通过限制器能够区分激波压缩和等熵压缩;方便应用在二维、三维,结构网格或者非结构网格上.     

守恒律方程组的大时间步长叠波格式

大时间步长叠波格式最初思想为LeVeque提出的大时间步长Godunov格式,通过叠加间断分解发出的强波来构造数值格式.原方法只给出了间断强波的穿越叠加方法,文章对其进行了完善,并推广到多维.针对膨胀波提出了一种网格单元分解法可以自动满足熵条件,避免出现非物理解.给出了格式的具体计算公式,并用单个守恒律方程、一维/多维Euler方程组进行了数值计算.计算结果表明,新格式除了可以采用大时间步长的优点外,在一定范围内随CFL数增加其耗散反而更低,因而对激波接触间断膨胀波的分辨率更高.     

非结构多边形网格解扭和重分方法

在多维流体动力学计算中，流体运动和计算网格的关系可以分为两种情况。一是Lagrangian方法，即网格跟随流体运动；二是Eulerian方法，即流体流过固定；下动的网格。一般计算网格的运动是任意的。这就对应于任意Lagrangian—Eulerian（ALE）方法。ALE方法的核心是通过调整网格运动，使得数值模拟的精度、效率有所提高。它的主要步骤是：显式Lagrangian步；网格重分，即得到新的计算网格；物理量重映，即将Lagrangian步的计算结果变换到新网格上。在这3步中，较少研究网格重分。数值模拟和网格重分的一个基本前提是网格是合理的，或者说网格不能发生翻转，网格应当是凸的。而Lagrangian步数值模拟会造成网格扭曲，因此在网格重分前进行网格解扭是十分必要的。文中描述了通用的网格解扭、重分算法，使得解扭、重分后的网格有较好的几何品质，同时尽可能接近Lagrangian网格。     

激波和单列水柱、气水固多界面相互作用的数值模拟

研究可压缩多介质流场的激波和多介质界面相互作用问题.在Descartes固定网格采用level-set方法追踪界面,气/气界面边界条件处理采用OGFM方法,采用修正的rGFM方法提高气/水和气/固界面处构造Riemann问题精度,将Riemann近似解得到的界面参数外推到两侧真实和虚拟流体,采用五阶WENO方法求解流场Euler方程和界面level-set方程,给出不同时刻流场数值纹影图像.结果表明：在可压缩流场嵌入固体和水、气体等目标,本文方法可较精确地分辨平面运动激波和单列水柱及包含气/气、气/水和气/固等界面作用后产生的复杂激波结构.和传统的分区与贴体变换方法不同,为Descartes网格包含多介质界面复杂流场计算提供新途径.     

考虑横向分离与超车期望的车辆跟驰模型

为了更加客观地描述实际的车辆跟驰行为, 在优化速度模型的基础上, 通过引入横向分离参数并提出超车期望和虚拟前车的概念, 建立了考虑横向分离与超车期望的车辆跟驰模型.对模型进行线性稳定性分析, 得到了模型稳定性条件, 发现车辆横向分离、超车期望和虚拟前车的位置的增加, 在车流密度较小、车速较快的情况下, 使得交通流稳定区域增大, 但在车流密度较大、车速较慢的情况下, 反而使得交通流稳定区域减小.数值模拟结果验证了模型稳定性分析的结果, 表明在交通瓶颈处等交通流密度较大、运行缓慢的区域, 为抑制交通拥堵, 应该限制车辆的横向偏移和超车行为的发生. 关键词： 交通流 车辆跟驰模型 横向分离 超车期望     

高压下黄铜中追赶稀疏波速度的辐射测量

 用光辐射法、以有机玻璃作窗口测量了143～182 GPa下H62黄铜的追赶稀疏波速度。给出了在对称碰撞下，用透明材料作窗口时，实验确定追赶稀疏波速度的公式。实验结果表明，H62黄铜的冲击熔化点在143～162 GPa之间。     

A cell-based study on pedestrian acceleration and overtaking in a transfer station corridor

Pedestrian speed in a transfer station corridor is faster than usual and sometimes running can be found among some of them. In this paper, pedestrians are divided into two categories. The first one is aggressive, and the other is conservative. Aggressive pedestrians weaving their way through crowd in the corridor are the study object of this paper. During recent decades, much attention has been paid to the pedestrians’ behavior, such as overtaking (also deceleration) and collision avoidance, and that continues in this paper. After sufficiently analyzing the characteristics of pedestrian flow in transfer station corridor, a cell-based model is presented in this paper, including the acceleration (also deceleration) and overtaking analysis. Acceleration (also deceleration) in a corridor is fixed according to Newton’s Law and then speed calculated with a kinematic formula is discretized into cells based on the fuzzy logic. After the speed is updated, overtaking is analyzed based on updated speed and force explicitly, compared to rule-based models, which herein we call implicit ones. During the analysis of overtaking, a threshold value to determine the overtaking direction is introduced. Actually, model in this paper is a two-step one. The first step is to update speed, which is the cells the pedestrian can move in one time interval and the other is to analyze the overtaking. Finally, a comparison between the rule-based cellular automata, the model in this paper and data in HCM 2000 is made to demonstrate our model can be used to achieve reasonable simulation of acceleration (also deceleration) and overtaking among pedestrians.     

轴对称环状静电模的漂移波湍流参量激发理论研究

本文所论述的轴对称环状静电模是指环形磁约束等离子体(如托卡马克)中环向模数为零的近理想静电流体模,它包含有测地声模和基频率与之较低的声模;也含有所谓的‘近零频带状流’.本文根据冷离子流体模型在圆形磁面构成的准环坐标系中的表示,对涉及以上三种模式的漂移波湍流参量激发理论,在一级环形效应近似下,进行了系统讨论,并证明了带状流的四个新命题.利用对漂移波能谱的参数化描写,注意到由漂移波能谱径向有限宽度所引发的特性,如波能传播量的双Landau奇点,揭示了有限宽度对径向δ谱所得结果的重要修正:如,对近零频带状流和测地声模的参量激发条件带来的严格限制.此外,还讨论了密度带状流在高q条件下被激发的可能性.本文选用合理的物理参数.采用图示方法详细地讨论了有关的数值结果.分析表明,测地声模和近零频带状流的参量激发不可能发生在同一小半径处;如果测地声模被参量激发,也应能观察到密度带状流.     

Nonlinear interactions between drift waves and zonal flows

Phase coherent interactions between drift waves and zonal flows are considered. For this purpose, mode coupling equations are derived by using a two-fluid model and the guiding center drifts. The equations are then Fourier analyzed to deduce the nonlinear dispersion relations. The latter depict the excitation of zonal flows due to the ponderomotive forces of drift waves. The flute-like zonal flows with insignificant density fluctuations have faster growth rates than those which have a finite wavelength along the magnetic field direction. The relevance of our investigation to drift wave driven zonal flows in computer simulations and laboratory plasmas is discussed. Received 5 April 2002 Published online 28 June 2002     

离子温度梯度模湍流的带状流最小自由度模型

在离子温度梯度模(ITG)湍流背景中,通过最小自由度模型中模耦合方式产生带状流,对此模型做了动力学稳定性分析及数值求解.并在此基础上初步探讨了湍流中漂移波与带状流的能量转移,以及雷诺协强与带状流的关系. 关键词： 等离子体 离子温度梯度模 湍流 带状流     

An optimization-based overtaking model for unidirectional pedestrian flow

We propose an optimization-based model for simulating the overtaking behaviour in the unidirectional pedestrian flow. A ‘visual area’ is introduced so that agents could receive the information regarding their surroundings and react by choosing one of three options: to move straight on, to dodge to the left, or to dodge to the right. And a side preference of each pedestrian for evading and overtaking is implemented based on traffic ‘social norms’. The model was validated by reproducing the experimentally obtained pedestrian flow patterns. The effects of the initial pedestrian formation on overtaking behaviour and the evacuation time have been analysed in different geometries. The results show that pedestrian flow patterns after overtaking are obviously influenced by both the initial positions and density of the slow pedestrians in the front. Phase changes of pedestrian formation are observed in both experiment and simulations. On the other hand, for sparse pedestrian crowds, the egress time of the fast individuals is mainly impacted by the horizontal distance between the initial positions of the slow pedestrians in the front, especially in the geometry with a bottleneck.     

Conservation laws,bilinear forms and solitons for a fifth-order nonlinear Schrödinger equation for the attosecond pulses in an optical fiber

Under investigation in this paper is a fifth-order nonlinear Schrödinger equation, which describes the propagation of attosecond pulses in an optical fiber. Based on the Lax pair, infinitely-many conservation laws are derived. With the aid of auxiliary functions, bilinear forms, one-, two- and three-soliton solutions in analytic forms are generated via the Hirota method and symbolic computation. Soliton velocity varies linearly with the coefficients of the high-order terms. Head-on interaction between the bidirectional two solitons and overtaking interaction between the unidirectional two solitons as well as the bound state are depicted. For the interactions among the three solitons, two head-on and one overtaking interactions, three overtaking interactions, an interaction between a bound state and a single soliton and the bound state are displayed. Graphical analysis shows that the interactions between the two solitons are elastic, and interactions among the three solitons are pairwise elastic. Stability analysis yields the modulation instability condition for the soliton solutions.     

一种可用于分析任意介电常数张量的各向异性波导导模的二维时域有限差分模型

为分析具有任意介电常数张量的各向异性波导的导模，本文通过把简化的二维时域有限差分（２－ＤＦＤＴＤ）法扩展至任意各向异性介质，提出了一种以Ｄ、Ｅ和Ｈ场为基础的统一的简化２－ＤＦＤＴＤ模型．利用该模型，研究了简化的复数２－ＤＦＤＴＤ方法与实变数２－ＤＦＤＴＤ方法之间的关系．文中还讨论了复变数方法和实变数方法的激励技术．     

Experiments on Ion-Acoustic Solitons in Plasmas Invited Review Article

Experiments on ion-acoustic solitons are reviewed. Theories and numerical simulations which are relevant to experimental results are also presented. The measured velocity and width of planar solitons are compared with the predictions of the Korteweg-deVries (KdV) equation which includes a finite ion temperature. The spatial evolution of compressive or rarefactive pulses is discussed. Cylindrical and spherical solitons are introduced together with a modified KdV equation and numerical results. Oblique collisions and their evolution, including overtaking and head-on collisions of two planar, cylindrical, or spherical solitons are described. Reflection, diffraction, and other topics related with ion-acoustic solitons are presented.