非定常欠膨胀射流的正格式数值模拟

 采用二阶正格式方法对非定常欠膨胀射流进行了数值模拟。将二维守恒方程的正格式方法推广到轴对称Euler方程组的求解,并对不同马赫数下的燃气射流进行了数值计算。计算结果表明,该方法能够较好地捕捉到包含膨胀波、入射激波、反射激波、马赫盘、射流边界以及三波点等复杂射流流场的波系结构,与实验照片反映的流动特征以及已有的数值结果相吻合。表明该方法对间断解具有较强的捕捉能力,在激波阵面上不会出现数值振荡。     

对运动斜激波绕射波纹壁面的研究

使用计算流体动力学的方法,对经典的运动激波绕射现象做数值模拟,研究了一类复杂激波反射问题一入射的运动斜激波绕射现象.给出一组运动斜激波绕射波纹壁面的非定常过程的模拟结果。计算结果显示出由运动斜激波绕射诱导的多波干扰产生的复杂流场结构。     

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

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

粘性物质中正激波稳定性分析

用线性稳定性理论,分析粘性物质中的正激波稳定性问题.粘性物质中任意强度的一维激波,其稳定性问题可归结为处理复数范围内的特征值问题,该特征值问题由两个一阶常微分方程及一个二阶常微分方程构成.这些常微分方程的系数依赖于流动的基本流场的物理量及其梯度.所获得的特征值问题由一个四阶精度的有限差分离散求解.分析考虑物质粘性的金属铝中的正激波稳定性,可以看出,正激波运动是稳定的,并且激波速度对波前和波后的小扰动量的衰减有相反的作用,而物质粘性有致稳的作用.     

一种超声速运动目标的被动声定位算法

超声速运动目标在空气中飞行时产生激波，本文利用激波到达时间，提出了一种基于波前方向矢量的算法，确定目标的速度和位置。最后，通过对三站均匀圆阵的估计精度的分析计算，讨论了该算法的有效性．     

计算物理第21卷（2004）总目录

第 1期研究论文基于几何区域分解的三维输运问题并行迭代算法阳述林　莫则尧　沈隆钧　 (1)…………………………变间距采样及准梯度下降法实现光束整形赵逸琼　伍　源　周秋桂　李永平　 (10 )………………………广义特征坐标系计算膨胀波与激波优越性的数值验证陈　正　石　静　吴子牛　 (15 )……………………格子Boltzmann方法中的边界条件聂德明　林建忠　 (2 1)……………………………………………………二维自适应结构网格的变分生成方法于　明　 (2 7)……………………………………………………………用直线稳定方法控制耦合哈密顿…     

激波和剪切层相互作用下的超音速射流

欠膨胀超音速射流处于螺旋模式下的中度欠膨胀时,其入射剪切层的激波具有很高强度,激波和剪切层发生了强烈的相互作用,远场辐射的拢动波出现了大间隔、交错的上下行类似螺旋锥面波形图像,该扰动波具有很强的向上游传播的指向性,导致上游噪声高于垂直喷嘴方向的声压级.而在相对压比较低的低度欠膨胀情况下,或高压比下的高度欠膨胀的情形,入射剪切层激波强度相对较弱,远场辐射没有大间隔、交错的上下行远场辐射 关键词： 超音速射流 啸叫 扰动波 激波 剪切层     

rGFM和level-set方法在水柱群和喷管流场计算中的应用

计算平面运动激波和水柱群相互作用以及喷管流场.在Descartes网格中利用level-set方法分别追踪气/水和气/固界面,采用rGFM方法处理气/水和气/固界面边界条件.将喷管内壁简化为气/固界面并施加固壁边界条件,内壁型线数据拟合采用三次样条插值.采用5阶WENO格式分别求解Euler方程、level-set方程和界面重新初始化方程.给出激波和水柱群相互作用流场密度纹影图和指定点p-t曲线以及喷管流场压力、密度云图和速度场.改进界面法线确定方法可提高Riemann问题构造精度.可分辨运动激波和水柱群作用产生的复杂激波波系,表明激波在各水柱界面的透射和反射、在列和行水柱界面的多次反射和透射.水柱群下游区域的激波波后压力下降,表明激波加热水柱群附近气流和反向运动的反射激波造成了激波衰减.喷管流场数值解和理论解相符.     

