Stream tube integration method for analysing subsurface fluid flow

A stream tube integration method is introduced to solve transient subsurface fluid flow problems. The method combines a geometry-embedded form of Darcy's Law and the notion of location of average. Two types of problems, transient radial flow to a well of finite radius in an areally infinite aquifer and in a double porosity system, are solved by the stream tube integration method and the integral finite difference method. Results of the solutions show that the stream tube integration method, with fixed coarse mesh, are more accurate and better behaved than the integral finite difference method, with fine mesh. The fixed mesh stream tube integration method is readily extended to the moving mesh method. With much coarse mesh, the moving mesh technique can obtain the same accurate results as the fixed mesh stream tube integration method. It is suggested that the stream tube integration method is a viable way to state, solve, interpret and verify numerical solutions. The method provides efficient computation and improved accuracy for analysing subsurface fluid flow. © 1998 John Wiley & Sons, Ltd.     

Stochastic and analytical approaches for sediment accumulation in river reservoirs

ABSTRACT

Sediment accumulation in a river reservoir is studied by stochastic time series models and analytical approach. The first-order moving average process is found the best for the suspended sediment discharge time series of the Juniata River at Newport, Pennsylvania, USA. Synthetic suspended sediment discharges are first generated with the chosen model after which analytical expressions are derived for the expected value and variance of sediment accumulation in the reservoir. The expected value and variance of the volume of sediment accumulation in the reservoir are calculated from a thousand synthetic time series each 38 years long and compared to the analytical approach. Stochastic and analytical approaches perfectly trace the observation in terms of the expected value and variability. Therefore, it is concluded that the expected value and variance of sediment accumulation in a reservoir could be estimated by analytical expressions without the cost of synthetic data generation mechanisms.     

一次强风暴三维结构的观测分析

2003年6月28日山东德州发生了一次降雹持续时间长、危害严重的强降雹天气过程。利用MICAPS、探空、NCEP资料,分析了降雹前的天气形势、能量场特征、含水量场。根据两部多普勒天气雷达资料,分析了雹云的演变过程及流场结构。研究表明,降雹前在降雹区上空西边低层存在高能舌,在中高层西边存在负能量平流或负能量平流中心。此次超级单体风暴过程由单个单体发展而来,具有强超级单体风暴的典型特征。该超级单体的移动方向在盛行风向右侧约20°左右,属于右移风暴。     

SOM拓扑保序机理的模拟弹性力分析

针对SOM内在原理和基本特性中一些未解决的重要问题,定义了输出阵列上神经元i承受的、正比于模式向量及权向量之间距离的2种模拟弹性力：输入模式向量对胜者i的拉力,以及阵列中胜者j对邻近i的牵动力;分析了弹性力的9个性质,对SOM学习过程中神经元权向量的演变过程及拓扑映射保序的机理给出了易于理解的、形象的解释。     

衬砌隧道中移动轴向激励作用下两相多孔介质动力响应

为了研究衬砌隧道中地铁振动在饱和地层中的传播情况,采用移动轴向激励模拟技术,建立了隧道-衬砌-两相介质的动力分析模型;利用波函数展开法、傅里叶变换法等,推导了频域内衬砌隧道移动轴向激励作用下两相多孔介质动力响应的解析解,并给出了两相多孔介质临界速度的经验公式;通过离散快速傅里叶逆变换得到了时-空域内两相介质的动力响应。结果表明：无衬砌隧道中移动的轴向常激励作用下,两相介质临界速度只与介质的剪切模量和密度有关,且数值接近其剪切波速的1.1倍;对于衬砌隧道,介质的临界速度随着衬砌剪切模量的增大而增大,随着衬砌密度的增大而减小,衬砌对振动向介质中的传播有一定的削弱作用,衬砌剪切模量与介质剪切模量相差越大,削弱越明显;动力响应频率越接近激振频率,其幅值越大,所对应的临界速度越小。     

Nonstationary stochastic response of structural systems equipped with nonlinear viscous dampers under seismic excitation

Nonlinear viscous dampers are supplemental devices widely used for enhancing the performance of structural systems exposed to seismic hazard. A rigorous evaluation of the effect of these damping devices on the seismic performance of a structural system should be based on a probabilistic approach and take into account the evolutionary characteristics of the earthquake input and of the corresponding system response. In this paper, an approximate analytical technique is proposed for studying the nonstationary stochastic response characteristics of hysteretic single degree of freedom systems equipped with viscous dampers subjected to a fully nonstationary random process representing the seismic input. In this regard, a stochastic averaging/linearization technique is utilized to cast the original nonlinear stochastic differential equation of motion into a simple first‐order nonlinear ordinary differential equation for the nonstationary system response variance. In comparison with standard linearization schemes, the herein proposed technique has the significant advantage that it allows to handle realistic seismic excitations with time‐varying frequency content. Further, it allows deriving a formula for determining the nonlinear system response evolutionary power spectrum. By this way, ‘moving resonance’ effects, related to both the evolutionary seismic excitation and the nonlinear system behavior, can be observed and quantified. Several applications involving various system and input properties are included. Furthermore, various response parameters of interest for the seismic performance assessment are considered as well. Comparisons with pertinent Monte Carlo simulations demonstrate the reliability of the proposed technique. Copyright © 2014 John Wiley & Sons, Ltd.     

