Analysis of heat transfer and fluid flow in a triangular enclosure in presence of an adiabatic fin

Natural convection thermal boundary layer adjacent to the heated inclined wall of a right angled triangle with an adiabatic fin attached to that surface is investigated by numerical simulations. The finite volume based unsteady numerical model is adopted for the simulation. It is revealed from the numerical results that the development of the boundary layer along the inclined surface is characterized by three distinct stages, i.e. a start-up stage, a transitional stage and a steady stage. These three stages can be clearly identified from the numerical simulations. Moreover, in presence of adiabatic fin, the thermal boundary layer adjacent to the inclined wall breaks initially. However, it is reattached with the downstream boundary layer next to the fin. More attention has been given to the boundary layer development near the fin area.

Preface to the proceedings of the 16th Computational Techniques and Applications Conference

CTAC2012 was the 16th biennial Computational Techniques and Applications Conference, and took place at Queensland University of Technology from 23 - 26 September, 2012. The ANZIAM Special Interest Group in Computational Techniques and Applications is responsible for the CTAC meetings, the first of which was held in 1981.

Travelling waves for a velocity–jump model of cell migration and proliferation

Cell invasion, characterised by moving fronts of cells, is an essential aspect of development, repair and disease. Typically, mathematical models of cell invasion are based on the Fisher–Kolmogorov equation. These traditional parabolic models can not be used to represent experimental measurements of individual cell velocities within the invading population since they imply that information propagates with infinite speed. To overcome this limitation we study combined cell motility and proliferation based on a velocity–jump process where information propagates with finite speed. The model treats the total population of cells as two interacting subpopulations: a subpopulation of left–moving cells, $L(x,t)$, and a subpopulation of right–moving cells, $R(x,t)$. This leads to a system of hyperbolic partial differential equations that includes a turning rate, $\Lambda \ge 0$, describing the rate at which individuals in the population change direction of movement. We present exact travelling wave solutions of the system of partial differential equations for the special case where $\Lambda = 0$ and in the limit that $\Lambda \to \infty$. For intermediate turning rates, $0 < \Lambda < \infty$, we analyse the travelling waves using the phase plane and we demonstrate a transition from smooth monotone travelling waves to smooth nonmonotone travelling waves as $\Lambda$ decreases through a critical value $\Lambda_{crit}$. We conclude by providing a qualitative comparison between the travelling wave solutions of our model and experimental observations of cell invasion. This comparison indicates that the small $\Lambda$ limit produces results that are consistent with experimental observations.

Improving ultrasound excitation systems using a flexible power supply with adjustable voltage and frequency to drive piezoelectric transducers

