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.     

三峡重庆库区深部地球物理特征与断裂构造

为了深入研究三峡重庆库区岩石圈动力学特征及其对断裂构造活动的控制和影响．并为之提供基础资料，为三峡重庆库区地震、地质灾害的监测与防治提供基础依据，在已有地球物理资料的基础上，从综合地球物理研究角度出发，通过实测地震测线资料的再解释．采用新的处理技术方法，对本区的东西向主剖面和南北向支测线的地震测深资料进行二维射线追踪处理、Pg波成像；选用场分离技术、位移数字成像技术重新处理了重力和航磁资料，通过联合反演来建立深部二维构造剖面，对剖面所揭示的基底构造特征和地壳结构特征、主要断裂构造特征以及莫霍面的起伏特征进行了精细分析和细致研究．从地球物理平面场特征出发建立了岩石圈构造三维框架。研究结果表明，这样的研究思路准确、方法得当．结论可靠；沿该剖面，把可解译的断裂分为Ⅲ级：Ⅰ级为超岩石圈断裂；Ⅱ级为壳断裂；Ⅲ级为盖层断裂。依据地球物理特征，准确揭示了库区的断裂构造特征，达到了预期效果。     

Sequences in the Hardness Ratio-Peak Energy Plane of Gamma-Ray Bursts

The narrowness of the distribution of the peak energy of the vFv spectrum of gamma-ray bursts (GRBs) and the unification of GRB populations are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB populations, the long GRBs, the short GRBs, and the X-ray flashes (XRFs), in the HR - Ep plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that the long GRBs and the XRFs observed by HETE-2 seem to follow the same sequence in the HR-Ep plane, with the XRFs at the low end of this sequence. We fit the sequence by a universal Band function, and find that this sequence is mainly defined by the low energy index α, and is insensitive to the high energy index,β. With fixed β=-5, a best fit is given by α=-1.00 with X^2min/dof=2.2. The long and short GRBs observed by BATSE follow significantly different sequences in the HR-Ep plane, with most of the short GRBs having a larger hardness ratio than the long GRBs at a given Ep. For the long GRBs a best-fit yields α=-0.30 and β=-2.05. For the short GRBs, a best fit gives α=-0.60 with x^2min=1.1 (withβfixed at-2.0 because it is numerically unstable). The α value for the short GRBs is significantly greater than that for the long GRBs. These results indicate that the global spectral behaviors of the long GRB sample and the XRF sample are similar, while that of the short GRBs is different. The short GRBs seem to be a unique subclass of GRBs, and they are not the higher energy extension of the long GRBs.     

CygX—1硬态高能光子的时延

硬X射线和软射线光子的时延时研究主能辐射的一个重要方面，傅里叶交叉谱常常被用于教育处硬软光子之间的时延，但是交叉说示能在高于30Hz的傅里叶频率上从CygX-1测量到有统计意义的快速光变，由李惕培在时域上发展的交叉相关函方法能有效测量不同时间尺度上的时延，利用此交叉相关分析方法，讨论了不同观测时期CygX-1硬态高能光子时延 的性质，分析结果表明，CygX-1硬态在短的尺度（＜0．1秒）上存存在有意义的时间延迟，并将短时间尺度上的观测结果和各种CygX-1硬态模型进行了比较讨论。     

Numerical study of the motions and drift force of a floating OWC device

The motion and the drift force of a floating OWC (oscillating water column) wave energy device in regular waves are studied taking account of the oscillating surface-pressure due to the pressure drop across the duct of the air chamber. The potential problem inside the chamber is formulated by making use of the Green integral equation associated with the Rankine-type Green function while the outer problem with the Kelvin-type Green function. The added mass, wave damping and excitation coefficients as well as the motion and drift force of the OWC device are calculated for various values of parameter related to the pressure drop.     

Prediction of local scour around circular piles under waves using a novel artificial intelligence approach

Abstract

The scour phenomena around vertical piles in oceans and under waves may influence the structure stability. Therefore, accurately predicting the scour depth is an important task in the design of piles. Empirical approaches often do not provide the required accuracy compared with data mining methods for modeling such complex processes. The main objective of this study is to develop three data-driven methods, locally weighted linear regression (LWLR), support vector machine (SVR), and multivariate linear regression (MLR) to predict the scour depth around vertical piles due to waves in a sand bed. It is the first effort to develop the LWLR to predict scour depth around vertical piles. The models simulate the scour depth mainly based on Shields parameter, pile Reynolds number, grain Reynolds number, Keulegan–Carpenter number, and sediment number. 111 laboratory datasets, derived from several experimental studies, were used for the modeling. The results indicated that the LWLR provided highly accurate predictions of the scour depths around piles (R?=?0.939 and RMSE = 0.075). Overall, this study demonstrated that the LWLR can be used as a valuable tool to predict the wave-induced scour around piles.     

地震折射波的虚拟射线理论及其应用

本文提出了一种地震折射液的虚拟射线理论,根据该理论,只要已知地表层的速度,即可直接由折射波信息提取地震参数,从而可实现折射界面的反演。 文中通过实际介质模型的计算机实验结果,验证了该理论的正确性。     

Wave run-up and overtopping on smooth and rough slopes of coastal structures

A series of hydraulic model tests has been carried out in a glass wave flume to investigate the influences of wave height, wave period, wave steepness, surf similarity parameter, roughness, layer thickness and porosity on wave run-up and overtopping of 1:2 sloped impermeable and permeable breakwaters fronted by a 1:10 gentle, smooth beach slope. The analysis of results involves the correlation between the overtopping energy transfer with the relative wall height and the relationship between wave run-up and overtopping rate. Further, measured wave run-up and overtopping rates are compared with the results given in the Shore Protection Manual (1984), Automated Coastal Engineering System (1992)and results of other investigators.     

Modeling of coupled tide–wave–surge process in the Yellow Sea

A coupled wave–tide–surge model has been developed in this study in order to investigate the effect of the interactions among tides, storm surges, and wind waves. The coupled model is based on the synchronous dynamic coupling of a third-generation wave model, WAM cycle 4, and the two-dimensional tide–surge model. The surface stress, which is generated by interactions between wind and wave, is calculated by using the WAM model directly based on an analytical approximation of the results using the quasi-linear theory of wave generation. The changes in bottom friction are created by the interactions between waves and currents and calculated by using simplified bottom boundary layer model. In consequence, the combined wave–current-induced bottom velocity and effective bottom drag coefficient were increased in the shallow waters during the strong storm conditions.     

Dynamic analysis of a vertical plate exposed to breaking wave impact

From the experimental studies in recent years, it has become known that when a wave breaks directly on a vertical faced coastal structure, high magnitude impact pressures are produced. The theoretical and experimental studies show that the dynamic response of such structures under wave impact loading is closely dependent on the magnitude and duration of the load history. The dynamic analysis and design of a coastal structure can be succeeded provided the design load history for the wave impact is available. Since these types of data are very scarce, it is much more convenient to follow a method which is based on static analysis for the dynamic design procedure. Therefore, to facilitate the dynamic design of a vertical plate that is exposed to breaking wave impact, a multiplication factor called “dynamic magnification factor” is herein presented which is defined as the ratio of the maximum value of the dynamic response to that found by static analysis. The computational results of the present study show that the dynamic magnification factor is a useful ratio to transfer the results of static analysis to the dynamic design of a coastal plate for the maximum impact pressure conditions of p_max/γH₀≤18.