A One-Dimensional Finite-Element Boundary-Layer Model with a Vertical Adaptive Grid

A one-dimensional atmospheric boundary-layer model is developed using the finite element method and a 1.5-order e-l turbulence closure scheme. A vertical adaptive strategy is implemented, based upon an heuristic error estimator that depends upon properties of the layer, such as the stratification. The model is used to simulate a moderately stratified stable boundary layer as described in the GABLS (GEWEX Atmospheric Boundary-Layer Study, where GEWEX is the Global Water and Energy Cycle Experiment) First Intercomparison Project, and then a more complicated diurnal cycle, as used in the GABLS Second Intercomparison Project. In the stable boundary-layer experiment, it is shown that including the adaptive strategy can improve the performance of the model such that the error in the model is significantly less (greater than an order of magnitude with an effective resolution of 8 m) than that of the model without adaptivity. The model’s turbulence closure scheme and the adaptivity strategy also successfully simulate the different stability regimes present in the diurnal cycle simulation, and represented all of the expected features.     

Effects of non‐uniform traction and specimen height in the direct shear test on stress and deformation in a rock fracture

In the direct shear test (DST), an internal moment is distributed within the rock specimen by non‐coaxial shear loads applied to the specimen, which cause non‐uniform distributions of both the traction on the loading planes and the stress and deformation in the specimen. To examine the validity of the DST for a rock fracture and to clarify the effect of specimen height, both the stress and deformation in a fracture in the DST were analyzed for specimens with three different heights using a three‐dimensional finite element method with quadratic joint elements for a fracture model. The constitutive law of the fracture considers the dependence of the non‐linear behavior of closure on shear displacement and that of shear stiffness on normal stress and was implemented in simulation code to give a conceptional fracture with uniform mechanical properties to extract only the effect of non‐uniform traction on the stress and deformation in the fracture. The results showed that both normal and shear stresses are concentrated near the end edges of the fracture, and these stress concentrations decrease with a decrease in the specimen height according to the magnitude of the moment produced by the non‐coaxial shear loads. Furthermore, although closure is greater near the end edges of the fracture, where normal stress is concentrated, this concentration of closure is not so significant within the range of this study because of the non‐linear behavior of closure, that is, closure does not significantly increase with an increase in normal stress at large normal stresses. Copyright © 2012 John Wiley & Sons, Ltd.     

Diurnal cycle of precipitation over the Maritime Continent simulated by a spectral cumulus parameterization

The validity of a spectral cumulus parameterization (spectral scheme) for simulating a diurnal cycle of precipitation over the Maritime Continent (MC) was examined using a regional atmospheric model. The impacts of entrainment parameterization and each type of convective closure, i.e., non-equilibrium (or equilibrium) closure for deep convection, mid-level, and shallow convective closures, were also examined. When vertically variable entrainment and appropriate convective closures were employed, the model adequately simulated a diurnal cycle of precipitation over both land and ocean as compared to the observation. Analysis regarding the entrainment parameterization revealed that variable entrainment parameterization was needed not only for simulating better mean patterns of precipitation, but also for more realistic phases of diurnal cycles. The impacts of convective closures appeared in the differences in the precipitation amplitude. Analysis on diurnal cycles of convective properties and tendencies revealed that the cycles between boundary layer forcing and convective heating determined convection strength and were affected by each type of convective closure. It can be concluded that the spectral scheme with appropriate convective closures is able to simulate a realistic diurnal cycle over the MC.     

大武水源地地下水环境模拟与石油污染控制研究

大武水源地是山东省内唯一的大型地下水水源地，由于相邻的堠皋地段地下水被石油污染，大武水源地面临着严重的石油污染威胁。本文通过对地下水环境模拟和地下水环境控制模拟的研究，提出在堠皋地段设置 地下水截流带，控制已污染的地下水向大武水源地径流。为达到保护大武水源地的目的，提出石油污染治理建议。     

Earthquake response of hysteretic mass-column using non-Gaussian closure

A non-Gaussian closure scheme based on the Edgeworth expansion of the probability density function is used to study the response of a hysteretic structure under random parametric excitation. The system considered consists of a weightless mass supporting a concentrated mass and it is subjected to the vertical and horizontal components of the ground acceleration modeled as nonstationary Gaussian white noise processes. The material of the structure exhibits bilinear hysteretic behaviour. The equation governing the motion of the system is transformed into an Itô stochastic differential equation. A set of ordinary differential equations governing the response statistics are obtained. These form an infinite hierarchy of equations which must be truncated in order to solve for moments of any order. The Edgeworth expansion of the joint density is used to truncate this infinite hierarchy. Such a closure scheme appears desirable since for hysteretic systems an explicit expression of the probability density is required. A frequently used closure scheme based on Gaussian assumption underestimates the response. The non-Gaussian density can be used in reliability studies.     

