Improved estimates of percolation and anisotropic permeability from 3-D X-ray microtomography using stochastic analyses and visualization

X-ray microtomography (micro-CT) with micron resolution enables new ways of characterizing microstructures and opens pathways for forward calculations of multiscale rock properties. A quantitative characterization of the microstructure is the first step in this challenge. We developed a new approach to extract scale-dependent characteristics of porosity, percolation, and anisotropic permeability from 3-D microstructural models of rocks. The Hoshen-Kopelman algorithm of percolation theory is employed for a standard percolation analysis. The anisotropy of permeability is calculated by means of the star volume distribution approach. The local porosity distribution and local percolation probability are obtained by using the local porosity theory. Additionally, the local anisotropy distribution is defined and analyzed through two empirical probability density functions, the isotropy index and the elongation index. For such a high-resolution data set, the typical data sizes of the CT images are on the order of gigabytes to tens of gigabytes; thus an extremely large number of calculations are required. To resolve this large memory problem parallelization in OpenMP was used to optimally harness the shared memory infrastructure on cache coherent Non-Uniform Memory Access architecture machines such as the iVEC SGI Altix 3700Bx2 Supercomputer. We see adequate visualization of the results as an important element in this first pioneering study.

Reconstructing 3D x-ray CT images of polymer gel dosimeters using the zero-scan method

In this study x-ray CT has been used to produce a 3D image of an irradiated PAGAT gel sample, with noise-reduction achieved using the ‘zero-scan’ method. The gel was repeatedly CT scanned and a linear fit to the varying Hounsfield unit of each pixel in the 3D volume was evaluated across the repeated scans, allowing a zero-scan extrapolation of the image to be obtained. To minimise heating of the CT scanner’s x-ray tube, this study used a large slice thickness (1 cm), to provide image slices across the irradiated region of the gel, and a relatively small number of CT scans (63), to extrapolate the zero-scan image. The resulting set of transverse images shows reduced noise compared to images from the initial CT scan of the gel, without being degraded by the additional radiation dose delivered to the gel during the repeated scanning. The full, 3D image of the gel has a low spatial resolution in the longitudinal direction, due to the selected scan parameters. Nonetheless, important features of the dose distribution are apparent in the 3D x-ray CT scan of the gel. The results of this study demonstrate that the zero-scan extrapolation method can be applied to the reconstruction of multiple x-ray CT slices, to provide useful 2D and 3D images of irradiated dosimetry gels.

A feasibility study of multislice x-ray CT imaging of gel dosimeters using the “zero scan” method

This study extends the ‘zero scan’ method for CT imaging of polymer gel dosimeters to include multi-slice acquisitions. Multi slice CT images consisting of 24 slices of 1.2 mm thickness were acquired of an irradiated polymer gel dosimeter, and processed with the zero scan technique. The results demonstrate that zero scan based gel readout can be successfully applied to generate a three dimensional image of the irradiated gel field. Compared to the raw CT images the processed figures and cross gel profiles demonstrated reduced noise and clear visibility of the penumbral region. Moreover these improved results further highlight the suitability of this method in volumetric reconstruction with reduced CT data acquisition per slice. This work shows that 3D volumes of irradiated polymer gel dosimeters can be acquired and processed with x-ray CT.

The first x-ray diffraction peak of water as a function of temperature in the range 4 – 63 °C

X-ray diffraction structure functions for water flowing in a 1.5 mm diameter siphon in the temperature range 4 – 63 °C were obtained using a 20 keV beam at the Australian Synchrotron. These functions were compared with structure functions obtained at the Advanced Light Source for a 0.5 mm thick sample of water in the temperature range 1 – 77 °C irradiated with an 11 keV beam. The two sets of structure functions are similar, but there are subtle differences in the shape and relative position of the two functions suggesting a possible differences between the structure of bulk and siphon water. In addition, the first structural peak (Q0) for water in a siphon, showed evidence of a step-wise increase in Q0 with increasing temperature rather than a smoothly varying increase. More experiments are required to investigate this apparent difference.

A comparison of surface doses for very small field size x-ray beams : Monte Carlo calculations and radiochromic film measurements

Stereotactic radiosurgery treatments involve the delivery of very high doses for a small number of fractions. To date, there is limited data in terms of the skin dose for the very small field sizes used in these treatments. In this work, we determine relative surface doses for small size circular collimators as used in stereotactic radiosurgery treatments. Monte Carlo calculations were performed using the BEAMnrc code with a model of the Novalis 15 Trilogy linear accelerator and the BrainLab circular collimators. The surface doses were calculated at the ICRU skin dose depth of 70 m all using the 6 MV SRS x-ray beam. The calculated surface doses varied between 15 – 12% with decreasing values as the field size increased from 4 to 30 mm. In comparison, surface doses were measured using Gafchromic EBT3 film positioned at the surface of a Virtual Water phantom. The absolute agreement between calculated and measured surface doses was better than 2.5% which is well within the 20 uncertainties of the Monte Carlo calculations and the film measurements. Based on these results, we have shown that the Gafchromic EBT3 film is suitable for surface dose estimates in very small size fields as used in SRS.

