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.     

Stochastic Modelling of the Impact of Flood Protection Measures Along the River Waal in the Netherlands

River flooding is a problem of international interest. In the past few years many countries suffered from severe floods. A large part of the Netherlands is below sea level and river levels. The Dutch flood defences along the river Rhine are designed for water levels with a probability of exceedance of 1/1250 per year. These water levels are computed with a hydrodynamic model using a deterministic bed level and a deterministic design discharge. Traditionally, the safety against flooding in the Netherlands is obtained by building and reinforcing dikes. Recently, a new policy was proposed to cope with increasing design discharges in the Rhine and Meuse rivers. This policy is known as the Room for the River (RfR) policy, in which a reduction of flood levels is achieved by measures creating space for the river, such as dike replacement, side channels and floodplain lowering. As compared with dike reinforcement, these measures may have a stronger impact on flow and sediment transport fields, probably leading to stronger morphological effects. As a result of the latter the flood conveyance capacity may decrease over time. An a priori judgement of safety against flooding on the basis of an increased conveyance capacity of the river can be quite misleading. Therefore, the determination of design water levels using a fixed-bed hydrodynamic model may not be justified and the use of a mobile-bed approach may be more appropriate. This problem is addressed in this paper, using a case study of the river Waal (one of the Rhine branches in the Netherlands). The morphological response of the river Waal to a flood protection measure (floodplain lowering in combination with summer levee removal) is analysed. The effect of this measure is subject to various sources of uncertainty. Monte Carlo simulations are applied to calculate the impact of uncertainties in the river discharge on the bed levels. The impact of the “uncertain” morphological response on design flood level predictions is analysed for three phenomena, viz. the impact of the spatial morphological variation over years, the impact of the seasonal morphological variation and the impact of the morphological variability around bifurcation points. The impact of seasonal morphological variations turns out to be negligible, but the other two phenomena appear to have each an appreciable impact (order of magnitude 0.05–0.1 m) on the computed design water levels. We have to note however, that other sources of uncertainty (e.g. uncertainty in hydraulic roughness predictor), which may be of influence, are not taken into consideration. In fact, the present investigation is limited to the sensitivity of the design water levels to uncertainties in the predicted bed level.     

Prediction of Flyrock in Mine Blasting: A New Computational Intelligence Approach

Natural Resources Research - Blasting is the predominant rock fragmentation technique in civil constructions, underground and surface mines. Flyrock is the unwanted throw of rock fragments during...     

Bathymetry enhancement by altimetry-derived gravity anomalies in the East Sea (Sea of Japan)

The gravity-geologic method (GGM) was used to enhance the bathymetry of the East Sea (Sea of Japan) with satellite altimetry-derived free-air gravity anomalies and shipborne depth measurements. By comparison with the bathymetry model of Smith and Sandwell’s (SAS) approach (1994), GGM was found to have an advantage with short wavelength (≤12 km) components, while SAS better predicts longer wavelength (≥25 km) components, despite its dependency on density contrast. To mitigate this limitation, a tuning density contrast of 10.25 g/cm³ between seawater and the seafloor was primarily estimated by the downward continuation method and then validated by the check points method with GGM. Similarly, SAS is limited by the “A” value in low-pass part of the Wiener filter, which defines the effective range of the wavelength components on bathymetry. As a final result, we present an enhanced GGM bathymetry model by integrating all available data.     

Astronomical calibration of the Matuyama-Brunhes boundary: Consequences for magnetic remanence acquisition in marine carbonates and the Asian loess sequences

We have compiled 19 records from marine carbonate cores in which the Matuyama-Brunhes boundary (MBB) has been reasonably well constrained within the astronomically forced stratigraphic framework using oxygen isotopes. By correlation of the δ¹⁸O data to a timescale based on astronomical forcing, we estimate astronomical ages for each of the MBB horizons. In all but one record the MBB occurs within Stage 19.

Most magnetostratigraphic sections in Asian Loess place the MBB within a loess interval. Since loess deposition is presumed to be associated with glacial intervals, loess horizons should correspond to even-numbered oxygen isotope stages. A glacial age for the MBB is at odds with the results presented here, which firmly place the MBB within interglacial Stage 19. Inconsistency among the many loess sections and between the loess and the marine records suggests that the magnetic interpretation of loess sections may be more complicated than hitherto supposed.

