The determination of three selected sites for the construction and operation of a low level radioactive waste storage and burial facility

Radioactive wastes are by-products of the use of radiation technologies. As with many technologies, the wastes are required to be disposed of in a safe manner so as to minimise risk to human health. This study examines the requirements for a hypothetical repository and develops techniques for decision making to permit the establishment of a shallow ground burial facility to receive an inventory of low-level radioactive wastes. Australia’s overall inventory is used as an example. Essential and desirable siting criteria are developed and applied to Australia's Northern Territory resulting in the selection of three candidate sites for laboratory investigations into soil behaviour. The essential quantifiable factors which govern radionuclide migration and ultimately influence radiation doses following facility closure are reviewed. Simplified batch and column procedures were developed to enable laboratory determination of distribution and retardation coefficient values for use in one-dimensional advection-dispersion transport equations. Batch and column experiments were conducted with Australian soils sampled from the three identified candidate sites using a radionuclide representative of the current national low-level radioactive waste inventory. The experimental results are discussed and site soil performance compared. The experimental results are subsequently used to compare the relative radiation health risks between each of the three sites investigated. A recommendation is made as to the preferred site to construct an engineered near-surface burial facility to receive the Australian low-level radioactive waste inventory.

Pressure filtration of Australian bagasse pulp

A one-dimensional pressure filtration model that can be used to predict the behaviour of bagasse pulp has been developed and verified in this study.The dynamic filtration model uses steady state compressibility parameters determined experimentally by uniaxial loading. The compressibility parameters M and N for depithed bagasse pulp were determined to be in the ranges 3000–8000kPa and 2.5–3.0 units, respectively. The model also incorporates experimentally determined steady state permeability data from separate experiments to predict the pulp concentration and fibre pressure throughout a pulp mat during dynamic filtration. Under steady state conditions, a variable Kozeny factor required different values for the permeability parameters when compared to a constant Kozeny factor. The specific surface area was 25–30% lower and the swelling factor was 20–25% higher when a variable Kozeny factor was used. Excellent agreement between experimental data and the dynamic filtration model was achieved when a variable Kozeny factor was used.

Hydrogen bonding in the anhydrous 1:1 proton-transfer compounds of isonipecotamide with nitro-substituted benzoic acids: the salts of the three isomeric mononitrobenzoic acids and of 3,5-dinitrobenzoic acid

The structures of the anhydrous 1:1 proton-transfer compounds of isonipecotamide (piperidine-4-carboxamide) with the three isomeric mononitro-substituted benzoic acids and 3,5-dinitrobenzoic acid, namely 4-carbamoylpiperidinium 2-nitrobenzoate (I), 4-carbamoylpiperidinium 3-nitrobenzoate (II), 4-carbamoylpiperidinium 4-nitrobenzoate (III), (C6H13N2O+ C7H4NO4-) and 4-carbamoylpiperidinium 3,5-dinitrobenzoate (IV) (C6H13N2O+ C7H5N2O6-)respectively, have been determined at 200 K. All salts form hydrogen-bonded structures: three-dimensional in (I), two-dimensional in (II) and (III) and one-dimensional in (IV). Featured in the hydrogen bonding of three of these [(I), (II) and (IV)] is the cyclic head-to-head amide--amide homodimer motif [graph set R2/2~(8)] through a duplex N---H...O association, the dimer then giving structure extension via either piperidinium or amide H-donors and carboxylate-O and in some examples [(II) and (IV)], nitro-O atom acceptors. In (I), the centrosymmetric amide-amide homodimers are expanded laterally through N-H...O hydrogen bonds via cyclic R2/4(8) interactions forming ribbons which extend along the c cell direction. These ribbons incorporate the 2-nitrobenzoate cations through centrosymmetric cyclic piperidine N-H...O(carboxyl) associations [graph set R4/4(12)], giving inter-connected sheets in the three-dimensional structure. In (II) in which no amide-amide homodimer is present, duplex piperidinium N-H...O(amide) hydrogen-bonding homomolecular associations [graph set R2/2(14)] give centrosymmetric head-to-tail dimers. Structure extension occurs through hydrogen-bonding associations between both the amide H-donors and carboxyl and nitro O-acceptors as well as a three-centre piperidinium N-H...O,O'(carboxyl) cyclic R2/1(4) association giving the two-dimensional network structure. In (III), the centrosymmetric amide-amide dimers are linked through the two carboxyl O-atom acceptors of the anions via bridging piperidinium and amide N-H...O,O'...H-N(amide) hydrogen bonds giving the two-dimensional sheet structure which features centrosymmetric cyclic R4/4(12) associations. In (IV), the amide-amide dimer is also centrosymmetric with the dimers linked to the anions through amide N-H...O(nitro) interactions. The piperidinium groups extend the structure into one-dimensional ribbons via N-H...O(carboxyl) hydrogen bonds. The structures reported here further demonstrate the utility of the isonipecotamide cation in molecular assembly and highlight the efficacy of the cyclic R2/2(8) amide-amide hydrogen-bonding homodimer motif in this process and provide an additional homodimer motif type in the head-to-tail R2/2(14) association.

