Fast iterative solution of large undrained soil‐structure interaction problems     

Kok‐Kwang Phoon  Swee‐Huat Chan  Kim‐Chuan Toh  Fook‐Hou Lee 《国际地质力学数值与分析法杂志》2003,27(3):159-181

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.  相似文献   

Analytical solutions of one‐dimensional infiltration before and after ponding     

Jiann‐Mou Chen  Yih‐Chi Tan  Chu‐Hui Chen 《水文研究》2003,17(4):815-822

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.  相似文献   

Sea‐surface temperature and salinity changes in the northwest Pacific since the Last Glacial Maximum     

Tadamichi Oba  Masafumi Murayama 《第四纪科学杂志》2004,19(4):335-346

The oxygen isotope records of both benthic and planktonic Foraminifera in five piston cores, collected from the region between the Oyashio and Kuroshio Currents near Japan, clearly show the marked latitudinal shifts of these two currents during the past 25 kyr. Under the present hydrographic condition, a clear relationship between the sea‐surface temperature (SST) and oxygen isotope differences from benthic to planktonic Foraminifera is observed in this region. Using this relationship, we find decreased SSTs of 12–13°C (maximum 15°C) in the southernmost core site at the Last Glacial Maximum (LGM), indicating the Oyashio Current shifted southward. The SSTs at the southern two core sites abruptly increased more than 10°C at 10–11 ka, suggesting the Kuroshio Current shifted northward over these sites at 10–11 ka. In contrast, the northern two core sites have remained under the influence of the cold Oyashio Current for the past 25 kyr. With the reasonable estimate of bottom‐water temperature decrease of 2.5°C at the LGM, the SSTs estimated by this new method give exactly the same SST values calculated from Mg/Ca ratio of planktonic Foraminifera, allowing palaeosea‐surface salinities to be reconstructed. The result suggests that the ice volume effect was 1.0 ± 0.1‰ at the LGM. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

REVIEWS OF GEOGRAPHIC SOFTWARE     

《The Professional geographer》1987,39(1):85-94

Reviews of geographic software in this article: DEMO-GRAPHICS: WORLD POPULATIONS AND PROJECTIONS. ESP GAUSS. CEMODEL S. Damus LIMDEP. William H. Greene MICROSTAT 4.1 OTIS PCIPS. (Personal Computer Image Processing System) . H.J. Meyers and R. Bernstein. REGRESSION ANALYSIS OF TIME SERIES (RATS) SPSS/PC+ URBAN DATA MANAGEMENT SOFTWARE (UDMS)  相似文献   

The importance of oxygen isotope provenance in relation to solute content of bulk meltwaters at Imersuaq Glacier,West Greenland     

Jacob C. Yde  N. Tvis Knudsen 《水文研究》2004,18(1):125-139

Stable oxygen isotope analysis and measurement of several dissolved cations and anions of bulk meltwater samples have provided information about the hydrochemical environment of the glacial hydrological system at Imersuaq Glacier, an outlet tongue from the Greenland ice‐sheet, West Greenland. The samples were collected at frequent intervals during the period 20–28 July 2000 in a small (<20 L s^?1) englacial meltwater outlet at the glacier margin. The results document the following findings: (i) a marked diurnal variation of δ¹⁸O is related to the composition of oxygen isotope provenances, mainly near‐marginal local superimposed ice and basal up‐sheared ice further up‐glacier; (ii) a relationship is seen between all base cations (Na⁺, K⁺, Ca²⁺, Mg²⁺), SO₄^2? and δ¹⁸O, indicating that solute acquisition is provided by solid–solution contact with the up‐sheared ice—as the relationship with Cl^? is weak the influence of seasalt‐derived solutes is small in the area; (iii) when the melt rate is high, two diurnal maxima of δ¹⁸O values and solute concentrations are measured, and it is suggested that a snow meltwater component is responsible for the second maximum of δ¹⁸O—a short residence time leads to a delayed decrease in ion concentrations. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

The Preboreal climate reversal and a subsequent solar‐forced climate shift     

J. van der Plicht  B. van Geel  S. J. P. Bohncke  J. A. A. Bos  M. Blaauw  A. O. M. Speranza  R. Muscheler  S. Bjrck 《第四纪科学杂志》2004,19(3):263-269

Accurate chronologies are essential for linking palaeoclimate archives. Carbon‐14 wiggle‐match dating was used to produce an accurate chronology for part of an early Holocene peat sequence from the Borchert (The Netherlands). Following the Younger Dryas–Preboreal transition, two climatic shifts could be inferred. Around 11 400 cal. yr BP the expansion of birch (Betula) forest was interrupted by a dry continental phase with dominantly open grassland vegetation, coeval with the PBO (Preboreal Oscillation), as observed in the GRIP ice core. At 11 250 cal. yr BP a sudden shift to a humid climate occurred. This second change appears to be contemporaneous with: (i) a sharp increase of atmospheric ¹⁴C; (ii) a temporary decline of atmospheric CO₂; and (iii) an increase in the GRIP ¹⁰Be flux. The close correspondence with excursions of cosmogenic nuclides points to a decline in solar activity, which may have forced the changes in climate and vegetation at around 11 250 cal. yr BP. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

