Hydrologic response of engineered media in living roofs and bioretention to large rainfalls: experiments and modeling

The hydrologic response of engineered media plays an important role in determining a stormwater control measure's ability to reduce runoff volume, flow rate, timing, and pollutant loads. Five engineered media, typical of living roof and bioretention stormwater control measures, were investigated in laboratory column experiments for their hydrologic responses to steady, large inflow rates. The inflow, medium water content response, and outflow were all measured. The water flow mechanism (uniform flow vs. preferential flow) was investigated by analyzing medium water content response in terms of timing, magnitude, and sequence with depth. Modeling the hydrologic process was conducted in the HYDRUS‐1D software, applying the Richards equation for uniform flow modeling, and a mobile–immobile model for preferential flow modeling. Uniform flow existed in most cases, including all initially dry living roof media with bimodal pore size distributions and one bioretention medium with unimodal pore size distribution. The Richards equation can predict the outflow hydrograph reasonably well for uniform flow conditions when medium hydraulic properties are adequately represented by appropriate functions. Preferential flow was found in two media with bimodal pore size distributions. The occurrence of preferential flow is more likely due to the interaction between the bimodal pore structure and the initial water content rather than the large inflow rate.     

Modeling nonlinear wave–wave interactions with the elliptic mild slope equation

An approach is developed to simulate wave–wave interactions using nonlinear elliptic mild-slope equation in domains where wave reflection, refraction, diffraction and breaking effects must also be considered. This involves the construction of an efficient solution procedure including effective boundary treatment, modification of the nonlinear equation to resolve convergence issues, and validation of the overall approach. For solving the second-order boundary-value problem, the Alternating Direction Implicit (ADI) scheme is employed, and the use of approximate boundary conditions is supplemented, for improved accuracy, with internal wave generation method and dissipative sponge layers. The performance of the nonlinear model is investigated for a range of practical wave conditions involving reflection, diffraction and shoaling in the presence of nonlinear wave–wave interactions. In addition, the transformation of a wave spectrum due to nonlinear shoaling and breaking, and nonlinear resonance inside a rectangular harbor are simulated. Numerical calculations are compared with the results from other relevant nonlinear models and experimental data available in literature. Results show that the approach developed here performs reasonably well, and has thus improved the applicability of this class of wave transformation models.     

Universal hydrofracturing algorithm for shear-thinning fluids: Particle velocity based simulation

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.     

天山南麓雷暴日短期预报研究

利用阿克苏及邻近地区12个气象站1980—2013年雷暴资料,以及同期高空资料,统计了各站年均雷暴日数,对发生区域雷暴天气的环流形势进行分类,归纳出各型的入型指标。通过逐步回归法,建立阿克苏及邻近地区区域雷暴概率回归预报模型,并对2013年进行试预报。结果表明:(1)阿克苏及邻近地区区域雷暴的影响系统主要分为4类:巴湖低槽型、急流型、西北气流型和温度槽型。(2)对2002—2012年5—9月(共1683 d)历史资料进行判别,满足入型条件的样本数为876 d,消空率为48%;对2013年5—9月(共153 d)历史资料进行判别,入型样本数为80 d,消空率为48%。(3)对2002—2012年5—9月所有入型样本进行回代检验,平均准确率为72.0%(平均TS评分为30.1%);对2013年5—9月所有入型样本进行试预报,平均准确率为63.2%(平均TS评分为28.2%)。     

1999年北京持续高温天气过程的诊断分析

利用北京观象台观测资料和NCEP/NCAR再分析数据,对1999年6月24日至7月2日北京一次持续性高温天气的演变和发展过程及非绝热加热作用对系统的影响进行诊断分析,结果表明：在此次高温天气发生前,欧亚大陆中高纬度环流经向度很大,欧洲北部和贝加尔湖以南为高压脊控制,中亚和我国东北地区则处于低压槽内。贝加尔湖南部的高压脊纬向延伸范围较广,在东移过程中长时间影响北京。随着贝加尔湖以南的高压脊逐渐东移,北京上空下沉增温与非绝热加热作用有所增强,北京逐渐受到高温天气影响。在高温天气发生的后半阶段,我国东北的低压槽入海后在120130°E附近维持并发展,槽前非绝热加热率很大。从垂直方向来看,加热率在500 hPa以下随高度迅速增加,根据全型涡度方程,强烈的非绝热加热率垂直分布不均作为一个明显的涡度源区,对入海低压槽的稳定维持有显著的作用。而入海低压槽的稳定维持,又阻碍了华北高压脊的东移,使其在北京地区长时间稳定少动,为北京带来多日的持续性高温天气。     