Factors affecting the extraction of absorbed dose information in a 3D polymer gel dosimeter by X ray computed tomography (IAEA-CN-96/122P)

Two sources of uncertainty in the X ray computed tomography imaging of polymer gel dosimeters are investigated in the paper.The first cause is a change in postirradiation density, which is proportional to the computed tomography signal and is associated with a volume change. The second cause of uncertainty is reconstruction noise.A simple technique that increases the residual signal to noise ratio by almost two orders of magnitude is examined.

Small angle X-ray scattering mapping and kinetics study of sub-critical CO2 sorption by two Australian coals

Time- and position-resolved synchrotron small angle X-ray scattering data were acquired from samples of two Australian coal seams: Bulli seam (Bulli 4, Ro=1.42%, Sydney Basin), which naturally contains CO2 and Baralaba seam (Ro=0.67%, Bowen Basin), a potential candidate for sequestering CO2. This experimental approach has provided unique, pore-size-specific insights into the kinetics of CO2 sorption in the micro- and small mesopores (diameter 5 to 175 Å) and the density of the sorbed CO2 at reservoir-like conditions of temperature and hydrostatic pressure. For both samples, at pressures above 5 bar, the density of CO2 confined in pores was found to be uniform, with no densification in near-wall regions. In the Bulli 4 sample, CO2 first flooded the slit pores between polyaromatic sheets. In the pore-size range analysed, the confined CO2 density was close to that of the free CO2. The kinetics data are too noisy for reliable quantitative analysis, but qualitatively indicate faster kinetics in mineral-matter-rich regions. In the Baralaba sample, CO2 preferentially invaded the smallest micropores and the confined CO2 density was up to five times that of the free CO2. Faster CO2 sorption kinetics was found to be correlated with higher mineral matter content but, the mineral-matter-rich regions had lower-density CO2 confined in their pores. Remarkably, the kinetics was pore-size dependent, being faster for smaller pores. These results suggest that injection into the permeable section of an interbedded coal-clastic sequence could provide a viable combination of reasonable injectivity and high sorption capacity.

Surface dosimetry for very small X-ray beams

Introduction The dose to skin surface is an important factor for many radiotherapy treatment techniques. It is known that TPS predicted surface doses can be significantly different from actual ICRP skin doses as defined at 70 lm. A number of methods have been implemented for the accurate determination of surface dose including use of specific dosimeters such as TLDs and radiochromic film as well as Monte Carlo calculations. Stereotactic radiosurgery involves delivering very high doses per treatment fraction using small X-ray fields. To date, there has been limited data on surface doses for these very small field sizes. The purpose of this work is to evaluate surface doses by both measurements and Monte Carlo calculations for very small field sizes. Methods All measurements were performed on a Novalis Tx linear accelerator which has a 6 MV SRS X-ray beam mode which uses a specially thin flattening filter. Beam collimation was achieved by circular cones with apertures that gave field sizes ranging from 4 to 30 mm at the isocentre. The relative surface doses were measured using Gafchromic EBT3 film which has the active layer at a depth similar to the ICRP skin dose depth. Monte Carlo calculations were performed using the BEAMnrc/EGSnrc Monte Carlo codes (V4 r225). The specifications of the linear accelerator, including the collimator, were provided by the manufacturer. Optimisation of the incident X-ray beam was achieved by an iterative adjustment of the energy, spatial distribution and radial spread of the incident electron beam striking the target. The energy cutoff parameters were PCUT = 0.01 MeV and ECUT = 0.700 - MeV. Directional bremsstrahlung splitting was switched on for all BEAMnrc calculations. Relative surface doses were determined in a layer defined in a water phantom of the same thickness and depth as compared to the active later in the film. Results Measured surface doses using the EBT3 film varied between 13 and 16 % for the different cones with an uncertainty of 3 %. Monte Carlo calculated surface doses were in agreement to better than 2 % to the measured doses for all the treatment cones. Discussion and conclusions This work has shown the consistency of surface dose measurements using EBT3 film with Monte Carlo predicted values within the uncertainty of the measurements. As such, EBT3 film is recommended for in vivo surface dose measurements.