The mean of the Stage 19 age estimates for the MBB is 777.9 ± 1.8 (N = 18). Inclusion of the single Stage 20 age results in a mean of 778.8 ± 2.5 (N = 19). The astronomical age estimate of the MBB compares favorably with an (unweighted) mean of 778.2 ± 3.5 (N = 10) from a compilation of ⁴⁰Ar/39^Ar results of transitional lava flows. Combining the two independent data sets yields a grand mean of 778.0 ± 1.7 (N = 28).

The new compilation shows virtually no trend in placement of the MBB within isotope Stage 19 as a function of sediment accumulation rate. We interpret this to mean that the average depth of remanence acquisition is within a few centimeters of the sediment-water interface.

Separating the cores into two geographic regions (an Indo-Pacific-Caribbean [IPC] Group and an Atlantic Group) results in a significant difference in the position of the mid-point of the reversal with respect to the astronomical time scale. The data presented here suggest a difference of several thousand years between the two regions. This observation could be caused by systematic differences between the two regions in sedimentation rate within the interval of interest, systematic differences in remanence acquisition, or by genuine differences in the timing of the directional changes between the two regions.     

Laboratory sulfurisation of the marine microalga Nannochloropsis salina

To understand more fully the mode of preservation of organic matter in marine sediments, laboratory sulfurisation of intact cells of the cultured microalga Nannochloropsis salina was performed using inorganic polysulfides in seawater at 50°C. Solvent extractable and non-extractable materials were analysed by GC–MS and Py–GC–MS, respectively, to study the incorporation of sulfur into the microalgal organic matter. No GC-amenable sulfur-containing compounds were found in the extracts apart from some minor thiophenes with a phytanyl carbon skeleton. The residue after extraction and hydrolysis contained abundant macromolecular sulfur-containing moieties as revealed by the presence of dominant C₂₈–C₃₂ thiols, thiophenes, thianes and thiolanes in the flash pyrolysates. These products are thought to be formed from moieties derived from sulfurisation of C₂₈–C₃₂ diols and alkenols, characteristic lipids of N. salina. C₁–C₂ alkylated thiophenes were also found in the pyrolysates and probably result from moieties formed upon sulfurisation of carbohydrates. The highly resistant biomacromolecule (algaenan) synthesised by N. salina remains unaffected by sulfurisation. The non-hydrolysable residue isolated from the sulfurised N. salina thus comprises algaenan and (poly)sulfide-bound long alkyl chains. The sulfurisation experiments show that both selective preservation of algaenans and lipid and carbohydrate “vulcanisation” can be involved in the preservation of algal organic matter in marine environments.     

Simulation of the extreme waves generated by typhoon Bolaven (1215) in the East China Sea and Yellow Sea

Record-breaking high waves occurred during the passage of the typhoon Bolaven (1215) (TYB) in the East China Sea (ECS) and Yellow Sea (YS) although its intensity did not reach the level of a super typhoon. Winds and directional wave measurements were made using a range of in-situ instruments mounted on an ocean tower and buoys. In order to understand how such high waves with long duration occurred, analyses have been made through measurement and numerical simulations. TYB winds were generated using the TC96 typhoon wind model with the best track data calibrated with the measurements. And then the wind fields were blended with the reanalyzed synoptic-scale wind fields for a wave model. Wave fields were simulated using WAM4.5 with adjustment of C_d for gust of winds and bottom friction for the study area. Thus the accuracy of simulations is considerably enhanced, and the computed results are also in better agreement with measured data than before. It is found that the extremely high waves evolved as a result of the superposition of distant large swells and high wind seas generated by strong winds from the front/right quadrant of the typhoon track. As the typhoon moved at a speed a little slower than the dominant wave group velocity in a consistent direction for two days, the wave growth was significantly enhanced by strong wind input in an extended fetch and non-linear interaction.     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.     

Shear strength and compressibility of oyster shell–sand mixtures

This paper investigates the fundamental characteristics of shear strength and deformation of crushed oyster shell–sand mixtures to stimulate recycling of waste oyster shells. Standard penetration tests (SPT) and large-scale direct-shear tests were carried out with different kinds of dry unit weight and mixing rate of oyster shell–sand mixture. Correlations between N-value, dry unit weight, and friction angle of mixtures were observed from the results of experimental tests, making it possible to estimate the in situ strength from SPT, and the coefficient of volume compressibility from the confined direct-shear compression test. These results also make it possible to compute the settlement of oyster shell–sand mixture when used in soft ground improvement.