Efficient DEM‐based overland flow routing using integrated recursive algorithms

Two primary concerns in performing watershed overland flow routing are the numerical instability and computational efficiency. The stability of executing an explicit scheme has to be maintained by observing the Courant–Friedrich–Lewy criterion, which is adopted to confirm that the numerical marching speed is larger than the wave celerity. Moreover, there is another criterion of time step devised in previous studies to avoid back‐and‐forth refluxing between adjacent grids. The situation of refluxing usually occurs on flat regions. In light of this, the selection of a small time increment to honor both restrictions simultaneously is believed to decrease the computational efficiency in performing overland flow routing. This study aims at creating a robust algorithm to relax both restrictions. The proposed algorithm was first implemented on a one‐dimensional overland plane to evaluate the accuracy of the numerical result by comparing it with an analytical solution. Then, the algorithm was further applied to a watershed for 2D runoff simulations. The results show that the proposed integrated algorithm can provide an accurate runoff simulation and achieve satisfactory performance in terms of computational speed.     

Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

The behavior of magneto-acoustic-gravity waves near the cusp resonance in a lossless,compressible, isothermal,stratified, electrically conducting and uniformly magnetized atmosphere. II: The valve effect

Abstract

In the present paper the behavior of (internal) magneto-acoustic-gravity waves near the cusp resonance in a lossless, compressible, isothermal, stratified, electrically conducting atmosphere that is permeated by a uniform, nearly horizontal magnetic field is re-addressed (Kamp, 1989). The previously analyzed linear conversion of long acoustic-gravity waves into short magneto-acoustic waves that carry off the energy from the resonance region along the magnetic field, is re-analyzed with boundary layer techniques that are based on the smallness of the vertical component of the magnetic field. More specifically the existence of the so-called valve effect for the generated magneto-acoustic mode near the critical level is explicitly demonstrated and shown to be governed by two rivalling effects.     

基于LLS的雷电流波头时间分布特征

利用湖北省2007—2013年雷电定位系统(Lightning Location System,LLS)监测的雷电流相关参数资料,采用数理统计方法,对雷电流波头时间的年、季、月、日变化以及累积概率和概率密度分布进行了统计分析。结果表明,正闪电平均波头时间为4.8μs,负闪和总闪电为3.6μs;一年中,正闪电7月波头时间最短,月变化特征大致呈"V"形;负闪和总闪电5—7月波头时间比其他月份短,月变化特征大致呈"U"形。正闪电波头时间季节、日变化明显,负闪和总闪电季节、日变化不明显。波头时间小于等于10μs时,正闪、负闪和总闪电的波头时间累积概率分别为95.5%、99.7%和99.6%,总闪电95%的波头时间大于1.5μs,50%的大于3.1μs,5%的大于6.0μs。波头时间在2~3μs时,正闪、负闪和总闪电的概率密度最大。气温较低月份,波头时间相对较长,反之,波头时间相对较短,其原因可能与雷电流幅值大小有关,具体原因有待进一步研究。     

强厄尔尼诺事件下2016年1月中国南方超级寒潮的动力学机制:瞬变波对大气长波异常的调制

本文采用Plumb三维波活动通量和局地Eliassen-Palm通量诊断方法,利用欧洲中期天气预报中心的逐日再分析资料ERA-Interim,分析了超强厄尔尼诺背景下2016年1月下旬中国南方超级寒潮的动力学机制：瞬变波对大气长波异常的调制。前期2015年12月的北大西洋海表热通量正异常,有助于后期大西洋阻塞形势的维持。大气长波能量沿大圆路径从大西洋阻高经乌拉尔地区向东亚中低纬度传播并在此辐合,导致了乌拉尔阻高和华北横槽的经向结构,更多强冷空气聚集在异常偏南的纬向槽线附近。寒潮爆发前夕,2支瞬变波列活跃在亚欧大陆。北支瞬变波列调制了北方的大气长波,使横槽转竖;南支瞬变波列协同调控了南方的大气长波,使南支印缅槽减弱;两者共同作用,促使冷平流大举南下,直达华南沿海,南方寒潮发生。     

天气尺度波列对长江中下游6月梅雨的影响

采用1979—2007年6月NCEP/NCAR2.5°×2.5°逐日再分析资料和中国743站逐日降水资料,利用相关分析、合成分析等方法,分析了天气尺度波列的特征及其对长江中下游6月梅雨的影响。结果表明:当长江中下游6月梅雨较少时,东亚及西太平洋区域存在一个天气尺度波列;该波列的延伸距离较短(从黄河河套地区经过长江中下游至南海、菲律宾海一带),维持时间也很短,且仅仅在500 h Pa以下较强。诊断及个例分析表明,当该波列异常显著时,长江中下游梅雨降水明显减少,而其南部区域降水则增多,说明该波列对预报长江中下游降水具有重要的指示意义。     

