Performance of Jacobi preconditioning in Krylov subspace solution of finite element equations

This paper examines the performance of the Jacobi preconditioner when used with two Krylov subspace iterative methods. The number of iterations needed for convergence was shown to be different for drained, undrained and consolidation problems, even for similar condition number. The differences were due to differences in the eigenvalue distribution, which cannot be completely described by the condition number alone. For drained problems involving large stiffness ratios between different material zones, ill‐conditioning is caused by these large stiffness ratios. Since Jacobi preconditioning operates on degrees‐of‐freedom, it effectively homogenizes the different spatial sub‐domains. The undrained problem, modelled as a nearly incompressible problem, is much more resistant to Jacobi preconditioning, because its ill‐conditioning arises from the large stiffness ratios between volumetric and distortional deformational modes, many of which involve the similar spatial domains or sub‐domains. The consolidation problem has two sets of degrees‐of‐freedom, namely displacement and pore pressure. Some of the eigenvalues are displacement dominated whereas others are excess pore pressure dominated. Jacobi preconditioning compresses the displacement‐dominated eigenvalues in a similar manner as the drained problem, but pore‐pressure‐dominated eigenvalues are often over‐scaled. Convergence can be accelerated if this over‐scaling is recognized and corrected for. Copyright © 2002 John Wiley & Sons, Ltd.     

Statistical significance of modal parameters of bridge systems identified from strong motion data

Modal parameters of structural systems have commonly been determined using system identification (SI) methods for damage detection and health monitoring. For determining the deterioration of the integrity of structural systems correctly, modal parameters of a healthy structure have to be obtained with adequate certainty so that these parameters can be used as reliable references for the healthy system to compare with those of the damaged system. In this study, the statistical significance of modal parameters identified using strong motion time histories recorded on two bridge structures is assessed. The confidence intervals of identified modal frequencies and damping ratios are obtained using Monte Carlo simulations and sensitivity analyses in conjunction with eigenrealization algorithm. The dependence of the statistical bounds on model parameters is examined. The effect of using different number of sensors on the statistical significance is evaluated using simulated time history data from a validated finite element model of a bridge. Copyright © 2005 John Wiley & Sons, Ltd.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.     

Hyporheic exchange flow around constructed in‐channel structures and implications for restoration design

In‐channel rock vane structures are widely used in stream restoration as a way to reduce stream channel erosion and create pool or riffle features. When these structures change hydraulic gradients they may affect ecological stream functions, such as hyporheic exchange flow (HEF) patterns. A study of constructed in‐channel structure controls on HEF was conducted in the third‐order Batavia Kill, New York using stream and hyporheic temperature amplitude analysis and computational fluid dynamics (CFD) hydraulic simulations. Temperature monitors were installed in the water column and channel bed at six locations around each of seven in‐channel restoration structures (three cross‐vanes and four J‐hooks) at baseflow in 2007. Elevation surveys of the structures were then used to simulate HEF using CFD. The results indicate a pattern of pronounced upwelling in the run section just below the structure, upwelling transitioning to downwelling within the pool, and pronounced downwelling in the glide out of the pool. This pattern is consistent with natural riffle pool sequences. The direction of HEF inferred from the temperature amplitude analysis agreed with the direction of flow simulated with CFD at 80% of the locations, and the few disagreements were expected due to model limitations. CFD simulation demonstrated that increasing stream flows result in changes in HEF spatial patterns and magnitude at each structure. This work illustrates how CFD simulations can guide design of in‐channel restoration structures for HEF function. Copyright © 2009 John Wiley & Sons, Ltd.     

三角网数字地面模型与农田平整设计

针对农田水利工程中农田纵横坡度设计与填挖土方计算问题，本文提出了三角网地面模型，推证了地面坐标转化为设计面坐标的坐标转化公式及三角形模型体积计算公式，给出了设计计算的逐步趋近算法。     

Testing solar forcing of pervasive Holocene climate cycles

The temporal and spatial extent of Holocene climate change is an area of considerable uncertainty, with solar forcing recently proposed to be the origin of cycles identified in the North Atlantic region. To address these issues we have developed an annually resolved record of changes in Irish bog tree populations over the last 7468 years which, together with radiocarbon‐dated bog and lake‐edge populations, extend the dataset back to ～9000 yr ago. The Irish trees underpin the internationally accepted radiocarbon calibration curve, used to derive a proxy of solar activity, and allow us to test solar forcing of Holocene climate change. Tree populations and age structures provide unambiguous evidence of major shifts in Holocene surface moisture, with a dominant cyclicity of 800 yr, similar to marine cycles in the North Atlantic, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. The cycles, however, are not coherent with changes in solar activity (both being on the same absolute timescale), indicating that Holocene North Atlantic climate variability at the millennial and centennial scale is not driven by a linear response to changes in solar activity. Copyright © 2005 John Wiley & Sons, Ltd.     

Detecting subtle hydrochemical anomalies with multivariate statistics: an example from ‘homogeneous’ groundwaters in the Great Artesian Basin,Australia