Contracting bubbles in Hele-Shaw cells with a power-law fluid

The problem of bubble contraction in a Hele-Shaw cell is studied for the case in which the surrounding fluid is of power-law type. A small perturbation of the radially symmetric problem is first considered, focussing on the behaviour just before the bubble vanishes, it being found that for shear-thinning fluids the radially symmetric solution is stable, while for shear-thickening fluids the aspect ratio of the bubble boundary increases. The borderline (Newtonian) case considered previously is neutrally stable, the bubble boundary becoming elliptic in shape with the eccentricity of the ellipse depending on the initial data. Further light is shed on the bubble contraction problem by considering a long thin Hele-Shaw cell: for early times the leading-order behaviour is one-dimensional in this limit; however, as the bubble contracts its evolution is ultimately determined by the solution of a Wiener-Hopf problem, the transition between the long-thin limit and the extinction limit in which the bubble vanishes being described by what is in effect a similarity solution of the second kind. This same solution describes the generic (slit-like) extinction behaviour for shear-thickening fluids, the interface profiles that generalise the ellipses that characterise the Newtonian case being constructed by the Wiener-Hopf calculation.

Lie group symmetry analysis for granular media stress equations

The Airy stress function, although frequently employed in classical linear elasticity, does not receive similar usage for granular media problems. For plane strain quasi-static deformations of a cohesionless Coulomb–Mohr granular solid, a single nonlinear partial differential equation is formulated for the Airy stress function by combining the equilibrium equations with the yield condition. This has certain advantages from the usual approach, in which two stress invariants and a stress angle are introduced, and a system of two partial differential equations is needed to describe the flow. In the present study, the symmetry analysis of differential equations is utilised for our single partial differential equation, and by computing an optimal system of one-dimensional Lie algebras, a complete set of group-invariant solutions is derived. By this it is meant that any group-invariant solution of the governing partial differential equation (provided it can be derived via the classical symmetries method) may be obtained as a member of this set by a suitable group transformation. For general values of the parameters (angle of internal friction and gravity g) it is found there are three distinct classes of solutions which correspond to granular flows considered previously in the literature. For the two limiting cases of high angle of internal friction and zero gravity, the governing partial differential equation admit larger families of Lie point symmetries, and from these symmetries, further solutions are derived, many of which are new. Furthermore, the majority of these solutions are exact, which is rare for granular flow, especially in the case of gravity driven flows.

Proton-transfer compounds of isnipecotamide with the aromatic dicarboxylic acids 4-nitrophthalic, 4,5-dichlorophthalic, 5-nitroisophthalic and terephthalic acid

The structures of the 1:1 proton-transfer compounds of isonipecotamide (4-piperidinecarboxamide) with 4-nitrophthalic acid, 4-carbamoylpiperidinium 2-carboxy-4-nitrobenzoate, C6H13N2O8+ C8H4O6- (I), 4,5-dichlorophthalic acid, 4-carbamoylpiperidinium 2-carboxy-4,5-dichlorobenzoate, C6H13N2O8+ C8H3Cl2O4- (II) and 5-nitroisophthalic acid, 4-carbamoylpiperidinium 3-carboxy-5-nitrobenzoate, C6H13N2O8+ C8H4O6- (III) as well as the 2:1 compound with terephthalic acid, bis(4-carbamoylpiperidinium)benzene-1,2-dicarboxylate dihydrate, 2(C6H13N2O8+) C8H4O42- . 2H2O (IV)have been determined at 200 K. All salts form hydrogen-bonded structures, one-dimensional in (II) and three-dimensional in (I), (III) and (IV). In (I) and (III) the centrosymmetric R2/2(8) cyclic amide-amide association is found while in (IV) several different types of water-bridged cyclic associations are present [graph sets R2/4(8), R3/4(10), R4/4(12), R3/3(18) and R4/6(22)]. The one-dimensional structure of (I), features the common 'planar' hydrogen 4,5-dichlorophthalate anion together with enlarged cyclic R3/3(13) and R3/4(17) associations. With the structures of (I) and (III) the presence of head-to-tail hydrogen phthalate chain substructures is found. In (IV) head-to-tail primary cation-anion associations are extended longitudinally into chains through the water-bridged cation associations and laterally by piperidinium N-H...O(carboxyl) and water O-H...O(carboxyl) hydrogen bonds. The structures reported here further demonstrate the utility of the isonipecotamide cation as a synthon for the generation of stable hydrogen-bonded structures. An additional example of cation--anion association with this cation is also shown in the asymmetric three-centre piperidinium N-H...O,O'(carboxyl) interaction in the first-reported structure of a 2:1 isonipecotamide-carboxylate salt.

