Numerical simulation of flow over urban-like topographies and evaluation of turbulence temporal attributes

We have used large-eddy simulation with an immersed boundary method to study turbulent flows over distributions of uniform height, staggered cubes. The computational domains were designed such that both the roughness sublayer and a region of the inertial layer are resolved. With this, we record vertical profiles of time series of fluctuating streamwise and vertical velocity at different locations throughout the domain. Contour images of these fluctuating quantities shown relative to elevation and time are studied; contour images of Reynolds shear stresses owing to ‘sweeps’ and ‘ejections’ are also studied. These images show that periods of momentum excess (deficit) in the inertial-layer precede excitation (subdual) of cube-scale coherent vortices in the roughness sublayer. We compute this time lag (termed advective lag) and demonstrate that it scales linearly with wall-normal elevation. The advective lag is attributed to coherent, low- and high-momentum regions in the aloft inertial layer. Vortex identification is used to illustrate the presence of hairpin packets encapsulating low-momentum regions. Based on this, the reported inclination angle associated with hairpin packets is used to guide the development of a model for prediction of advective lag with height. The model captures the advective lag profiles reasonably well. In the interest of generality, additional cases of flow over homogeneous roughness (aerodynamic drag imposed with the equilibrium logarithmic law) are considered. We again observe that advective lag scales linearly with wall-normal elevation. Advective lag predictions from the aforementioned model agree well with results for these cases.     

湍流煤粉火焰中多相燃烧的跨尺度模拟方法及应用

通过“湍流涡团尺度”与“边界层反应尺度”的关联,建立了湍流煤粉火焰中多相燃烧的跨尺度模拟方法。该方法能够预报湍流脉动宏观规律对静止颗粒边界层内气相反应（如挥发分火焰、CO火焰）的影响。将该方法用于煤粉旋流燃烧数值模拟中,结果显示：与完全忽略边界层气相反应的单膜模型相比,跨尺度模拟的预报结果与Lockwood实验数据有更好的符合。     

气动光学效应湍流脉动模型系数修正

以光学窗口的气动光学效应研究为背景,研究影响光学传输预测的流场湍流脉动量预测.选取平板、压缩拐角、凹腔流动模型,采用大涡模拟（LES）方法研究超声速湍流脉动工程预测模型的系数修正.结果表明：LES方法能够获得无激波干扰、强激波干扰及底部大分离条件下,湍流密度脉动的定量分布,并据此给出脉动工程模型的系数修正,结果已经应用于型号飞行器的光学窗口气动光学效应预测.     

New expressions for the surface roughness length and displacement height in the atmospheric boundary layer

An alternative model for the prediction of surface roughness length is developed. In the model a new factor is introduced to compensate for the effects of wake diffusion and interactions between the wake and roughness obstacles. The experiments are carried out by the use of the hot wire anemometry in the simulated atmospheric boundary layer in a wind tunnel. Based on the experimental data, a new expression for the zero-plane displacement height is proposed for the square arrays of roughness elements, which highlights the influence of free-stream speed on the roughness length. It appears that the displacement height increases with the wind speed while the surface roughness length decreases with Reynolds number increasing. It is shown that the calculation results based on the new expressions are in reasonable agreement with the experimental data.     

Large eddy simulations and parameterisation of roughness element orientation and flow direction effects in rough wall boundary layers