The ability of a piezoelectric transducer in energy conversion is rapidly expanding in several applications. Some of the industrial applications for which a high power ultrasound transducer can be used are surface cleaning, water treatment, plastic welding and food sterilization. Also, a high power ultrasound transducer plays a great role in biomedical applications such as diagnostic and therapeutic applications. An ultrasound transducer is usually applied to convert electrical energy to mechanical energy and vice versa. In some high power ultrasound system, ultrasound transducers are applied as a transmitter, as a receiver or both. As a transmitter, it converts electrical energy to mechanical energy while a receiver converts mechanical energy to electrical energy as a sensor for control system. Once a piezoelectric transducer is excited by electrical signal, piezoelectric material starts to vibrate and generates ultrasound waves. A portion of the ultrasound waves which passes through the medium will be sensed by the receiver and converted to electrical energy. To drive an ultrasound transducer, an excitation signal should be properly designed otherwise undesired signal (low quality) can deteriorate the performance of the transducer (energy conversion) and increase power consumption in the system. For instance, some portion of generated power may be delivered in unwanted frequency which is not acceptable for some applications especially for biomedical applications. To achieve better performance of the transducer, along with the quality of the excitation signal, the characteristics of the high power ultrasound transducer should be taken into consideration as well. In this regard, several simulation and experimental tests are carried out in this research to model high power ultrasound transducers and systems. During these experiments, high power ultrasound transducers are excited by several excitation signals with different amplitudes and frequencies, using a network analyser, a signal generator, a high power amplifier and a multilevel converter. Also, to analyse the behaviour of the ultrasound system, the voltage ratio of the system is measured in different tests. The voltage across transmitter is measured as an input voltage then divided by the output voltage which is measured across receiver. The results of the transducer characteristics and the ultrasound system behaviour are discussed in chapter 4 and 5 of this thesis. Each piezoelectric transducer has several resonance frequencies in which its impedance has lower magnitude as compared to non-resonance frequencies. Among these resonance frequencies, just at one of those frequencies, the magnitude of the impedance is minimum. This resonance frequency is known as the main resonance frequency of the transducer. To attain higher efficiency and deliver more power to the ultrasound system, the transducer is usually excited at the main resonance frequency. Therefore, it is important to find out this frequency and other resonance frequencies. Hereof, a frequency detection method is proposed in this research which is discussed in chapter 2. An extended electrical model of the ultrasound transducer with multiple resonance frequencies consists of several RLC legs in parallel with a capacitor. Each RLC leg represents one of the resonance frequencies of the ultrasound transducer. At resonance frequency the inductor reactance and capacitor reactance cancel out each other and the resistor of this leg represents power conversion of the system at that frequency. This concept is shown in simulation and test results presented in chapter 4. To excite a high power ultrasound transducer, a high power signal is required. Multilevel converters are usually applied to generate a high power signal but the drawback of this signal is low quality in comparison with a sinusoidal signal. In some applications like ultrasound, it is extensively important to generate a high quality signal. Several control and modulation techniques are introduced in different papers to control the output voltage of the multilevel converters. One of those techniques is harmonic elimination technique. In this technique, switching angles are chosen in such way to reduce harmonic contents in the output side. It is undeniable that increasing the number of the switching angles results in more harmonic reduction. But to have more switching angles, more output voltage levels are required which increase the number of components and cost of the converter. To improve the quality of the output voltage signal with no more components, a new harmonic elimination technique is proposed in this research. Based on this new technique, more variables (DC voltage levels and switching angles) are chosen to eliminate more low order harmonics compared to conventional harmonic elimination techniques. In conventional harmonic elimination method, DC voltage levels are same and only switching angles are calculated to eliminate harmonics. Therefore, the number of eliminated harmonic is limited by the number of switching cycles. In the proposed modulation technique, the switching angles and the DC voltage levels are calculated off-line to eliminate more harmonics. Therefore, the DC voltage levels are not equal and should be regulated. To achieve this aim, a DC/DC converter is applied to adjust the DC link voltages with several capacitors. The effect of the new harmonic elimination technique on the output quality of several single phase multilevel converters is explained in chapter 3 and 6 of this thesis. According to the electrical model of high power ultrasound transducer, this device can be modelled as parallel combinations of RLC legs with a main capacitor. The impedance diagram of the transducer in frequency domain shows it has capacitive characteristics in almost all frequencies. Therefore, using a voltage source converter to drive a high power ultrasound transducer can create significant leakage current through the transducer. It happens due to significant voltage stress (dv/dt) across the transducer. To remedy this problem, LC filters are applied in some applications. For some applications such as ultrasound, using a LC filter can deteriorate the performance of the transducer by changing its characteristics and displacing the resonance frequency of the transducer. For such a case a current source converter could be a suitable choice to overcome this problem. In this regard, a current source converter is implemented and applied to excite the high power ultrasound transducer. To control the output current and voltage, a hysteresis control and unipolar modulation are used respectively. The results of this test are explained in chapter 7.

Fractional models for the migration of biological cells in complex spatial domains

Travelling wave phenomena are observed in many biological applications. Mathematical theory of standard reaction-diffusion problems shows that simple partial differential equations exhibit travelling wave solutions with constant wavespeed and such models are used to describe, for example, waves of chemical concentrations, electrical signals, cell migration, waves of epidemics and population dynamics. However, as in the study of cell motion in complex spatial geometries, experimental data are often not consistent with constant wavespeed. Non-local spatial models have successfully been used to model anomalous diffusion and spatial heterogeneity in different physical contexts. In this paper, we develop a fractional model based on the Fisher-Kolmogoroff equation and analyse it for its wavespeed properties, attempting to relate the numerical results obtained from our simulations to experimental data describing enteric neural crest-derived cells migrating along the intact gut of mouse embryos. The model proposed essentially combines fractional and standard diffusion in different regions of the spatial domain and qualitatively reproduces the behaviour of neural crest-derived cells observed in the caecum and the hindgut of mouse embryos during in vivo experiments.

