Novel numerical methods for time-space fractional reaction diffusion equations in two dimensions

We consider time-space fractional reaction diffusion equations in two dimensions. This equation is obtained from the standard reaction diffusion equation by replacing the first order time derivative with the Caputo fractional derivative, and the second order space derivatives with the fractional Laplacian. Using the matrix transfer technique proposed by Ilic, Liu, Turner and Anh [Fract. Calc. Appl. Anal., 9:333--349, 2006] and the numerical solution strategy used by Yang, Turner, Liu, and Ilic [SIAM J. Scientific Computing, 33:1159--1180, 2011], the solution of the time-space fractional reaction diffusion equations in two dimensions can be written in terms of a matrix function vector product $f(A)b$ at each time step, where $A$ is an approximate matrix representation of the standard Laplacian. We use the finite volume method over unstructured triangular meshes to generate the matrix $A$, which is therefore non-symmetric. However, the standard Lanczos method for approximating $f(A)b$ requires that $A$ is symmetric. We propose a simple and novel transformation in which the standard Lanczos method is still applicable to find $f(A)b$, despite the loss of symmetry. Numerical results are presented to verify the accuracy and efficiency of our newly proposed numerical solution strategy.

Ocular characteristics in myopic anisometropia

PURPOSE: To investigate the interocular symmetry of ocular optical, biometric and biomechanical characteristics between the more and less ametropic eyes of myopic anisometropes. METHODS: Thirty-four young, healthy myopic anisometropic adults (≥ 1 D spherical equivalent difference between eyes) without amblyopia or strabismus were recruited. A range of biometric and optical parameters were measured in the more and less ametropic eye of each subject including; axial length, ocular aberrations, intraocular pressure and corneal topography, thickness and biomechanics. Morphology of the anterior eye in primary and downward gaze was examined using custom software analysis of high resolution digital images. Ocular sighting dominance was assessed using the hole-in-the-card test. RESULTS: Mean absolute spherical equivalent anisometropia was 1.74 ± 0.74 D. There was a strong correlation between the degree of anisometropia and the interocular difference in axial length (r = 0.81, p < 0.001). The more and less ametropic fellow eyes displayed a high degree of interocular symmetry for the majority of biometric, biomechanical and optical parameters measured. When the level of anisometropia exceeded 1.75 D (n = 10), the more myopic eye was the dominant sighting eye in nine of these ten subjects. Subjects with greater levels of anisometropia (> 1.75 D) also showed high levels of correlation between the dominant and non-dominant eyes in their biometric, biomechanical and optical characteristics. CONCLUSIONS: Although significantly different in axial length, anisometropic eyes display a high degree of interocular symmetry for a range of anterior eye biometrics and optical parameters. For higher levels of anisometropia, the more myopic eye tends to be the dominant sighting eye.

Ethane-1,2-diaminium 4,4'-sulfonyldibenzoate

In the title compound 2(C2H10N2+) C14H8O6S2-, both the ethylenediaminium cations and the 4,4'-sulfonyldibenzoate dianions have crystallographic two-fold rotational symmetry and are interlinked by aminium N-H...O(carboxyl) hydrogen-bonding associations giving sheets which extend along (101) and are further linked down b forming a three-dimensional structure.

The flammability of metallic materials in normal and reduced gravity