Crystal growth from undercooling melts under strong gradients of temperatures generating in short-duration microgravity

Crystal growth of bulk CdTe in short-duration microgravity is performed by the unidirectional cooling method. The largest growth grains in microgravity samples are 4X2mm. The cooling profiles indicate undercooling melts in microgravity. Cooling melt samples in microgravity generate strong gradient of temperature due to stop thermal convections. Temperature distribution in the melt is calculated by the one-dimensional equation of heat conduction, and about 100 K-undercooling is considered to occur at the cooling surface.

Predicting residue-wise contact orders in proteins by support vector regression

Background The residue-wise contact order (RWCO) describes the sequence separations between the residues of interest and its contacting residues in a protein sequence. It is a new kind of one-dimensional protein structure that represents the extent of long-range contacts and is considered as a generalization of contact order. Together with secondary structure, accessible surface area, the B factor, and contact number, RWCO provides comprehensive and indispensable important information to reconstructing the protein three-dimensional structure from a set of one-dimensional structural properties. Accurately predicting RWCO values could have many important applications in protein three-dimensional structure prediction and protein folding rate prediction, and give deep insights into protein sequence-structure relationships. Results We developed a novel approach to predict residue-wise contact order values in proteins based on support vector regression (SVR), starting from primary amino acid sequences. We explored seven different sequence encoding schemes to examine their effects on the prediction performance, including local sequence in the form of PSI-BLAST profiles, local sequence plus amino acid composition, local sequence plus molecular weight, local sequence plus secondary structure predicted by PSIPRED, local sequence plus molecular weight and amino acid composition, local sequence plus molecular weight and predicted secondary structure, and local sequence plus molecular weight, amino acid composition and predicted secondary structure. When using local sequences with multiple sequence alignments in the form of PSI-BLAST profiles, we could predict the RWCO distribution with a Pearson correlation coefficient (CC) between the predicted and observed RWCO values of 0.55, and root mean square error (RMSE) of 0.82, based on a well-defined dataset with 680 protein sequences. Moreover, by incorporating global features such as molecular weight and amino acid composition we could further improve the prediction performance with the CC to 0.57 and an RMSE of 0.79. In addition, combining the predicted secondary structure by PSIPRED was found to significantly improve the prediction performance and could yield the best prediction accuracy with a CC of 0.60 and RMSE of 0.78, which provided at least comparable performance compared with the other existing methods. Conclusion The SVR method shows a prediction performance competitive with or at least comparable to the previously developed linear regression-based methods for predicting RWCO values. In contrast to support vector classification (SVC), SVR is very good at estimating the raw value profiles of the samples. The successful application of the SVR approach in this study reinforces the fact that support vector regression is a powerful tool in extracting the protein sequence-structure relationship and in estimating the protein structural profiles from amino acid sequences.

Shape discrimination using invariants defined from higher-order spectra

An approach to pattern recognition using invariant parameters based on higher-order spectra is presented. In particular, bispectral invariants are used to classify one-dimensional shapes. The bispectrum, which is translation invariant, is integrated along straight lines passing through the origin in bifrequency space. The phase of the integrated bispectrum is shown to be scale- and amplification-invariant. A minimal set of these invariants is selected as the feature vector for pattern classification. Pattern recognition using higher-order spectral invariants is fast, suited for parallel implementation, and works for signals corrupted by Gaussian noise. The classification technique is shown to distinguish two similar but different bolts given their one-dimensional profiles

Synthesis and characterisation of metal oxyhydroxide and oxide nanomaterials

