Composite materials dynamic fracture studies by generalized Shmuely difference algorithm

The generalized Shmuely Difference Algorithm (GSDA) is presented here to analyze the dynamic fracture performance of orthogonal-anisotropic composite materials, such as glass fibre reinforced phenolplast. The difference recurrence Formulae and boundary condition difference extrapolation formulae are derived and programmed. The dynamic stress intensity factors (DSIF) of the isotropic and anisotropic centrally cracked plates are computed respectively using GSDA and compared with that published previously. GSDA is proved effective and reliable. Copyright (C) 1996 Elsevier Science Ltd.

Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Stress-distribution around a rigid line in dissimilar media

The elastic plane problem of a rigid line inclusion between two dissimilar media was considered. By solving the Riemann-Hilbert problem, the closed-form solution was obtained and the stress distribution around the rigid line was investigated. It was found that the modulus of the singular behavior of the stress remains proportional to the inverse square root of the distance from the rigid line end, but the stresses possess a pronounced oscillatory character as in the case of an interfacial crack tip.

Coflow jets with large velocity difference - a new aerodynamic principle of flame stabilization and combustion intensification

A new aerodynamic principle of flame stabilization and combustion intensification, the coflow jets with large velocity difference, is described. One or more small high-velocity jets of air or steam, injected off the axis and in the same direction as the low-velocity main fuel-air flow into the combustor, create a large recirculation zone of high turbulence intensity in which the combustibles and high temperature gases are effectively mixed, so that stable and intensive combustion can be maintained even for fuels with poor ignition. A pulverized coal combustor based on the principle mentioned above is shown to be characteristic of excellent combustoom and a simple structure. A number of precombustors of this type are in operation at some power stations and industrial boilers of China. Using such precombustor, successtul startups and part-load operation of the boilers have become available under conditions of unpreheated air and low-grade coal with volatiles as low as 15% and ash content as high as 30%. This principle shows good promise as an attractive new technology of combustion.

Spectral-line profile of turbulent gas

phases should be specified when the particle Reynolds number is higher than

Kinetic model of continuous-wave flow chemical lasers

In this paper an analysis of the kinetic theory of the continuous-wave flow chemical lasers(CWFCL) is presented with emphasis being laid on the effects of inhomogeneous broadeningon CWFCL's performance. The results obtained are applicable to the case where laser fre-quency is either coincident or incoincident with that of the eenter of the line shape. This rela-tion has been,compared with that of the rate model in common use. These two models are almostidentical as the broadening parameter η is larger than 1. The smaller the value of η, thegreater the difference between the results of these two models will be. For fixed η, the dif-ferences between fhe results of the two models increase with the increase of the frequencyshift parameter ξ. When η is about less than 0.2. the kinetic model can predict exactly the in-homogeneous broadening effects,while the rate model cannot.

沉降对不同粒径聚苯乙烯胶乳球异向聚集的影响

R.Folkersma等报道了在微重力环境下2μm聚苯乙烯（PS）胶乳球的异向聚集速率有明显增大的结果,本文作者之一孙祉伟等的实验结果却与此有很大差异.为此作者在孙祉伟等的实验基础上对1,2,3μm以及1μm+2μm PS胶乳球混和体系的相对聚集速率进行了研究.作者对原有实验装置进行了改进,并验证了改进后的实验装置的可行性.用密度匹配法实现模拟微重力条件,用快聚集过程中浊度随时间的变化表示相对聚集速率.结果表明,重力引起的沉降对所研究体系聚集速率的影响是很小的,在实验误差范围内可以忽略.作者认为与Stein等结果显著不同的原因是二者使用的样品表面性质不同以及实验方法的差异.

Flow similariy applied in convection-based tilt sensor

The environment temperature has inevitable effects on property of the convect ion-based tilt sensors. It not only redefines the application, but also prevents the improvement of the sensor performance. Numerical simulation of the fluid flow in the chamber of a sensor was performed and the influence of the environment temperature was studied in this paper. At zero tilt angle, the temperature distribution along the perpendicular line cross the heat source at various environment temperatures was presented. It was found that the flow varied dramatically at different environment temperatures, which would cause the output signal vary accordingly, even when the tilt angle was kept at a constant, because this device works by sensing the change of flow. At the same condition, we present the numerical results when the temperature difference across the heat source and the environment was kept at the same, in those results, it was found that the temperature difference at every point along the perpendicular line cross the heat source keep the same, this result confirms the similarity principle of nature convection. Second, A method of eliminating environment temperature infect on property of convect ion-based tilt sensor, which is based on the theory of flow similarity, is proposed. It was found that a thermal transistance can be piped on the circuit of heat source to compensate the temperature of the heat source. A compensative circuit was specially designed which can keep flow similarity by changing heat source temperature in order to eliminate the influence of environment temperature. The experiment results show that above 70% temperature drift can be eliminated by this compensative circuit.

可压缩流气动声场的数值模拟

