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.     

对月日潮汐摩擦耗散的估计

月日潮汐摩擦和地球惯量矩变化是日长长期变化的主要原因．在本文中，利用最新的地球物理和古生物钟数据，对过去15亿年以来的月日潮汐摩擦、地球惯量矩变化和日长长期变化等作了数值对比研究．由此得到二个重要结论：一是仅利用地球的自转形变不能解释J2的变化，这说明地球的重力分异现象至今仍存在着；其二是在几亿年前的潮汐摩擦比现在大得多，若取潮汐耗散与距离的立方成反比时，理论结果与由古生物钟得到的回归年日数和朔望月日数数据较为符合。     

Influence of the Atmospheric Tides on the Earth Rotation

Spectral analysis of the components of the relative atmospheric angular momentum vector is performed based on the series of these components for the 6 h intervals within the period of 1958–2000. These series have been computed in the Subbureau of the Atmospheric Angular Momentum of the International Earth Rotation Service using the NCEP/NCAR reanalysis of atmospheric observations. The basic harmonics of diurnal tides are determined. New results on the fortnight's and week's duration oscillations of the equatorial components of the atmospheric angular momentum are obtained. The zonal tides transformation mechanisms in the atmosphere are discussed. It is shown that the main mechanism of the zonal tides effect on the atmospheric variability is the amplitude modulation of daily oscillations of the relative atmospheric angular momentum. The effects of the atmospheric tides on the Earth rotation are discussed.     

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.     

低风速下海面附近气温场对海气动量通量的影响

实用测的海上和层平均风剖面数据和温度剖面数据，通过数据回归和迭代方法计算出了在不同大气稳定情况下的海面阻力系数。得到了与前人理论计算一致的结论：海面阻力系数随海面大气稳定度的增加而减小，另外，我们还发现：在海面风速小于１３ｍ／ｓ时，不能认为气温剖面外推到海面的值与海面水温的是一致的。这样若用海气温差作为衡量海面上方大气的稳定程度，难于得到上面给出了理结论。这一点同前人的理论计算结果是不相同的。     

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.     

Gravity torques for surface waves

The effects of the gravity torques acting on the angular momentum of surface gravity waves are calculated theoretically. For short crested waves the gravity torque is caused by the force of gravity on the orbiting fluid particles acting down the slopes of the crests and troughs and in the direction parallel to the crests and troughs. The gravity torque tries to rotate the angular momentum vectors, and thus the waves themselves, counterclockwise in the horizontal plane, as viewed from above, in both hemispheres. The amount of rotation per unit time is computed to be significant assuming reasonable values for the along-crest and trough slopes for waves in a storm area. The gravity torque has a frequency which is double the frequency of the waves. For long crested waves the gravity torque acts in the vertical plane of the orbit and tries to decelerate the particles when they rise and accelerate them when they fall. By disrupting the horizontal cyclostrophic balance of forces on the fluid particles (centrifugal force versus pressure force) the gravity torque accounts qualitatively for the three characteristics of breaking waves: that they break at the surface, that they break at the crest, and that the crest breaks in the direction of wave propagation.