Phantom Scatter and the A-SI EPID

Introduction: The use of amorphous-silicon electronic portal imaging devices (a-Si EPIDs) for dosimetry is complicated by the effects of scattered radiation. In photon radiotherapy, primary signal at the detector can be accompanied by photons scattered from linear accelerator components, detector materials, intervening air, treatment room surfaces (floor, walls, etc) and from the patient/phantom being irradiated. Consequently, EPID measurements which presume to take scatter into account are highly sensitive to the identification of these contributions. One example of this susceptibility is the process of calibrating an EPID for use as a gauge of (radiological) thickness, where specific allowance must be made for the effect of phantom-scatter on the intensity of radiation measured through different thicknesses of phantom. This is usually done via a theoretical calculation which assumes that phantom scatter is linearly related to thickness and field-size. We have, however, undertaken a more detailed study of the scattering effects of fields of different dimensions when applied to phantoms of various thicknesses in order to derive scattered-primary ratios (SPRs) directly from simulation results. This allows us to make a more-accurate calibration of the EPID, and to qualify the appositeness of the theoretical SPR calculations. Methods: This study uses a full MC model of the entire linac-phantom-detector system simulated using EGSnrc/BEAMnrc codes. The Elekta linac and EPID are modelled according to specifications from the manufacturer and the intervening phantoms are modelled as rectilinear blocks of water or plastic, with their densities set to a range of physically realistic and unrealistic values. Transmissions through these various phantoms are calculated using the dose detected in the model EPID and used in an evaluation of the field-size-dependence of SPR, in different media, applying a method suggested for experimental systems by Swindell and Evans [1]. These results are compared firstly with SPRs calculated using the theoretical, linear relationship between SPR and irradiated volume, and secondly with SPRs evaluated from our own experimental data. An alternate evaluation of the SPR in each simulated system is also made by modifying the BEAMnrc user code READPHSP, to identify and count those particles in a given plane of the system that have undergone a scattering event. In addition to these simulations, which are designed to closely replicate the experimental setup, we also used MC models to examine the effects of varying the setup in experimentally challenging ways (changing the size of the air gap between the phantom and the EPID, changing the longitudinal position of the EPID itself). Experimental measurements used in this study were made using an Elekta Precise linear accelerator, operating at 6MV, with an Elekta iView GT a-Si EPID. Results and Discussion: 1. Comparison with theory: With the Elekta iView EPID fixed at 160 cm from the photon source, the phantoms, when positioned isocentrically, are located 41 to 55 cm from the surface of the panel. At this geometry, a close but imperfect agreement (differing by up to 5%) can be identified between the results of the simulations and the theoretical calculations. However, this agreement can be totally disrupted by shifting the phantom out of the isocentric position. Evidently, the allowance made for source-phantom-detector geometry by the theoretical expression for SPR is inadequate to describe the effect that phantom proximity can have on measurements made using an (infamously low-energy sensitive) a-Si EPID. 2. Comparison with experiment: For various square field sizes and across the range of phantom thicknesses, there is good agreement between simulation data and experimental measurements of the transmissions and the derived values of the primary intensities. However, the values of SPR obtained through these simulations and measurements seem to be much more sensitive to slight differences between the simulated and real systems, leading to difficulties in producing a simulated system which adequately replicates the experimental data. (For instance, small changes to simulated phantom density make large differences to resulting SPR.) 3. Comparison with direct calculation: By developing a method for directly counting the number scattered particles reaching the detector after passing through the various isocentric phantom thicknesses, we show that the experimental method discussed above is providing a good measure of the actual degree of scattering produced by the phantom. This calculation also permits the analysis of the scattering sources/sinks within the linac and EPID, as well as the phantom and intervening air. Conclusions: This work challenges the assumption that scatter to and within an EPID can be accounted for using a simple, linear model. Simulations discussed here are intended to contribute to a fuller understanding of the contribution of scattered radiation to the EPID images that are used in dosimetry calculations. Acknowledgements: This work is funded by the NHMRC, through a project grant, and supported by the Queensland University of Technology (QUT) and the Royal Brisbane and Women's Hospital, Brisbane, Australia. The authors are also grateful to Elekta for the provision of manufacturing specifications which permitted the detailed simulation of their linear accelerators and amorphous-silicon electronic portal imaging devices. Computational resources and services used in this work were provided by the HPC and Research Support Group, QUT, Brisbane, Australia.