Metallic materials exposed to oxygen-enriched atmospheres – as commonly used in the medical, aerospace, aviation and numerous chemical processing industries – represent a significant fire hazard which must be addressed during design, maintenance and operation. Hence, accurate knowledge of metallic materials flammability is required. Reduced gravity (i.e. space-based) operations present additional unique concerns, where the absence of gravity must also be taken into account. The flammability of metallic materials has historically been quantified using three standardised test methods developed by NASA, ASTM and ISO. These tests typically involve the forceful (promoted) ignition of a test sample (typically a 3.2 mm diameter cylindrical rod) in pressurised oxygen. A test sample is defined as flammable when it undergoes burning that is independent of the ignition process utilised. In the standardised tests, this is indicated by the propagation of burning further than a defined amount, or „burn criterion.. The burn criterion in use at the onset of this project was arbitrarily selected, and did not accurately reflect the length a sample must burn in order to be burning independent of the ignition event and, in some cases, required complete consumption of the test sample for a metallic material to be considered flammable. It has been demonstrated that a) a metallic material.s propensity to support burning is altered by any increase in test sample temperature greater than ~250-300 oC and b) promoted ignition causes an increase in temperature of the test sample in the region closest to the igniter, a region referred to as the Heat Affected Zone (HAZ). If a test sample continues to burn past the HAZ (where the HAZ is defined as the region of the test sample above the igniter that undergoes an increase in temperature of greater than or equal to 250 oC by the end of the ignition event), it is burning independent of the igniter, and should be considered flammable. The extent of the HAZ, therefore, can be used to justify the selection of the burn criterion. A two dimensional mathematical model was developed in order to predict the extent of the HAZ created in a standard test sample by a typical igniter. The model was validated against previous theoretical and experimental work performed in collaboration with NASA, and then used to predict the extent of the HAZ for different metallic materials in several configurations. The extent of HAZ predicted varied significantly, ranging from ~2-27 mm depending on the test sample thermal properties and test conditions (i.e. pressure). The magnitude of the HAZ was found to increase with increasing thermal diffusivity, and decreasing pressure (due to slower ignition times). Based upon the findings of this work, a new burn criterion requiring 30 mm of the test sample to be consumed (from the top of the ignition promoter) was recommended and validated. This new burn criterion was subsequently included in the latest revision of the ASTM G124 and NASA 6001B international test standards that are used to evaluate metallic material flammability in oxygen. These revisions also have the added benefit of enabling the conduct of reduced gravity metallic material flammability testing in strict accordance with the ASTM G124 standard, allowing measurement and comparison of the relative flammability (i.e. Lowest Burn Pressure (LBP), Highest No-Burn Pressure (HNBP) and average Regression Rate of the Melting Interface(RRMI)) of metallic materials in normal and reduced gravity, as well as determination of the applicability of normal gravity test results to reduced gravity use environments. This is important, as currently most space-based applications will typically use normal gravity information in order to qualify systems and/or components for reduced gravity use. This is shown here to be non-conservative for metallic materials which are more flammable in reduced gravity. The flammability of two metallic materials, Inconel® 718 and 316 stainless steel (both commonly used to manufacture components for oxygen service in both terrestrial and space-based systems) was evaluated in normal and reduced gravity using the new ASTM G124-10 test standard. This allowed direct comparison of the flammability of the two metallic materials in normal gravity and reduced gravity respectively. The results of this work clearly show, for the first time, that metallic materials are more flammable in reduced gravity than in normal gravity when testing is conducted as described in the ASTM G124-10 test standard. This was shown to be the case in terms of both higher regression rates (i.e. faster consumption of the test sample – fuel), and burning at lower pressures in reduced gravity. Specifically, it was found that the LBP for 3.2 mm diameter Inconel® 718 and 316 stainless steel test samples decreased by 50% from 3.45 MPa (500 psia) in normal gravity to 1.72 MPa (250 psia) in reduced gravity for the Inconel® 718, and 25% from 3.45 MPa (500 psia) in normal gravity to 2.76 MPa (400 psia) in reduced gravity for the 316 stainless steel. The average RRMI increased by factors of 2.2 (27.2 mm/s in 2.24 MPa (325 psia) oxygen in reduced gravity compared to 12.8 mm/s in 4.48 MPa (650 psia) oxygen in normal gravity) for the Inconel® 718 and 1.6 (15.0 mm/s in 2.76 MPa (400 psia) oxygen in reduced gravity compared to 9.5 mm/s in 5.17 MPa (750 psia) oxygen in normal gravity) for the 316 stainless steel. Reasons for the increased flammability of metallic materials in reduced gravity compared to normal gravity are discussed, based upon the observations made during reduced gravity testing and previous work. Finally, the implications (for fire safety and engineering applications) of these results are presented and discussed, in particular, examining methods for mitigating the risk of a fire in reduced gravity.

Biomaterial scaffolds in cartilage-subchondral bone defects influencing the repair of autologous articular cartilage transplants