In this work, a range of nanomaterials have been synthesised based on metal oxyhydroxides MO(OH), where M=Al, Co, Cr, etc. Through a self-assembly hydrothermal route, metal oxyhydroxide nanomaterials with various morphologies were successfully synthesised: one dimensional boehmite (AlO(OH)) nanofibres, zero dimensional indium hydroxide (In(OH)3) nanocubes and chromium oxyhydroxide (CrO(OH)) nanoparticles, as well as two dimensional cobalt hydroxide and oxyhydroxide (Co(OH)2 & CoO(OH)) nanodiscs. In order to control the synthetic nanomaterial morphology and growth, several factors were investigated including cation concentration, temperature, hydrothermal treatment time, and pH. Metal ion doping is a promising technique to modify and control the properties of materials by intentionally introducing impurities or defects into the material. Chromium was successfully applied as a dopant for fabricating doped boehmite nanofibres. The thermal stability of the boehmite nanofibres was enhanced by chromium doping, and the photoluminescence property was introduced to the chromium doped alumina nanofibres. Doping proved to be an efficient method to modify and functionalize nanomaterials. The synthesised nanomaterials were fully characterised by X-ray diffraction (XRD), transmission electron microscopy (TEM) combined with selected area electron diffraction (SAED), scanning electron microscopy (SEM), BET specific surface area analysis, X-ray photoelectron spectroscopy (XPS) and thermo gravimetric analysis (TGA). Hot-stage Raman and infrared emission spectroscopy were applied to study the chemical reactions during dehydration and dehydroxylation. The advantage of these techniques is that the changes in molecular structure can be followed in situ and at the elevated temperatures.

Kinetics of electron recombination of dye-senitized solar cells based on TiO2 nanorod arrays sensitized with different dyes

The performance and electron recombination kinetics of dye-sensitized solar cells based on TiO2 films consisting of one-dimensional nanorod arrays (NR-DSSCs) which are sensitized with dye N719, C218 and D205 respectively have been studied. It has been found that the best efficiency is obtained with the dye C218 based NR-DSSCs, benefiting from a 40% higher short-circuit photocurrent density. However, the open circuit photovoltage of the N719 based cell is 40 mV higher than that of the organic dye C218 and D205 based devices. Investigation of the electron recombination kinetics of the NR-DSSCs has revealed that the effective electron lifetime, τn, of the N719 based NR-DSSC is the lowest whereas the τn of the C218 based NR-DSSC is the highest among the three dyes. The higher Voc with the N719 based NR-DSSC is originated from the more negative energy level of the conduction band of the TiO2 film. In addition, in comparison to the DSSCs with conventional nanocrystalline particles based TiO2 films, the NR-DSSCs have shown over two orders of magnitude higher τn when employing N719 as the sensitizer. Nevertheless, the τn of the DSSCs with the C218 based nanorod arrays is only ten-fold higher than the that of the nanoparticles based devices. The remarkable characteristic of the dye C218 in suppressing the electron recombination of DSSCs is discussed.

Structural Systematics of the Anhydrous 1:1 Proton-Transfer Compounds of 3,5-Dinitrosalicylic Acid with Aniline and Monosubstituted Anilines

The crystal structures of the proton-transfer compounds of 3,5-dinitrosalicylic acid (DNSA) with a series of aniline-type Lewis bases [aniline, 2-hydroxyaniline, 2-methoxyaniline, 3-methoxyaniline, 4-fluoroaniline, 4-chloroaniline and 2-aminoaniline] have been determined and their hydrogen-bonding systems analysed. All are anhydrous 1:1 salts: [(C6H8N)+(C7H3N2O7)-], (1), [(C6H8NO)+(C7H3N2O7)-], (2), [(C7H10NO)+(C7H3N2O7)-], (3), [(C7H10NO)+(C7H3N2O7)-], (4), [(C6H7FN)+(C7H3N2O7)-], (5), [(C6H7ClN)+(C7H3N2O7)-], (6), and [(C6H9N2)+(C7H3N2O7)-], (7) respectively. Crystals of 1 and 6 are triclinic, space group P-1 while the remainder are monoclinic with space group either P21/n (2, 4, 5 and 7) or P21 (3). Unit cell dimensions and contents are: for 1, a = 7.2027(17), b = 7.5699(17), c = 12.9615(16) Å, α = 84.464(14), β = 86.387(15), γ = 75.580(14)o, Z = 2; for 2, a = 7.407(3), b = 6.987(3), c = 27.653(11) Å, β = 94.906(7)o, Z = 4; for 3, a = 8.2816(18), b = 23.151(6), c = 3.9338(10), β = 95.255(19)o, Z = 2; for 4, a = 11.209(2), b = 8.7858(19), c = 15.171(3) Å, β = 93.717(4)o, Z = 4; for 5, a = 26.377(3), b = 10.1602(12), c = 5.1384(10) Å, β = 91.996(13)o, Z = 4; for 6, a = 11.217(3), b = 14.156(5), c = 4.860(3) Å, α = 99.10(4), β = 96.99(4), γ = 76.35(2)o, Z = 2; for 7, a = 12.830(4), b = 8.145(3), c = 14.302(4) Å, β = 102.631(6)o, Z = 4. In all compounds at least one primary linear intermolecular N+-H…O(carboxyl) hydrogen-bonding interaction is present which, together with secondary hydrogen bonding results in the formation of mostly two-dimensional network structures, exceptions being with compounds 4 and 5 (one-dimensional) and compound 6 (three-dimensional). In only two cases [compounds 1 and 4], are weak cation-anion or cation-cation π-π interactions found while weak aromatic C-H…O interactions are insignificant. The study shows that all compounds fit the previously formulated classification scheme for primary and secondary interactive modes for proton-transfer compounds of 3,5-dinitrosalicylic acid but there are some unusual variants.

