再入尾迹电磁特性的湍流效应研究

本文旨意在于通过探讨高超声速再入尾迹中的湍流等离子体与电磁波相互作用的机理，以及建立能反映此机理的应用性理论模型，从而提供一套可进行目标特性分析的方法，以便为工程部门的突防技术服务。本题目在再入气动物理现象研究中具有重要意义。综合分析指出，地面雷达观测到的非相干散射信号主要来源于再入尾迹的亚密湍流区产生的体积散射。因此，电磁散射特性分析主要针对尾迹亚密湍流等离子体。并且，这里所有的分析都是根据在工程应用中最成熟的一阶畸变波Born近似理论模型。再入尾迹电磁特性的湍流效应研究，着眼点就在于湍流等离子体场的研究。对湍流等离子体场理论模型，本文试图通过模式理论来表达，即求解平均化的全Navier-Stokes方程及其封闭方程k-ε-g模型，从而准确获得流动平均场和脉动场信息。这种表达方式较以前有了较大改进。注意到高超声速流动具有强烈可压缩性的特点，故使用的N-S平均方程由质量加权平均过程产生，湍流模型方程也经过可压缩性修正。方程的离散求解方法，都是运用带矢通量分裂的二阶TVD格式的有限体积法。再入尾迹湍流场的初始条件由近尾迹（底部）流动经N-S方程求解给定，初始值更加准确可靠。尾迹从层流到湍流的转捩过程采用相对成熟的半经验公式确定。飞行器的高超声速再入过程必然导致它周围的空气温度升高，使得流动表现出真实气体效应。对重点考察的湍流流动而言，真实气体效应主要表现为气体处于热化学平衡状态。就工程部门面临的实际问题，把一阶畸变波Born近似的解算方法做些改进，使其能够处理的范围从轴对称尾迹扩展到三维湍流等离子体场是必要的。这为深入的理论分析提供了有力的保障。在能够准确模拟湍流流动的刻划雷达散射截面的基础上，考察亚密湍流等离子体对电磁散射的影响。通过选择的几个有代表性的因素进行讨论，初步结果表明：湍流转捩方式、湍流尺度对尾迹雷达散射截面值计算影响不大，而电子组份脉动能初始值影响较明显，且在特定条件下湍流模型的影响亦不大。但由于湍流模型涉及脉动初始值，其影响需进一步确定。同时，一些今后开展继续此项研究工作的有益建议也提了出来。

障碍物扰动对半开口管道中火焰传播特性的影响

自从1926年Chanman和Wheeler率先开创有障碍物管道中的火焰传播研究工作以来，管道中障碍物扰动引起的火焰加速现象引起了广泛的关注。由于这类现象在燃烧科学上的学术意义及其在生产中诸如安全问题等方面的实际意义，人们相继进行了一系列的研究工作。归纳起来，可以分为两大类：1）封闭管道的火焰传播。人们在不同形状的管道或容器中研究了障碍物对火焰加速的影响，在理论分析和数值计算方面也作了一些有益的工作；2）开口管道的火焰传播。相对于闭口管道，开口管道中的火焰传播研究则逊色得多。这方面尚缺乏系统的实验数据，理论方面的分析也欠缺得多。但当前在实际中，并不乏开口体系的应用，如，现已广泛应用于电力系统的一种燃气除灰装置，是一燃烧气体燃料的半开口系统。因此，研究半开口管道中非稳态燃烧的加速机制具有重要意义。在本论文工作中，通过大量的实验研究，比较系统地研究了障碍物的扰动对预混火焰传播特性的影响。实验在一长L=5m、内径D＝80mm的一端封闭、一端开口的火焰传播管内进行，管内均匀布置障碍物，通过改变障碍物的形状、间距、阻塞比大小，同时选用五种不同的可燃气体，探索了障碍物结构对预混湍流火焰加速和管内压力上升的影响。实验表明，对于敏感气体如氢气和乙炔，由于障碍物扰动产生的影响，火焰不断加速，并最终达到一准稳定状态；在适当的条件下，火焰传播状态可由爆燃向爆轰转变，此时火焰速度发生跃变；而对干不敏感气体如甲烷，则爆燃转爆轰现象不容易发生。在不同的火焰传播状态，障碍物结构特性对火焰速度和压力产生的影响各不相同。在缓燃态，随着阻塞比的变化，最大火焰速度先上升后下降，在BR=0.3～0.4之间存在一最大值；在銮塞态，最大火焰速度受阻塞比变化的影响不明显，略低于燃烧产物的声速；在阻塞比BR＝0.5附近，压力达到最大值。而在爆轰态，随着阻塞比的增加，最大火焰速度和压力逐渐降低，爆燃转爆轰的浓度范围变小。由于在有障碍物的管道中，火焰速度很容易达到声速（变塞态）或超声速（爆轰态），必须考虑流体的马赫数效应。本文在前人研究成果的基础上，给出了湍流马赫数修正的可压缩性两方程湍流模型，模拟了半开口狭长管道中重复布置的障碍物引起的湍流火焰加速现象。最大火焰速度和管内压力的计算结果与实验测量值吻合良好，这表明用修正后的湍流模型能够比较真实地模拟障碍物管内预混火焰的发展过程。通过对管道内障碍物扰动引起的燃烧波加速的机理和技术研究，对流场扰动对燃烧波产生和发展的影响规律有了比较全面的了解，研究结果对气脉冲除灰技术的完善具有直接的指导意义。

Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle- resolved simulations are briefly discussed.

