Soil moisture monitoring at the field scale using neutron probe

Measurement of the moisture variation in soils is required for geotechnical design and research because soil properties and behavior can vary as moisture content changes. The neutron probe, which was developed more than 40 years ago, is commonly used to monitor soil moisture variation in the field. This study reports a full-scale field monitoring of soil moisture using a neutron moisture probe for a period of more than 2 years in the Melbourne (Australia) region. On the basis of soil types available in the Melbourne region, 23 sites were chosen for moisture monitoring down to a depth of 1500 mm. The field calibration method was used to develop correlations relating the volumetric moisture content and neutron counts. Observed results showed that the deepest “wetting front” during the wet season was limited to the top 800 to 1000 mm of soil whilst the top soil layer down to about 550mmresponded almost immediately to the rainfall events. At greater depths (550 to 800mmand below 800 mm), the moisture variations were relatively low and displayed predominantly periodic fluctuations. This periodic nature was captured with Fourier analysis to develop a cyclic moisture model on the basis of an analytical solution of a one-dimensional moisture flow equation for homogeneous soils. It is argued that the model developed can be used to predict the soil moisture variations as applicable to buried structures such as pipes.

A simple method for EPID-based in vivo dosimetry for radiotherapy treatments of head-and-neck cancers

This study used a homogeneous water-equivalent model of an electronic portal imaging device (EPID), contoured as a structure in a radiotherapy treatment plan, to produce reference dose images for comparison with in vivo EPID dosimetry images. Head and neck treatments were chosen as the focus of this study, due to the heterogeneous anatomies involved and the consequent difficulty of rapidly obtaining reliable reference dose images by other means. A phantom approximating the size and heterogeneity of a typical neck, with a maximum radiological thickness of 8.5 cm, was constructed for use in this study. This phantom was CT scanned and a simple treatment including five square test fields and one off-axis IMRT field was planned. In order to allow the treatment planning system to calculate dose in a model EPID positioned a distance downstream from the phantom to achieve a source-to-detector distance (SDD) of 150 cm, the CT images were padded with air and the phantom’s “body” contour was extended to encompass the EPID contour. Comparison of dose images obtained from treatment planning calculations and experimental irradiations showed good agreement, with more than 90% of points in all fields passing a gamma evaluation, at γ (3%, 3mm )Similar agreement was achieved when the phantom was over-written with air in the treatment plan and removed from the experimental beam, suggesting that water EPID model at 150 cm SDD is capable of providing accurate reference images for comparison with clinical IMRT treatment images, for patient anatomies with radiological thicknesses ranging from 0 up to approximately 9 cm. This methodology therefore has the potential to be used for in vivo dosimetry during treatments to tissues in the neck as well as the oral and nasal cavities, in the head-and-neck region.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Computer simulation of scaffold degradation

Scaffolds are porous biocompatible materials with suitable microarchitectures that are designed to allow for cell adhesion, growth and proliferation. They are used in combination with cells in regenerative medicine to promote tissue regeneration by means of a controlled deposition of natural extracellular matrix by the hosted cells therein. This healing process is in many cases accompanied by scaffold degradation up to its total disappearance when the scaffold is made of a biodegradable material. This work presents a computational model that simulates the degradation of scaffolds. The model works with three-dimensional microstructures, which have been previously discretised into small cubic homogeneous elements, called voxels. The model simulates the evolution of the degradation of the scaffold using a Monte Carlo algorithm, which takes into account the curvature of the surface of the fibres. The simulation results obtained in this study are in good agreement with empirical degradation measurements performed by mass loss on scaffolds after exposure to an etching alkaline solution.

A review on the production of levulinic acid and furanics from sugars

Sugarcane products represent an abundant and relatively low cost carbon resource that can be utilised to produce chemical intermediates such as levulinic acid and furanics. These chemicals can be easily upgraded to commodity and specialty chemicals and biofuels by high yielding and well established technologies. However, there are challenges and technical hurdles that need to be overcome before these chemical intermediates can be cost-effectively produced in commercial quantities. The paper reviews production of levulinic acid and furanics from sugars by homogeneous mineral acid catalysts, and reports on preliminary studies on the production of these compounds with environmentally friendly biodegradable sulfonic acids. The yields (>50% of theoretical) of levulinic acid, formic acid and furfural obtained with these organic acids are comparable to that of sulphuric acid currently used for their production.

Automatic classification of human epithelial type 2 cell indirect immunofluorescence images using cell pyramid matching

This paper describes a novel system for automatic classification of images obtained from Anti-Nuclear Antibody (ANA) pathology tests on Human Epithelial type 2 (HEp-2) cells using the Indirect Immunofluorescence (IIF) protocol. The IIF protocol on HEp-2 cells has been the hallmark method to identify the presence of ANAs, due to its high sensitivity and the large range of antigens that can be detected. However, it suffers from numerous shortcomings, such as being subjective as well as time and labour intensive. Computer Aided Diagnostic (CAD) systems have been developed to address these problems, which automatically classify a HEp-2 cell image into one of its known patterns (eg. speckled, homogeneous). Most of the existing CAD systems use handpicked features to represent a HEp-2 cell image, which may only work in limited scenarios. We propose a novel automatic cell image classification method termed Cell Pyramid Matching (CPM), which is comprised of regional histograms of visual words coupled with the Multiple Kernel Learning framework. We present a study of several variations of generating histograms and show the efficacy of the system on two publicly available datasets: the ICPR HEp-2 cell classification contest dataset and the SNPHEp-2 dataset.

