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.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

含有集约束向量极值问题的最优性条件

在序线性拓扑空间里研究了含有集约束向量极值问题的最优性条件,并建立了充分性和必要性条件.     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

海浪的条件特征周期

在海浪波面高度为正态分布的假定下,导出一种以给定波高为条件的条件周期概率密度函数.与风浪槽中测量数据比较,结果表明,在窄谱情况下此概率密度函数与实验室风浪的实际符合良好.根据此密度函数定义了3种条件特征周期,并导出它们与平均周期的关系式.根据这些关系对有关海洋工程上的一些问题作了解释和讨论.     

Internal generation of waves: Delta source function method and source term addition method

In this study, we investigate two internal wave generation methods in numerical modeling of time-dependent equations for water wave propagation, i.e., delta source function method and source term addition method, the latter of which has been called the line source method in literatures. We derive delta source functions for the Boussinesq-type equations and extended mild-slope equations. By applying the fractional step splitting method, we show that the delta source function method is equivalent to the source term addition method employing the energy velocity. This suggests that the energy velocity should be used rather than the phase velocity for the transport of incident wave energy in the source term addition method. Finally, the performance of the delta source function method is verified by accurately generating nonlinear cnoidal waves as well as linear waves for horizontally one-dimensional cases.     

A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

Modelling of periodic wave transformation in the inner surf zone

In this paper, we analyse the ability of the nonlinear shallow-water (NSW) equations to predict wave distortion and energy dissipation of periodic broken waves in the inner surf zone. This analysis is based on the weak-solution theory for conservative equations. We derive a new one-way model, which applies to the transformation of non-reflective periodic broken waves on gently sloping beaches. This model can be useful to develop breaking-wave parameterizations (in particular broken-wave celerity expression) in both time-averaged wave models and time-dependent Boussinesq-type models. We also derive a new wave set-up equation which provides a simple and explicit relation between wave set-up and energy dissipation. Finally, we compare numerical simulations of both, the NSW model and the simplified one-way model, with spilling wave breaking experiments and we find a good agreement.     

浙江西沪港重金属铜的配位容量和形态分析

根据2000年6月10～14日在浙江省西沪港采集的海水样品,利用AA-800石墨炉原子吸收分光光度技术和阳极溶出伏安法测定样品中重金属铜的含量,获得铜在海水中受不同的有机配体控制.不同粒级的铜表观配位容量表明西沪港海水过孔径1.00μm微孔膜的(ACuCC)较高,为144.4nmol/dm3;过0.40和0.20μm滤膜的(ACuCC)分别为103.0和102nmol/dm3;铜的有机配体条件稳定常数的对数值在7.25～9.14之间.铜的总量为21.72nmol/dm3.铜全部为稳定溶解态,其中pH2酸溶态占95.0%,强有机结合态占5.0%.溶解态铜中有机结合态占过滤海水中总铜的61.6%.     

关于波浪Boussinesq方程的研究

对有关波浪 Boussinesq方程的研究成果进行了系统的归纳总结和评述 ,以期对本学科的发展起到一定的引导和促进作用