Interactions between kinetic and radiative processes inside moving absorbing fluids

The aim of this paper is to present several features of the couplings occurring between radiative transfer and the kinetics of a moving dielectric. After determining how the velocity field affects the apparent thermo-optical properties of matter, the energy transport problem is investigated in instationary regime and the general form of transient radiative transfer equation inside a moving medium is built. Then, the model is applied to the particular case of turbulent flows: a system of two equations for mean and fluctuating radiative energies is presented, and the resolution of this system is finally carried out.     

研究两相流中固粒对流体湍动特性影响的新方法

本文提出了计算两相流中固粒对流体湍动特性影响的一种新方法，得到不同情况下固粒对流体端动特性的影响。将该方法用于槽流湍流场的求解，说明该方法是可行的。     

Analysis of a turbulent jet mixing flow by using a PIV-PLIF combined system

In the present paper, the development of a high-resolution PIV-PLIF combined system for the simultaneous measurements of velocity and passive scalar concentration fields is described. The high-resolution PIV-PLIF combined system is used to perform the simultaneous whole-field measurements of velocity and concentration in the near field of a turbulent jet mixing flow. The characteristics of the mass transfer process and momentum transfer process in the near field of the jet mixing flows are discussed in the terms of ensemble-averaged velocity and concentration, turbulent intensity, concentration standard deviation and the correlation terms between the fluctuating velocities and concentration.     

Electromagnetic scattering model of the Kelvin wake and turbulent wake by a moving ship

Taking the attenuation character of the Kelvin wake and the limitation of the traditional two-scale method into account, the practical electromagnetic (EM) scattering model of the Kelvin wake is obtained by using a facet-based model; and for a turbulent ship wake, it is produced by dealing with the wave energy loss rate due to turbulence with the width of turbulent wake closely following a moving ship. The volume scattering theory of foam or bubbles, by using the vector radiation transfer equation, is introduced to modify the traditional two-scale method, a special method for the EM scattering of a turbulent wake, which brings the scattering of ship wakes into better agreement with the real situation.     

钝体后湍流预混燃烧的PDF模拟

本文采用PDF方法对矩形燃烧室内钝体后的湍流预混火焰进行了数值模拟。脉动速度-频率-标量联合的PDF输运方程用Monte Carlo方法求解,质量、动量和能量的平均值由基于无结构网格的有限体积法求解,压力通过状态方程获得。PDF方程中所需的平均密度、平均速度和压力由有限体积法提供,并将用Monte Carlo方法求出的雷诺应力、化学反应源项和比热比传递给有限体积法。本文对丙烷和空气燃烧的不同简化化学反应机理进行了研究,并与实验结果进行比较,获得满意的结果。     

Generation of a magnetic field by dynamo action in a turbulent flow of liquid sodium

We report the observation of dynamo action in the von Kármán sodium experiment, i.e., the generation of a magnetic field by a strongly turbulent swirling flow of liquid sodium. Both mean and fluctuating parts of the field are studied. The dynamo threshold corresponds to a magnetic Reynolds number R(m) approximately 30. A mean magnetic field of the order of 40 G is observed 30% above threshold at the flow lateral boundary. The rms fluctuations are larger than the corresponding mean value for two of the components. The scaling of the mean square magnetic field is compared to a prediction previously made for high Reynolds number flows.     

A Modeling Study on Particle Dispersion in Wall-Bounded Turbulent Flows

Three physical mechanisms which may affect dispersion of particle's motion in wall-bounded turbulent flows, including the effects of turbulence, wall roughness in particle-wall collisions, and inter-particle collisions, are numerically investigated in this study. Parametric studies with different wall roughness extents and with different mass loading ratios of particles are performed in fully developed channel flows with the Eulerian-Lagrangian approach. A low-Reynolds-number $k-\epsilon$ turbulence model is applied for the solution of the carrier-flow field, while the deterministic Lagrangian method together with binary-collision hard-sphere model is applied for the solution of particle motion. It is shown that the mechanism of inter-particle collisions should be taken into account in the modeling except for the flows laden with sufficiently low mass loading ratios of particles. Influences of wall roughness on particle dispersion due to particle-wall collisions are found to be considerable in the bounded particle-laden flow. Since the investigated particles are associated with large Stokes numbers, i.e., larger than $\mathcal{O}(1)$, in the test problem, the effects of turbulence on particle dispersion are much less considerable, as expected, in comparison with another two physical mechanisms investigated in the study.     