Comparison of mechanical and ultrasound elastic modulus of ovine tibial cortical bone

Finite element models of bones can be created by deriving geometry from anx-ray CT scan. Material properties such as the elastic modulus can then be applied using either a single or set of homogeneous values, or individual elements can have local values mapped onto them. Values for the elastic modulus can be derived from the CT density values using an elasticityversus density relationship. Many elasticity–density relationships have been reported in the literature for human bone. However, while ovine in vivo models are common in orthopaedic research, no work has been done to date on creating FE models of ovine bones. To create these models and apply relevant material properties, an ovine elasticity-density relationship needs to be determined. Using fresh frozen ovine tibias the apparent density of regions of interest was determined from a clinical CT scan. The bones were the sectioned into cuboid samples of cortical bone from the regions of interest. Ultrasound was used to determine the elastic modulus in each of three directions – longitudinally, radially and tangentially. Samples then underwent traditional compression testing in each direction. The relationships between apparent density and both ultrasound, and compression modulus in each directionwere determined. Ultrasound testing was found to be a highly repeatable non-destructive method of calculating the elastic modulus, particularly suited to samples of this size. The elasticity-density relationships determined in the longitudinal direction were very similar between the compression and ultrasound data over the density range examined.A clear difference was seen in the elastic modulus between the longitudinal and transverse directions of the bone samples, and a transverse elasticity-density relationship is also reported.

A light and electron microscopic study of NADPH-diaphorase-,calretinin- and parvalbumin-containing neurons in the rat nucleus accumbens

The rat nucleus accumbens contains medium-sized, spiny projection neurons and intrinsic, local circuit neurons, or interneurons. Sub-classes of interneurons, revealed by calretinin (CR) or parvalbumin (PV) immunoreactivity or reduced nicotinamide adenine dinucleotide phosphate (NADPH)-diaphorase histochemistry, were compared in the nucleus accumbens core, shell and rostral pole. CR, PV and NADPH-diaphorase-containing neurons are shown to form three non-co-localising populations in these three areas. No significant differences in neuronal population densities were found between the subterritories. NADPH-diaphorase-containing neurons could be further separated morphologically into three sub-groups, but CR- and PV-immunoreactive neurons form homogeneous populations. Ultrastructurally, NADPH-diaphorase-, CR- and PV-containing neurons in the nucleus accumbens all possess nuclear indentations. These are deeper and fewer in neurons immunoreactive for PV than in CR- and NADPH-diaphorase-containing neurons. CR-immunoreactive boutons form asymmetrical and symmetrical synaptic specialisations on spines, dendrites and somata, while PV-immunoreactive boutons make only symmetrical synaptic specialisations. Both CR- and PV-immunoreactive boutons form symmetrical synaptic specialisations with medium-sized spiny neurons and contact other CR- and PV-immunoreactive somata, respectively. A novel non-carcinogenic substrate for the peroxidase reaction (Vector Slate Grey, SG) was found to be characteristically electron-dense and may be distinguishable from the diaminobenzidine reaction product. We conclude that the three markers used in this study are localised in distinct populations of nucleus accumbens interneurons. Our studies of their synaptic connections contribute to an increased understanding of the intrinsic circuitry of this area.

Vibrational spectroscopy of the sulphate mineral sturmanite from Kuruman manganese deposits, South Africa

The mineral sturmanite is a hydrated calcium iron aluminium manganese sulphate tetrahydroxoborate hydroxide of formula Ca6(Fe, Al, Mn)2(SO4)2(B(OH)4)(OH)12•26H2O. We have studied the mineral sturmanite using a number of techniques, including SEM with EPMA and vibrational spectroscopy. Chemical analysis shows a homogeneous phase, composed by Ca, Fe, Mn, S, Al and Si. B is not determined in this EPMA technique. An intense Raman band at 990 cm−1 is assigned to the SO42− symmetric stretching mode. Raman spectroscopy identifies multiple sulphate symmetric stretching modes in line with the three sulphate crystallographically different sites. Raman spectroscopy also identifies a band at 1069 cm−1 which may be attributed to a carbonate symmetric stretching mode, indicating the presence of thaumasite. Infrared spectra display two bands at 1080 and 1107 cm−1 assigned to the SO42− antisymmetric stretching modes. The observation of multiple bands in this ν4 spectral region offers evidence for the reduction in symmetry of the sulphate anion from Td to C2v or even lower symmetry. The Raman band at 3622 cm−1 is assigned to the OH unit stretching vibration and the broad feature at around 3479 cm−1 to water stretching bands. Infrared spectroscopy shows a set of broad overlapping bands in the OH stretching region. Vibrational spectroscopy enables an assessment of the molecular structure of sturmanite to be made.