基于特征线理论的辨识法在激波装配中的应用

准确地给出激波位置信息对于激波装配极为重要.但是,在使用计算流体力学（computational fluid dyna-mics,CFD）方法模拟复杂流动时很难准确地给出激波的位置.根据激波捕捉得到的流场信息确定的激波位置往往带有极大误差,在定常问题的模拟中,这种误差可以随着迭代逐渐消除,然而在非定常问题的模拟中,这种误差往往会积累甚至导致计算崩溃.文章将基于特征线理论的激波辨识技术应用到激波装配中,根据已有流场信息准确判断激波的位置.对于定常问题,该方法的应用加速了收敛速度;对于非定常问题,该方法的应用可以极大地避免初始误差的产生.      

Lagrange坐标系下二维气动方程组的RKDG有限元方法

构造Lagrange坐标系下二维可压缩气动方程组的RKDG(Runge-Kutta Discontinuous Galerkin)有限元方法.将流体力学方程组和几何守恒律统-求解,所有计算都在固定的网格上进行,计算过程中不需要网格节点的速度信息.对几个数值算例进行数值模拟,得到较好的数值模拟结果.     

Particle paths in Stokes’ edge wave

We treat the particle motion in Stokes’ linear edge wave along a uniformly sloping beach. By a rotation of the coordinate frame, we show that there is no particle motion in the direction orthogonal to the sloping beach, and conclude that particles have a longshore drift in the direction of wave propagation which decreases with depth and distance from the shoreline. We discuss the application of this rotated coordinate frame to higher mode (Ursell) and weakly nonlinear (Whitham) edge waves, and show that the weakly nonlinear case is identical to that for two-dimensional deep-water Stokes waves.     

Roles of positively charged dust in the formation of ion-acoustic subsonic solitary waves in electron-ion-positively charged dust plasmas

The electron-ion-positively charged dust plasma system containing Boltzmann distributed electron species, cold inertial ion species, and stationary positively charged dust (pcd) species are considered. The roles of pcd species in the formation of ion-acoustic (IA) subsonic solitary waves (SWs) are investigated by the pseudo-potential approach, which is valid for arbitrary amplitude time-independent subsonic SWs, as well as by the reductive perturbation method, which is valid for the time-dependent small amplitude subsonic SWs. It is observed that the presence of the pcd species reduces the phase speed of the IA waves, and consequently supports the IA subsonic SWs with the positive wave potential in such electron-ion-pcd plasmas. This is due to the reduction of the space charge electric field by the presence of the pcd species. The applications of the work in space environments (viz. Earth's mesosphere, cometary tails, Jupiter's magnetosphere, etc.), where pcd species have been detected, are briefly discussed.     

Nonlinear fast growth of water waves under wind forcing

In the wind-driven wave regime, the Miles mechanism gives an estimate of the growth rate of the waves under the effect of wind. We consider the case where this growth rate, normalised with respect to the frequency of the carrier wave, is of the order of the wave steepness. Using the method of multiple scales, we calculate the terms which appear in the nonlinear Schrödinger (NLS) equation in this regime of fast-growing waves. We define a coordinate transformation which maps the forced NLS equation into the standard NLS with constant coefficients, that has a number of known analytical soliton solutions. Among these solutions, the Peregrine and the Akhmediev solitons show an enhancement of both their lifetime and maximum amplitude which is in qualitative agreement with the results of tank experiments and numerical simulations of dispersive focusing under the action of wind.     

Effects of the dust size distribution in one-dimensional quantum dusty plasma

The effects of the dust size distribution in ultracold quantum dusty plasmas are investigated. The amplitude φ_m and width ω of quantum dust acoustic waves are studied with different dust size distributions in the system. The φ_m and ω of the quantum dust acoustic waves are found to increase as the total number density increases. The φ_m and ω are greater for unusual dusty plasmas than for typical dusty plasmas. Moreover, as the Fermi temperature of the dust grains increases, the φ_m of the wave decreases. The ω of quantum dust acoustic waves increases as the speed u₀ of the wave increases.     

A numerical algorithm for viscous incompressible interfacial flows

