A mass-conservative control volume-finite element method for solving Richards’ equation in heterogeneous porous media

We present a mass-conservative vertex-centred finite volume method for efficiently solving the mixed form of Richards’ equation in heterogeneous porous media. The spatial discretisation is particularly well-suited to heterogeneous media because it produces consistent flux approximations at quadrature points where material properties are continuous. Combined with the method of lines, the spatial discretisation gives a set of differential algebraic equations amenable to solution using higher-order implicit solvers. We investigate the solution of the mixed form using a Jacobian-free inexact Newton solver, which requires the solution of an extra variable for each node in the mesh compared to the pressure-head form. By exploiting the structure of the Jacobian for the mixed form, the size of the preconditioner is reduced to that for the pressure-head form, and there is minimal computational overhead for solving the mixed form. The proposed formulation is tested on two challenging test problems. The solutions from the new formulation offer conservation of mass at least one order of magnitude more accurate than a pressure head formulation, and the higher-order temporal integration significantly improves both the mass balance and computational efficiency of the solution.

A new control-volume finite-element scheme for heterogeneous porous media : application to the drying of softwood

An improved mesoscopic model is presented for simulating the drying of porous media. The aim of this model is to account for two scales simultaneously: the scale of the whole product and the scale of the heterogeneities of the porous medium. The innovation of this method is the utilization of a new mass-conservative scheme based on the Control-Volume Finite-Element (CV-FE) method that partitions the moisture content field over the individual sub-control volumes surrounding each node within the mesh. Although the new formulation has potential for application across a wide range of transport processes in heterogeneous porous media, the focus here is on applying the model to the drying of small sections of softwood consisting of several growth rings. The results conclude that, when compared to a previously published scheme, only the new mass-conservative formulation correctly captures the true moisture content evolution in the earlywood and latewood components of the growth rings during drying.

Experimental and finite element analysis of a double strap joint between steel plates and normal modulus CFRP

Strengthening of steel structures using externally-bonded carbon fibre reinforced polymers ‘CFRP’ is a rapidly developing technique. This paper describes the behaviour of axially loaded flat steel plates strengthened using carbon fibre reinforced polymer sheets. Two steel plates were joined together with adhesive and followed by the application of carbon fibre sheet double strap joint with different bond lengths. The behaviour of the specimens was further investigated by using nonlinear finite element analysis to predict the failure modes and load capacity. In this study, bond failure is the dominant failure mode for normal modulus (240 GPa) CFRP bonding which closely matched the results of finite elements. The predicted ultimate loads from the FE analysis are found to be in good agreement with experimental values.

Volcanology and petrology of submarine volcanoes of the New Hebrides island arc

The New Hebrides Island Arc, an intra-oceanic island chain in the southwest Pacific, is formed by subduction of the Indo-Australian Plate beneath the Pacific Plate. The southern end of the New Hebrides Island Arc is an ideal location to study the magmatic and tectonic interaction of an emerging island arc as this part of the island chain is less than 3 million years old. A tectonically complex island arc, it exhibits a change in relative subduction rate from ~12cm/yr to 6 cm/yr before transitioning to a left-lateral strike slip zone at its southern end. Two submarine volcanic fields, Gemini-Oscostar and Volsmar, occur at this transition from normal arc subduction to sinistral strike slip movement. Multi-beam bathymetry and dredge samples collected during the 2004 CoTroVE cruise onboard the RV Southern Surveyor help define the relationship between magmatism and tectonics, and the source for these two submarine volcanic fields. Gemini-Oscostar volcanic field (GOVF), dominated by northwest-oriented normal faults, has mature polygenetic stratovolcanoes with evidence for explosive subaqueous eruptions and homogeneous monogenetic scoria cones. Volsmar volcanic field (VVF), located 30 km south of GOVF, exhibits a conjugate set of northwest and eastwest-oriented normal faults, with two polygenetic stratovolcanoes and numerous monogenetic scoria cones. A deep water caldera provides evidence for explosive eruptions at 1500m below sea level in the VVF. Both volcanic fields are dominated by low-K island arc tholeiites and basaltic andesites with calcalkalic andesite and dacite being found only in the GOVF. Geochemical signatures of both volcanic fields continue the along-arc trend of decreasing K2O with both volcanic fields being similar to the New Hebrides central chain lavas. Lavas from both fields display a slight depletion in high field strength elements and heavy rare earth elements, and slight enrichments in large-ion lithophile elements and light rare earth elements with respect to N-MORB mantle. Sr and Nd isotope data correlate with heavy rare earth and high field strength element data to show that both fields are derived from depleted mantle. Pb isotopes define Pacific MORB mantle sources and are consistent with isotopic variation along the New Hebrides Island Arc. Pb isotopes show no evidence for sediment contamination; the subduction component enrichment is therefore a slab-derived enrichment. There is a subtle spatial variation in source chemistry which sees a northerly trend of decreasing enrichment of slab-derived fluids.