Time-averaging strategies in the finite-volume/particle hybrid algorithm for the joint PDF equation of turbulent reactive flows

The influence of time-averaging on bias is investigated in the finite-volume/particle hybrid algorithm for the joint PDF equation for statistically-stationary turbulent reactive flows. It is found that the time-averaging of the mean fluctuating velocity (TAu) leads to the same variances of the fluctuating velocity before and after the velocity correction, whereas without TAu the estimates are different, and an additional numerical dissipation rate is introduced for the turbulent kinetic energy (TKE). When 100 particles per cell are used without TAu, a large bias error is found to be involved in the unconditional statistics of the statistically-stationary solutions of two tested turbulent flames, the Cabra H₂/N₂ lifted flame and the Sandia piloted flame E. The use of TAu reduces this bias dramatically for the same number of particles per cell. The conditional statistics in these flames, however, are hardly affected by TAu. To a large extent, the effect of the bias error on the unconditional statistics is similar to the effect of increasing the model constant C _{ω 1} in the stochastic turbulence frequency model.     

A mixed-time-scale low-Reynolds-number one-equation turbulence model

In this paper, a new low-Reynolds-number (LRN) one-equation turbulence model for eddy viscosity is proposed. A mixed time scale, representing a combination of three time scales: two time scales made of strain-rate parameter S and vorticity parameter Ω and the turbulent time scale k/?, is introduced into this model. The proposed model is derived from an LRN k?? two-equation model where the mixed time scale has been proved to be very effective for predicting local flows over complex terrains. In the transport equation of the model, the mixed time scale is included in the production and the dissipation terms. The new model is evaluated in channel flows at various Reynolds numbers, boundary layer flows with or without pressure gradient and backward-facing step flows with different expansion ratios and Reynolds numbers. Then the grid convergence of the model is investigated. Finally, the model performance for different values of the weighting constant C_s in the mixed time scale is assessed. The results show that the proposed model reproduces the correct wall-limiting behaviour of turbulent quantities and performs well in the near-wall region of turbulent flows. The model could be expected to be adopted in hybrid Reynolds averaged Navier–Stokes/large eddy simulation methodology for complex wall-bounded flows at high Reynolds numbers.     

各向异性浮力紊流的显式代数应力模型及其验证

在Wallin & Johansson的显式代数雷诺应力模型(EARSM)基础上,引入浮力的紊动效应,使该模型能够反映浮力对紊流运动的影响.利用同位网格上的SIMPLEC算法,模拟了T-型连接管内的温度混合流动.与实验资料和DNS计算结果比较表明,提出的模型能够反映各向异性的浮力紊流特征,而且模型在计算进程中显示出良好的数值收敛性能,与涡粘性两方程模型的计算花费相当.     

A modified PANS model for computations of unsteady turbulence cavitating flows

A modification to the PANS (partially averaged Navier-Stokes) model is proposed to simulate unsteady cavitating flows. In the model, the parameter f _k is modified to vary as a function of the ratios between the water density and the mixture density in the local flows. The objective of this study is to validate the modified model and further understand the interaction between turbulence and cavitation around a Clark-Y hydrofoil. The comparisons between the numerical and experiment results show that the modified model can be improved to predict the cavity evolution, vortex shedding frequency and the lift force fluctuating in time fairly well, as it can effectively modulate the eddy viscosity in the cavitating region and various levels of physical turbulent fluctuations are resolved. In addition, from the computational results, it is proved that cavitation phenomenon physically influences the turbulent level, especially by the vortex shedding behaviors. Also, the mean u-velocity profiles demonstrate that the attached cavity thickness can alter the local turbulent shear layer.     

