Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.     

集合卡尔曼滤波同化中雷达位置的敏感性研究

针对对流尺度集合卡尔曼滤波（EnKF）雷达资料同化中雷达位置对同化的影响进行研究。为了考察强对流出现在雷达不同方位时集合卡尔曼滤波同化雷达资料的能力,以一个理想风暴为例,设计了8个均匀分布在模拟区域周围的模拟雷达进行试验。单雷达同化试验中,初期同化对雷达位置较敏感,而十几个循环后对雷达方位的敏感性降低。造成初期同化效果较差的雷达观测位于模拟区域正南和正北方向,这两部雷达与模拟区域中心的连线垂直于风暴移动方向（即环境气流的方向）。双雷达试验的结果表明,正东、正南、正西和正北方向的雷达组合观测会使同化初期误差较大,这说明并不是所有与风暴连线成90°的雷达组合都能在短时同化中得到合理的分析结果,还需要都处于模拟区域对角线上（即与环境气流成45°夹角）,同化效果才较好。短时同化后的确定性预报结果表明,较大分析误差也会导致较大预报误差。这些分析误差主要是由于同化初期不准确的集合平均场驱动出的不合理的背景误差协方差造成的。当背景场随着同化循环得到改进后,驱动出的合理的背景误差协方差使得不同位置雷达同化造成的差异逐步减小。基于上述结果,引入迭代集合均方根滤波（iEnSRF）算法,结果显示使用该算法后,雷达位置对同化效果的影响减小,同化不同位置的雷达资料均能有效降低分析和预报误差。      

A TWO-LAYER PRIMITIVE EQUATION OCEAN CURRENT MODEL AND ITS SIMULATION

A two-layer primitive equation model is developed in this paper. The capabilities of this model aretested by the use of multiyearly averaged January and July sea surface level pressure fields and windfields which can be diagnosed from the pressure fields. The results show that the ocean surface currentsand undercurrents in the second layer driven by the sea surface wind and the sea surface pressure areclose to the observation. The results are also compared with that of the IAP OGCM and the OSUOGCM.     

早新生代温室气候及冰期气候转型的模拟研究

早新生代是地质史上最后一个温室气候期,随后南极冰盖形成,地球进入到晚新生代冰期。温室气候的成因和冰期气候转型的机制一直是国际相关学界关注的问题。评述国际上对此开展的古气候模拟,反映了早新生代温室气候受到了海洋和大陆的地理位置、暖海洋温盐环流和海洋热输送、太阳辐射和大气CO2浓度变化的作用和影响。古气候模拟还反映了早新生代温室气候转向冰期气候,受到了大洋通道改变和高原构造隆起、大气成分变化以及海陆生态系相互的作用和反馈。这些古气候模拟试验锁定在气候变化的关键时段和重要驱动因子,对测试地球内外驱动力和地球各圈层反馈作用提供了重要的科学依据;温室气候以及趋向冰期气候的模拟研究对探讨气候变化内在机制、预测未来气候具有重要意义。      

Numerical modeling investigation on turbulent oscillatory flow over a plane rough bed composed by randomly arrayed particles

A three-dimensional numerical model is established to simulate the turbulent oscillatory boundary layer over a fixed and rough bed composed by randomly arrayed solid spheres based on the lattice Boltzmann method and the large eddy simulation model.The equivalent roughness height,the location of the theoretical bed and the time variation of the friction velocity are investigated using the log-fit method.The time series of turbulent intensity and Reynolds stress are also investigated.The equivalent roughness height of cases with Reynolds numbers of 1×10~4–6×10~4 is approximately 2.81 d(grain size).The time variation of the friction velocity in an oscillatory cycle exhibits sinusoidal-like behavior.The friction factor depends on the relative roughness in the rough turbulent regime,and the pattern of solid particles arrayed as the rough bed in the numerical simulations has no obvious effect on the friction factor.     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.     

Probabilistic stability analyses of undrained slopes by 3D random fields and finite element methods