城市大气边界层的数值研究

本文利用能量闭合的二维非线性、非定常模式,结合地面热量平衡方程,研究了城市热岛和城市热岛环流以及其它相关的城市边界层参量的发展演变规律;同时探讨了在城市污染面源的情况下,污染物的扩散、输送问题.模拟结果得到了天津市城市热岛观测资料的证实.     

复杂下垫面地域边界层结构的三维细网格数值模拟

建立了一个可供复杂下垫面地域使用的非静力的三维细网格边界层模式，就复杂下垫面条件下的边界层结构和湍流特征作了以实例为对照的数值模拟试验，模式采用能量闭合方案，舍弃了静力近似。以实测资料为初台输入，同时还做了一些数值试验，分别获得了采用静力与非静力模式和不同闭合方案对ＰＢＬ模型的结果。     

A model of the wave and current boundary‐layer structure

A model of the wave and current boundary‐layer structure was developed using the k–ε turbulent closure model. The finite‐difference method was used to solve the governing equations. Vertical logarithmic grids and equal time steps were adopted. The following modelled simulations were obtained: (1) vertical profiles of wave velocity amplitude, eddy viscosity coefficient and turbulent kinetic energy with waves only; (2) vertical profiles of wave velocity amplitude, mean current velocity, eddy viscosity coefficient and turbulent kinetic energy with waves having a following current. To test the validity and the rationality of the present model, vertical profiles of modelled wave velocity amplitude and mean velocity were compared with corresponding experimental results available in the literature. Copyright © 2006 John Wiley & Sons, Ltd.     

Fish Canyon Tuff榍石He扩散及(U- Th)/He年代学研究

榍石富含U和Th,是(U-Th)/He定年的理想矿物之一。本文以Fish Canyon Tuff榍石为例,开展了榍石He扩散行为和榍石(U-Th)/He定年实验方法研究。榍石分步加热扩散实验结果表明He扩散系数ln(D/a²)与温度倒数呈负相关,与期望的热活化扩散过程一致。测试Fish Canyon Tuff榍石(U-Th)/He年龄分布在28.3～24.6 Ma之间,平均值为26.7±1.2 Ma (1σ),Th/U分布在4.6～5.5之间,平均值为5.2±0.2,在误差范围内与国际上已出版数据一致,表明建立的榍石(U-Th)/He定年实验方法可靠。本次测试15粒榍石碎片外表层(～20μm)存在不同程度的磨蚀(即不完整晶体),且榍石表层磨蚀厚度随着等效半径的增加而增大。榍石碎片(U-Th)/He年龄介于完整晶体(U-Th)/He年龄和真实(U-Th)/He年龄之间,且随着榍石等效半径及表层磨蚀厚度(<20μm)的增大,(U-Th)/He年龄更接近真实年龄,这表明榍石(U-Th)/He年龄不确定度与等效半径大小和表层磨蚀厚度有关。     

Process-based model for nearshore hydrodynamics, sediment transport and morphological evolution in the surf and swash zones

Hydrodynamics and sediment transport in the nearshore zone were modeled numerically taking into account turbulent unsteady flow. The flow field was computed using the Reynolds Averaged Navier–Stokes equations with a k–ε turbulence closure model, while the free surface was tracked using the Volume-Of-Fluid technique. This hydrodynamical model was supplemented with a cross-shore sediment transport formula to calculate profile changes and sediment transport in the surf and swash zones. Based on the numerical solutions, flow characteristics and the effects of breaking waves on sediment transport were studied. The main characteristic of breaking waves, i.e. the instantaneous sediment transport rate, was investigated numerically, as was the spatial distribution of time-averaged sediment transport rates for different grain sizes. The analysis included an evaluation of different values of the wave friction factor and an empirical constant characterizing the uprush and backwash. It was found that the uprush induces a larger instantaneous transport rate than the backwash, indicating that the uprush is more important for sediment transport than the backwash. The results of the present model are in reasonable agreement with other numerical and physical models of nearshore hydrodynamics. The model was found to predict well cross-shore sediment transport and thus it provides a tool for predicting beach morphology change.