Characteristic boundary conditions for simulations of compressible reacting flows with multi-dimensional,viscous and reaction effects

A generalized formulation of the characteristic boundary conditions for compressible reacting flows is proposed. The new and improved approach resolves a number of lingering issues of spurious solution behaviour encountered in turbulent reacting flow simulations in the past. This is accomplished (a) by accounting for all the relevant terms in the determination of the characteristic wave amplitudes and (b) by accommodating a relaxation treatment for the transverse gradient terms with the relaxation coefficient properly determined by the low Mach number asymptotic expansion. The new boundary conditions are applied to a comprehensive set of test problems including: vortex-convection; turbulent inflow; ignition front propagation; non-reacting and reacting Poiseuille flows; and counterflow cases. It is demonstrated that the improved boundary conditions perform consistently superior to existing approaches, and result in robust and accurate solutions with minimal acoustic wave interactions at the boundary in hostile turbulent combustion simulation conditions.     

A multi-scale simulation method for high Reynolds number wall-bounded turbulent flows

We present an assessment and enhancement of the hybrid two-level large-eddy simulation method (A.G. Gungor and S. Menon, A new two-scale model for large eddy simulation of wall-bounded flows, Prog. Aerosp. Sci. 46 (2010), pp. 28–45), a multi-scale formulation for simulation of high Reynolds number wall-bounded turbulent flows. The assessment of the method is performed by examining role of static and dynamic blending functions used to perform hybridisation of two-level simulation (K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows, J. Comput. Phys. 220 (2006), pp. 290–311; K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows. Part 2: Non-homogeneous flows, J. Comput. Phys. 222 (2007), pp. 673–701) and large-eddy simulation methods. The sensitivity of first- and second-order turbulence statistics to the type of blending functions is investigated by simulating a fully developed turbulent flow in a channel at a friction Reynolds number Re_τ = 395 and comparing the results with those obtained using a direct numerical simulation. The first-order statistics do not show any significant differences for different blending functions, but the second-order statistics show some minor differences. The dynamic evaluation of the hybrid region and the blending function is necessary for non-equilibrium and complex flows where use of a static blending function can lead to inaccurate results. We propose two criteria for the dynamic evaluation; first evaluates extent of the hybrid region based on the subgrid turbulent kinetic energy and the second estimates the blending function based on a characteristic length scale. The computational efficiency of the method is enhanced by incorporating a hybrid programming paradigm where a standard domain decomposition by the message-passing-interface library is combined with the open multi-processing based parallelisation. A further enhancement of the method is achieved by incorporating a closure model for the unclosed hybrid terms in the governing equations, which appear due to hybridisation of two-level- and large-eddy-simulation methods. The model is based on an order of magnitude approximation and a preliminary assessment of the model shows improvement of turbulence statistics when used to simulate turbulent flow in a periodic channel. The assessment and improvements to the multi-scale method make it more suitable for simulation of practical wall-bounded turbulent flows at higher Reynolds number than a conventional large-eddy simulation. This is demonstrated by simulating two representative cases; turbulent flow at high Reynolds number in a periodic channel and flow over a bump placed on the lower surface of a channel, where a relatively coarser computational grid is found to be sufficient for reasonably accurate results.     

An Improved Method to Determine Particle Dispersion Width for Efficient Modeling of Turbulent Two‐Phase Flows

An improved approach is presented for the hybrid Eulerian‐Lagrangian modeling of turbulent two‐phase flows. The hybrid model consists of a nonlinear k–ε model for the fluid flow and an efficient Lagrangian trajectory model for the particulate flow. The improved approach avoids an empirical correlation required to determine the dispersion width for the existing Stochastic‐Probabilistic Efficiency Enhanced Dispersion (SPEED) model. The improved SPEED model is validated using experimental data for a poly‐dispersed water spray interacting with a turbulent annular air jet behind a bluff‐body. Numerical results for the number‐mean and Sauter‐mean droplet diameters, as well as mean and fluctuating droplet velocities are compared with the experimental data and with the predictions of other dispersion models. It is demonstrated that higher computational efficiency and smoother profiles of Sauter‐mean diameter can be obtained with the improved stochastic‐probabilistic model than with the eddy‐interaction model.     