A long slope consisting of spatially random soils is a common geographical feature. This paper examined the necessity of three-dimensional(3 D) analysis when dealing with slope with full randomness in soil properties. Although 3 D random finite element analysis can well reflect the spatial variability of soil properties, it is often time-consuming for probabilistic stability analysis. For this reason, we also examined the least advantageous(or most pessimistic) cross-section of the studied slope. The concept of"most pessimistic" refers to the minimal cross-sectional average of undrained shear strength. The selection of the most pessimistic section is achievable by simulating the undrained shear strength as a 3 D random field. Random finite element analysis results suggest that two-dimensional(2 D) plane strain analysis based the most pessimistic cross-section generally provides a more conservative result than the corresponding full 3 D analysis. The level of conservativeness is around 15% on average. This result may have engineering implications for slope design where computationally tractable 2 D analyses based on the procedure proposed in this study could ensure conservative results.     

Comparison of chaos detection methods in the circular restricted three‐body problem

In this study, a sample of orbits is considered in the framework of the planar circular restricted three‐body problem. In order to separate ordered from chaotic orbits three numerical methods are compared: the Largest Lyapunov Characteristic Exponent (LLCE) and the Smaller Alignment Index (SALI) provide a fairly good characterization of the chaotic motions, while the computational time required is of the same order; the Correlation Dimension (CD) has the advantage of correctly classifying sticky orbits, but at the expense of a longer computational time. In order to classify a given orbit, any pair of the three methods can be considered, but LLCE and SALI are recommended due to their speed. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

蒙特卡洛模型在现代黄河三角洲浅层沉积物压实速率模拟中的应用

在进行黄河三角洲浅层沉积物固结压实研究中,过去传统的方法通常基于钻孔和剖面,进行局部计算分析。为克服钻孔分布不均的局限性需研究整个现代黄河三角洲地区浅层沉积物的压实沉降特征,基于1987年黄河三角洲地区综合工程地质勘察获取的钻孔资料,采用了数学地质中的蒙特卡洛方法,结合土力学的分层总和法对整个现代黄河三角洲浅层沉积物的压实速率进行模拟,并分析了压实模拟的具体影响因素。结果表明当前浅层沉积物的90%累积概率(P90)的压实速率变化范围为0.18~9.07 mm/a,影响模拟沉积层的压实速率主要因素为沉积物的初始孔隙比、压缩系数和平均沉积速率,上述要素均与压实速率呈显著正相关,软土层应是浅沉积地层压实沉降的主要贡献层。提出按照时间序列分别对不同流路时期叶瓣进行模拟,更符合现代黄河三角洲地区的沉积环境和历史,并对比23个地面沉降水准监测数据,定量分离了浅层沉积物固结压实分量,其约占沉降总量的13%。     

Numerical simulations of a highway bridge structure employing passive negative stiffness device for seismic protection

A new passive seismic response control device has been developed, fabricated, and tested by the authors and shown to be capable of producing negative stiffness via a purely mechanical mechanism, thus representing a new generation of seismic protection devices. Although the concept of negative stiffness may appear to be a reversal on the desired relationship between the force and displacement in structures (the desired relationship being that the product of restoring force and displacement is nonnegative), when implemented in parallel with a structure having positive stiffness, the combined system appears to have substantially reduced stiffness while remaining stable. Thus, there is an ‘apparent weakening and softening’ of the structure that results in reduced forces and increased displacements (where the weakening and softening is of a non‐damaging nature in that it occurs in a seismic protection device rather than within the structural framing system). Any excessive displacement response can then be limited by incorporating a damping device in parallel with the negative stiffness device. The combination of negative stiffness and passive damping provides a large degree of control over the expected performance of the structure. In this paper, a numerical study is presented on the performance of a seismically isolated highway bridge model that is subjected to various strong earthquake ground motions. The Negative Stiffness Devices (NSDs) are described along with their hysteretic behavior as obtained from a series of cyclic tests wherein the tests were conducted using a modified design of the NSDs (modified for testing within the bridge model). Using the results from the cyclic tests, numerical simulations of the seismic response of the isolated bridge model were conducted for various configurations (with/without negative stiffness devices and/or viscous dampers). The results demonstrate that the addition of negative stiffness devices reduces the base shear substantially, while the deck displacement is limited to acceptable values. This assessment was conducted as part of a NEES (Network for Earthquake Engineering Simulation) project which included shaking table tests of a quarter‐scale highway bridge model. Copyright © 2014 John Wiley & Sons, Ltd.