A vibrational spectroscopic study of the phosphate mineral rimkorolgite (Mg,Mn2+)5(Ba, Sr)(PO4)4·8H2O from Kovdor massif, Kola Peninsula, Russia

We have studied aspect of the molecular structure of the phosphate mineral rimkorolgite from Zheleznyi iron mine, Kovdor massif, Kola Peninsula, Russia, using SEM with EDX and vibrational spectroscopy. Qualitative chemical analysis shows a homogeneous phase, composed by P, Mg, Ba, Mn and Ca. Small amounts of Si were also observed. An intense Raman peak at 975 cm−1 is assigned to the PO43− ν1 symmetric stretching mode. The Raman band at 964 cm−1 is attributed to the HPO42− ν1 symmetric stretching vibration. Raman bands observed at 1016, 1035, 1052, 1073, 1105 and 1135 cm−1 are attributed to the ν3 antisymmetric stretching vibrations of the HPO42− and PO43− units. Complexity in the spectra of the phosphate bending region is observed. The broad Raman band at 3272 cm−1 is assigned to the water stretching vibration. Vibrational spectroscopy enables aspects on the molecular structure of rimkorolgite to be undertaken.

A vibrational spectroscopic study of the borate mineral takedaite Ca3(BO3)2

We have studied the mineral takedaite Ca3(BO3)2, a borate mineral of calcium using SEM with EDX and vibrational spectroscopy. Chemical analysis shows a homogeneous phase, composed of Ca. Boron was not detected. A very intense Raman band at 1087 cm−1 is assigned to the BO stretching vibration of BO3 units. Additional Raman bands may be due to isotopic splitting. In the infrared spectrum, bands at 1218 cm−1 and at 1163, 1262 and 1295 cm−1 are assigned to the trigonal borate stretching modes. Raman bands at 712 and 715 cm−1 are assigned to the in-plane bending modes of the BO3 units. Vibrational spectroscopy enables aspects of the molecular structure of takedaite to be assessed.

A comparison of stability and bifurcation criteria for inflated spherical elastic shells

The nonlinear stability analysis introduced by Chen and Haughton [1] is employed to study the full nonlinear stability of the non-homogeneous spherically symmetric deformation of an elastic thick-walled sphere. The shell is composed of an arbitrary homogeneous, incompressible elastic material. The stability criterion ultimately requires the solution of a third-order nonlinear ordinary differential equation. Numerical calculations performed for a wide variety of well-known incompressible materials are then compared with existing bifurcation results and are found to be identical. Further analysis and comparison between stability and bifurcation are conducted for the case of thin shells and we prove by direct calculation that the two criteria are identical for all modes and all materials.

Plasma–surface interactions at nanoscales : a combinatorial theoretical, process diagnostics and surface microanalysis approach

This paper introduces an integral approach to the study of plasma-surface interactions during the catalytic growth of selected nanostructures (NSs). This approach involves basic understanding of the plasma-specific effects in NS nucleation and growth, theoretical modelling, numerical simulations, plasma diagnostics, and surface microanalysis. Using an example of plasma-assisted growth of surface-supported single-walled carbon nanotubes, we discuss how the combination of these techniques may help improve the outcomes of the growth process. A specific focus here is on the effects of nanoscale plasma-surface interactions on the NS growth and how the available techniques may be used, both in situ and ex situ to optimize the growth process and structural parameters of NSs.

Crystalline Si nanoparticles below crystallization threshold : effects of collisional heating in non-thermal atmospheric-pressure microplasmas

Nucleation and growth of highly crystalline silicon nanoparticles in atmospheric-pressure low-temperature microplasmas at gas temperatures well below the Si crystallization threshold and within a short (100 μs) period of time are demonstrated and explained. The modeling reveals that collision-enhanced ion fluxes can effectively increase the heat flux on the nanoparticle surface and this heating is controlled by the ion density. It is shown that nanoparticles can be heated to temperatures above the crystallization threshold. These combined experimental and theoretical results confirm the effective heating and structure control of Si nanoparticles at atmospheric pressure and low gas temperatures.

Ion impact distribution over plasma exposed nanocone arrays

The effect of an ordered array of nanocones on a conducting substrate immersed in the plasma on the transport of the plasma ions is investigated. The real conical shape of the cones is rigorously incorporated into the model. The movement of 10^5 CH3+ ions in the plasma sheath modified by the nanocone array is simulated. The ions are driven by the electric fields produced by the sheath and the nanostructures. The surface charge density and the total charge on the nanotips with different aspect ratios are computed. The ion transport simulation provides important characteristics of the displacement and velocity of the ions. The relative ion distribution along the lateral surfaces of the carbon nanotips is computed as well. It is shown that a rigorous account of the realistic nanostructure shape leads to very different distribution of the ion fluxes on the nanostructured surfaces compared to the previously reported works. The ion flux distribution is a critical factor in the nucleation process on the substrate and determines the nanostructure growth patterns.