Prevalence and predictors of refractive error in a genetically isolated population : the Norfolk Island Eye Study

BACKGROUND: We aimed to determine the prevalence and associations of refractive error on Norfolk Island. DESIGN: Population-based study on Norfolk Island, South Pacific. PARTICIPANTS: All permanent residents on Norfolk Island aged ≥ 15 years were invited to participate. METHODS: Patients underwent non-cycloplegic autorefraction, slit-lamp biomicroscope examination and biometry assessment. Only phakic eyes were analysed. MAIN OUTCOME MEASURES: Prevalence and multivariate associations of refractive error and myopia. RESULTS: There were 677 people (645 right phakic eyes, 648 left phakic eyes) aged ≥ 15 years were included in this study. Mean age of participants was 51.1 (standard deviation 15.7; range 15-81). Three hundred and seventy-six people (55.5%) were female. Adjusted to the 2006 Norfolk Island population, prevalence estimates of refractive error were as follows: myopia (mean spherical equivalent ≥ -1.0 D) 10.1%, hypermetropia (mean spherical equivalent ≥ 1.0 D) 36.6%, and astigmatism 17.7%. Significant independent predictors of myopia in the multivariate model were lower age (P < 0.001), longer axial length (P < 0.001), shallower anterior chamber depth (P = 0.031) and increased corneal curvature (P < 0.001). Significant independent predictors of refractive error were increasing age (P < 0.001), male gender (P = 0.009), Pitcairn ancestry (P = 0.041), cataract (P < 0.001), longer axial length (P < 0.001) and decreased corneal curvature (P < 0.001). CONCLUSIONS: The prevalence of myopia on Norfolk Island is lower than on mainland Australia, and the Norfolk Island population demonstrates ethnic differences in the prevalence estimates. Given the significant associations between refractive error and several ocular biometry characteristics, Norfolk Island may be a useful population in which to find the genetic basis of refractive error.

Ultimate load behaviour of box-columns under combined loading of axial compression and torsion

This paper describes a study of the theoretical and experimental behaviour of box-columns of varying b/t ratios under loadings of axial compression and torsion and their combinations. Details of the testing rigs and the testing methods, the results obtained such as the load-deflection curves and the interaction diagrams, and experimental observations regarding the behaviour of box-models and the types of local plastic mechanisms associated with each type of loading are presented. A simplified rigid-plastic analysis is carried out to study the collapse behaviour of box-columns under these loadings, based on the observed plastic mechanisms, and the results are compared with those of experiments.

Examination of the role of nitric oxide synthase and renal kallikrein as candidate genes for essential hypertension

Nitric oxide synthase and renal kallikrein are both involved in blood pressure regulation. Genes for these enzymes may, therefore, be considered candidates for hypertension pathogenesis. 2. In the present study, genotypes for nitric oxide synthase and renal kallikrein microsatellite markers were determined in a cross-sectional association analysis of hypertensive patients and normotensive control subjects. 3. Results from this study did not indicate an association of either of the candidate gene polymorphisms with essential hypertension. Hence, findings for this study do not support a role for these genes in human hypertension.

Acetylcholine receptor antibody in the diagnosis and management of myasthenia gravis

The relationship of acetylcholine receptor (AchR) antibodies to disease activity in myasthenia gravis (MG) is controversial. Some authors claim a direct correlation with disease activity and treatment, in particular plasmapheresis therapy, whereas others have commented on the poor overall correlation of antibody levels with clinical state. Antibody levels were examined in a population of MG patients and correlated with disease activity and response to treatment. Antibodies to skeletal muscle AchR were found in most patients with generalised MG (24/25) and in about half of the patients with purely ocular MG (6/10) and in neither of 2 patients with congenital MG. There was scant correlation with disease activity or response to treatment. It is concluded that the assay is more useful for diagnosis than for management of MG.