We present a new algorithm to numerically simulate two-dimensional viscous incompressible flows with moving interfaces. The motion is updated in time by using the backward difference formula through an iterative procedure. At each iteration, the pseudo-spectral technique is applied in the horizontal direction. The resulting semi-discretized equations constitute a boundary value problem in the vertical coordinate which is solved by decoupling growing and decaying solutions. Numerical tests justify that this method achieves fully second-order accuracy in both the temporal variable and vertical coordinate. As an application of this algorithm, we study the motion of Stokes waves in the presence of viscosity. Our numerical results are consistent with the recently published asymptotic solution for Stokes waves in slightly viscous fluids.     

Interfacial waves with free-surface boundary conditions: an approach via a model equation

In a two-fluid system where the lower fluid is bounded below by a rigid bottom and the upper fluid is bounded above by a free surface, two kinds of solitary waves can propagate along the interface and the free surface: classical solitary waves characterized by a solitary pulse or generalized solitary waves with nondecaying oscillations in their tails in addition to the solitary pulse. The classical solitary waves move faster than the generalized solitary waves. The origin of the nonlocal solitary waves can be understood from a physical point of view. The dispersion relation for the above system shows that short waves can propagate at the same speed as a “slow” solitary wave. The interaction between the solitary wave and the short waves creates a nonlocal solitary wave. In this paper, the interfacial-wave problem is reduced to a system of ordinary differential equations by using a classical perturbation method, which takes into consideration the possible resonance between short waves and “slow” solitary waves. In the past, classical Korteweg–de Vries type models have been derived but cannot deal with the resonance. All solutions of the new system of model equations, including classical as well as generalized solitary waves, are constructed. The domain of validity of the model is discussed as well. It is also shown that fronts connecting two conjugate states cannot occur for “fast” waves. For “slow” waves, fronts exist but they have ripples in their tails.     

Brillouin-zone mapping of the spin waves in a ferromagnetic bilayer system

The eigenproblems of spin waves in a heterogeneous ferromagnetic bilayer system with periodic boundary conditions are solved using the interface-rescaling approach. Brillouin zone mapping and the eigenmodes of the system are investigated. We find three types of spin waves may exist in the system: the bulk mode, the interface mode, and the perfect confined mode. The fine structure of the energy band in the heterogeneous bilayer system is first given for the whole two-dimensional Brillouin zone. Conditions for the existence of the interface mode are discussed. Finally, we analyze the resonant-confined spin waves in bulk modes and their oscillating behavior.     

Exact plasma wave solutions for isothermal electron fluid plasma

Quasistationary electron plasma waves of arbitrary amplitude and speed that are exact solutions of the isothermal electron fluid equations are shown to exist.     

在光纤零群速色散区传输的光孤子波

通过对超短光脉冲在单模光纤中传输方程的分析研究,给出了在零群速色散传输方程的亮,反波解。结果表明,超短光脉中在光纤的零群群速色散仍能以亮,暗孤波的形式传输,且不存在孤子自频移现象。     

拖曳体激发内波时空特性实验及其理论模型

在具有密度跃层的分层流体中,采用沿水槽中纵剖面对称布置电导率探头阵列的方法,对1个球体和2个不同长径比细长体在拖曳运动下激发内波的时空特性进行了系列实验.结果表明:存在一个与长径比近似为线性关系的临界Froude数Frc,当FrFr_c时,内波相关速度均与物体运动速度一致,体积效应内波为主控内波,内波波高均随拖曳速度增大而先增大后减小,Lee波峰值对应速度随长径比增大而增大;当FrFr_c时,内波相关速度均小于物体运动速度,其相关速度Froude数Fr_(iw)均在0.43—1.18之间的一个条带内变化,尾迹效应内波为主控内波,内波波高均随拖曳速度增大而近似线性增大.此外,从波形结构上看,体积效应内波关于水槽中纵剖面是对称的,而尾迹效应内波关于水槽中纵剖面是不对称的.结合上述实验结果,在已有针对拖曳球产生内波的等效源理论模型基础上,针对体积效应内波,提出了不同长径比模型的等效源移动速度和体积的设置方法;针对尾迹效应内波正对称和反对称这一特性,提出了正对称组合源和反对称组合源理论模型及其参数设置方法.所得计算结果在波高、波形结构和波系分布上与实验结果符合良好,表明了所提出的理论模型及其参数设置方法的合理性和有效性.