Small-slope expansion for thermal and reflected radiation from a rough surface

Small-perturbation and small-slope expansions are applied and compared for calculation of thermal emission from a rough sea surface and reflected sky radiation. The comparison shows that the expansions are identical. This permits us to use the small-perturbation expansion to calculate microwave thermal radiation from the ocean instead of the two-scale model and the small-slope expansion. At grazing observation angles multiple scattering and shadowing become crucial and high-order terms of the expansion must be taken into account. This imposes a limitation on the applicability of the small-slope approximation.     

圆柱尾流温度标量场的小波分析$lt;font color="#ff0000"$gt;$lt;font color="#ff0000"$gt;（已撤稿）$lt;/font$gt;$lt;/font$gt;

通过使用一排16根冷线探头排在多个空间点同时测量微加热圆柱的尾流温度场, 用小波分析技术对瞬时温度场的时间序列信号进行多尺度分析, 目的是研究不同尺度脉动温度对总体温度场的贡献.直径为d = 12.7 mm 的圆柱产生了被测的尾流, 对应的雷诺数为5500, 测量区域位于下游距离为2d 和 20d 之间. 基于小波多尺度分辨技术, 尾流温度场被分解为不同温度脉动特征尺度的小波分量. 通过分析这些小波分量的瞬时温度等值线图, 能够直接观测到不同特征尺度的涡结构运动特征和湍流间歇过程. 特别地, 我们在近场区从原始信号分解获得的高频区域中发现了K-H涡的存在. 不同尺度的温度方差沿流向的变化表明, 在下游距离为x=3d和 20d之间, 中等尺度的结构比大尺度和小尺度结构对总的温度均方根的贡献更大. 不同尺度的自相关函数表明, 大尺度和中等尺度的结构显示出较大的相关性, 而高频的小波分量则更快地失去了原有的拟序性. 关键词： 湍流尾流 被动标量 小波分析     

Response of plasma turbulence against externally-controlled perturbations

A theory of the dynamic response of the turbulent plasma against the externally-controlled perturbations is reported. Based on Mori’s method [Prog. Theor. Phys. 33 423 (1965)], the nonlinear force is assumed to be separated into the memory function and the nonlinear fluctuating force. The former corresponds to the damping term, and the latter is categorized into the noise term. The response of the turbulent plasma against the externally-controlled source is formulated. The response kernel, which connects the externally-controlled source and the response of the turbulent field, is shown to have both the nonlocal property (in space) and the non-Markovian response (in time). A discussion is made on the nonlocal and non-Markovian response, including the case of disparate-scale interactions. A new method is proposed to observe experimentally the nonlocal interaction in the drift wave turbulence via the zonal flows.     

Determinations of turbulent velocity fluctuations and mean particle radii in flames using scattered laser-power spectra

We describe a procedure for simultaneous determinations of mean particle radii and mean values of fluctuating velocity components in turbulent flows containing monodisperse particles. The results of in situ particle-size measurements are compared with estimates derived after withdrawals of samples.     

模型加力燃烧室燃烧流场的并行计算

本文采用基于MPI的并行算法,采用动态内存分配、分区算法和多点重合交错网格系统,在贴体网格下对带V形槽稳定器模型加力燃烧室素流化学反应流场进行数值模拟,湍流模型采用k方程亚网格尺度模型,燃烧模型采用亚网格EBU模型,采用热通量辐射模型估算辐射通量。计算结果表明并行计算对复杂形状的化学反应流动计算效率很高,是模拟大规模的燃烧问题的有力工具。