An investigation of goethite-seeded Al(OH)3 precipitation using in situ X-ray diffraction and Rietveld-based quantitative phase analysis

An in situ X-ray diffraction investigation of goethite-seeded Al(OH)3 precipitation from synthetic Bayer liquor at 343 K has been performed. The presence of iron oxides and oxyhydroxides in the Bayer process has implications for alumina reversion, which causes significant process losses through unwanted gibbsite precipitation, and is also relevant for the nucleation and growth of scale on mild steel process equipment. The gibbsite, bayerite and nordstrandite polymorphs of Al(OH)3 precipitated from the liquor; gibbsite appeared to precipitate first, with subsequent formation of bayerite and nordstrandite. A Rietveld-based approach to quantitative phase analysis was implemented for the determination of absolute phase abundances as a function of time, from which kinetic information for the formation of the Al(OH)3 phases was determined.

Assessment of body composition in Indian adults: comparison between dual-energy X-ray absorptiometry and isotope dilution technique

Dual-energy X-ray absorptiometry (DXA) and isotope dilution technique have been used as reference methods to validate the estimates of body composition by simple field techniques; however, very few studies have compared these two methods. We compared the estimates of body composition by DXA and isotope dilution (18O) technique in apparently healthy Indian men and women (aged 19–70 years, n 152, 48 % men) with a wide range of BMI (14–40 kg/m2). Isotopic enrichment was assessed by isotope ratio mass spectroscopy. The agreement between the estimates of body composition measured by the two techniques was assessed by the Bland–Altman method. The mean age and BMI were 37 (SD 15) years and 23·3 (SD 5·1) kg/m2, respectively, for men and 37 (SD 14) years and 24·1 (SD 5·8) kg/m2, respectively, for women. The estimates of fat-free mass were higher by about 7 (95 % CI 6, 9) %, those of fat mass were lower by about 21 (95 % CI 218,223) %, and those of body fat percentage (BF%) were lower by about 7·4 (95 % CI 28·2, 26·6) % as obtained by DXA compared with the isotope dilution technique. The Bland–Altman analysis showed wide limits of agreement that indicated poor agreement between the methods. The bias in the estimates of BF% was higher at the lower values of BF%. Thus, the two commonly used reference methods showed substantial differences in the estimates of body composition with wide limits of agreement. As the estimates of body composition are method-dependent, the two methods cannot be used interchangeably

Analysis of ancient pottery and ceramic objects using X-ray fluorescence spectrometry

Archaeology has been called 'the science of the artefact' and nothing demonstrates this point better than the current interest displayed in provenance studies of archaeological objects. In theory, every vessel carries a chemical compositional pattern or 'fingerprint' identical with the clay from which it was made and this relationship is basic to provenance studies. The reasoning behind provenance or sourcing studies is to probe into this past and attempt to re-create prehistory by obtaining information on exchange and social interaction. This paper discusses the use of XRF spectrometry for the analysis of ancient pottery and ceramics to examine whether it is possible to predict prehictoric cultural exchanges.

Assessment of the effects of CT dose in averaged x-ray CT images of a dose-sensitive polymer gel

The signal-to-noise ratio achievable in x-ray computed tomography (CT) images of polymer gels can be increased by averaging over multiple scans of each sample. However, repeated scanning delivers a small additional dose to the gel which may compromise the accuracy of the dose measurement. In this study, a NIPAM-based polymer gel was irradiated and then CT scanned 25 times, with the resulting data used to derive an averaged image and a "zero-scan" image of the gel. Comparison between these two results and the first scan of the gel showed that the averaged and zero-scan images provided better contrast, higher contrast-to- noise and higher signal-to-noise than the initial scan. The pixel values (Hounsfield units, HU) in the averaged image were not noticeably elevated, compared to the zero-scan result and the gradients used in the linear extrapolation of the zero-scan images were small and symmetrically distributed around zero. These results indicate that the averaged image was not artificially lightened by the small, additional dose delivered during CT scanning. This work demonstrates the broader usefulness of the zero-scan method as a means to verify the dosimetric accuracy of gel images derived from averaged x-ray CT data.

Crystal structures and hydrogen-bonding in the co-crystalline adducts of 3,5-dinitrobenzoic acid with 4-aminosalicylic acid and 2-hydroxy-3-(1H-indol-3-yl)propenoic acid

