两相湍流统一色噪声法概率密度函数模型

统一色噪声近似方法对简单一维色噪声问题研究较为充分, 本文将统一色噪声法应用到高度复杂的多维气固两相湍流系统之中.首先从颗粒运动Langevin方程出发, 利用统一色噪声法获得两相湍流Fokker-Planck方程, 然后以此为基础建立颗粒轨道两阶矩模型.文中建立的新模型成功应用于后台阶两相湍流流场的数值模拟, 预报合理正确.研究表明, 对于多维两相湍流系统, 统一色噪声法仍然行之有效.     

The generating functional for the probability density functions of Navier-Stokes turbulence

A generating functional for the equal-time spatial probability density functions which represent the ensemble of turbulent incompressible Navier-Stokes fluids is introduced. By formally solving the linear evolution equation satisfied by this functional, the probability densities are represented as functional integrals. It is shown that the generating functional can be regarded as the space characteristic functional of a generalized random field defined on the phase space spanned by the material position and velocity fields of a fluid particle. The interpretation of this random field, which satisfies a dynamical equation similar to Vlasov's, is clarified through the formal analogies between the statistics of molecules and fluid particles at the functional level. A class of statistically realizable and solvable models is also considered within the context of the present formalism.     

随机系统的概率密度函数形状调节

针对受高斯白噪声激励的非线性随机系统,提出了使状态响应的概率密度函数形状跟踪期望形状的调节方法.首先,确立了非线性随机系统的多项式反馈机制,同时对系统中的非线性部分进行多项式展开;然后,以Fokker-Planck-Kolmogorov方程为工具,导出了与控制增益相关的各阶矩递推方程,并根据跟踪问题的要求,构造了矩逼近优化问题,用梯度搜索法求解该优化问题,获得了调节函数;再依据特征函数与概率密度函数构成Fourier对的关系,对状态响应的概率密度函数进行重构;最后,通过两个例子仿真,验证了本文方法的有效性.     

谐振子的概率密度、概率流密度与时间的关系

以一维线性谐振子为例,对非定态情况,通过数值计算给出了不同时间的概率密度和概率流密度分布,并且讨论了概率密度和概率流密度随时间的变化的基本特征.     

氢原子电子云密度分布分析

通过数值计算形象地分析了氢原子中电子概率密度的分布情况,对此做了一些有别于传统认识的讨论,并由此引出了对粒子运动和波函数概率解释的认识。     

A probability density function of liftoff velocities in mixed-size wind sand flux

With the discrete element method (DEM), employing the diameter distribution of natural sands sampled from the Tengger Desert, a mixed-size sand bed was produced and the particle-bed collision was simulated in the mixed-size wind sand movement. In the simulation, the shear wind velocity, particle diameter, incident velocity and incident angle of the impact sand particle were given the same values as the experimental results. After the particle-bed collision, we collected all the initial velocities of rising sand particles, including the liftoff angular velocities, liftoff linear velocities and their horizontal and vertical components. By the statistical analysis on the velocity sample for each velocity component, its probability density functions were obtained, and they are the functions of the shear wind velocity. The liftoff velocities and their horizontal and vertical components are distributed as an exponential density function, while the angular velocities are distributed as a normal density function. Supported by the Key Project of the National Natural Science Foundation of China (Grant No. 10532040)     

Statistics and structures of pressure and density in compressible isotropic turbulence

We study statistics and structures of pressure and density in the presence of large-scale shock waves in a forced compressible isotropic turbulence using high-resolution numerical simulation. The spectra for pressure and density exhibit a ?2 scaling over an operational definition of the inertial range. Both the numerical simulation and a heuristic PDF model reveal that the PDFs of pressure increment exhibit a ?2 power law region for the separation in the operational definition of inertial range, quantitatively similar to the PDF of pressure gradient, which also displays a ?2 power law region. Moreover, the statistical relation between density increment and pressure increment has been investigated through a shock-relation model. There is a positive correlation between the vorticity magnitude and pressure, which is different from the case of incompressible turbulence. We argue that this difference is due to large-scale shock waves, another type of intermittent structures in addition to vortex structures in incompressible turbulence.     

Correlation dimension analysis and capillary wave turbulence in Dragon-Wash phenomena

