Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

分组密码的安全性与强化设计

分组密码是数据通讯中最常用的数据加密方式，以DES为例分析现有分组加密算法的安全隐患，并提出了可变密钥加密和变长密文输出两个新思路，可应用于所有现有分组加密算法以提高安全性，并就该方法的安全性、效率、具体应用做出了分析。     

About Darcy's law in non-Galilean frame

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.     

Cokriging nonstationary data

Universal cokriging is used to obtain predictions when dealing with multivariate random functions. An important type of nonstationarity is defined in terms of multivariate random functions with increments which are stationary of orderk. The covariance between increments of different variables is modeled by means of the pseudo-cross-covariance function. Criteria are formulated to which the parameters of pseudo-cross-covariance functions must comply so as to ensure positive-definiteness. Cokriging equations and the induced cokriging equations are given. The study is illustrated by an example from soil science.     

低频扩频中伪随机序列产生方法的研究

伪随机序列的产生是实现BPC低频时码发播台扩展频谱技术升级的一个重要环节。针对BPC低频时码发播台的现有系统，提出了伪随机序列的产生方案，并对其性能作了分析。     

Spectral Analysis of Nonlinear Drag Forces

The spectral properties of nonlinear drag forces of random waves on vertical circular cylinders are analyzed in this paper by means of nonlinear spectral analysis. The analysis provides basic parameters for estimation of the characteristic drag forces. Numerical computation is also performed for the investigation of the effects of nonlinearity of the drag forces.The results indicate that the wave drag forces calculated by linear wave theory are larger than those calculated by the third order Stokes wave theory for given waves. The difference between them increases with wave height. The wave drag forces calculated by use of hnear approximation are about 5% smaller than their actual values when measured in the peak values of spectral densities. This will result in a safety problem for the design of offshore structures. Therefore, the nonlinear effect of wave drag forces should be taken into comidemtion in design and application of important offshore structures.     

声波在水-多孔介质海底界面上的反射与透射

声波在海底界面的反射和透射是海底散射、海底混响、海底目标探测的重要问题。利用Biot多孔介质声传播理论对声波在水-多孔海底界面上的反射和透射进行了分析,具体给出反射声波的反射系数,3种透射声波的透射系数以及声能透射系数随入射波入射角和频率(10~40 kHz)的变化关系,分析了各种透射波对透射声能的贡献。多孔海底介质参数分别采用Stoll和Chotiros给出的2组参数进行理论计算。     

随机海浪的列率谱分析研究

首次将列率谱分析的方法引入到随机海浪的研究。根据不同环境条件下实验室资料，对列率谱与频率谱进行了比较。结果表明，列率频对频率谱高频部分有明显改善，特别对二倍频含能段的分辨具有明显的优越性。另外，列率谱具有计算速度快、简单易行、分辨率高等优点。