The influence of macrophyte growth,typical of eutrophic conditions,on river flow velocities and turbulence production     

Pamela Naden  Ponnambalam Rameshwaran  Owen Mountford  Coralie Robertson 《水文研究》2006,20(18):3915-3938

The influence of emergent and submerged macrophytes on flow velocity and turbulence production is demonstrated in a 140 m reach of the River Blackwater in Farnborough, Hampshire, UK. Macrophyte growth occurs in patches and is dominated by Sparganium erectum and Sparganium emersum. In May 2001, patches of S. erectum were already established and occupied 18% of the channel area. The flow adjusted to these (predominantly lateral) patches by being channelled through a narrower cross‐section. The measured velocity profiles showed a logarithmic form, with deviations attributable to topographic control. The channel bed was the main source of turbulence. In September 2001, in‐stream macrophytes occupied 27% of the channel, and overhanging bank vegetation affected 32% of the area. Overall flow resistance, described by Manning's n, showed a threefold increase that could be attributed to the growth of S. emersum in the middle of the channel. Velocity profiles showed different characteristic forms depending on their position relative to plant stems and leaves. The overall velocity field had a three‐dimensional structure. Turbulence intensities were generally higher and turbulence profiles tended to mirror the velocity profiles. Evidence for the generation of coherent eddies was provided by ratios of the root mean square velocities. Spectral analysis identified deviations from the Kolmogorov ?5/3 power law and provided statistical evidence for a spectral short‐cut, indicative of additional turbulence production. This was most marked for the submerged vegetation and, in some instances, the overhanging bank vegetation. The long strap‐like leaves of S. emersum being aligned approximately parallel to the flow and the highly variable velocity field created by the patch arrangement of macrophytes suggest that the dominant mechanism for turbulence production is vortex shedding along shear zones. Wake production around individual stems of S. emersum close to the bed may also be important locally. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)     

MARC J.P. GOUW 《Sedimentology》2008,55(5):1487-1516

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.  相似文献   

The Lower Permian Wasp Head Formation,Sydney Basin: high‐latitude,shallow marine sedimentation following the late Asselian to early Sakmarian glacial event in eastern Australia     

MICHAEL C. RYGEL  CHRISTOPHER R. FIELDING  KERRIE L. BANN  TRACY D. FRANK  LAUREN BIRGENHEIER  STUART C. TYE 《Sedimentology》2008,55(5):1517-1540

The Lower Permian Wasp Head Formation (early to middle Sakmarian) is a ～95 m thick unit that was deposited during the transition to a non‐glacial period following the late Asselian to early Sakmarian glacial event in eastern Australia. This shallow marine, sandstone‐dominated unit can be subdivided into six facies associations. (i) The marine sediment gravity flow facies association consists of breccias and conglomerates deposited in upper shoreface water depths. (ii) Upper shoreface deposits consist of cross‐stratified, conglomeratic sandstones with an impoverished expression of the Skolithos Ichnofacies. (iii) Middle shoreface deposits consist of hummocky cross‐stratified sandstones with a trace fossil assemblage that represents the Skolithos Ichnofacies. (iv) Lower shoreface deposits are similar to middle shoreface deposits, but contain more pervasive bioturbation and a distal expression of the Skolithos Ichnofacies to a proximal expression of the Cruziana Ichnofacies. (v) Delta‐influenced, lower shoreface‐offshore transition deposits are distinguished by sparsely bioturbated carbonaceous mudstone drapes within a variety of shoreface and offshore deposits. Trace fossil assemblages represent distal expressions of the Skolithos Ichnofacies to stressed, proximal expressions of the Cruziana Ichnofacies. Impoverished trace fossil assemblages record variable and episodic environmental stresses possibly caused by fluctuations in sedimentation rates, substrate consistencies, salinity, oxygen levels, turbidity and other physio‐chemical stresses characteristic of deltaic conditions. (vi) The offshore transition‐offshore facies association consists of mudstone and admixed sandstone and mudstone with pervasive bioturbation and an archetypal to distal expression of the Cruziana Ichnofacies. The lowermost ～50 m of the formation consists of a single deepening upward cycle formed as the basin transitioned from glacioisostatic rebound following the Asselian to early Sakmarian glacial to a regime dominated by regional extensional subsidence without significant glacial influence. The upper ～45 m of the formation can be subdivided into three shallowing upward cycles (parasequences) that formed in the aftermath of rapid, possibly glacioeustatic, rises in relative sea‐level or due to autocyclic progradation patterns. The shift to a parasequence‐dominated architecture and progressive decrease in ice‐rafted debris upwards through the succession records the release from glacioisostatic rebound and amelioration of climate that accompanied the transition to broadly non‐glacial conditions.  相似文献   

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume     

G. POSTMA  M. G. KLEINHANS  P. TH. MEIJER  J. T. EGGENHUISEN 《Sedimentology》2008,55(6):1541-1557

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.  相似文献