This paper describes the evolution of surface capillary waves of deep water excited by gradually increasing the lateral external force at a single frequency. The vertical velocities of the water surface are measured by using a Polytec Laser Vibrometer with a thin layer of aluminium powder scattering on the surface to reflect the laser beam. Nonlinear interaction processes result in a stationary Fourier spectrum of the vertical surface velocities （the same as the surface elevation）, i.e. Iω -ω^-3-5. The observed spectrum can be interpreted as a wave-turbulent Kolmogorov spectrum for the case of ‘narrowband pumping＇ for a direct cascade of energy. Correlation dimension analysis of the whole development process reveals four distinct stages during the wave structure development and identifies the wave turbulence stage.     

Single point velocity distributions in homogeneous isotropic turbulence

The starting point for this paper lies in the results obtained by Tatsumi (2004) for isotropic turbulence with the self-preserving hypothesis. A careful consideration of the mathematical structure of the one-point velocity distribution function equation obtained by Tatsumi (2004) leads to an exact analysis of all possible cases and to all admissible solutions of the problem. This paper revisits this interesting problem from a new point of view, and obtains a new complete set of solutions. Based on these exact solutions, some physically significant consequences of recent advances in the theory of homogenous statistical solution of the Navier--Stokes equations are presented. The comparison with former theory was also made. The origin of non--Gaussian character could be deduced from the above exact solutions.     

海洋湍流中光波特征参量和短期光束扩展的研究

Nikishov等建立的海洋湍流功率谱模型中,假设了海水有着稳定的分层.但是,实际海水通常不是稳定分层的,温度与盐度的涡流扩散率是不相等的.2017年,Elamassie等建立了考虑这些因素的更合理的海洋湍流功率谱模型.湍流介质中光波空间相干长度等基本特征参量在表征湍流强度和光传输相位校正技术等方面起着重要作用.本文基于Elamassie海洋湍流功率谱模型,重新推导出了海洋湍流中光波结构函数、光波空间相干长度和Fried参数的解析公式,并校验了所得公式的正确性.研究发现：当温度变化引起的光学湍流占主导地位时,Nikishov海洋湍流功率谱模型把湍流强度低估了;当盐度变化引起的光学湍流占主导地位时,Nikishov海洋湍流功率谱模型把湍流强度高估了.基于Elamassie海洋湍流功率谱模型,本文推导出了高斯光束短期光束扩展的半解析公式,并验证了其正确性.研究还表明：海水稳定分层与否,短期光束扩展差异很大.本文研究结果对水下湍流环境中的光通信、成像和传感等应用具有重要意义.     

Multi-environment probability density function method for modelling turbulent combustion using realistic chemical kinetics

A computational fluid dynamics (CFD) tool for performing turbulent combustion simulations that require finite-rate chemistry is developed and tested by modelling a series of bluff-body stabilized flames that exhibit different levels of finite-rate chemistry effects ranging from near equilibrium to near global extinction. The new modelling tool is based on the multi-environment probability density function (MEPDF) methodology and combines the following: the direct quadrature method of moments (DQMOM); the interaction-by-exchange-with-the-mean (IEM) mixing model; and realistic combustion chemistry. Using DQMOM, the MEPDF model can be derived from the transport PDF equation by depicting the joint composition PDF as a weighted summation of a finite number of multi-dimensional Dirac delta functions in the composition space. The MEPDF method with multiple reactive scalars retains the unique property of the joint PDF method of treating chemical reactions exactly. However, unlike the joint PDF methods that typically must resort to particle-based Monte-Carlo solution schemes, the MEPDF equations (i.e. the transport equations of the weighted delta-peaks) can be solved by traditional Eulerian grid-based techniques. In the current study, a pseudo time-splitting scheme is adopted to solve the MEPDF equations; the reaction source terms are computed with a highly efficient and accurate in-situ adaptive tabulation (ISAT) algorithm. A 19-species reduced mechanism based on quasi-steady state assumptions is used in the simulations of the bluff-body flames. The modelling results are compared with the experimental data, including mixing, temperature, major species and important minor species such as CO and NO. Compared with simulations using a Monte-Carlo joint PDF method, the new approach shows comparable accuracy.     

时变反应扩散系统中螺旋波和湍流的控制

研究了一类参数时变的反应扩散系统中螺旋波和湍流对外电场的响应问题.在数值模拟中,以一类改进的Fitzhugh-Nagumo模型为研究对象(在恰当参数值下可分别描述激发介质和振荡介质),考虑随机和不确定因素(如内外噪声、气压、温度梯度分布和介质形变等)所引起的系统参数涨落对斑图演化的影响,在模拟中选取的参数涨落范围确保系统可以观测到稳定旋转的螺旋波、漫游的螺旋波和湍流,经历一定的暂态过程后,对介质施加极化电场,研究螺旋波和湍流在外电场中的演化.数值计算结果表明：在系统参数发生涨落和外电场强度比较小情况下,主 关键词： 螺旋波 湍流 时变系统 Fitzhugh-Nagumo模型     