We conduct a series of large eddy simulations (LES) of turbulent boundary layers over arrays of cuboidal roughness elements at arbitrary orientation angles (non-frontal orientations with the incident flow). Flow response to changing roughness orientation is systematically studied at two ground coverage densities, λ_p = 0.06 and 0.11. As expected, the effective roughness heights z_o measured from LES are higher for λ_p = 0.11 than for λ_p = 0.06, although appreciable changes both in z_o and wall shear stress (friction velocity) are observed at both ground coverage densities as the roughness orientation angle changes. This suggests the necessity of accounting for detailed rough wall topology (including more information than just λ_p, λ_f) when relating rough wall morphology to its aerodynamic properties. To this end, a recently developed analytical rough wall parameterisation is used to predict the aerodynamic properties of the simulated rough surfaces. In this rough wall model, wake interactions among roughness elements are explicitly modelled using the concept of sheltering height and exponential attenuation coefficient. As a result, the parameterisation is responsive to detailed ground roughness arrangements and flow conditions, including roughness height variations, element orientation, incident flow direction, transverse displacements, etc. Model-predicted effective roughness heights, wall stress, mean velocity at the height of the roughness, and in some cases displacement height, are compared against the LES measurements from this study as well as numerical/experiment measurements from other authors. The predictions from the model are found to agree well with the measurements both in trends and in absolute values, thus extending the applicability of the analytical rough wall model to more general surfaces than those previously tested.     

Acoustic force model for the fluid flow under standing waves

An acoustic Strouhal number is introduced to demonstrate that the viscosity of fluid can be ignored in the process of sound propagation and acoustic streaming is independent on the frequency of the acoustic wave. Furthermore, acoustic force based on the periodic velocity fluctuation caused by standing acoustic wave is introduced into Navier–Stokes equation in order to describe the fluid flow in the acoustic boundary layer. The numerical results show that the predicted results are consistent with the analytic solution. And the effect of the nonlinear terms cannot be ignored so the analytic solution derived by boundary-velocity condition is only an approximation for acoustic streaming.     

Numerical simulation of the flow within and over an intersection model with Reynolds-averaged Navier--Stokes method

In this study, the Reynolds-averaged Navier--Stokes (RANS) method is employed to simulate the flow within and over an intersection model with three kinds of $k$--$\varepsilon$ turbulence closure schemes, namely, standard model, renormalization group (RNG) model and realizable $k$--$\varepsilon$model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizable $k$--$\varepsilon$ model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable $k$--$\varepsilon$ model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.     

A Web services accessible database of turbulent channel flow and its use for testing a new integral wall model for LES

The output from a direct numerical simulation (DNS) of turbulent channel flow at Re_τ ≈ 1000 is used to construct a publicly and Web services accessible, spatio-temporal database for this flow. The simulated channel has a size of 8πh × 2h × 3πh, where h is the channel half-height. Data are stored at 2048 × 512 × 1536 spatial grid points for a total of 4000 time samples every 5 time steps of the DNS. These cover an entire channel flow-through time, i.e. the time it takes to traverse the entire channel length 8πh at the mean velocity of the bulk flow. Users can access the database through an interface that is based on the Web services model and perform numerical experiments on the slightly over 100 terabytes (TB) DNS data on their remote platforms, such as laptops or local desktops. Additional technical details about the pressure calculation, database interpolation, and differentiation tools are provided in several appendices. As a sample application of the channel flow database, we use it to conduct an a-priori test of a recently introduced integral wall model for large eddy simulation of wall-bounded turbulent flow. The results are compared with those of the equilibrium wall model, showing the strengths of the integral wall model as compared to the equilibrium model.     

Dual solutions in boundary layer flow of Maxwell fluid over a porous shrinking sheet

An analysis is carried out for dual solutions of the boundary layer flow of Maxwell fluid over a permeable shrinking sheet. In the investigation, a constant wall mass transfer is considered. With the help of similarity transformations, the governing partial differential equations（PDEs） are converted into a nonlinear self-similar ordinary differential equation（ODE）. For the numerical solution of transformed self-similar ODE, the shooting method is applied. The study reveals that the steady flow of Maxwell fluid is possible with a smaller amount of imposed mass suction compared with the viscous fluid flow. Dual solutions for the velocity distribution are obtained. Also, the increase of Deborah number reduces the boundary layer thickness for both solutions.     

Turbulence modulation model for gas–particle flow based on probability density function approach