The structures of the cocrystalline adducts of 3,5-dinitrobenzoic acid (3,5-DNBA) with 4-aminosalicylic acid (PASA), the 1:1 partial hydrate, C7H4N2O6 .C7H7NO3 . 2H2O, (I) and 2-hydroxy-3-(1H-indol-3-yl)propenoic acid (HIPA) and the 1:1:1 d6-dimethylsulfoxide solvate, C7H4N2O6 . C11H9NO3 . C2D6OS, (II) are reported. The crystal substructure of (I) comprises two centrosymmetric hydrogen-bonded R2/2(8) homodimers, one with 3,5-DNBA, the other with PASA, and an R2/2(8) 3,5-DNBA-PASA heterodimer. In the crystal, inter-unit amine N-H...O and water O-H...O hydrogen bonds generate a three-dimensional supramolecular structure. In (II), the asymmetric unit consists of the three constituent molecules which form an essentially planar cyclic hydrogen-bonded heterotrimer unit [graph set R2/3(17)] through carboxyl, hydroxy and amino groups. These units associate across a crystallographic inversion centre through the HIPA carboxylic acid group in an R2/2~(8) hydrogen-bonding association, giving a zero-dimensional structure lying parallel to (100). In both structures, pi--pi interactions are present [minimum ring centroid separations: 3.6471(18)A in (I) and 3.5819(10)A in (II)].

Crystal structures of the co-crystalline adduct 5-(4-bromophenyl)-1,3,4-thiadiazol-2-amine--4-nitrobenzoic acid (1/1) and the salt 2-amino-5-(4-bromophenyl)-1,3,4-thiadiazol-3-ium 2-carboxy-4,6-dinitrophenolate

The structures of the 1:1 co-crystalline adduct C8H6BrN3S . C7H5NO4 (I) and the salt C8H7BrN3S+ C7H3N2O7- (II) from the interaction of 5-(4-bromophenyl)-1,3,4-thiadiazol-2-amine with 4-nitrobenzoic acid and 3,5-dinitrosalicylic acid, respectively, have been determined. The primary inter-species association in both (I) and (II) is through duplex R2/2(8) (N-H...O/O-H...O) or (N-H...O/N-H...O) hydrogen bonds, respectively, giving heterodimers. In (II), these are close to planar [dihedral angles between the thiadiazole ring and the two phenyl rings are 2.1(3)deg. (intra) and 9.8(2)deg. (inter)], while in (I) these angles are 22.11(15) and 26.08(18)deg., respectively. In the crystal of (I), the heterodimers are extended into a one-dimensional chain along b through an amine N-...N(thiadiazole) hydrogen bond but in (II), a centrosymmetric cyclic heterotetramer structure is generated through N-H...O hydrogen bonds to phenol and nitro O-atom acceptors and features, together with the primary R2/2(8) interaction, conjoined R4/6(12), R2/1(6) and S(6) ring motifs. Also present in (I) are pi--pi interactions between thiadiazole rings [minimum ring centroid separation, 3.4624(16)deg.] as well as short Br...O(nitro) interactions in both (I) and (II) [3.296(3)A and 3.104(3)A, respectively].

Crystal structures and hydrogen bonding in the proton transfer salts of nicotine with 3,5-dinitrosalicylic acid and 5-sulfosalicylic acid

The structures of the 1:1 anhydrous salts of nicotine (NIC) with 3,5-dinitrosalicylic acid (DNSA) and 5-sulfosalicylic acid (5-SSA), namely (1R,2S)-1-methyl-2-(3-pyridyl)-1H-pyrrolidin-1-ium 2-carboxy-4,6-dinitrophenolate, C10H15N2+ C7H3N2O7-, (I) and (1R,2S)-1-methyl-2-(3-pyridyl)-1H-pyrrolidin-1-ium 3-carboxy-4-hydroxybenzenesulfonate, C10H15N2+ C7H5O6S-, (II) are reported. The asymmetric units of both (I) and (II) comprise two independent nicotinium cations (C and D) and either two DNSA or two 5-SSA anions (A and B), respectively. One of the DNSA anions shows a 25% rotational disorder in the benzene ring system. In the crystal of (I), inter-unit pyrrolidinium N-H...N(pyridine) hydrogen bonds generate zigzag NIC cation chains which extend along a while the DNSA anions are not involved in any formal inter-species hydrogen bonding but instead form pi--pi associated stacks which parallel the NIC chains along a [ring centroid separation, 3.857(2)A]. Weak C-H...O interactions between chain substructures give an overall three-dimensional structure. With (II), A and B anions form independent zigzag chains with C and D cations, respectively, through carboxylic acid O-H...N(pyridine) hydrogen bonds. These chains, which extend along b are pseudo-centrosymmetrically related and give pi--pi interactions between the benzene rings of anions A and B and the pyridine rings of the NIC cations C and D, respectively [ring centroid separations, 3.6422(19) and 3.7117(19)A]. Present also are weak intermolecular C-H...O hydrogen-bonding interactions between the chains, giving an overall three-dimensional structure.