The repair of articular cartilage typically involves the repair of cartilage-subchondral bone tissue defects. Although various bioactive materials have been used to repair bone defects, how these bioactive materials in subchondral bone defects influence the repair of autologous cartilage transplant remains unclear. The aim of this study was to investigate the effects of different subchondral biomaterial scaffolds on the repair of autologous cartilage transplant in a sheep model. Cylindrical cartilage-subchondral bone defects were created in the right femoral knee joint of each sheep. The subchondral bone defects were implanted with hydroxyapatite-β-tricalcium phosphate (HA-TCP), poly lactic-glycolic acid (PLGA)-HA-TCP dual-layered composite scaffolds (PLGA/HA-TCP scaffolds), or autologous bone chips. The autologous cartilage layer was placed on top of the subchondral materials. After three months, the effect of different subchondral scaffolds on the repair of autologous cartilage transplant was systematically studied by investigating the mechanical strength, structural integration and histological responses. The results showed that the transplanted cartilage layer supported by HA-TCP scaffolds had better structural integration and higher mechanical strength than that supported by PLGA/HA-TCP scaffolds. Furthermore, HA-TCP supported cartilage showed higher expression of acid mucosubstances and glycol-amino-glycan (GAG) contents than that supported by PLGA/HA-TCP scaffolds. Our results suggested that the physicochemical properties, including the inherent mechanical strength and material chemistry of the scaffolds, play important roles in influencing the repair of autologous cartilage transplants. The study may provide useful information for the design and selection of proper subchondral biomaterials to support the repair of both subchondral bone and cartilage defects.

Ocular characteristics of anisometropia

Animal models of refractive error development have demonstrated that visual experience influences ocular growth. In a variety of species, axial anisometropia (i.e. a difference in the length of the two eyes) can be induced through unilateral occlusion, image degradation or optical manipulation. In humans, anisometropia may occur in isolation or in association with amblyopia, strabismus or unilateral pathology. Non-amblyopic myopic anisometropia represents an interesting anomaly of ocular growth, since the two eyes within one visual system have grown to different endpoints. These experiments have investigated a range of biometric, optical and mechanical properties of anisometropic eyes (with and without amblyopia) with the aim of improving our current understanding of asymmetric refractive error development. In the first experiment, the interocular symmetry in 34 non-amblyopic myopic anisometropes (31 Asian, 3 Caucasian) was examined during relaxed accommodation. A high degree of symmetry was observed between the fellow eyes for a range of optical, biometric and biomechanical measurements. When the magnitude of anisometropia exceeded 1.75 D, the more myopic eye was almost always the sighting dominant eye. Further analysis of the optical and biometric properties of the dominant and non-dominant eyes was conducted to determine any related factors but no significant interocular differences were observed with respect to best-corrected visual acuity, corneal or total ocular aberrations during relaxed accommodation. Given the high degree of symmetry observed between the fellow eyes during distance viewing in the first experiment and the strong association previously reported between near work and myopia development, the aim of the second experiment was to investigate the symmetry between the fellow eyes of the same 34 myopic anisometropes following a period of near work. Symmetrical changes in corneal and total ocular aberrations were observed following a short reading task (10 minutes, 2.5 D accommodation demand) which was attributed to the high degree of interocular symmetry for measures of anterior eye morphology, and corneal biomechanics. These changes were related to eyelid shape and position during downward gaze, but gave no clear indication of factors associated with near work that might cause asymmetric eye growth within an individual. Since the influence of near work on eye growth is likely to be most obvious during, rather than following near tasks, in the third experiment the interocular symmetry of the optical and biometric changes was examined during accommodation for 11 myopic anisometropes. The changes in anterior eye biometrics associated with accommodation were again similar between the eyes, resulting in symmetrical changes in the optical characteristics. However, the more myopic eyes exhibited slightly greater amounts of axial elongation during accommodation which may be related to the force exerted by the ciliary muscle. This small asymmetry in axial elongation we observed between the eyes may be due to interocular differences in posterior eye structure, given that the accommodative response was equal between eyes. Using ocular coherence tomography a reduced average choroidal thickness was observed in the more myopic eyes compared to the less myopic eyes of these subjects. The interocular difference in choroidal thickness was correlated with the magnitude of spherical equivalent and axial anisometropia. The symmetry in optics and biometrics between fellow eyes which have undergone significantly different visual development (i.e. anisometropic subjects with amblyopia) is also of interest with respect to refractive error development. In the final experiment the influence of altered visual experience upon corneal and ocular higher-order aberrations was investigated in 21 amblyopic subjects (8 refractive, 11 strabismic and 2 form deprivation). Significant differences in aberrations were observed between the fellow eyes, which varied according to the type of amblyopia. Refractive amblyopes displayed significantly higher levels of 4th order corneal aberrations (spherical aberration and secondary astigmatism) in the amblyopic eye compared to the fellow non-amblyopic eye. Strabismic amblyopes exhibited significantly higher levels of trefoil, a third order aberration, in the amblyopic eye for both corneal and total ocular aberrations. The results of this experiment suggest that asymmetric visual experience during development is associated with asymmetries in higher-order aberrations, proportional to the magnitude of anisometropia and dependent upon the amblyogenic factor. This suggests a direct link between the development of higher-order optical characteristics of the human eye and visual feedback. The results from these experiments have shown that a high degree of symmetry exists between the fellow eyes of non-amblyopic myopic anisometropes for a range of biomechanical, biometric and optical parameters for different levels of accommodation and following near work. While a single specific optical or biomechanical factor that is consistently associated with asymmetric refractive error development has not been identified, the findings from these studies suggest that further research into the association between ocular dominance, choroidal thickness and higher-order aberrations with anisometropia may improve our understanding of refractive error development.