The major ion composition of Great Artesian Basin groundwater in the lower Namoi River valley is relatively homogeneous in chemical composition. Traditional graphical techniques have been combined with multivariate statistical methods to determine whether subtle differences in the chemical composition of these waters can be delineated. Hierarchical cluster analysis and principal components analysis were successful in delineating minor variations within the groundwaters of the study area that were not visually identified in the graphical techniques applied. Hydrochemical interpretation allowed geochemical processes to be identified in each statistically defined water type and illustrated how these groundwaters differ from one another. Three main geochemical processes were identified in the groundwaters: ion exchange, precipitation, and mixing between waters from different sources. Both statistical methods delineated an anomalous sample suspected of being influenced by magmatic CO₂ input. The use of statistical methods to complement traditional graphical techniques for waters appearing homogeneous is emphasized for all investigations of this type. Copyright © 2006 John Wiley & Sons, Ltd.     

Relative sea‐level changes,glacial isostatic modelling and ice‐sheet reconstructions from the British Isles since the Last Glacial Maximum

While contributing <1 m equivalent eustatic sea‐level rise the British Isles ice sheet produced glacio‐isostatic rebound in northern Britain of similar magnitude to eustatic sea‐level change, or global meltwater influx, over the last 18 000 years. The resulting spatially variable relative sea‐level changes combine with observations from far‐field locations to produce a rigorous test for quantitative models of glacial isostatic adjustment, local ice‐sheet history and global meltwater influx. After a review of the attributes of relative sea‐level observations significant for constraining large‐scale models of the isostatic adjustment process we summarise long records of relative sea‐level change from the British Isles and far‐field locations. We give an overview of different global theoretical models of the isostatic adjustment process before presenting intercomparisons of observed and predicted relative sea levels at sites in the British Isles and far‐field for a range of Earth and ice model parameters in order to demonstrate model sensitivity and the resolving power available from using evidence from the British Isles. For the first time we show a good degree of fit between relative sea‐level observations and predictions that are based upon global Earth and ice model parameters, independently derived from analysis of far‐field data, with a terrain‐corrected model of the British Isles ice sheet that includes extensive glaciation of the North Sea and western continental shelf, that does not assume isostatic equilibrium at the Last Glacial Maximum and keeps to trimline constraints of ice surface elevation. We do not attempt to identify a unique solution for the model lithosphere thickness parameter or the local‐scale detail of the ice model in order to provide a fit for all sites, but argue that the next stage should be to incorporate an ice‐sheet model that is based on quantitative, glaciological model simulations. We hope that this paper will stimulate this debate and help to integrate research in glacial geomorphology, glaciology, sea‐level change, Earth rheology and quantitative modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

Variations in dissolved CO2 and CH4 in a first‐order stream and catchment: an investigation of soil–stream linkages

Spatial and seasonal variations in CO₂ and CH₄ concentrations in streamwater and adjacent soils were studied at three sites on Brocky Burn, a headwater stream draining a peatland catchment in upland Britain. Concentrations of both gases in the soil atmosphere were significantly higher in peat and riparian soils than in mineral soils. Peat and riparian soil CO₂ concentrations varied seasonally, showing a positive correlation with air and soil temperature. Streamwater CO₂ concentrations at the upper sampling site, which mostly drained deep peats, varied from 2·8 to 9·8 mg l^?1 (2·5 to 11·9 times atmospheric saturation) and decreased markedly downstream. Temperature‐related seasonal variations in peat and riparian soil CO₂ were reflected in the stream at the upper site, where 77% of biweekly variation was explained by an autoregressive model based on: (i) a negative log‐linear relationship with stream flow; (ii) a positive linear relationship with soil CO₂ concentrations in the shallow riparian wells; and (iii) a negative linear relationship with soil CO₂ concentrations in the shallow peat wells, with a significant 2‐week lag term. These relationships changed markedly downstream, with an apparent decrease in the soil–stream linkage and a switch to a positive relationship between stream flow and stream CO₂. Streamwater CH₄ concentrations also declined sharply downstream, but were much lower (<0·01 to 0·12 mg l^?1) than those of CO₂ and showed no seasonal variation, nor any relationship with soil atmospheric CH₄ concentrations. However, stream CH₄ was significantly correlated with stream flow at the upper site, which explained 57% of biweekly variations in dissolved concentrations. We conclude that stream CO₂ can be a useful integrative measure of whole catchment respiration, but only at sites where the soil–stream linkage is strong. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of soil interaction on the performance of liquid column dampers for seismic applications

The effects of soil–structure interaction (SSI) while designing the liquid column damper (LCD) for seismic vibration control of structures have been presented in this study. The formulation for the input–output relation of a flexible‐base structure with attached LCD has been presented. The superstructure has been modelled by a single‐degree‐of‐freedom (SDOF) system. The non‐linearity in the orifice damping of the LCD has been replaced by equivalent linear viscous damping by using equivalent linearization technique. The force–deformation relationships and damping characteristics of the foundation have been described by complex valued impedance functions. Through a numerical stochastic study in the frequency domain, the various aspects of SSI on the functioning of the LCD have been illustrated. A simpler approach for studying the LCD performance considering SSI, using an equivalent SDOF model for the soil–structure system available in literature by Wolf (Dynamic Soil–Structure Interaction. International Series in Civil Engineering and Engineering Mechanics. Prentice‐Hall: Englewood Cliffs, NJ, 1985) has also been presented. Copyright © 2005 John Wiley & Sons, Ltd.