The development of Monte Carlo techniques for the verification of radiotherapy treatments

Radiotherapy is a cancer treatment modality in which a dose of ionising radiation is delivered to a tumour. The accurate calculation of the dose to the patient is very important in the design of an effective therapeutic strategy. This study aimed to systematically examine the accuracy of the radiotherapy dose calculations performed by clinical treatment planning systems by comparison againstMonte Carlo simulations of the treatment delivery. A suite of software tools known as MCDTK (Monte Carlo DICOM ToolKit) was developed for this purpose, and is capable of: • Importing DICOM-format radiotherapy treatment plans and producing Monte Carlo simulation input files (allowing simple simulation of complex treatments), and calibrating the results; • Analysing the predicted doses of and deviations between the Monte Carlo simulation results and treatment planning system calculations in regions of interest (tumours and organs-at-risk) and generating dose-volume histograms, so that conformity with dose prescriptions can be evaluated. The code has been tested against various treatment planning systems, linear acceleratormodels and treatment complexities. Six clinical head and neck cancer treatments were simulated and the results analysed using this software. The deviations were greatest where the treatment volume encompassed tissues on both sides of an air cavity. This was likely due to the method the planning system used to model low density media.

Effect of differentially heating and consequent transient behaviour for natural convection within an enclosure of sinusoidal corrugated side walls

Unsteady natural convection due to differentially heating of the sinusoidal corrugated side walls of a modified square enclosure has been numerically investigated. The fluid inside the enclosure is air, initially as quiescent. The flat top and bottom surfaces are considered as adiabatic. The numerical scheme is based on the finite element method adapted to triangular non-uniform mesh element by a non-linear parametric solution algorithm. The results are obtained for the Rayleigh number, Ra ranging from 1e+05 to 1e+08 for different corrugation amplitude and frequency with constant physical properties for the fluid medium considered. The streamlines, isotherms and average Nusselt numbers are presented to observe the effect of sudden heating and its consequent transient behavior on fluid flow and heat transfer characteristics for the range of governing parameters. The present results show that the transient phenomena are greatly influenced by the variation of the aforementioned parameters.

An improved boundary layer scaling with ramp heating on a sloping plate

A scaling analysis for the natural convection boundary layer adjacent to an inclined semi-infinite plate subject to a non-instantaneous heating in the form of an imposed wall temperature which increases linearly up to a prescribed steady value over a prescribed time is reported. The development of the boundary layer flow from start-up to a steady-state has been described based on scaling analyses and verified by numerical simulations. The analysis reveals that, if the period of temperature growth on the wall is sufficiently long, the boundary layer reaches a quasi-steady mode before the growth of the temperature is completed. In this mode the thermal boundary layer at first grows in thickness and then contracts with increasing time. However, if the imposed wall temperature growth period is sufficiently short, the boundary layer develops differently, but after the wall temperature growth is completed, the boundary layer develops as though the startup had been instantaneous. The steady state values of the boundary layer for both cases are ultimately the same.

Monte Carlo simulations of dynamic radiotherapy treatments