Vortex dislocations in wake-type flow induced by spanwise disturbances

Vortex dislocations in wake-type flow induced by three types of spanwise disturbances superimposed on an upstream velocity profile are investigated by direct numerical simulations. Three distinct modes of vortex dislocations and flow transitions have been found. A local spanwise exponential decay disturbance leads to the appearance of a twisted chainlike mode of vortex dislocation. A stepped spanwise disturbance causes a streamwise periodic spotlike mode of vortex dislocation. A spanwise sinusoidal wavy disturbance with a moderate waviness causes a strong unsteadiness of wake behavior. This unsteadiness starts with a systematic periodic mode of vortex dislocation in the spanwise direction followed by the spanwise vortex shedding suppressed completely with increased time and the near wake becoming a steady shear flow. Characteristics of these modes of vortex dislocation and complex vortex linkages over the dislocation, as well as the corresponding dynamic processes related to the appearance of dislocations, are described by examining the variations of vortex lines and vorticity distribution. The nature of the vortex dislocation is demonstrated by the substantial vorticity modification of the spanwise vortex from the original spanwise direction to streamwise and vertical directions, accompanied by the appearance of noticeable vortex branching and complex vortex linking, all of which are produced at the locations with the biggest phase difference or with a frequency discontinuity between shedding cells. The effect of vortex dislocation on flow transition, either to an unsteady irregular vortex flow or suppression of the Kaacutermaacuten vortex shedding making the wake flow steady state, is analyzed. Distinct similarities are found in the mechanism and main flow phenomena between the present numerical results obtained in wake-type flows and the experimental-numerical results of cylinder wakes reported in previous studies.

Numerical Investigation On The Flowfield Of "Swallowtail" Cavity For Supersonic Mixing Enhancement

A "swallowtail" cavity for the supersonic combustor was proposed to serve as an efficient flame holder for scramjets by enhancing the mass exchange between the cavity and the main flow. A numerical study on the "swallowtail" cavity was conducted by solving the three-dimensional Reynolds-averaged Navier-Stokes equations implemented with a k-epsilon turbulence model in a multi-block mesh. Turbulence model and numerical algorithms were validated first, and then test cases were calculated to investigate into the mechanism of cavity flows. Numerical results demonstrated that the certain mass in the supersonic main flow was sucked into the cavity and moved spirally toward the combustor walls. After that, the flow went out of the cavity at its lateral end, and finally was efficiently mixed with the main flow. The comparison between the "swallowtail" cavity and the conventional one showed that the mass exchanged between the cavity and the main flow was enhanced by the lateral flow that was induced due to the pressure gradient inside the cavity and was driven by the three-dimensional vortex ring generated from the "swallowtail" cavity structure.

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.

基于畸变相位波前分形特征产生矩形湍流相屏

在包含时间进程的光波大气传输及其自适应光学相位校正的数值模拟研究中,如长曝光成像和自适应光学系统的动态控制过程,矩形湍流相屏的产生和应用尤为重要.而现在通常使用的功率谱反演法产生的是正方形的湍流相屏,只采用其中的矩形部分显然造成计算机资源的浪费;并且谱反演法产生的湍流相屏需要进行低频补偿,从而明显地增加计算量.基于大气湍流所造成的畸变相位波前的分形特征,提出了一种产生矩形湍流相屏的新方法,并与解析理论结果进行对比,验证了这种矩形相屏产生方法的正确性.与已有的方法相比,此算法具有两个明显的优点:算法简单、计算效率高,节省计算机资源;与大气湍流介质统计特性无论在高频部分还是在低频部分均符合得较好.

Performance analysis of bit error rate for free space optical communication with tip-tilt compensation based on gamma-gamma distribution

In this paper, the gamma-gamma probability distribution is used to model turbulent channels. The bit error rate (BER) performance of free space optical (FSO) communication systems employing on-off keying (OOK) or subcarrier binary phase-shift keying (BPSK) modulation format is derived. A tip-tilt adaptive optics system is also incorporated with a FSO system using the above modulation formats. The tip-tilt compensation can alleviate effects of atmospheric turbulence and thereby improve the BER performance. The improvement is different for different turbulence strengths and modulation formats. In addition, the BER performance of communication systems employing subcarrier BPSK modulation is much better than that of compatible systems employing OOK modulation with or without tip-tilt compensation.

Quantifying chaos: practical estimation of the correlation dimension

Be it a physical object or a mathematical model, a nonlinear dynamical system can display complicated aperiodic behavior, or "chaos." In many cases, this chaos is associated with motion on a strange attractor in the system's phase space. And the dimension of the strange attractor indicates the effective number of degrees of freedom in the dynamical system.

In this thesis, we investigate numerical issues involved with estimating the dimension of a strange attractor from a finite time series of measurements on the dynamical system.

Of the various definitions of dimension, we argue that the correlation dimension is the most efficiently calculable and we remark further that it is the most commonly calculated. We are concerned with the practical problems that arise in attempting to compute the correlation dimension. We deal with geometrical effects (due to the inexact self-similarity of the attractor), dynamical effects (due to the nonindependence of points generated by the dynamical system that defines the attractor), and statistical effects (due to the finite number of points that sample the attractor). We propose a modification of the standard algorithm, which eliminates a specific effect due to autocorrelation, and a new implementation of the correlation algorithm, which is computationally efficient.

Finally, we apply the algorithm to chaotic data from the Caltech tokamak and the Texas tokamak (TEXT); we conclude that plasma turbulence is not a low- dimensional phenomenon.