Adaptive element free Galerkin method applied to analysis of earthquake induced liquefaction

An automatically adaptive element free method is presented to analyze the seismic response of liquefiable soils.The method is based on the element free Galerkin method (EFGM) and the fission procedure that is part of h-refinement,indicated by error estimation. In the proposed method, a posteriori error estimate procedure that depends on the energy normof stress and the T-Belytschko (TB) stress recovery scheme is incorporated. The effective cyclic elasto-plastic constitutivemodel is used to describe the nonlinear behavior of the saturated soil. The governing equations are established by u-pformulation. The proposed method can effectively avoid the volumetric locking due to large deformation that usually occursin numerical computations using the finite element method (FEM). The efficiency of the proposed method is demonstratedby evaluating the seismic response of an embankment and comparing it to results obtained through FEM. It is shown that theproposed method provides an accurate seismic analysis of saturated soil that includes the effects of liquefaction     

济南边界层风廓线雷达与L波段雷达大风探空测风对比

基于相同站点的济南边界层风廓线雷达(WPR)和L波段雷达大风探空测风(以下简称LW)资料,采用相关、拟合、廓线分析等方法,系统地对比两者大风(≥10.8 m·s~(-1))数据的各统计特征、相关性、时空变化规律的异同。结果表明:(1)两者u分量、v分量、风向、风速的相关系数分别为0.973、0.965、0.994、0.665,标准偏差分别为2.04 m·s~(-1)、2.88 m·s~(-1)、10.82°、2.53 m·s~(-1),大风数据相关性总体较高,且相关系数风向优于风速的、降水样本优于非降水的,表明WPR观测大风数据可信、降水期间WPR水平大风数据可以使用。(2)两者大风风向、风速基本一致,但在低层差异较高层的大,尤其是340 m及以下高度层更显著,降水、非降水、全部样本在2980、1900、2740 m以下WPR风向值均略小于LW,而在1300、2740、1660 m以下WPR风速值均略大于LW。(3)两者资料互补性总体较好,大风风向、风速、v分量差值随高度变化符合对数律递增或递减规律,但u分量需分段拟合;在进行回归方程求算和资料互补时,应考虑差异随高度等的变化。     

Performance Evaluation of the GRAPES Model in Wind Simulations Over South China

In the present study, the performance of the GRAPES model in wind simulation over south China was assessed. The simulations were evaluated by using surface observations and two sounding stations in south China. The results show that the GRAPES model could provide a reliable simulation of the distribution and diurnal variation of the wind. It showed a generally overestimated southerly wind speed especially over the Pearl River Delta region and the south of Jiangxi Province as well as the coastal region over south China. GRAPES also exhibited a large number of stations with the opposite surface wind directions over the east of Guangxi and the south of Jiangxi during the nocturnal- to-morning period, as well as an overall overestimation of surface wind over the coastal regions during the afternoon. Although GRAPES could simulate the general evolutional characteristics of vertical wind profile, it underestimated wind speed above 900 hPa and overestimated wind speed below 900 hPa. Though the parameterization scheme of gravity wave drag proved to be an effective method to alleviate the systematic deviation of wind simulation, GRAPES still exhibited large errors in wind simulation, especially in the lower and upper troposphere.     

建筑物对大气污染物扩散影响的大涡模拟

利用一个大涡模式对一个方形建筑物周围的气流场进行了模拟,并与相应的风洞实验结果进行比较,对比结果表明大涡模拟方法可以精细地反映建筑物周围的流场特征.在此基础上,将拉格朗日随机游动扩散模式与大涡模式相结合,对在受建筑物影响的气流场中的大气污染物扩散进行模拟,模拟结果表明该方法可以很好地模拟出在建筑物影响下的气流变形所引起的各种污染物散布情况.建筑物周围的气流结构特性使得建筑物顶部污染源位置的细小变化可能造成建筑物周围污染物分布形势的很大不同,特别是对建筑物背风侧的空腔区内地面污染物浓度有着很大的影响.当排放     

TURBULENCE CHARACTERISTICS OF THE ATMOSPHERIC SURFACE LAYER OVER AN OPEN GRASSLAND*