The effects of tumour motion during radiation therapy delivery have been widely investigated. Motion effects have become increasingly important with the introduction of dynamic radiotherapy delivery modalities such as enhanced dynamic wedges (EDWs) and intensity modulated radiation therapy (IMRT) where a dynamically collimated radiation beam is delivered to the moving target, resulting in dose blurring and interplay effects which are a consequence of the combined tumor and beam motion. Prior to this work, reported studies on the EDW based interplay effects have been restricted to the use of experimental methods for assessing single-field non-fractionated treatments. In this work, the interplay effects have been investigated for EDW treatments. Single and multiple field treatments have been studied using experimental and Monte Carlo (MC) methods. Initially this work experimentally studies interplay effects for single-field non-fractionated EDW treatments, using radiation dosimetry systems placed on a sinusoidaly moving platform. A number of wedge angles (60º, 45º and 15º), field sizes (20 × 20, 10 × 10 and 5 × 5 cm2), amplitudes (10-40 mm in step of 10 mm) and periods (2 s, 3 s, 4.5 s and 6 s) of tumor motion are analysed (using gamma analysis) for parallel and perpendicular motions (where the tumor and jaw motions are either parallel or perpendicular to each other). For parallel motion it was found that both the amplitude and period of tumor motion affect the interplay, this becomes more prominent where the collimator tumor speeds become identical. For perpendicular motion the amplitude of tumor motion is the dominant factor where as varying the period of tumor motion has no observable effect on the dose distribution. The wedge angle results suggest that the use of a large wedge angle generates greater dose variation for both parallel and perpendicular motions. The use of small field size with a large tumor motion results in the loss of wedged dose distribution for both parallel and perpendicular motion. From these single field measurements a motion amplitude and period have been identified which show the poorest agreement between the target motion and dynamic delivery and these are used as the „worst case motion parameters.. The experimental work is then extended to multiple-field fractionated treatments. Here a number of pre-existing, multiple–field, wedged lung plans are delivered to the radiation dosimetry systems, employing the worst case motion parameters. Moreover a four field EDW lung plan (using a 4D CT data set) is delivered to the IMRT quality control phantom with dummy tumor insert over four fractions using the worst case parameters i.e. 40 mm amplitude and 6 s period values. The analysis of the film doses using gamma analysis at 3%-3mm indicate the non averaging of the interplay effects for this particular study with a gamma pass rate of 49%. To enable Monte Carlo modelling of the problem, the DYNJAWS component module (CM) of the BEAMnrc user code is validated and automated. DYNJAWS has been recently introduced to model the dynamic wedges. DYNJAWS is therefore commissioned for 6 MV and 10 MV photon energies. It is shown that this CM can accurately model the EDWs for a number of wedge angles and field sizes. The dynamic and step and shoot modes of the CM are compared for their accuracy in modelling the EDW. It is shown that dynamic mode is more accurate. An automation of the DYNJAWS specific input file has been carried out. This file specifies the probability of selection of a subfield and the respective jaw coordinates. This automation simplifies the generation of the BEAMnrc input files for DYNJAWS. The DYNJAWS commissioned model is then used to study multiple field EDW treatments using MC methods. The 4D CT data of an IMRT phantom with the dummy tumor is used to produce a set of Monte Carlo simulation phantoms, onto which the delivery of single field and multiple field EDW treatments is simulated. A number of static and motion multiple field EDW plans have been simulated. The comparison of dose volume histograms (DVHs) and gamma volume histograms (GVHs) for four field EDW treatments (where the collimator and patient motion is in the same direction) using small (15º) and large wedge angles (60º) indicates a greater mismatch between the static and motion cases for the large wedge angle. Finally, to use gel dosimetry as a validation tool, a new technique called the „zero-scan method. is developed for reading the gel dosimeters with x-ray computed tomography (CT). It has been shown that multiple scans of a gel dosimeter (in this case 360 scans) can be used to reconstruct a zero scan image. This zero scan image has a similar precision to an image obtained by averaging the CT images, without the additional dose delivered by the CT scans. In this investigation the interplay effects have been studied for single and multiple field fractionated EDW treatments using experimental and Monte Carlo methods. For using the Monte Carlo methods the DYNJAWS component module of the BEAMnrc code has been validated and automated and further used to study the interplay for multiple field EDW treatments. Zero-scan method, a new gel dosimetry readout technique has been developed for reading the gel images using x-ray CT without losing the precision and accuracy.

Heating effects in nanofocusing metal wedges

This paper investigates theoretically and numerically local heating effects in plasmon nanofocusing structures with a particular focus on the sharp free-standing metal wedges. The developed model separates plasmon propagation in the wedge from the resultant heating effects. Therefore, this model is only applicable where the temperature increments in a nanofocusing structure are sufficiently small not to result in significant variations of the metal permittivity in the wedge. The problem is reduced to a one-dimensional heating model with a distributed heat source resulting from plasmon dissipation in the metal wedge. A simple heat conduction equation governing the local heating effects in a nanofocusing structure is derived and solved numerically for plasmonic pulses of different lengths and reasonable energies. Both the possibility of achieving substantial local temperature increments in the wedge (with a significant self-influence of the heating plasmonic pulses), and the possibility of relatively weak heating (to ensure the validity of the previously developed nanofocusing theory) are demonstrated and discussed, including the future applications of the obtained results. Applicability conditions for the developed model are also derived and discussed.

Sensory imagery in craving : from cognitive psychology to new treatments for addiction

Sensory imagery is a powerful tool for inducing craving because it is a key component of the cognitive system that underpins human motivation. The role of sensory imagery in motivation is explained by Elaborated Intrusion (EI) theory. Imagery plays an important role in motivation because it conveys the emotional qualities of the desired event, mimicking anticipated pleasure or relief, and continual elaboration of the imagery ensures that the target stays in mind. We argue that craving is a conscious state, intervening between unconscious triggers and consumption, and summarise evidence that interfering with sensory imagery can weaken cravings. We argue that treatments for addiction can be enhanced by the application of EI theory to maintain motivation, and assist in the management of craving in high-risk situations.

Finite volume approach of natural convection in a triangular enclosure with localized heating from below

In this study, natural convection heat transfer and buoyancy driven flows have been investigated in a right angled triangular enclosure. The heater located on the bottom wall while the inclined wall is colder and the remaining walls are maintained as adiabatic. Governing equations of natural convection are solved through the finite volume approach, in which buoyancy is modeled via the Boussinesq approximation. Effects of different parameters such as Rayleigh number, aspect ratio, prantdl number and heater location are considered. Results show that heat transfer increases when the heater is moved toward the right corner of the enclosure. It is also revealed that increasing the Rayleigh number, increases the strength of free convection regime and consequently increases the value of heat transfer rate. Moreover, larger aspect ratio enclosure has larger Nusselt number value. In order to have better insight, streamline and isotherms are shown.

Thermo-magnetic convection of paramagnetic fluids with non-instantaneous heating