2007年夏季江淮强降水过程中10~30d低频变化及其与对流层上层波包活动的联系

利用NCEP/NCAR再分析和全国740站逐日降水资料,运用一点滞后相关等方法,对2007年夏季江淮流域强降水期间低频振荡的波动活动特征及其与降水低频变化的联系进行了分析。结果表明,在2007年夏季降水中,降水低频分量起着重要作用。降水的低频振荡主周期为10~30d,降水距平时间序列与10~30d低频分量具有较好的对应关系。低频扰动在对流层上层和低层都呈现波列状分布,且在降水活跃位相时,低频环流在高、低层具有斜压结构。在对流层上层,低频扰动有缓慢的东移倾向,相速度为每天2~3个经度。西风带中存在多次移动性波列向下游的传播,且在120°E以西以每天14经度的群速度向下游频散能量,表明10~30d低频波动具有明显的下游发展特征。在强降水开始5d前,低频波动与能量可起源于高纬的乌拉尔山附近,沿着西北-东南向的路径向下游传播。下游发展的低频波动为江淮流域带来了能量,为强降水的发生提供了条件。这些结果加深了人们对低频波动在江淮流域强降水过程中所起作用的认识,可为寻找江淮流域强降水过程预报线索提供科学依据。     

基于势流理论和OpenFOAM的耦合模型对多向不规则波浪的模拟

基于势流理论的数值水池可以快速计算波浪的传播及其对建筑物的作用,但是势流理论是基于波浪的无黏性假设的,而在工程中,通常需要在固体边界及波浪破碎的区域考虑黏性效应。针对基于求解Navier-Stokes方程的黏性水池计算量较大、速度较慢的缺点,采用耦合的方法模拟多向不规则波浪的传播,即在外域通过基于势流理论的数值水池产生多向不规则波浪,内域采用求解Navier-Stokes方程和流体体积法(VOF方法)对自由表面进行追踪,通过外域所提供的边界波浪,内域计算可以在较小区域进行计算,从而达到减少计算量、提高计算效率的目的。     

腾冲火山构造区马站岩浆囊地球物理特征的再探讨

云南腾冲火山盆地深部存在岩浆囊是人们关注的问题,多方面进行了大量的研究且存在不同意见。作者此前依据低电阻异常推断出位于小空山、大空山、黑空山等火山口的下部,深13~30 km范围内东西向25 km,南北向30km的低电阻率的岩浆囊,其东侧为大盈江断裂。文章利用反射地震、重力、航空磁测资料进一步探讨了岩浆囊的其他地球物理特征。由于火山喷发,不仅在火山口附近形成了局部物质亏损,而且,岩浆囊本身也亏损了岩浆物质,在重力上延图中滤去了表层的火山局部负重力异常后,呈现出大空山北侧的与低阻体相对应的负重力异常,这与岩浆囊高温熔融状态含水低密度物质有密切关系。整体的重力负异常印证了低阻岩浆囊的存在。腾冲马站岩浆囊的温度为397~651°C,平均温度为524°C,其主体中心部位不可能存在磁性体。对航磁资料的处理,消除表层火山岩影响后,仍然发现有较强磁异常,其平面投影范围与低阻体相近,推断磁异常是岩浆囊顶层,在15 km深度以上范围内,随着岩浆囊体温度逐步下降,当低于480°C以后会形成新的铁磁性矿物,因此,在囊体上层出现了范围与囊体相近的较强的磁性体。反射地震仅仅在浅部较清楚地观测到较连续的反射波界面,是浅部火山岩、含水层、以及花岗岩顶界面的反映。向深部反射界面很不连续推测是由于多次岩浆上涌,其喷发时间有先后,成分有差异,故岩浆囊物质的不均匀性,虽然由于含水、矿化、熔融体构成了低电阻的共性,可温度的不均匀,却又显示了不同期次岩浆成分有变化,由岩浆囊向地表的通道上,东部花岗岩下为高速,相对于较低温度的物质,在小空山、大空山和黑空山深部沿大盈江断裂带有东西两个大的隐伏花岗岩体之间是最新的火山喷发通道,具有相对低速的通道。     

基于FLUENT的内孤立波质量源造波方法

为研究内孤立波质量源数值造波方法,本文采用两个点源形式的质量源,分别放置于两层流体的上下层中作为内孤立波激发源。推导源项表达式,从不可压缩流体的Navier-Stokes方程出发,结合内孤立波Kd V、e Kd V理论,基于商业软件FLUENT发展了一种内孤立波质量源数值造波方法。通过数值模拟,分析了质量源造波过程中内孤立波的生成过程,并将数值模拟结果同理论及实验作对比。结果表明:基于此方法生成的内孤立波波形、波高及波致水平速度与理论及实验吻合度较好,该方法是可行的,并且耗时短、效率高。