移动荷载作用下土体动力反应数值模拟

在研究轨道振动荷载的产生机理基础上,归纳出一个能够涵盖速度、线路不平顺、轮重、轨枕对轮载的分散作用等因素在内的轨道荷载解析表达式,以此作为输入施加到土体模型上,利用ANSYS模拟地基土在轨道移动荷载作用下的三维瞬态弹塑性反应。通过计算得到土体不同位置位移和加速度时程曲线,并得出如下结论:荷载刚施加时地面位移会有一个突变,但是随着荷载的移动会达到一个稳定值,并随着水平面内远离轨道,该值逐渐衰减;荷载刚刚施加时加速度会有一个突变值,随着荷载的移动逐渐趋于零。轨道机车速度对竖向方向的位移和加速度影响不大。     

移动荷载作用下非均匀地基的动力响应分析

基于线弹性动力学理论,结合坐标变换,建立了移动荷载作用下非均匀弹性半平面地基的动力控制方程,利用半解析法研究了移动荷载作用下二维非均匀地基的动力响应问题。采用傅里叶（Fourier）级数展开,假设了响应函数的级数形式,通过理论推导获得了剪切模量随深度任意变化的非均匀地基在移动荷载作用下各物理量的解析表达式。考虑土体的剪切模量沿厚度方向按幂函数梯度变化,通过数值算例分析并讨论了地基非均匀参数、荷载移动速度以及地基表面的剪切模量等对地基力学响应的影响规律,并与均质地基的计算结果进行了比较。数值结果表明：地基中各点的竖向位移随着土体表面剪切模量和表征土体非均匀性的梯度因子的增大而减小,随着荷载移动速度的增大而增大。在移动荷载作用下,非均匀地基与均匀地基的动力响应有着显著的区别。     

Sensitivity of interfacial hydraulics to the microtopographic roughness of water‐lain gravels

Flow within the interfacial layer of gravel‐bed rivers is poorly understood, but this zone is important because the hydraulics here transport sediment, generate flow structures and interact with benthic organisms. We hypothesized that different gravel‐bed microtopographies generate measurable differences in hydraulic characteristics within the interfacial layer. This was tested using a high density of spatially and vertically distributed, velocity time series measured in the interfacial layers above three surfaces of contrasting microtopography. These surfaces had natural water‐worked textures, captured in the field using a casting procedure. Analysis was repeated for three discharges, with Reynolds numbers between 165000 and 287000, to evaluate whether discharge affected the impact of microtopography on interfacial flows. Relative submergence varied over a small range (3.5 to 8.1) characteristic of upland gravel‐bed rivers. Between‐surface differences in the median and variance of several time‐averaged and turbulent flow parameters were tested using non‐parametric statistics. Across all discharges, microtopographic differences did not affect spatially averaged (median) values of streamwise velocity, but were associated with significant differences in its spatial variance, and did affect spatially averaged (median) turbulent kinetic energy. Sweep and ejection events dominated the interfacial region above all surfaces at all flows, but there was a microtopographic effect, with Q2 and Q4 events less dominant and structures less persistent above the surface with the widest relief distribution, especially at the highest Reynolds number flow. Results are broadly consistent with earlier work, although this analysis is unique because of the focus on interfacial hydraulics, spatially averaged ‘patch scale’ metrics and a statistical approach to data analysis. An important implication is that observable differences in microtopography do not necessarily produce differences in interfacial hydraulics. An important observation is that appropriate roughness parameterizations for gravel‐bed rivers remain elusive, partly because the relative contributions to flow resistance of different aspects of bed microtopography are poorly constrained. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

改进的DisTrad模型在地形起伏区地表温度空间分辨率提升的应用

DisTrad（Disaggregation procedure for radiometric surface temperature）模型是常用于遥感地表温度空间分辨率提升的主要模型之一。DisTrad模型常面向空间范围有限、地形相对平坦的研究区域,且常选用植被参数（如植被指数或植被覆盖度等）作为关键参数。然而在空间范围较大、地形起伏地区,地表温度的空间变异可能无法完全通过植被参数解释。本研究选取四川盆地及毗邻地区为研究区,通过模拟数据研究DisTrad模型在地形起伏区地表温度空间分辨率提升中的适用性。数字高程模型（Digital Elevation Model,DEM）、归一化差值植被指数（Normalized Difference Vegetation Index,NDVI）等参数,采用滑动窗口逐步回归,将空间分辨率为6km的地表温度提升至空间分辨率为1km。研究结果表明,改进的模型在平原及海拔较低的高原地区提升获得的地表温度空间分辨率具有较高精度,均方根误差（Root Mean Square Error,RMSE）为0.5K左右;在地形起伏较大的地区,RMSE为4K,验证了改进的模型提升的可行性。