中国东部城市群发展对南海夏季风爆发影响的模拟研究

本文利用NCAR开发的CAM5.1(Community Atmosphere Model Version 5.1)模式,针对我国东部大规模城市下垫面发展对南海夏季风爆发的影响进行了数值模拟研究。结果表明我国东部大规模城市群发展可能使得南海夏季风提前1候爆发;机理分析表明:在南海夏季风爆发之前,中国东部城市群发展引起的陆面增温,使得南海及其附近地区南北温差提前逆转、中国东部区域海平面气压降低,导致中南半岛到南海地区西南气流加强,中南半岛到南海地区降水增加,而凝结潜热垂直变化强迫出的异常环流,促进了南亚高压的加强及提前北跳,相伴随的高层抽吸作用有助于季风对流的建立和西太平洋副高的减弱东撤,从而形成了有利于南海夏季风爆发的高低层环流条件,导致南海夏季风提前爆发。另外,观测结果表明1993年之后南海夏季风爆发的日期相对上一个年代明显提前约2候,城市化快速发展阶段与南海夏季风爆发的年代际变化存在时间段的吻合,表明城市下垫面发展可能是南海夏季风提前爆发的原因之一。     

利用压密方程计算分析西安市地面沉降特征

利用西安市地下水位监测资料,基于水文地质三维结构模型,在不考虑粘土层滞后压缩变形的理想情况下,采用压密方程计算获得西安市抽取地下水可引起的理论地面沉降量及沉降分布特征,结果表明:计算所得沉降显著区位于西安市西南部的高新区及东南部的曲江新区,沉降量呈由北向南递增的特征,与InSAR监测结果整体趋势具有较好的一致性;地面沉降分布特征受到活动地裂缝影响,沉降曲线呈近NE向偏转展布。研究结果可为沉降灾害预防研究提供宏观的参考信息。     

一种裸露土壤湿度反演方法

针对目前土壤湿度反演方法研究较少且缺少实时性的现状,该文提出一种土壤湿度反演方法——最小二乘支持向量机技术。以积分方程模型为正向算法,数值模拟不同雷达参数(频率、入射角及极化)下后向散射系数随土壤含水量和地表粗糙度的变化情况。经过数据敏感性分析,选取C-波段和X-波段、小入射角下的同极化后向散射系数作为支持向量回归的训练样本信息;经过适当的训练,利用支持向量回归技术对土壤含水量进行了反演研究;并考虑通过多频率、多极化、多入射角数据的组合,消除地表粗糙度的影响,提高反演精度。模拟结果表明,该方法反演土壤湿度具有较高的精度和较好的实时性;同时,与人工神经网络方法的结果比较,证明了该方法的有效性,为土壤湿度的反演研究提供了一种方法。     

垂荡式波浪能装置运动模型数值分析

垂荡式波浪能装置在海洋可再生能源开发中被广泛应用,通过浮子与摇杆在垂荡方向的相对运动吸收波浪能。在以往相关的运动预报数值分析中,通常基于微幅波假设,仅考虑浮子与摇杆在垂荡方向的运动,忽略摇杆其他自由度运动。建立并求解了垂荡式波浪能装置的非线性联合运动方程组,分析垂荡式波浪能装置的波浪载荷、浮子与平台连接处的受力情况。数值计算系统的运动响应,并将计算结果与已有的试验数据进行比较验证,结果表明数值模拟的垂荡式波能装置的运动响应与试验结果相符合。最后,应用本计算方法分析PTO(power take-off)参数对波能装置发电性能的影响。     

TTI介质拟声波方程数值模拟

TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。