Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Screening of polymers on selected Hawaii soils for erosion reduction and particle settling

In recent years, high‐molecular‐weight anionic polyacrylamides (PAMs) have been tested on a variety of soils, primarily in temperate climates. However, little information is available regarding the effectiveness of PAM for preventing soil loss through runoff in tropical settings. Screening tests were performed using three negatively charged PAMs and one positively charged PAM on five Hawaii soils (two Oxisols, one Vertisol, and two Aridisols) to determine erosion loss, sediment settling, and aggregate stability. A laboratory‐scale rainfall simulator was used to apply erosive rainfall at intensities from 5 to 8·5 cm h^?1 at various PAM doses applied in both dry and solution forms. Soil detachment due to splash and runoff, as well as the runoff and percolate water volumes, were measured for initial and successive storms. The impact of PAM on particle settling and aggregate stability was also evaluated for selected soil‐treatment combinations. Among the PAMs, Superfloc A‐836 was most effective, and significantly reduced runoff and splash sediment loss for the Wahiawa Oxisol and Pakini Andisol at rates varying between 10 and 50 kg ha^?1. Reduced runoff and splash sediment loss were also noted for PAM Aerotil‐D when applied in solution form to the Wahiawa Oxisol. Significant reductions in soil loss were not noted for either the Lualualei Vertisol or the Holomua Oxisol. It is believed that the high montmorillonite content of the Lualualei Vertisol and the low cation‐exchange capacity of the Holomua Oxisol diminished the effectiveness of the various PAMs tested. The polymers were also found to enhance sediment settling of all soils and helped improve their aggregate stability. This screening study shows the potential use of PAM for tropical soils for applications such as infiltration enhancement, runoff reduction, and enhanced sedimentation of detention ponds. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

冬小麦遥感冠层温度监测土壤含水量的试验研究

在冬小麦主要生育期(2002年4月初到5月底),对不灌溉的冬小麦测定了冠层温度、地温、气温以及土壤含水量,计算了冠气温差且分析了冠层温度和冠气温差与不同土层厚度的土壤含水量相关关系。结果表明:14:00的冠层温度能较好地反映20cm土层的土壤含水量变化,但与其它各土层相关性有较大的波动性;14:00的冠气温差能较好地反映40cm以上土层的土壤含水量变化,二者的相关性很高,在20cm、40cm土层,两者相关系数R2分别为0.98866、0.99389,这为用区域遥感数据反演主要生育期冬小麦的冠气温差进而监测区域40cm土壤含水量提供了实验性的依据;拔节期和灌浆期,用14:00冠气温差来拟合各土壤层的土壤含水量有较高的精度,从而为用区域遥感数据监测区域土壤含水量提供了经验性的模型。     

中国兴凯湖北岸平原晚全新世花粉记录及泥炭沼泽形成

通过兴凯湖北岸平原泥炭剖面高分辩率花粉分析研究，对晚全新世花粉划分4个组合带。XKH-4组合时期(1857～1746aB．P．)为沼泽发育前期，这一时期花粉浓度小，陆生草本植物占优势，气候干冷。XKH-3组合时期(1746～1287aB．P．)为沼泽发育早期，这一时期花粉浓度较大，且水生植物花粉含量为剖面最高，喜温落叶阔叶植物大发展时期，气候温和湿润。XKH-2组合时期(1287～602aB．P．)为沼泽发育中期，这一时期花粉浓度最小，但以陆生草本植物为主，木本植物为辅，水生植物急剧减少，针叶植物出现两次高峰，气候向冷干方向发展。XKH-1组合时期(602aB．P．至今)为沼泽发育盛期，这一时期花粉浓度最大，陆生草本植物大发展时期，气候波动较大。     

Post-depositional modification of NO-3 in snow layers at EastAntarctica and at the headwater of ürumqi River

Acidic species, such as Nitrate, in polar snow and firn layers are “reversibly” deposited, and are sufficiently volatile to undergo significant postdepositional exchange between snow/firn and the atmosphere. Through comparison of the snowpit and snowpack nitrate concentrations from central East Antarctica and the headwater of ürumqi River, we conclude that the nitrate peaks in the uppermost surface snow layers in central Antarctica are not related to an atmospheric signal and must account for post-depositional effects. Such effects, however, are not found in the surface snowpack nitrate profiles from the headwater of ürumqi River. Two reasons may account for the post-depositional difference. At first, nitrate in the polar snow and firn layers appears to be hydrated ion, which can be taken up by the atmosphere, while at the headwater of ürumqi River it seems mainly as mineral ion, which assembles the behavior of aerosol-derived species that are “irreversibly” deposited and do not undergo significant post-depositional exchange with the atmosphere. Secondly, the chemical features of the snow and ice on the Antarctica are mainly determined by wet deposition, to the contrary, dry deposition is more significant at the headwater of lUrumqi River than that on the East Antarctic Plateau.     