气动声学是一门流动力学和声学之间的交叉学科，主要研究流动及其与物体相互作用产生噪声的机理。动用计算技术研究气动声学问题的手段称为计算气动声学。本文的目的是，基于高精度数值算法的研究，分别运用Lighthill比拟理论、Kirchhoff积分和直接数值模拟等方法，针对翼型绕流、激波-涡干扰和轴对称射流，研究了物面非定常脉动压力、涡脱落、激波-涡干扰以及涡对并等产生噪声的机理。首先针对声场与主流场在能级和特征尺度等方面的差异，从空间离散角度分析了几种差分格式，表明迎风紧致格式/对称紧致格式有较小的数值色散、耗散和各向异性误差，因而适用于气动噪声的计算。以Runge-Kutta格式为例，对时间离散带来的误差进行了分析。指出对声波计算来说，仅考虑格式稳定性是不够的，时间步长还受到允许色散误差和耗散误差的限制。基于保色戎关系的思想，构造了优化Runge-Kutta格式。处例显示优化Runge-Kutta格式相对于经典格式有更高的计算效率。采用3阶迎风紧致格式和3阶Runge-Kutta格式数值模拟了NACA0012翼型的可压缩非定常绕流流场，并将此流场作为近场声源，运用声学比拟理论对偶极子声和四极子声进行研究。结果指出，主流速度对远场声压有决定性影响，在来流马赫数较大时，四极子噪声和偶极子噪声具有相同量级，不能被忽略，表明了可压缩效应对声场的影响。采用5阶迎风紧致格式和4阶Runge-Kutta格式求解非定常可压缩Navier-Stokes方程，对激波-单涡/双涡干扰导致的声场进行了直接数值模拟。详细研究了激波-涡干扰产生噪声的机理，指出噪声的产生及其性质和激波变形密切相关。研究了近场噪声衰减和传播距离r的关系，发现噪声衰减大致和r~(4/5)而不是r~(1/2)成反比关系，提出这种差异是由流场的非线性效应引起的。构造了Kirchhoff积分和非定常流动计算相结合的算法。采用5阶迎风紧致格式和3阶Runge-Kutta格式对亚声速轴对称射流进行直接数值模拟。将射流流场作为近场声源，结合Kirchhoff方法求解远场 气动噪声。数值结果表明远场噪声具有方向性，噪声声压在离开对称轴20°处达到最大值。随着传播距离增大，噪声方向性逐渐减弱。

若干流体力学问题的数学性质研究

该文利用高智的扩散抛物化方程组理论及流体力学基本方程组的特征次特征理论,流体大小尺度（LSS）方程组理论以及摄动有限差分（PFD）方法,研究若干流体力学问题的数学性质.该文得到的主要结论有：1.利用湍流大小尺度（LSS）方程组推导出湍流大小尺度涡量（LSSV）方程组,并证明两个关于湍流大小尺度涡量的命题,从而得到湍流封闭大小尺度涡量（CLSSV）方程组,并对已有的近程相互作用命题进行推广.2.根据扩散抛物化方程组理论和流体力学层次结构方程组的特征和次特征方法,研究了抛物化稳定性方程组（PSE）的特征和次特征以及消除PSE的剩余椭圆特性的问题.3.利用摄动有限差分（PFD）方法得到对流扩散反应方程的变步长摄动有限差分格式,是等步长摄动有限差分格式的推广.

槽道湍流的直接数值模拟

通过直接数值模拟(DNS)研究槽道湍流的性质和机理。包含五个部分：1）湍流直接数值模拟的差分方法研究。2）求解不可压N－S方程的高效算法和不可压槽道湍流的直接数值模拟。3）可压缩槽道湍流的直接数值模拟和压缩性机理分析。4）“二维湍流”的机理分析。5）槽道湍流的标度律分析。1．针对壁湍流计算网格变化剧烈的特点，构造了基于非等距网格的的迎风紧致格式。该方法直接针对计算网格构造格式中的系数，克服了传统方法采用 Jacobian 变换因网格变化剧烈而带来的误差。针对湍流场的多尺度特性分析了差分格式的精度、网格尺度与数值模拟能分辨的最小尺度的关系，给出不同差分格式对计算网格步长的限制。同时分析了计算中混淆误差的来源和控制方法，指出了迎风型紧致格式能很好地控制混淆误差。2．将上述格式与三阶精度的Adams半隐格式相结合，构造了不可压槽道湍流直接数值模拟的高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项，避免了压力边界条件处理的困难。利用FFT对方程中的隐式部分进行解耦，解耦后的方程采用追赶法（LU分解法）求解，大大减少了计算量。为了检验该方法，进行了三维不可压槽道湍流的直接数值模拟，得到了Re＝2800的充分发展不可压槽道湍流，并对该湍流场进行了统计分析。包括脉动速度偏斜因子在内的各阶统计量与实验结果及Kim等人的计算结果吻合十分理想，说明本方法是行之有效的。3．进行了三维充分发展的可压缩槽道湍流的直接数值模拟。得到了 Re=3300，Ma=0.8的充分发展可压槽道湍流的数据库。流场的统计特征（如等效平均速度分布，“半局部”尺度无量纲化的脉动速度均方根）和他人的数值计算结果吻合。得到了可压槽道湍流的各阶统计量，其中脉动速度的偏斜因子和平坦因子等高阶统计量尚未见其他文献报道。同时还分析了压缩性效应对壁湍流影响的机理，指出近壁处的压力－膨胀项将部分湍流脉动的动能转换成内能，使得可压湍流近壁速度条带结构更加平整。4．模拟了二维不可压槽道流动的饱和态（所谓“二维湍流”），分析了“二维槽道湍流”的非线性行为特征。分析了流场中的上抛－下扫和间歇现象，研究了“二维湍流”与三维湍流的区别。指出“二维湍流”反映了三维湍流的部分特征，同时指出了展向扰动对于湍流核心区发展的重要性。5．首次对可压缩槽道湍流及“二维槽道湍流”标度律进行了分析，得出了以下结论：a）槽道湍流中，在槽道中心线附近较宽的区域，存在标度律。b）该区域流场存在扩展自相似性（ESS）。c）在Mach数不是很高时，压缩性对标度指数影响不大。本文结果同SL标度律的理论值吻合较好，有效支持了该理论。对“二维槽道湍流”也有相似的结论，但与三维湍流不同的是，“二维槽道湍流”存在标度律的区域更宽，近壁处的标度指数比中心处有所升高。