The paper focuses on the turbulence modulation problem in gas–particle flow with the use of probability density function(PDF) approach. By means of the PDF method, a general statistical moment turbulence modulation model without considering the trajectory difference between two phases is derived from the Navier–Stokes equations. A new turbulence production term induced by the dispersed-phase is analyzed and considered. Furthermore, the trajectory difference between two media is taken into account. Subsequently, a new k–ε turbulence modulation model in dilute particle-laden flow is successfully set up. Then, the changes to several terms, including the turbulence production, dissipation, and diffusion terms, are well described consequently. The promoted model provides a more probable explanation for the modification of particles on the turbulence. Finally, we applied the model to simulate a gas–particle turbulence flow case in a wall jet, and found that the simulation results agree well with the experimental data.     

Approximate solution of the flow over a nonlinear magneto-hydrodynamic stretching sheet

The approximate solution of the magneto-hydrodynamic(MHD) boundary layer flow over a nonlinear stretching sheet is obtained by combining the Lie symmetry method with the homotopy perturbation method.The approximate solution is tabulated,plotted for the values of various parameters and compared with the known solutions.It is found that the approximate solution agrees very well with the known numerical solutions,showing the reliability and validity of the present work.     

Analytic solution for magnetohydrodynamic boundary layer flow of Casson fluid over a stretching/shrinking sheet with wall mass transfer

In this analysis,the magnetohydrodynamic boundary layer flow of Casson fluid over a permeable stretching/shrinking sheet in the presence of wall mass transfer is studied.Using similarity transformations,the governing equations are converted to an ordinary differential equation and then solved analytically.The introduction of a magnetic field changes the behavior of the entire flow dynamics in the shrinking sheet case and also has a major impact in the stretching sheet case.The similarity solution is always unique in the stretching case,and in the shrinking case the solution shows dual nature for certain values of the parameters.For stronger magnetic field,the similarity solution for the shrinking sheet case becomes unique.     

On the effect of an anisotropy-resolving subgrid-scale model on large eddy simulation predictions of turbulent open channel flow with wall roughness

A large eddy simulation (LES) was conducted of turbulent flow in a channel with a rough wall on one side and a free surface on the other by adopting an anisotropy-resolving subgrid-scale (SGS) model. A shear Reynolds number of Re_τ = 395 was used based on the mean friction velocity and channel height. To investigate the grid dependency of the LES results caused by the SGS model, three grid resolutions were tested under the same definition of a roughness shape by using the immersed boundary method. The results obtained were compared with direct numerical simulation data with and without the wall roughness and those without the extra anisotropic term. The primary focus was on how the present anisotropic SGS model with coarser grid resolutions can properly provide the effects of roughness on the mean velocity and turbulent stresses, leading to a considerable reduction of the computational cost of LES.     

ENSO事件随机动力学模型的渐近分析

利用边界层型函数，研究了ENSO事件随机动力学的某一模型，给出了这一问题的n阶渐 近展开式，将相关结论应用于特殊的ENSO事件，并得到了零阶渐近解，为分析ENSO事件的变 化状态提供了依据. 关键词： ENSO事件 边界层型函数 渐近展开式     

大气边界层白天温度测量用转动拉曼激光雷达

设计了一个转动拉曼激光雷达系统,对大气边界层温度进行全天候高精度测量。系统选用波长355 nm的紫外激光作光源,采用高光谱分辨力光栅,将雷达接收到的散射信号以及太阳背景光从空间上分离,配合边缘反射镜,反射转动拉曼信号,去除大部分米氏-瑞利散射和太阳背景光噪声信号,并用两个窄带干涉滤光片,分离中心波长为353.9 nm和353.1 nm转动拉曼散射信号,同时对噪声信号进行2次高精度剔除,以保证白天测量的需要。对系统进行了分析及数值计算,结果表明在激光脉冲能量300 mJ,望远镜有效口径25 cm,测量时间10 min的条件下,可以在白天太阳背景光辐射度为3×108W/(m2.sr.nm)的边界层内测量高度2.5 km以下的大气温度分布,并在大气散射比低于2.5的条件下,温度测量精度可达到1 K。     