冬春季切变类冰雹发生条件的对比分析

本文以 1996年 12月 31日和 1981年 5月 1日为例 ,对冬、春季节发生在江苏的较大范围的切变类冰雹天气过程作了对比分析。结果指出 ,无论冬季或春季当高原东部有深槽东移 ,冷暖空气在江淮地区交汇 ,地面抬升系统为暖切 ,并有大气层结不稳定 (Δθse( 50 0 - 850 ) <0℃ =中心和较强的风向和风速垂直切变、85 0hPa西南急流轴、85 0hPa最大水汽通量轴线、5 0 0和 85 0hPa正涡度中心等相配置时 ,就可能导致江苏地区较大范围强对流天气的发生。     

Geochemical correlations between effusive and explosive silicic volcanics in the Saraykent region (Yozgat), central Anatolia,Turkey

Two main volcanic events are distinguished between Saraykent and Akçakışla in the Yozgat province of central Anatolia: (1) early Late Cretaceous–Palaeocene effusive activity, that produced a sequence of intermediate to felsic ‘basal lavas’; and (2) marginally later Palaeocene explosive activity that formed a series of covering ignimbrite flows. Due to their close temporal and spatial relation, geochemical comparisons were made between the silicic members of the lavas and ignimbrites, to identify chemical groups and their relative petrogenesis. The basal lavas range from calc‐alkaline basaltic andesites to dominant rhyolites. Based on trace element correlations three main geochemical groups were identified: the Akçakışla rhyolites (present as domes); Akçakışla rhyodacites‐dacites (lava flows); and Ozan‐Saraykent rhyolites (lava flows). Large‐ion lithophile elements have been mobile in all the groups, but mainly in the Akçakışla rhyolites. Rare earth element (REE) patterns show marked similarity between the Ozan and Saraykent basal lavas. The Akçakışla dome rhyolites are more fractionated with lower La_N/Yb_N ratios (c.10), whereas the Akçakışla basal lavas have much higher La_N/Yb_N ratios (c.30). The chemical coherence and petrographic similarities between the Saraykent and Ozan lavas suggest a single suite related via fractionation. Three geochemical groups were also established for the ignimbrites: Saraykent ignimbrite; Bağlıca ignimbrite‐Toklu‐Kızıldağ crystal tuffs; and Keklikpınar ignimbrite. The ignimbrites, like the basal lavas, display a pronounced depletion in Ba on ORG‐normalized plots. Relative to the basal lavas, chondrite‐normalized patterns for the ignimbrites are different in displaying negative Eu anomalies that indicate feldspar fractionation. The lack of geochemical overlap or coherence between any of the lava and ignimbrite groups suggests that they represent distinct eruptive events and are not related in any simple volcanic development and cogenetic sense. Two geochemical features are common to all the volcanic rock groups: (1) the presence of a Nb‐Ta anomaly, which is generally accepted as a crustal signature; and (2) the relatively low Y abundances which appear characteristic for the region as a whole. These fundamental features of the local silicic volcanism largely reflect source composition and effects. Copyright © 2002 John Wiley & Sons, Ltd.     

Finite element analysis of the effect of corrosion on the behavior of reinforced earth walls

A simple model for the corrosion‐induced loss of stiffness and strength of the steel strips of earth‐reinforced walls was introduced in a finite element simulation of the long‐term behavior of the wall, in which the backfill‐strips interactions are taken into account by means of a generalized homogenization procedure (called a multiphase model). The results show an initial phase of slow displacements induced by the loss of stiffness, followed after a few decades by a steep acceleration of the displacements, leading to wall failure. The influences of the parameter controlling corrosion, the backfill cohesion and the heterogeneity of the corrosion process are discussed. Results are used to discuss a strategy for reinforced earth wall surveillance. Copyright © 2011 John Wiley & Sons, Ltd.