Scrape-off layer width set by non-linear streamer flows in drift wave turbulence

In this work, we discuss the physics behind the excitation of non-linear streamer flows in drift wave turbulence and how to selectively excite these flows. Streamer flows are modelled as a non-linear, radially elongated convective cell in drift wave turbulence. It is shown that density modulation is key for exciting streamer flows. We show that streamer flows have a finite frequency, albeit smaller than that of drift waves. Streamers propagate in the ion direction. These theoretical predictions are compared against experimental data, which shows reasonable agreement. Finally, the scrape-off layer width set by streamer flows is calculated, and a scaling law against macroscopic plasma parameters is obtained.     

谐振子在含时均匀外场作用下跃迁概率的精确解

采用尝试波函数的方法求解了含时受迫振子的严格波函数，并给出了跃迁概率的精确解。     

Stabilization of spiral wave and turbulence in the excitable media using parameter perturbation scheme

This paper proposes a scheme of parameter perturbation to suppress the stable rotating spiral wave, meandering spiral wave and turbulence in the excitable media, which is described by the modified Fitzhug-Nagumo （MFHN） model. The controllable parameter in the MFHN model is perturbed with a weak pulse and the pulse period is decided by the rotating period of the spiral wave approximatively. It is confirmed that the spiral wave and spiral turbulence can be suppressed greatly. Drift and instability of spiral wave can be observed in the numerical simulation tests before the whole media become homogeneous finally.     

Large eddy simulation of turbulent premixed combustion using tabulated detailed chemistry and presumed probability density function

A method of chemistry tabulation combined with presumed probability density function (PDF) is applied to simulate piloted premixed jet burner flames with high Karlovitz number using large eddy simulation. Thermo-chemistry states are tabulated by the combination of auto-ignition and extended auto-ignition model. To evaluate the predictive capability of the proposed tabulation method to represent the thermo-chemistry states under the condition of different fresh gases temperature, a-priori study is conducted by performing idealised transient one-dimensional premixed flame simulations. Presumed PDF is used to involve the interaction of turbulence and flame with beta PDF to model the reaction progress variable distribution. Two presumed PDF models, Dirichlet distribution and independent beta distribution, respectively, are applied for representing the interaction between two mixture fractions that are associated with three inlet streams. Comparisons of statistical results show that two presumed PDF models for the two mixture fractions are both capable of predicting temperature and major species profiles, however, they are shown to have a significant effect on the predictions for intermediate species. An analysis of the thermo-chemical state-space representation of the sub-grid scale (SGS) combustion model is performed by comparing correlations between the carbon monoxide mass fraction and temperature. The SGS combustion model based on the proposed chemistry tabulation can reasonably capture the peak value and change trend of intermediate species. Aspects regarding model extensions to adequately predict the peak location of intermediate species are discussed.     

在HT-7托卡马克上密度扫描实验中的边界湍流和输运

在HT-7托卡马克上密度扫描实验中使用一个快速扫描气动朗缪尔探针来研究边界湍流和输运。随着中心弦平均等离子体密度的增加，在等离子体边界观察到径向电场剪切增强，径向电场的增强能够解释边界输运的减小和全局粒子约束的改善。粒子约束时间的增加和垂直扩散系数的减小证明粒子输运降低了。在等离子体边界区域发现电离和辐射驱动湍流的特征，表明电离和辐射在湍流驱动过程中具有重要性。     

圆柱热尾流中温度的概率密度函数

本文通过实验研究圆柱热尾流中温度的概率密度函数在不同雷诺数下沿中心线的演变以及其与湍流混合程度的关系.实验中雷诺数的取值范围是1200-8600,温度是由直径为0.63 μm的冷线探针测量的.实验结果表明,温度概率密度函数在尾流中场区随空间位置变化显著.雷诺数的增加加快了这个变化过程,特别加速了尾流中心线温度的概率密度函数从非高斯向近似高斯分布的演变.     

关于双曲函数方法求孤波解的注记

针对文献[18]提出的求解非线性波动方程孤波解的双曲函数方法和文献[19]的分析和改进,给出一个注记.并进一步讨论了它的应用.表明这种方法确实是一种简单而实用的方法. 关键词： 双曲函数方法 孤波解 非线性波动方程