天文图像空域重建新方法：迭代位移叠加法

在天文斑点成像中 ,基于斑点图中一个新的位移基准点———目标的自相关极大值点 ,提出了一种在空间域中重建天文图像的迭代位移叠加法 (ISA法 ) ,它能有效地消除大气湍流影响 ,实现望远镜衍射受限分辨率天文成像。在阐述了该方法的基本原理之后 ,给出了一些天文目标的像复原实验     

Marangoni flow and mass transfer of power-law non-Newtonian fluids over a disk with suction and injection

We scrutinize the approximate analytical solutions by the optimal homotopy analysis method (OHAM) for the flow and mass transfer within the Marangoni boundary layer of power-law fluids over a disk with suction and injection in the present paper. Concentration distribution on the surface of a disk varies in a power-law form. The non-Newtonian fluid flow is due to the surface concentration gradient without considering gravity and buoyancy. According to the conservation of mass, momentum and concentration, the governing partial differential equations are established, and the appropriate generalized Kármán transformation is found to reduce them to the nonlinear ordinary differential equations. OHAM is used to access the approximate analytical solution. The influences of Marangoni the number, suction/injection parameters and power-law exponent on the flow and mass transfer are examined.     

A kinetic model for droplet growth in the transition regime

A variant of the moment expansion method, used in an earlier paper to describe the flow of a gas toward an absorbing sphere, is applied to a more realistic model of a droplet condensing from a supersaturated vapor. In the simplest version a spherical droplet absorbs all incoming vapor molecules, but spontaneously emits molecules with a Maxwellian distribution at the droplet temperature and with the corresponding saturated vapor density. From a solution of the stationary linearized Boltzmann equation with these boundary conditions we obtain expressions for the heat and mass currents toward the sphere as a function of the supersaturation and the temperature difference between the droplet and the vapor at infinity. For small droplet radii the known free flow limit is obtained in a natural way. From the calculated expressions for the heat and mass current we derive evolution equations for the radius and temperature of the droplet. The temperature evolves more rapidly and can thus be eliminated adiabatically; the resulting growth curve for the radius shows a sharp transition from a kinetically controlled regime for small radii to a regime dominated by heat conduction for large radii. The effect of incomplete absorption at the surface is also studied. The actual calculations are carried out for Maxwell molecules, with parameters corresponding to argon at 0.65T _c and 100% supersaturation.     

基于辐射传输模型的环境一号卫星CCD相机的水体大气校正方法研究

水体大气校正问题是开展我国环境一号卫星水色遥感定量化应用的关键。针对环境卫星CCD相机的特点,以水气耦合的辐射传输模型构建大气校正参数查找表,研究以地面气象数据辅助的逐像元水体大气校正方法,实现水体离水反射率和遥感反射比的反演。以现场测量数据和MODIS数据为参考进行水体大气校正效果验证,研究发现CCD相机的反演结果在蓝、绿波段的精度较高而红、近红的反演结果系统偏大。研究结果还表明气溶胶模型是影响水体大气校正精度的重要因素。     

Transfer matrix spectrum for the finite-width Ising model with adjustable boundary conditions: Exact solution

Using the spinor approach, we calculate exactly the complete spectrum of the transfer matrix for the finite-width, planar Ising model with adjustable boundary conditions. Specifically, in order to control the boundary conditions, we consider an Ising model wrapped around the cylinder, and introduce along the axis a seam of defect bonds of variable strength. Depending on the boundary conditions used, the mass gap is found to vanish algebraically or exponentially with the size of the system. These results are compared with recent numerical simulations, and with random-walk and capillary-wave arguments.