A linkage and cross-sectional study of hypertension and obesity using a poly(A) Alu-repeat polymorphism at the glucagon receptor gene locus (17q25)

1. Previous glucagon receptor gene (GCGR) studies have shown a Gly40Ser mutation to be more prevalent in essential hypertension and to affect glucagon binding affinity to its receptor. An Alu-repeat poly(A) polymorphism colocalized to GCGR was used in the present study to test for association and linkage in hypertension as well as association in obesity development. 2. Using a cross-sectional approach, 85 hypertensives and 95 normotensives were genotyped using polymerase chain reaction primers flanking the Alu-repeat. Both hypertensive and normotensive populations were subdivided into lean and obese categories based on body mass index (BMI) to determine involvement of this variant in obesity. For the linkage study, 89 Australian Caucasian hypertension affected sibships (174 sibpairs) were genotyped and the results were analysed using GENE-HUNTER, Mapmaker Sibs, ERPA and SPLINK (all freely available from http://linlkage.rockefeller. edu/soft/list.html). 3. Cross-sectional results for both hypertension and obesity were analysed using Chi-squared and Monte Carlo analyses. Results did not show an association of this variant with either hypertension (χ2 = 6.9, P = 0.14; Monte Carlo χ2 = 7.0, P = 0.11; n = 5000) or obesity (χ2 = 3.3, P = 0.35; Monte Carlo χ2 = 3.26, P = 0.34; n = 5000). In addition, results from the linkage study using hypertensive sib-pairs did not indicate linkage of the poly(A) repent with hypertension. Hence, results did not indicate a role far the Alu-repeat in either hypertension or obesity. However, as the heterozygosity of this poly(A) repeat is low (35%), a larger number of hypertensive sib-pairs may be required to draw definitive conclusions.

Molecular genetic analyses of RFLPs for PCR-amplified growth hormone gene, renal kallikrein gene and atrial natriuretic factor gene in essential hypertension

The present study examined polymorphisms of genes that might be involved in the onset of essential hypertension (HT). These included the (i) growth hormone gene (GH1), whose locus has recently been linked to elevated blood pressure (BP) in the stroke-prone SHR, although recent sib-pair analysis of a polymorphism near the human chorionic somatomammotropin gene (a member of the GH cluster) was unable to show linkage with HT; (ii) renal kallikrein gene (KLK1); and (iii) atrial natriuretic factor gene (ANF), where a primary defect in production or activity of kallikrein or ANF could cause NaCl retention and vasoconstriction. Association analyses were conducted to compare restriction fragment length polymorphisms (RFLPs) of each gene in 85 HT and 95 normotensive (NT) Caucasian subjects whose parents had a similar BP status at age ≥50 years. The frequency of the minor allele of (i) a RsaI RFLP in the promoter of GH1, amplified from leukocyte DNA by the polymerase chain reaction, was 0.15 in the HT group and 0.14 in the NT group (χ1=0.34, P=0.55); (ii) a TaqI RFLP for KLK1 was 0.035 in the HT group and 0.015 in the NT group (χ2=1.5, P=0.21); and (iii) a XhoI RFLP for ANF was 0.50 in HTs and 0.46 in NTs (χ2=0.20, P=0.65). Studies of HT pedigrees found one family in which the ANF locus and HT were not linked, owing to an obligate recombinant. The present data thus provide no evidence for involvement of the growth hormone, renal kallikrein, nor ANF gene in the causation of essential hypertension.

Incidence of hypotony and sympathetic ophthalmia following trans scleral cyclophotocoagulation for glaucoma and a report of risk factors

Background: To report the incidence and risk factors for hypotony and estimate the risk of sympathetic ophthalmia following diode laser trans-scleral cyclophotocoagulation (TSCPC). Design: Retrospective study using data from a private tertiary glaucoma clinic and review of the literature. Participants: Seventy eyes of 70 patients with refractory glaucoma who received TSCPC treatment. Methods: Review of the records of consecutive patients who underwent TSCPC by a single ophthalmic surgeon and review of the literature. Main Outcome Measures: Hypotony (including phthisis bulbi), sympathetic ophthalmia. Results: Seven eyes (10%; CI 5-19%) developed hypotony and included 4 eyes that developed phthisis. Higher total energy delivered during TSCPC treatment was associated with an increased risk of hypotony: eyes that developed hypotony received a mean total energy of 192.5 ± 73.2 joules, compared to a mean of 152.9 ± 83.2 joules in hypotony-free cases. The difference in mean energy delivered between the hypotony and non-hypotony group was 38.53 (95% CI: -27.57 to 104.63). The risk of sympathetic ophthalmia estimated from a review of the published literature and current series was one in 1512, or 0.07% (CI 0.03% - 0.17%). Conclusions: Total laser energy is one of several risk factors that act in a sufficient component cause-model to produce hypotony in an individual patient. The small sample size precluded inference for other individual putative risk factors but titrating laser energy may help decrease the occurrence of hypotony. The risk of sympathetic ophthalmia calculated from the literature is likely an overestimate caused by publication bias.

Cyclooxygenase-2-linked attenuation of hypoxia-induced pulmonary hypertension and intravascular thrombosis

Exogenous prostacyclin is effective in reducing pulmonary vascular resistance in some forms of human pulmonary hypertension (PH). To explore whether endogenous prostaglandins played a similar role in pulmonary hypertension, we examined the effect of deleting cyclooxygenase (COX)-gene isoforms in a chronic hypoxia model of PH. Pulmonary hypertension, examined by direct measurement of right ventricular end systolic pressure (RVESP), right ventricular hypertrophy (n = 8), and hematocrit (n = 3), was induced by 3 weeks of hypobarichypoxia in wild-type and COX-knockout (KO) mice. RVESP was increased in wild-type hypoxic mice compared with normoxic controls (24.4 ± 1.4 versus 13.8 ± 1.9 mm Hg; n = 8; p < 0.05). COX-2 KO mice showed a greater increase in RVESP following hypoxia (36.8 ± 2.7 mm Hg; p < 0.05). Urinary thromboxane (TX)B2 excretion increased following hypoxia (44.6 ± 11.1 versus 14.7 ± 1.8 ng/ml; n = 6; p < 0.05), an effect that was exacerbated by COX-2 gene disruption (54.5 ± 10.8 ng/ml; n = 6). In contrast, the increase in 6-keto-prostacyclin1α excretion following hypoxia was reduced by COX-2 gene disruption (29 ± 3 versus 52 ± 4.6 ng/ml; p < 0.01). Tail cut bleed times were lower following hypoxia, and there was evidence of intravascular thrombosis in lung vessels that was exacerbated by disruption of COX-2 and reduced by deletion of COX-1. The TXA2/endoperoxide receptor antagonist ifetroban (50 mg/kg/day) offset the effect of deleting the COX-2 gene, attenuating the hypoxia-induced rise in RVESP and intravascular thrombosis. COX-2 gene deletion exacerbates pulmonary hypertension, enhances sensitivity to TXA2, and induces intravascular thrombosis in response to hypoxia. The data provide evidence that endogenous prostaglandins modulate the pulmonary response to hypoxia. Copyright © 2008 by The American Society for Pharmacology and Experimental Therapeutics.

Experimental arts practice incorporating motion capture

My practice-led research explores and maps workflows for generating experimental creative work involving inertia based motion capture technology. Motion capture has often been used as a way to bridge animation and dance resulting in abstracted visuals outcomes. In early works this process was largely done by rotoscoping, reference footage and mechanical forms of motion capture. With the evolution of technology, optical and inertial forms of motion capture are now more accessible and able to accurately capture a larger range of complex movements. Made by Motion is a collaboration between digital artist Paul Van Opdenbosch and performer and choreographer Elise May; a series of studies on captured motion data used to generate experimental visual forms that reverberate in space and time. The project investigates the invisible forces generated by and influencing the movement of a dancer. Along with how the forces can be captured and applied to generating visual outcomes that surpass simple data visualisation, projecting the intent of the performer’s movements. The source or ‘seed’ comes from using an Xsens MVN – Inertial Motion Capture system to capture spontaneous dance movements, with the visual generation conducted through a customised dynamics simulation. In my presentation I will be displaying and discussing a selected creative works from the project along with the process and considerations behind the work.