Described in this paper is an experiment on atmosphere-surface turbulent exchange and boundary layer turbulence properties conducted in July 1994 over the Kerqin Grassland,Jilin,China.The characteristics of the turbulent spectrum,and the relationships of the standard deviation of the turbulent velocity components and sensible heat flux with the atmospheric stability are studied using data from a sonic anemometer and a fast-response platinum resistance thermometer mounted on a 100 m tower.The results show that in the surface layer over a flat,uniform and open grassland,for a broad stability range(-22.12≤Z_iL≤17.98),the velocity spectra obey the-2/3 power law in the inertial subrange,and 1 power law at low frequencies.Under near neutral stratification,σ_u/u_*=1.20,σ_v/u_*=1.23 and σ_w/u_*=1.02.For Z_i/L≤-0.2,the standard deviations of the turbulent velocity components follow a 1/3 power law.For Z_i/L≤-0.1,the standard deviation of the temperature fluctuations follows a-1/3 power law,and as Z_i/L≤-0.08,it exhibits nonlinear behavior.Sensible heat flux is well correlated to the stability parameter.     

近50年我国持续性暴雨的统计分析及其大尺度环流背景

根据1951～2005年中国大陆730个台站的逐日降水资料和持续性暴雨发生的特点, 提出了采用局地持续性暴雨和区域持续性暴雨事件的两种客观定义来描述和统计近50年我国持续性暴雨时空特征及其变化。单站逐日降水量连续三天或三天以上均大于等于50 mm为一次局地持续性暴雨；区域持续性暴雨指在一定区域范围内连续三天降水量总和大于等于100 mm且每天降水量大于等于25 mm的面积超过某一阈值 (详细定义见正文）。根据局地持续性暴雨定义, 指出: 近50年中国局地持续性暴雨事件主要发生在江南和华南地区, 发生季节以6月为最多；根据区域持续性暴雨定义, 统计分析了四类典型的区域持续性暴雨类型, 分别是: 渤海辽西型、北方经向型、南方锋面型和华南低压型。其中南方锋面型又可根据持续性暴雨易发生的地理位置分为江淮型、江南型和华南型三种。对这六类区域持续性暴雨的历史个例进行同比分析, 研究了不同类型持续性暴雨发生的季节性和年际变化,以及它们在大尺度环流背景方面的共性特征。     

风预报的数值产品释用研究

根据地转风和偏差风原理，建立了风预报的数值产品释用方案。通过一年试验预报和３年业务预报检验，证明此方案是一个有效的实用的方案。     

北半球500 hPa高度场定常波不平稳性分析

提出并阐明Lorenz环流分解意义下的定常波不平稳性概念,它是月平均网上纬向波动分量气候变率与定常波强度相对大小的表征.根据Lorenz环流分解,定义全域(局域)定常波不平稳度Ius(I1us),分析了北半球500 hPa位势高度场定常波强度较大的30°-60°N纬带的定常波不平稳性,结果表明:(1)全城定常波不平稳带位置存在季节件北进、南退过程.平稳的定常波出现在冬季的35°-55°N的中纬度带和夏季的副热带地区(35°N以南),分别与冬季的东亚大槽、北美槽和较弱的欧洲槽,以及夏季的副热带高压等系统相联系.不平稳度的高值中心出现在春季的35°N和夏季的50°N,这与定常波强度季节变化和月平均图上槽脊位置、强度年际异常有关.(2)局域定常波不平稳度存在着明显的纬向不对称性.平稳带通常位于定常波的强槽强脊所控制的区域,而不平稳带通常位于定常波强度较弱的区域.副热带(35°N及以南)局域定常波不平稳度冬强于夏,中纬度(35°N及以北)则夏强于冬.夏季局域定常波不平稳度地理分布具有复杂的结构.但无论冬夏,北欧是定常波最不平稳的地区,北美大陆附近的定常波则相对平稳.(3)夏季,从华北经东北至北太平洋存在一个定常波不平稳度高值带,其高值中心位于中国黑龙江省东部(45°N,130°E),主要影响中国北方(东北、华北、西北),可能是该区夏季气候脆弱带的环流成因.     

大涡模拟模式研究对流层湍流

利用三维大涡模拟模式研究不同地面热力状况和风切变对湍流产生的影响，在对流层湍流发展过程中，这两种作用是同时存在的、同时起作用，水平风速的瞬时垂直切变可以在自由大气中激发湍流，但它不能维持很长时间，对流热泡也可以在自由大气中激发湍流，其湍流强度与地面的热量强度有关。     

华南前汛期大范围暴雨的合成分析

利用多时间尺度 NCEP/ NCAR再分析资料、TBB资料 ,通过合成分析、动力诊断等方法 ,对华南大范围暴雨高低空共同特征进行了探讨 ,得出 :华南暴雨区位于 2 0 0 h Pa高空急流入口的右侧 ,高空辐散中心区下方 ,同时位于 2 0 0 h Pa南亚高压北缘和中纬度脊前辐散气流中。对于低层而言 ,暴雨区位于 85 0 h Pa低空急流轴的左侧。