Velocity-jump models with crowding effects

Velocity jump processes are discrete random walk models that have many applications including the study of biological and ecological collective motion. In particular, velocity jump models are often used to represent a type of persistent motion, known as a “run and tumble”, which is exhibited by some isolated bacteria cells. All previous velocity jump processes are non-interacting, which means that crowding effects and agent-to-agent interactions are neglected. By neglecting these agent-to-agent interactions, traditional velocity jump models are only applicable to very dilute systems. Our work is motivated by the fact that many applications in cell biology, such as wound healing, cancer invasion and development, often involve tissues that are densely packed with cells where cell-to-cell contact and crowding effects can be important. To describe these kinds of high cell density problems using a velocity jump process we introduce three different classes of crowding interactions into a one-dimensional model. Simulation data and averaging arguments lead to a suite of continuum descriptions of the interacting velocity jump processes. We show that the resulting systems of hyperbolic partial differential equations predict the mean behavior of the stochastic simulations very well.

Exploring Rhythmic Automata

The use of Cellular Automata (CA) for musical purposes has a rich history. In general the mapping of CA states to note-level music representations has focused on pitch mapping and downplayed rhythm. This paper reports experiments in the application of one-dimensional cellular automata to the generation and evolution of rhythmic patterns. A selection of CA tendencies are identified that can be used as compositional tools to control the rhythmic coherence of monophonic passages and the polyphonic texture of musical works in broad-brush, rather than precisely deterministic, ways. This will provide the composer and researcher with a clearer understanding of the useful application of CAs for generative music.

Mathematical modelling of tumour-induced angiogenesis

The growth of solid tumours beyond a critical size is dependent upon angiogenesis, the formation of new blood vessels from an existing vasculature. Tumours may remain dormant at microscopic sizes for some years before switching to a mode in which growth of a supportive vasculature is initiated. The new blood vessels supply nutrients, oxygen, and access to routes by which tumour cells may travel to other sites within the host (metastasize). In recent decades an abundance of biological research has focused on tumour-induced angiogenesis in the hope that treatments targeted at the vasculature may result in a stabilisation or regression of the disease: a tantalizing prospect. The complex and fascinating process of angiogenesis has also attracted the interest of researchers in the field of mathematical biology, a discipline that is, for mathematics, relatively new. The challenge in mathematical biology is to produce a model that captures the essential elements and critical dependencies of a biological system. Such a model may ultimately be used as a predictive tool. In this thesis we examine a number of aspects of tumour-induced angiogenesis, focusing on growth of the neovasculature external to the tumour. Firstly we present a one-dimensional continuum model of tumour-induced angiogenesis in which elements of the immune system or other tumour-cytotoxins are delivered via the newly formed vessels. This model, based on observations from experiments by Judah Folkman et al., is able to show regression of the tumour for some parameter regimes. The modelling highlights a number of interesting aspects of the process that may be characterised further in the laboratory. The next model we present examines the initiation positions of blood vessel sprouts on an existing vessel, in a two-dimensional domain. This model hypothesises that a simple feedback inhibition mechanism may be used to describe the spacing of these sprouts with the inhibitor being produced by breakdown of the existing vessel's basement membrane. Finally, we have developed a stochastic model of blood vessel growth and anastomosis in three dimensions. The model has been implemented in C++, includes an openGL interface, and uses a novel algorithm for calculating proximity of the line segments representing a growing vessel. This choice of programming language and graphics interface allows for near-simultaneous calculation and visualisation of blood vessel networks using a contemporary personal computer. In addition the visualised results may be transformed interactively, and drop-down menus facilitate changes in the parameter values. Visualisation of results is of vital importance in the communication of mathematical information to a wide audience, and we aim to incorporate this philosophy in the thesis. As biological research further uncovers the intriguing processes involved in tumourinduced angiogenesis, we conclude with a comment from mathematical biologist Jim Murray, Mathematical biology is : : : the most exciting modern application of mathematics.