Effect of geometry on the melting rates of iron rods burning in high pressure oxygen

The effect of sample geometry on the melting rates of burning iron rods was assessed. Promoted-ignition tests were conducted with rods having cylindrical, rectangular, and triangular cross-sectional shapes over a range of cross-sectional areas. The regression rate of the melting interface (RRMI) was assessed using a statistical approach which enabled the quantification of confidence levels for the observed differences in RRMI. Statistically significant differences in RRMI were observed for rods with the same cross-sectional area but different cross-sectional shape. The magnitude of the proportional difference in RRMI increased with the cross-sectional area. Triangular rods had the highest RRMI, followed by rectangular rods, and then cylindrical rods. The dependence of RRMI on rod shape is shown to relate to the action of molten metal at corners. The corners of the rectangular and triangular rods melted faster than the faces due to their locally higher surface area to volume ratios. This phenomenon altered the attachment geometry between liquid and solid phases, increasing the surface area available for heat transfer, causing faster melting. Findings relating to the application of standard flammability test results in industrial situations are also presented.

Single crystal raman spectroscopy of natural finnemanite and comparison with the synthesised analogue

The single crystal Raman spectra of natural mineral finnemanite Pb5(AsO3)3Cl from the Långban locality, Filipstad district, Värmland province, Sweden are presented for the first time. It is a hexagonal mineral belonging to the ortho arsenite group, where the [AsO3]3- ion is isolated. The spectra of finnemanite are characterized by a strong band at 734 cm-1 overlying a shoulder at 726 cm-1, and broad overlapping bands in the lower wavenumber with the strongest band positioned at 174 cm-1. Band assignments were made based on band symmetry, experimental band positions from literature and DFT calculated Raman spectrum, and spectral comparison with other ortho arsenite minerals reinerite, cafarsite, and nealite and synthetic lead arsenite compounds Pb2(AsO2)3Cl, Pb2As2O5, and PbAs2O4 . The band at 734 cm-1 was assigned to υ1(AsO3), bands at 726 and 640 cm-1 assigned to υ3, 372 and 357 cm-1 to υ2, and 244, 239 and 207 cm-1 to υ4. The single crystal spectra of finnemanite showed good mode separation, allowing bands to be assigned a symmetry species of Ag, E1g, or E2g.

Single crystal raman spectroscopy of natural paulmooreite Pb2As2O5 and in comparison with the synthesised analogue

