Certified solutions for hydraulic structures using the node‐based smoothed point interpolation method (NS‐PIM)

A meshfree node‐based smoothed point interpolation method (NS‐PIM), which has been recently developed for solid mechanics problems, is applied to obtain certified solutions with bounds for hydraulic structure designs. In this approach, shape functions for displacements are constructed using the point interpolation method (PIM), and the shape functions possess the Kronecker delta property and permit the straightforward enforcement of essential boundary conditions. The generalized smoothed Galerkin weak form is then applied to construct discretized system equations using the node‐based smoothed strains. As a very novel and important property, the approach can obtain the upper bound solution in energy norm for hydraulic structures. A 2D gravity dam problem and a 3D arch dam problem are solved, respectively, using the NS‐PIM and the simulation results of NS‐PIM are found to be the upper bounds. Together with standard fully compatible FEM results as a lower bound, we have successfully determined the solution bounds to certify the accuracy of numerical solutions. This confirms that the NS‐PIM is very useful for producing certified solutions for the analysis of huge hydraulic structures. Copyright © 2009 John Wiley & Sons, Ltd.     

Dynamics of evapotranspiration from a riparian pond complex in the Western Boreal Forest,Alberta, Canada

The relative contributions to total actual evapotranspiration (AET) from pond and riparian areas in a pond‐wetland complex in the Western Boreal Plain (WBP) of northern Alberta are measured using the Bowen ratio energy balance technique. Measurements show that a pond typical of the WBP evaporates at a rate more than twice that of the adjacent riparian peatland. Relating the actual to potential evapotranspiration over both surfaces yields Priestley–Taylor α coefficients of 0·69 and 1·11 for the peatland and pond respectively. Further results demonstrate that the sheltering and turbulent influences of the adjacent forested areas must be considered in the processes governing the permanence of WBP ponds. That is, forestry practices may inadvertently enhance the evaporative losses from the ponds, over and above the controls exerted by the regional climate. Copyright © 2006 John Wiley & Sons, Ltd.     

Deconvolution with wavelets and vaguelettes

The use of wavelets for the solution of convolution equations is studied as a possible alternative to the well-established Fast Fourier Transform (FFT) technique. Two possible solution strategies are investigated: (1) The use of wavelets for the representation of both the given data and the unknown solution. This leads to an algorithm with good de-noising and data-compression properties. In terms of computational efficiency this algorithm is inferior to FFT. (2) The use of wavelets for the representation of the unknown solution and of so-called vaguelettes for the representations of the given data. This leads to an algorithm which is even faster than FFT. Received: 14 October 1998 / Accepted: 30 November 1999     

Discontinuous Galerkin methods for advective transport in single-continuum models of fractured media

Accurate simulation of flow and transport processes in fractured rocks requires that flow in fractures and shear zones to be coupled with flow in the porous rock matrix. To this end, we will herein consider a single-continuum approach in which both fractures and the porous rock are represented as volumetric objects, i.e., as cells in an unstructured triangular grid with a permeability and a porosity value associated with each cell. Hence, from a numerical point of view, there is no distinction between flow in the fractures and the rock matrix. This enables modelling of realistic cases with very complex structures. To compute single-phase advective transport in such a model, we propose to use a family of higher-order discontinuous Galerkin methods. Single-phase transport equations are hyperbolic and have an inherent causality in the sense that information propagates along streamlines. This causality is preserved in our discontinuous Galerkin discretization. We can therefore use a simple topological sort of the graph of discrete fluxes to reorder the degrees-of-freedom such that the discretized linear system gets a lower block-triangular form, from which the solution can be computed very efficiently using a single-pass forward block substitution. The accuracy and utility of the resulting transport solver is illustrated through several numerical experiments.     

Comparison of three models to estimate evapotranspiration for a temperate mixed forest

Forest evapotranspiration is one of the main components in the regional water budget. A comparison between measured and estimated eddy covariance (EC) data, considering the Katerji–Perrier (KP), Todorovic (TD) and Priestley–Taylor (PT) actual evapotranspiration methods, was carried out. These models, relying on more easily obtainable data, are valuable when long‐term direct measurements are not available. The objective of this paper is to compare the effectivity of these three models. In this paper, experimental data were obtained within the temperate mixed forest of broad‐leaved and coniferous trees of the Changbai Mountains in northeastern China during the growing seasons of 2003 to 2005. The KP method gave the most effective values for half‐hourly and daily evapotranspiration computed by summing up half‐hourly estimates, and the TD method overestimated evapotranspiration by about 30%. The diurnal courses of estimated and measured evapotranspiration showed bell curves, similar to that of net radiation, except for a slight increase at about 14:30 solar time due to a peak value of vapour pressure deficit (VPD). For the case of daily evapotranspiration using daily mean micrometeorological variables, the PT method presented the closest values to the measurements. Accuracy of estimation related to VPD negatively (especially for VPD > 1·5 kPa). The KP parameters, considered to be vegetation dependent, were a = 0·545 and b = 1·31 at the experimental site. A constant PT parameter (α = 1·18) was applied to estimated evapotranspiration. Daily values of α responded to VPD (negatively) more strongly than to soil moisture (positively) in this forest. The experiment showed the inherent limits and advantages of the three methods. The KP method, a semi‐empirical approach, was preferred to estimate half‐hourly evapotranspiration. The TD method was a mechanistic approach to estimate reference evapotranspiration and always overestimated actual evapotranspiration. The PT method, being site dependent and the simplest approach, was effective enough to estimate large time‐scale (at least daily) evapotranspiration. Copyright © 2008 John Wiley & Sons, Ltd.     

土层参数的随机性对场地传递函数的影响

本文采用秦勒级数展开估计土层参数随机性对场地传递函数的影响。对单层土模型在土层参数服从独立高斯分布假定下进行了公式推导和实例计算,比较了不同的泰勒级数阶次、土层参数变异系数以及阻尼比对传递系数随机性的综合效应。结果表明,即使在小参数摄动范围内采用一阶泰勒级数也是不够的;一般情况下采用二阶泰勒级数即可给出满意的估计。     

20个CMIP5模式对中亚地区年平均气温模拟能力评估

利用CMIP5中20个模式的历史模拟结果和英国东英格利亚大学CRU观测数据,采用泰勒图、趋势分析、滑动平均和EOF等方法从气候态和气候变率两方面检验各个模式对中亚年平均气温的模拟能力。结果表明：各模式能较好地模拟1951—2005年中亚地区显著增温趋势和年平均气温的空间分布特点,尤其是高、低值中心和等值线数值分布。泰勒图分析显示,大部分模式的均方根误差在0.5左右,空间相关系数在0.85～0.90之间,标准差在0.5～1.0之间。EOF分析结果表明,模式集合平均能够较好地表现出中亚气温一致升高以及南北反位相波动这两个主要模态的时空变化特征。     

���ڸ߽�̩�ռ���չ��ʽ�Ķ������ӳ����ڲ��о�

???????μ???巽??????????????????????????ZTD?????÷????????????????????????????????????????????幫?????????????????????????????????????????Ч????ZTD????徫???????SCIGN?????GPS?????????????????????????Χ?????????????????????????????????????ZTD???????????????????????????????????????????????????????????????????????????????????????????????????????????????????