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.     

STUDY ON GIS FOR YIELD ESTIMATION BY REMOTE SENSING IN JILIN MAIZE BELT

ＳＴＵＤＹ　ＯＮ　ＧＩＳ　ＦＯＲ　ＹＩＥＬＤ　ＥＳＴＩＭＡＴＩＯＮ　ＢＹ　ＲＥＭＯＴＥ　ＳＥＮＳＩＮＧ　ＩＮ　ＪＩＬＩＮ　ＭＡＩＺＥ　ＢＥＬＴＳＴＵＤＹＯＮＧＩＳＦＯＲＹＩＥＬＤＥＳＴＩＭＡＴＩＯＮＢＹＲＥＭＯＴＥＳＥＮＳＩＮＧＩＮＪＩＬＩＮＭＡＩＺ...     

Multiple bipolar structure of normal magnetic field at the magnetopause

We describe how a local plasma structure can be changed by a transverse shear flow using numerical simulation to investigate the disturbance process near the magnetopause. The results show that magnetic field lines are bent by transverse shear flow disturbance near the current sheet region. There are multiple bipolar structures of the normal magnetic field in the numerical simulation. We term this new feature as K-point magnetic reconnection, realistic for discussing space observations.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 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.     

Diffuse radiation,twilight, and photochemistry — I

A photochemical scheme which includes a detailed treatment of multiple scattering up to solar zenith angles of 96° (developed for use in a GCM) has been used to study partitioning within chemical families. Attention is drawn to the different zenith angle dependence of diffuse radiation for the two spectral regions <310 nm and >310 nm. The effect that this has on the so-called 40 km ozone problem is discussed. The importance of correctly including multiple scattering for polar ozone studies is emphasised.     

Multiple discriminant analysis of the Drava River alluvial plain sediments

Three discriminant function models are raised and cross-compared in order to distinguish geochemical patterns characteristic for the Drava River floodplain sediments. Based on data representing total element concentrations in samples collected from alluvium (A), terrace (T), and unconsolidated bedrock (B) at the border of a floodplain, four element clusters emerged accounting for discrimination between the referred groups of sediments. The most prominent is contaminant/carbonate cluster characteristic for alluvium. The other two are: silicate cluster typical for unconsolidated geological substrate (Neogene sedimentary rocks); and naturally dispersed heavy metal cluster separating terrace from the former two groups. Models introducing depth intervals and single profiles as grouping criteria reveal identical sediment-heavy metal matrices. The second important issue of this paper is possibility of reclassification of samples originally assigned to one of the a priori defined groups of sediments, based on established geochemical pattern. The mapped geological units can be reconsidered by the post hoc assignments to a different group if geological border between alluvium and terrace or between terrace and bedrock can not be established geologically with absolute certainty.     

引水工程秦巴段隧洞地应力模拟及工程地质问题

利用Ansys有限元软件对三峡引水工程秦巴段线路不同深度、不同截面形态隧洞围岩的应力重分布情况进行模拟计算,得到圆形隧洞、城门形隧洞和马蹄形隧洞围岩的应力数值和等值线图.利用库仑一纳维尔强度准则,对花岗岩区及灰岩区隧洞围岩的剪切破坏进行分析;利用三种岩爆应力判别指标,对岩爆进行初步预测;进而对不同深度、不同截面形态的隧洞围岩的稳定性进行分析;最后对引水工程隧洞设计中截面形态的选择给出了初步建议.     

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.