The single crystal Raman spectra of natural mineral paulmooreite Pb2As2O5 from the Långban locality, Filipstad district, Värmland province, Sweden are presented for the first time. It is a monoclinic mineral containing an isolated [As2O5]4-. Depolarised and single crystal spectra of the natural and synthetic sample compare favorably and are characterized by strong bands around 186 and 140 cm-1 and three medium bands at 800 – 700 cm-1. Band assignments were made based on band symmetry and spectral comparison between experimental band positions and those resulting from Hartree-Fock calculation of an isolated [As2O5]4- ion. Spectral comparison was also made with lead arsenites such as synthetic PbAs2O4 and Pb2(AsO2)3Cl and natural finnemanite in order to determine the contribution of the terminal and bridging O in paulmooreite. Bands at 760 – 733 cm-1 were assigned to terminal As-O vibrations, whereas stretches of the bridging O occur at 562 and 503 cm-1. The single crystal spectra showed good mode separation, allowing bands to be assigned a symmetry species of Ag or Bg.

Monte Carlo verification of gel dosimetry measurements for stereotactic radiotherapy

The quality assurance of stereotactic radiotherapy and radiosurgery treatments requires the use of small-field dose measurements that can be experimentally challenging. This study used Monte Carlo simulations to establish that PAGAT dosimetry gel can be used to provide accurate, high resolution, three-dimensional dose measurements of stereotactic radiotherapy fields. A small cylindrical container (4 cm height, 4.2 cm diameter) was filled with PAGAT gel, placed in the parietal region inside a CIRS head phantom, and irradiated with a 12 field stereotactic radiotherapy plan. The resulting three-dimensional dose measurement was read out using an optical CT scanner and compared with the treatment planning prediction of the dose delivered to the gel during the treatment. A BEAMnrc DOSXYZnrc simulation of this treatment was completed, to provide a standard against which the accuracy of the gel measurement could be gauged. The three dimensional dose distributions obtained from Monte Carlo and from the gel measurement were found to be in better agreement with each other than with the dose distribution provided by the treatment planning system's pencil beam calculation. Both sets of data showed close agreement with the treatment planning system's dose distribution through the centre of the irradiated volume and substantial disagreement with the treatment planning system at the penumbrae. The Monte Carlo calculations and gel measurements both indicated that the treated volume was up to 3 mm narrower, with steeper penumbrae and more variable out-of-field dose, than predicted by the treatment planning system. The Monte Carlo simulations allowed the accuracy of the PAGAT gel dosimeter to be verified in this case, allowing PAGAT gel to be utilised in the measurement of dose from stereotactic and other radiotherapy treatments, with greater confidence in the future.

Determination of stress-strain characteristics of railhead steel using image analysis

True stress-strain curve of railhead steel is required to investigate the behaviour of railhead under wheel loading through elasto-plastic Finite Element (FE) analysis. To reduce the rate of wear, the railhead material is hardened through annealing and quenching. The Australian standard rail sections are not fully hardened and hence suffer from non-uniform distribution of the material property; usage of average properties in the FE modelling can potentially induce error in the predicted plastic strains. Coupons obtained at varying depths of the railhead were, therefore, tested under axial tension and the strains were measured using strain gauges as well as an image analysis technique, known as the Particle Image Velocimetry (PIV). The head hardened steel exhibit existence of three distinct zones of yield strength; the yield strength as the ratio of the average yield strength provided in the standard (σyr=780MPa) and the corresponding depth as the ratio of the head hardened zone along the axis of symmetry are as follows: (1.17 σyr, 20%), (1.06 σyr, 20%- 80%) and (0.71 σyr, > 80%). The stress-strain curves exhibit limited plastic zone with fracture occurring at strain less than 0.1.

Raman spectroscopy of the multianion mineral gartrellite-PbCu(Fe3+,Cu)(AsO4)2(OH,H2O)2

The multianion mineral gartrellite PbCu(Fe3+,Cu)(AsO4)2(OH,H2O)2 has been studied by a combination of Raman and infrared spectroscopy. The vibrational spectra of two gartrellite samples from Durango and Ashburton Downs were compared. Gartrellite is one of the tsumcorite mineral group based upon arsenate and sulphate anions. Crystal symmetry is either triclinic in the case of an ordered occupation of two cationic sites, triclinic due to ordering of the H bonds in the case of species with 2 water molecules per formula unit, or monoclinic in the other cases. Characteristic Raman spectra of the minerals enable the assignment of the bands to specific vibrational modes. These spectra are related to the structure of gartrellite. The position of the hydroxyl and water stretching vibrations are related to the strength of the hydrogen bond formed between the OH unit and the AsO4 anion.