The unsteady boundary-layer development for thermomagnetic convection of paramagnetic fluids inside a square cavity has been considered in this study. The cavity is placed in a microgravity condition (no gravitation acceleration) and under a uniform magnetic field which acts vertically. A ramp temperature boundary condition is applied on left vertical side wall of the cavity where the temperature initially increases with time up to some specific time and maintain constant thereafter. A distinct magnetic convection boundary layer is developed adjacent to the left vertical wall due to the effect of the magnetic body force generated on the paramagnetic fluid. An improved scaling analysis has been performed using triple-layer integral method and verified by numerical simulations. The Prandtl number has been chosen greater than unity varied over 5-100. Moreover, the effect of various values of the magnetic parameter and magnetic Rayleigh number on the fluid flow and heat transfer has been shown.

The formation of YBa2Cu3O7-x in melt-texture heat treatments

A study of the bulk formation of YBa2Cu3O7-x from the Y2BaCuO5 plus liquid regime reveals that phase formation occurs at appreciable rates below 950°C in air. This result has been observed for phase-pure YBa2Cu3O7-x starting material given two types of heat treatment: held at 1100°C and slow-cooled from 1030°C at 6°C/h or heat-treated isothermally. Differential thermal analysis, with a cooling rate of 10°C/min indicates that the degree of undercooling for the peritectic formation of YBa2Cu3O7-x is greater than 100°C. © 1994.

A safe system-based approach to selection of clear zones, safety barriers and other roadside treatments

On average, 560 fatal run-off-road crashes occur annually in Australia and 135 in New Zealand. In addition, there are more than 14,000 run-off-road crashes causing injuries each year across both countries. In rural areas, run-off-road casualty crashes constitute 50-60% of all casualty crashes. Their severity is particularly high with more than half of those involved sustaining fatal or serious injuries. This paper reviews the existing approach to roadside hazard risk assessment, selection of clear zones and hazard treatments. It proposes a modified approach to roadside safety evaluation and management. It is a methodology based on statistical modelling of run-off-road casualty crashes, and application of locally developed crash modification factors and severity indices. Clear zones, safety barriers and other roadside design/treatment options are evaluated with a view to minimise fatal and serious injuries – the key Safe System objective. The paper concludes with a practical demonstration of the proposed approach. The paper is based on findings from a four-year Austroads research project into improving roadside safety in the Safe System context.

Enhanced photocatalytic activity of titanium oxide nanotubes after heating treatment

Titanium oxide nanotubes were obtained by an electrochemical anodization method. Scanning electron microscope results demonstrate that the diameter of the tubes is about 120 nm and the length of the tubes is around 13 μm. Transmission electron microscope results indicate that the nanotubes are assembled by numerous nanoparticles and tube-like structure remains well after heat treatment at 400-600 °C. The photocatalysis performance of the nanotubes was evaluated in terms of the decomposition rate of methyl orange under UV irradiation. The results show that the photocatalytic activity was enhanced through the heating treatment of the nanotubes, and the nanotubes heated at 600 °C exhibits the best photocatalytic activity. X-ray diffraction patterns indicate that there is no phase transformation during the heat treatment. Therefore, the enhanced activity can be attributed to the improvement of nanotubes crystallinity, which may provide more insights about the effect of the crystallinity on the photocatalytic performance.

Development of a practical design guide for balcony acoustic treatments

Balcony acoustic treatments can be demonstrated to provide important benefits in reducing road traffic noise within the balcony space and consequently internally for any adjacent room. The actual effect on road traffic noise is derived from a multitude of variables that can be broadly categorized into (a) acoustical and (b) geometrical for two distinct propagation volumes being (i) the street space, and (ii) the balcony space. A series of recent research activities in this area has incorporated the use of a combined image and diffuse source model, which can be used to predict the effect of balconies on road traffic noise for large number of scenarios. This paper investigates and presents a method and capability to summarize predictive data into user friendly guidelines aimed for use by acoustical professionals and architects and possible implementation in building design policies for environmental noise. The paper concludes with a presentation of the likely format of a potential design guide.

The significance of balcony acoustic treatments in mitigating urban road traffic noise

Urban road traffic noise in cities is an ongoing and increasing problem across much of the world. Consequently a large amount of effort is expended in attempts to address this problem, especially in the area of acoustic design of buildings. Acoustic design policies developed by government authorities will typically focus on required transport noise reductions through a building façade to meet a specified internal noise levels. The significance of balcony acoustic treatments has been highlighted in recent decades yet this area has potentially been considered less important than the need for acoustic isolation of building facades. This paper outlines recent research that has been conducted in determining the significance of balcony acoustic treatments in mitigating urban road traffic noise. It summarizes recent literature, some of which focuses on technological advances in the knowledge of balcony acoustic design and some literature discusses the overall aims and benefits of balcony acoustic design. The aim of this paper is to promote the use of balcony acoustic design as a significant element in the overall solution towards mitigating road traffic noise in modern cities.