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.     

Image processing of 2D resistivity data for imaging faults

A methodology to locate automatically limits or boundaries between different geological bodies in 2D electrical tomography is proposed, using a crest line extraction process in gradient images. This method is applied on several synthetic models and on field data set acquired on three experimental sites during the European project PALEOSIS where trenches were dug. The results presented in this work are valid for electrical tomographies data collected with a Wenner-alpha array and computed with an l₁ norm (blocky inversion) as optimization method. For the synthetic cases, three geometric contexts are modelled: a vertical and a dipping fault juxtaposing two different geological formations and a step-like structure. A superficial layer can cover each geological structure. In these three situations, the method locates the synthetic faults and layer boundaries, and determines fault displacement but with several limitations. The estimated fault positions correlate exactly with the synthetic ones if a conductive (or no superficial) layer overlies the studied structure. When a resistive layer with a thickness of 6 m covers the model, faults are positioned with a maximum error of 1 m. Moreover, when a resistive and/or a thick top layer is present, the resolution significantly decreases for the fault displacement estimation (error up to 150%). The tests with the synthetic models for surveys using the Wenner-alpha array indicate that the proposed methodology is best suited to vertical and horizontal contacts. Application of the methodology to real data sets shows that a lateral resistivity contrast of 1:5–1:10 leads to exact faults location. A fault contact with a resistivity contrast of 1:0.75 and overlaid by a resistive layer with a thickness of 1 m gives an error location ranging from 1 to 3 m. Moreover, no result is obtained for a contact with very low contrasts (1:0.85) overlaid by a resistive soil. The method shows poor results when vertical gradients are greater than horizontal ones. This kind of image processing technique should be systematically used for improving the objectiveness of tomography interpretation when looking for limits between geological objects.     

Oxygen compounds in Titan's stratosphere as observed by Cassini CIRS

We have investigated the abundances of Titan's stratospheric oxygen compounds using 0.5 cm⁻¹ resolution spectra from the Composite Infrared Spectrometer on the Cassini orbiter. The CO abundance was derived for several observations of far-infrared nadir spectra, taken at a range of latitudes (75° S-35° N) and emission angles (0°-60°), using rotational lines that have not been analysed before the arrival of Cassini at Saturn. The derived volume mixing ratios for the different observations are mutually consistent regardless of latitude. The weighted mean CO volume mixing ratio is 47±8 ppm if CO is assumed to be uniform with latitude. H₂O could not be detected and an upper limit of 0.9 ppb was determined. CO₂ abundances derived from mid-infrared nadir spectra show no significant latitudinal variations, with typical values of 16±2 ppb. Mid-infrared limb spectra at 55° S were used to constrain the vertical profile of CO₂ for the first time. A vertical CO₂ profile that is constant above the condensation level at a volume mixing ratio of 15 ppb reproduces the limb spectra very well below 200 km. This is consistent with the long chemical lifetime of CO₂ in Titan's stratosphere. Above 200 km the CO₂ volume mixing ratio is not well constrained and an increase with altitude cannot be ruled out there.     

Isolation of protein from Chlorella sorokiniana CY1 using liquid biphasic flotation assisted with sonication through sugaring-out effect

Microalgae, a sustainable source of multi beneficial components has been discovered and could be utilised in pharmaceutical, bioenergy and food applications. This study aims to investigate the sugaring-out effect on the recovery of protein from wet green microalga, Chlorella sorokiniana CY1 which was assisted with sonication. A comparison of monosaccharides and disaccharides as one of the phaseforming constituents shows that the monosaccharides, glucose was the most suitable sugar in forming the phases with acetonitrile to enhance the production of protein(52% of protein). The protein productivity of microalgae was found to be significantly influenced by the volume ratio of both phases, as the yield of protein increased to 77%. The interval time between the sonication as well as the sonication modes were influencing the protein productivity as well. The optimum protein productivity was obtained with 10 s of resting time in between sonication. Pulse mode of sonication was suitable to break down the cell wall of microalgae compared to continuous mode as a lower protein yield was obtained with the application of continuous mode. The optimum condition for protein extraction were found as followed: 200 g/L glucose as bottom phase with volume ratio of I:1.25, 10 s of resting time for ultrasonication, 5 s of ultrasonication in pulse mode and 0.25 g of biomass weight. The high yield of protein about 81% could be obtained from microalgae which demonstrates the potential of this source and expected to play an important role in the future.     

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.     

Pilot‐scale field validation of the long electrode electrical resistivity tomography method

A validation experiment, carried out in a scaled field setting, was attempted for the long electrode electrical resistivity tomography method in order to demonstrate the performance of the technique in imaging a simple buried target. The experiment was an approximately 1/17 scale mock‐up of a region encompassing a buried nuclear waste tank on the Hanford site. The target of focus was constructed by manually forming a simulated plume within the vadose zone using a tank waste simulant. The long electrode results were compared to results from conventional point electrodes on the surface and buried within the survey domain. Using a pole‐pole array, both point and long electrode imaging techniques identified the lateral extents of the pre‐formed plume with reasonable fidelity but the long electrode method was handicapped in reconstructing vertical boundaries. The pole‐dipole and dipole‐dipole arrays were also tested with the long electrode method and were shown to have the least favourable target properties, including the position of the reconstructed plume relative to the known plume and the intensity of false positive targets. The poor performance of the pole‐dipole and dipole‐dipole arrays was attributed to an inexhaustive and non‐optimal coverage of data at key electrodes, as well as an increased noise for electrode combinations with high geometric factors. However, when comparing the model resolution matrix among the different acquisition strategies, the pole‐dipole and dipole‐dipole arrays using long electrodes were shown to have significantly higher average and maximum values within the matrix than any pole‐pole array. The model resolution describes how well the inversion model resolves the subsurface. Given the model resolution performance of the pole‐dipole and dipole‐dipole arrays, it may be worth investing in tools to understand the optimum subset of randomly distributed electrode pairs to produce maximum performance from the inversion model.     

Isolation of protein from Chlorella sorokiniana CY1 using liquid biphasic flotation assisted with sonication through sugaring-out effect

Microalgae,a sustainable source of multi beneficial components has been discovered and could be utilised in pharmaceutical,bioenergy and food applications.This study aims to investigate the sugaring-out effect on the recovery of protein from wet green microalga,Chlorella sorokiniana CY1 which was assisted with sonication.A comparison of monosaccharides and disaccharides as one of the phaseforming constituents shows that the monosaccharides,glucose was the most suitable sugar in forming the phases with acetonitrile to enhance the production of protein(52% of protein).The protein productivity of microalgae was found to be significantly influenced by the volume ratio of both phases,as the yield of protein increased to 77%.The interval time between the sonication as well as the sonication modes were influencing the protein productivity as well.The optimum protein productivity was obtained with 10s of resting time in between sonication.Pulse mode of sonication was suitable to break down the cell wall of microalgae compared to continuous mode as a lower protein yield was obtained with the application of continuous mode.The optimum condition for protein extraction were found as followed:200g/L glucose as bottom phase with volume ratio of 1:1.25,10s of resting time for ultrasonication,5s of ultrasonication in pulse mode and 0.25g of biomass weight.The high yield of protein about 81% could be obtained from microalgae which demonstrates the potential of this source and expected to play an important role in the future.     

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.