Averaging methods for numerical simulations of flows through heterogeneous porous media

Many natural rock systems contain small patches of different permeability which affect the flow of fluids through them. As these heterogeneities become smaller and more numerous, they become harder to model numerically. We consider how to reduce the computational effort required in simulations by incorporating their effects in the boundary conditions at the edges of each grid block. This is in contrast with current methods which involve often arbitrary changes in the fluid properties. The method is restricted to the case of widely-spaced patches, which simplifies interaction effects. The system then reduces to an array of dipoles, and two averaging methods are proposed for finite grid blocks. Several infinite systems, including vertical and horizontal bands, are also considered as further approximations. There is a great wealth of existing results from different fields which lead to identical mathematical problems and which can be used in these cases. Finally, we consider how to use these techniques when the precise configuration of the grid block is not known, but only its statistical properties. This can lead to results which are very different from the deterministic case.     

Analysis of dynamic features in intersecting pedestrian flows

This paper focuses on the simulation analysis of stripe formation and dynamic features of intersecting pedestrian flows.The intersecting flows consist of two streams of pedestrians and each pedestrian stream has a desired walking direction.The model adopted in the simulations is the social force model, which can reproduce the self-organization phenomena successfully. Three scenarios of different cross angles are established. The simulations confirm the empirical observations that there is a stripe formation when two streams of pedestrians intersect and the direction of the stripes is perpendicular to the sum of the directional vectors of the two streams. It can be concluded from the numerical simulation results that smaller cross angle results in higher mean speed and lower level of speed fluctuation. Moreover, the detailed pictures of pedestrians' moving behavior at intersections are given as well.     

HL-2A 装置的边缘参数测量

HL-2A 装置中平面边缘的等离子体特性通过磁力传动的马赫/ 雷诺协强/ 朗缪尔10 探针组进行了研究。10 探针组安装在可径向向里和向外移动, 并可绕轴旋转360^o 的传动杆上, 用于测量主等离子体边缘的温度、密度、悬浮电位、空间电位、径向和极向电场、湍流的雷诺协强、径向和极向等离子体流速及其径向分布。HL- 2A 装置的实验结果表明, 边缘等离子体扰动诱发的雷诺协强产生了边缘极向流; 雷诺协强的径向梯度驱动带状流抑制了湍流输运。     

Momentum transport and torque scaling in Taylor-Couette flow from an analogy with turbulent convection

We generalize an analogy between rotating and stratified shear flows. This analogy is summarized in Table 1. We use this analogy in the unstable case (centrifugally unstable flow vs. convection) to compute the torque in Taylor-Couette configuration, as a function of the Reynolds number. At low Reynolds numbers, when most of the dissipation comes from the mean flow, we predict that the non-dimensional torque G = T/ν² L, where L is the cylinder length, scales with Reynolds number R and gap width η, G = 1.46η^3/2(1 - η)^-7/4 R ^3/2. At larger Reynolds number, velocity fluctuations become non-negligible in the dissipation. In these regimes, there is no exact power law dependence the torque versus Reynolds. Instead, we obtain logarithmic corrections to the classical ultra-hard (exponent 2) regimes: G = 0.50 . These predictions are found to be in excellent agreement with avail-able experimental data. Predictions for scaling of velocity fluctuations are also provided. Received 7 June 2001 and Received in final form 7 December 2001     

等厚薄板翼型的气动力特性研究

本文从吹风实验及数据分析两方面来研究等厚薄板翼型的气动力特性。在实验中对单园弧、双圆弧、抛物线三种类型的薄板翼型进行了孤立翼型的吹风试验，得出了各翼型升力系数，阻力系数随冲角变化的结果。然后利用非交错网格下的ＳＩＭＰＬＥ方法，计算了等厚薄板翼型流场，计算结果和实验结果吻合较好。     

平板混合层流动中两相脉动关联矩的大涡模拟

在气粒两相平板混合层流动中,对气相流动采用大涡模拟,对颗粒相流动采用轨道模拟,研究了两相脉动关联矩。由两相各自瞬时速度出发可以直接获得两相脉动关联矩的统计结果。模拟获得的颗粒相一阶矩、二阶矩以及两相脉动关联矩与实验结果定量符合,表明基于细观模拟是不通过模型获得和研究两相脉动关联矩的可行途径。     

A pair of stream functions for three-dimensional vortex flows

Introducing a vector potential, that is based on a pair of stream functions, and a velocity potential, antisymmetric equations for the stream functions are derived with the help of a variational principle. It is found that the equations are in a suitable form to investigate flows with helical symmetry, and, for example, to connect upstream axisymmetric flows with downstream helical flows. The special case of a transition from an upstream solid-body vortex to a downstream helical flow is investigated in detail. Furthermore, the stream-function equations are particularly useful to investigate general small-amplitude inertia waves on vortex flows. Time-dependent helical flows that are time-independent in a suitably rotating frame of reference can also be discussed with the proposed method.     

Determination of shock wave velocity from photographs obtained using a polarization interferometer

A polarization shear interferometer is described in which one of the two images of the flow being studied is produced outside the interference field and serves as a shadow pattern. This image is recorded again using an additional light source — a pulsed laser. In this way one can determine the shock wave velocity and establish, on the basis of gas-dynamic relationships, the proper numeration of the interference fringes on the interferogram.     

Dynamics of the thermocline in the equatorial region of the Pacific ocean

This paper is devoted to the subsurface current dynamics in equatorial regions, where the hallmark of a strong stratification is a sharp interface (thermocline), separating two layers of different density, and whose depth is dependent upon the strength of the winds above the ocean's surface. We give here a few monotonicity results concerning the dynamics of the thermocline in the equatorial region. The most important one asserts that the level of the thermocline decreases as the strength of the wind at ten meters above the ocean surface, denoted |U_w|, increases. Moreover, the strength of the current at the thermocline decreases as |U_w| increases.     

用统一二阶矩模型模拟强旋湍流气粒两相流动的初探

本文初步探索用统一二阶矩（USM）两相湍流模型及K-ε-kp模型预报旋流数为1.5的轴向-切向进风、轴向供粉的强旋气粒两相流动，预报结果与实验的对比表明，USM模型在预报强族流动两相时平均切向速度场上优于K-ε-Kp模型，并且能够合理地给出与实验定性一致的两相湍流脉动的各向异性．