Binocular summation improves performance to defocus-induced blur

PURPOSE. To assess whether there are any advantages of binocular over monocular vision under blur conditions. METHODS. We measured the effect of defocus, induced by positive lenses, on the pattern reversal Visual Evoked Potential (VEP) and on visual acuity (VA). Monocular (dominant eye) and binocular VEPs were recorded from thirteen volunteers (average age: 28±5 years, average spherical equivalent: -0.25±0.73 D) for defocus up to 2.00 D using positive powered lenses. VEPs were elicited using reversing 10 arcmin checks at a rate of 4 reversals/second. The stimulus subtended a circular field of 7 degrees with 100% contrast and mean luminance 30 cd/m2. VA was measured under the same conditions using ETDRS charts. All measurements were performed at 1m viewing distance with best spectacle sphero-cylindrical correction and natural pupils. RESULTS. With binocular stimulation, amplitudes and implicit times of the P100 component of the VEPs were greater and shorter, respectively, in all cases than for monocular stimulation. Mean binocular enhancement ratio in the P100 amplitude was 2.1 in-focus, increasing linearly with defocus to be 3.1 at +2.00 D defocus. Mean peak latency was 2.9 ms shorter in-focus with binocular than for monocular stimulation, with the difference increasing with defocus to 8.8 ms at +2.00 D. As for the VEP amplitude, VA was always better with binocular than with monocular vision, with the difference being greater for higher retinal blur. CONCLUSIONS. Both subjective and electrophysiological results show that binocular vision ameliorates the effect of defocus. The increased binocular facilitation observed with retinal blur may be due to the activation of a larger population of neurons at close-to-threshold detection under binocular stimulation.

Analysis of conical-cylinder shells fundamental characteristics for structural health diagnostics

This paper describes the formulation for the free vibration of joined conical-cylindrical shells with uniform thickness using the transfer of influence coefficient for identification of structural characteristics. These characteristics are importance for structural health monitoring to develop model. This method was developed based on successive transmission of dynamic influence coefficients, which were defined as the relationships between the displacement and the force vectors at arbitrary nodal circles of the system. The two edges of the shell having arbitrary boundary conditions are supported by several elastic springs with meridional/axial, circumferential, radial and rotational stiffness, respectively. The governing equations of vibration of a conical shell, including a cylindrical shell, are written as a coupled set of first order differential equations by using the transfer matrix of the shell. Once the transfer matrix of a single component has been determined, the entire structure matrix is obtained by the product of each component matrix and the joining matrix. The natural frequencies and the modes of vibration were calculated numerically for joined conical-cylindrical shells. The validity of the present method is demonstrated through simple numerical examples, and through comparison with the results of previous researchers.

A vibrational spectroscopic study of the phosphate mineral Wardite NaAl3(PO4)2(OH)4•2(H2O)

Three wardite mineral samples from different origins have been analysed by vibrational spectroscopy. The mineral is unusual in that it belongs to a unique symmetry class, namely the tetragonal-trapezohedral group. The structure of wardite contains layers of corner-linked –OH bridged MO6 octahedra stacked along the tetragonal C-axis in a four-layer sequence and linked by PO4 groups. Consequentially not all phosphate units are identical. Thus, two intense Raman bands observed at 995 and 1051 cm-1 are assigned to the ν1 PO43- symmetric stretching mode. Intense Raman bands are observed at 605 and 618 cm-1 with shoulders at 578 and 589 cm-1 are assigned to the ν4 out of plane bending modes of the PO43-. The observation of multiple bands supports the concept of non-equivalent phosphate units in the structure. Sharp infrared bands are observed at 3544 and 3611 cm-1 are attributed to the OH stretching vibrations of the hydroxyl units. Vibrational spectroscopy enables subtle details of the molecular structure of wardite to be determined.