Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

Bernstein instability driven by thermal ring distribution

The classic Bernstein waves may be intimately related to banded emissions detected in laboratory plasmas, terrestrial, and other planetary magnetospheres. However, the customary discussion of the Bernstein wave is based upon isotropic thermal velocity distribution function. In order to understand how such waves can be excited, one needs an emission mechanism, i.e., an instability. In non-relativistic collision-less plasmas, the only known Bernstein wave instability is that associated with a cold perpendicular velocity ring distribution function. However, cold ring distribution is highly idealized. The present Brief Communication generalizes the cold ring distribution model to include thermal spread, so that the Bernstein-ring instability is described by a more realistic electron distribution function, with which the stabilization by thermal spread associated with the ring distribution is demonstrated. The present findings imply that the excitation of Bernstein waves requires a sufficiently high perpendicular velocity gradient associated with the electron distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dust charging processes with a Cairns-Tsallis distribution function with negative ions

Dust grain charging processes are presented in a non-Maxwellian dusty plasma following the Cairns-Tsallis (q, α)–distribution, whose constituents are the electrons, as well as the positive/negative ions and negatively charged dust grains. For this purpose, we have solved the current balance equation for a negatively charged dust grain to achieve an equilibrium state value (viz., q{sub d} = constant) in the presence of Cairns-Tsallis (q, α)–distribution. In fact, the current balance equation becomes modified due to the Boltzmannian/streaming distributed negative ions. It is numerically found that the relevant plasma parameters, such as the spectral indexes q and α, the positive ion-to-electron temperature ratio, and the negative ion streaming speed (U{sub 0}) significantly affect the dust grain surface potential. It is also shown that in the limit q → 1 the Cairns-Tsallis reduces to the Cairns distribution; for α = 0 the Cairns-Tsallis distribution reduces to pure Tsallis distribution and the latter reduces to Maxwellian distribution for q → 1 and α = 0.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model-centric distribution automation: Capacity, reliability, and efficiency

A series of analyses along with field validations that evaluate efficiency, reliability, and capacity improvements of model-centric distribution automation are presented. With model-centric distribution automation, the same model is used from design to real-time control calculations. A 14-feeder system with 7 substations is considered. The analyses involve hourly time-varying loads and annual load growth factors. Phase balancing and capacitor redesign modifications are used to better prepare the system for distribution automation, where the designs are performed considering time-varying loads. Coordinated control of load tap changing transformers, line regulators, and switched capacitor banks is considered. In evaluating distribution automation versus traditional system design and operation, quasi-steady-state power flow analysis is used. In evaluating distribution automation performance for substation transformer failures, reconfiguration for restoration analysis is performed. In evaluating distribution automation for storm conditions, Monte Carlo simulations coupled with reconfiguration for restoration calculations are used. As a result, the evaluations demonstrate that model-centric distribution automation has positive effects on system efficiency, capacity, and reliability.

25 ENERGY STORAGE↗

THE INFLUENCE OF PICK-UP IONS ON SPACE PLASMA DISTRIBUTIONS

The addition and incorporation of highly ordered distributions of pickup ions can increase the ordering of space plasmas, decreasing their entropy, and driving them away from equilibrium. In this study, we present a model for describing the entropy decrease of heliospheric plasmas caused by the pick-up ions and their highly ordered arrangement. We consider a plasma flow consisting of two proton populations: (1) a solar wind distribution that is in quasi-equilibrium and (2) a spherical shell distribution of pick-up protons. Because all protons are indistinguishable, they equally share a hybrid state, where the hybrid probability distribution is given by the normalized sum of the distributions of the solar wind and pick-up protons. We derive an analytical formulation of the thickness of the pick-up proton shell distribution. Then we show that the entropy of the hybrid state depends on a dimensionless fluctuation number, {Theta}, which is a measure of the organizing role of pick-up protons. This number characterizes and compares two types of fluctuations, the thermal motion of solar wind protons and the bounded motion of pick-up protons that forms the spherical shell distribution. We find that the entropy variation between the hybrid and original states is proportional to the negative logarithm of {Theta}. When the two competing fluctuations balance each other, {Theta} {approx} 1, the pick-up protons have no effect on the entropy of solar wind protons. Remarkably, while sunward of the termination shock pick-up ions act to slightly increase entropy, throughout the inner heliosheath their effect is to decrease the entropy dramatically.

79 ASTRONOMY AND ASTROPHYSICS↗

PICK-UP ION DISTRIBUTIONS AND THEIR INFLUENCE ON ENERGETIC NEUTRAL ATOM SPECTRAL CURVATURE

This paper focuses on the analysis and significance of the spectral curvature of energetic neutral atoms (ENAs) detected by the Interstellar Boundary Explorer. The flux versus energy spectrum is analytically expressed in terms of the source proton distributions, namely: (1) the solar wind kappa distribution of protons and (2) the coexisting filled spherical shell distribution of pick-up ions (PUIs). The influence of PUIs on the spectral index and curvature is modeled and investigated in detail. It is analytically shown that (1) the PUI speed upper limit is restricted by the Earthward PUI velocity vector, (2) the PUI distribution causes a positive spectral curvature, and (3) the exact expressions of the spectral index and curvature can be used to extract information about the governing parameters of the parent proton distributions. The sky maps of the spectral curvature reveal a possible band-like configuration of positive spectral curvature that is missing in the original flux sky maps. This band can be roughly separated into the north/south polar regions and two ecliptic meridional 'columns' located around the ecliptic longitudes {approx}5 Degree-Sign and {approx}150 Degree-Sign . The geometric locus between the two cones with noseward axis, and apertures {approx}60 Degree-Sign and {approx}120 Degree-Sign , configures the band-like region of (1) the positive curvature and (2) the maximum values of PUI distribution. Indeed, the observed curvature band is highly correlated with PUI distributions, and is possibly caused by the influence of PUIs on bending the spectrum from linear (log-log scale) to concave upward, thus increasing its spectral curvature.

79 ASTRONOMY AND ASTROPHYSICS↗

THE DENSITY DISTRIBUTION IN TURBULENT BISTABLE FLOWS

We numerically study the volume density probability distribution function (n-PDF) and the column density probability distribution function ({Sigma}-PDF) resulting from thermally bistable turbulent flows. We analyze three-dimensional hydrodynamic models in periodic boxes of 100 pc by side, where turbulence is driven in the Fourier space at a wavenumber corresponding to 50 pc. At low densities (n {approx}< 0.6 cm{sup -3}), the n-PDF is well described by a lognormal distribution for an average local Mach number ranging from {approx}0.2 to {approx}5.5. As a consequence of the nonlinear development of thermal instability (TI), the logarithmic variance of the distribution of the diffuse gas increases with M faster than in the well-known isothermal case. The average local Mach number for the dense gas (n {approx}> 7.1 cm{sup -3}) goes from {approx}1.1 to {approx}16.9 and the shape of the high-density zone of the n-PDF changes from a power law at low Mach numbers to a lognormal at high M values. In the latter case, the width of the distribution is smaller than in the isothermal case and grows slower with M. At high column densities, the {Sigma}-PDF is well described by a lognormal for all of the Mach numbers we consider and, due to the presence of TI, the width of the distribution is systematically larger than in the isothermal case but follows a qualitatively similar behavior as M increases. Although a relationship between the width of the distribution and M can be found for each one of the cases mentioned above, these relations are different from those of the isothermal case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling the distribution of Mg II absorbers around galaxies using background galaxies and quasars

We present joint constraints on the distribution of Mg II absorption around high redshift galaxies obtained by combining two orthogonal probes, the integrated Mg II absorption seen in stacked background galaxy spectra and the distribution of parent galaxies of individual strong Mg II systems as seen in the spectra of background quasars. We present a suite of models that can be used to predict, for different two- and three-dimensional distributions, how the projected Mg II absorption will depend on a galaxy's apparent inclination, the impact parameter b and the azimuthal angle between the projected vector to the line of sight and the projected minor axis. In general, we find that variations in the absorption strength with azimuthal angles provide much stronger constraints on the intrinsic geometry of the Mg II absorption than the dependence on the inclination of the galaxies. In addition to the clear azimuthal dependence in the integrated Mg II absorption that we reported earlier in Bordoloi et al., we show that strong equivalent width Mg II absorbers (W{sub r} (2796) ≥ 0.3 Å) are also asymmetrically distributed in azimuth around their host galaxies: 72% of the absorbers in Kacprzak et al., and 100% of the close-in absorbers within 35 kpc of the center of their host galaxies, are located within 50° of the host galaxy's projected semi minor axis. It is shown that either composite models consisting of a simple bipolar component plus a spherical or disk component, or a single highly softened bipolar distribution, can well represent the azimuthal dependencies observed in both the stacked spectrum and quasar absorption-line data sets within 40 kpc. Simultaneously fitting both data sets, we find that in the composite model the bipolar cone has an opening angle of ∼100° (i.e., confined to within 50° of the disk axis) and contains about two-thirds of the total Mg II absorption in the system. The single softened cone model has an exponential fall off with azimuthal angle with an exponential scale length in opening angle of about 45°. We conclude that the distribution of Mg II gas at low impact parameters is not the same as that found at high impact parameters. At larger impact parameters beyond 40 kpc, there is evidence for a much more symmetric distribution, significantly different from that closer in to the galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-parametric projections of national income distribution consistent with the Shared Socioeconomic Pathways

Understanding and projecting income distributions within countries and regions is important to understanding consumption trends and the distributional consequences of climate impacts and responses. Several global, country-level projections of income distribution are available but most project only the Gini coefficient (a summary statistic of the distribution) or utilize the Gini along with the assumption of a lognormal distribution. We test the lognormal assumption and find that it typically underestimates income in the highest deciles and over-estimates it in others. We find that a new model based on two principal components of national time series data for income distribution provides a better fit to the data for all deciles, especially for the highest and lowest. We also construct a projection model in which the first principal component is driven by the Gini coefficient and the second captures deviations from this relationship. We use the model to project income distribution by decile for all countries for the five shared socioeconomic pathways. We find that inequality is consistently higher than projections based on the Gini and the lognormal functional form, with some countries reaching ratios of the highest to lowest income deciles that are almost three times their value using the lognormal assumption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Narrowing Frequency Probability Density Function for Achieving Minimized Uncertainties in Power Systems Operation – a Stochastic Distribution Control Perspective

With the increased penetration of renewables and other distributed energy resources (DERs), the frequency response in power systems exhibits stochastic nature, where the standard deterministic swing equation should be either represented or transferred into an Ito stochastic differential equation whose solution should be the non-stationary transient behavior of the probability density function (PDF) of the frequency. This transfers the frequency quality control for power systems into a stochastic distribution control problem, where the shape of the frequency distribution, namely the frequency PDF, needs to be directly controlled. In this context, this requires the development of control strategies via controlling the power balance so that the PDF of the frequency is made as narrow as possible whilst centered at its required mean value simple because such a sharp distribution shape means less uncertainties or randomness for the frequency response. In this paper, the formulation of the stochastic swing equation will be made first taking into account of DERs. This will then be followed by the development of stochastic distribution control model that links the power sources with the PDF of the frequency using Fokker Planck Kolmogorov (FPK) equations. A generic constrained optimization problem will be formulated where the cost function is composed of a functional distance between the actual and the desired PDFs of the frequency. A feasible solution using B-spine Neural Networks based stochastic distribution control model will be described. Using the obtained stochastic distribution control model, a set of control algorithms will be described that uses available power sources to shape the PDF of the frequency or to minimize the randomness of the frequency via minimized entropy approach. Future directions will be discussed in the later part of the paper

Swing equation, Power Systems, stchastic distribut↗

Enhancing Distribution Grid Resilience Through Model Predictive Controller Enabled Prioritized Load Restoration Strategy

Effective resilience improvement strategies enable the power grid to cope with disruptive extreme events. Most power grid outages are caused by disruptions in distribution grids. Motivated by the urgent need for power system resilience research, this paper proposes a priority-weighted optimal load restoration technique to enhance the resilience of distribution grids against extreme events. The proposed technique is based on smart distribution technology and framed as sequential multi-step decision process (MDP) and mixed integer linear program (MILP). It is formulated as optimal control problem with a model predictive control (MPC) approach. We applied the devised MILP-MPC-based load restoration technique to a simplified single-bus version of the IEEE 13-bus distribution system with integrated distributed energy resources (DERs) such as wind turbine, photovoitaic array, microturbine, and energy storage device. The technique executes a reducing and rolling horizon optimization in each control step in real-time using the forecasted information of the renewables, the fuel status of the microturbine and the state of charge of the energy storage device. We consider an extreme event which triggered outage of the upstream utility grid and caused islanded operation of the distribution grid. We demonstrated the effectiveness of the proposed MPC approach in restoring the distribution grid loads based on their priority during the main grid outage-caused islanded operation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Sequential Mitigation Solutions to Enable Distributed PV Grid Integration

This paper presents a method to assess the technical impacts of a series of distribution system mitigations for integrating increasing penetrations of distributed generation from solar photovoltaics (DGPV). Solutions considered include distribution grid upgrades and advanced control technologies. This study helps pave the way to estimate the distribution integration costs needed to reach various PV penetration levels. The approach uses three stream-lined, bounding scenarios for increased PV penetration patterns. It begins by studying the distribution impact of PV power and quantifying the PV hosting capacity. Then, it studies multiple mitigation solutions to increase PV penetration on distribution feeders. These mitigation solutions are implemented in a sequential manner according to their typical costs. The proposed technique is demonstrated on two representative distribution feeders.

14 SOLAR ENERGY↗

Synthetic, Realistic Transmission and Distribution Co-Simulation for Voltage Control Benchmarking

As distributed energy resources (DERs) play an increasing role in generation and grid services, control and market paradigms will need to be tested with models that incorporate transmission and distribution systems. The test systems will need to be realistic and accessible for researchers to form conclusions about realistic behaviors and maintain comparability and repeatability of results. These test systems also must have high resolution models of sufficient size for realistic distribution and transmission interaction. This paper demonstrates the tightly coupled HELICS-based co-simulation of realistic synthetic distribution and transmission models for the Austin metropolitan area and how this framework can be used to analyze resultant behaviors from distributed resources and their controls. These results are compared to the same distribution system run without co-simulation using a fixed feed-in voltage. Results show that even without advanced voltage controls, the feed-in voltage response when coupled to the transmission system results a lower level of voltage excursions, suggesting the importance of using this type of realistic, tightly coupled, transmission and distribution cosimulation when assessing DERs.

27 ARPA - Advanced Research Projects Agency-Energy↗

Building Large-Scale U.S. Synthetic Electric Distribution System Models

Rapid increases in distributed energy resources on distribution systems are prompting research efforts to improve and evaluate electric power distribution algorithms; however, there is a shortage of realistic, large-scale, U.S.-style test systems for the evaluation of such advanced algorithms. Some available tools to build large-scale test systems are of European style, and their application to studies focused on the United States might not be desirable given very different characteristics between the European and U.S. distribution designs. Motivated by this need, this paper develops detailed algorithms to build large-scale U.S. distribution systems and incorporates them in a new Reference Network Model, RNM-US. The approach starts with information from street maps and a catalog with electric equipment that includes power lines, transformers, voltage regulators, capacitors, and switching devices. The paper presents the algorithms through an illustrative case study of the different steps that comprise the process of building a synthetic distribution grid. Finally, the paper presents a medium- and a large-scale data set covering 10 million electrical nodes and 120,000 km of power lines, demonstrating the applicability of the proposed method to build very large-scale synthetic distribution systems.

27 ARPA - Advanced Research Projects Agency-Energy↗

Fast Tuning-Free Distributed Algorithm for Solving the Network-Constrained Economic Dispatch

With the increasing penetration of distributed energy resources (DERs) and their participation in the electricity market, it becomes more desirable to apply distributed algorithms for resource allocation in order to address the resulting computational and communicational challenges. Most of the existing distributed algorithms for solving the network-constrained economic dispatch (NCED) problem require the tuning of certain auxiliary parameters. As a result, the robustness of these algorithms against the varieties in DERs is greatly undermined. In this paper, a new distributed algorithm, optimality condition consensus (OCC), is proposed to solve the NCED problem by using distributed power flow (DPF) and ratio consensus as fundamental tools. It inherits the advantages of existing distributed algorithms for the NCED problem but removes the need for parameter tuning to improve performance in practice. In conclusion, the effectiveness of the proposed distributed algorithm in terms of efficiency, scalability, and robustness is demonstrated through detailed case studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Maximizing the biological effect of proton dose delivered with scanned beams via inhomogeneous daily dose distributions

Purpose: Biological effect of radiation can be enhanced with hypofractionation, localized dose escalation, and, in particle therapy, with optimized distribution of linear energy transfer (LET). The authors describe a method to construct inhomogeneous fractional dose (IFD) distributions, and evaluate the potential gain in the therapeutic effect from their delivery in proton therapy delivered by pencil beam scanning. Methods: For 13 cases of prostate cancer, the authors considered hypofractionated courses of 60 Gy delivered in 20 fractions. (All doses denoted in Gy include the proton's mean relative biological effectiveness (RBE) of 1.1.) Two types of plans were optimized using two opposed lateral beams to deliver a uniform dose of 3 Gy per fraction to the target by scanning: (1) in conventional full-target plans (FTP), each beam irradiated the entire gland, (2) in split-target plans (STP), beams irradiated only the respective proximal hemispheres (prostate split sagittally). Inverse planning yielded intensity maps, in which discrete position control points of the scanned beam (spots) were assigned optimized intensity values. FTP plans preferentially required a higher intensity of spots in the distal part of the target, while STP, by design, employed proximal spots. To evaluate the utility of IFD delivery, IFD plans were generated by rearranging the spot intensities from FTP or STP intensity maps, separately as well as combined using a variety of mixing weights. IFD courses were designed so that, in alternating fractions, one of the hemispheres of the prostate would receive a dose boost and the other receive a lower dose, while the total physical dose from the IFD course was roughly uniform across the prostate. IFD plans were normalized so that the equivalent uniform dose (EUD) of rectum and bladder did not increase, compared to the baseline FTP plan, which irradiated the prostate uniformly in every fraction. An EUD-based model was then applied to estimate tumor control probability (TCP) and normal tissue complication probability (NTCP). To assess potential local RBE variations, LET distributions were calculated with Monte Carlo, and compared for different plans. The results were assessed in terms of their sensitivity to uncertainties in model parameters and delivery. Results: IFD courses included equal number of fractions boosting either hemisphere, thus, the combined physical dose was close to uniform throughout the prostate. However, for the entire course, the prostate EUD in IFD was higher than in conventional FTP by up to 14%, corresponding to the estimated increase in TCP to 96% from 88%. The extent of gain depended on the mixing factor, i.e., relative weights used to combine FTP and STP spot weights. Increased weighting of STP typically yielded a higher target EUD, but also led to increased sensitivity of dose to variations in the proton's range. Rectal and bladder EUD were same or lower (per normalization), and the NTCP for both remained below 1%. The LET distributions in IFD also depended strongly on the mixing weights: plans using higher weight of STP spots yielded higher LET, indicating a potentially higher local RBE. Conclusions: In proton therapy delivered by pencil beam scanning, improved therapeutic outcome can potentially be expected with delivery of IFD distributions, while administering the prescribed quasi-uniform dose to the target over the entire course. The biological effectiveness of IFD may be further enhanced by optimizing the LET distributions. IFD distributions are characterized by a dose gradient located in proximity of the prostate's midplane, thus, the fidelity of delivery would depend crucially on the precision with which the proton range could be controlled.

60 APPLIED LIFE SCIENCES↗

Preliminary investigations on the determination of three-dimensional dose distributions using scintillator blocks and optical tomography

Purpose: Clinical QA in teletherapy as well as the characterization of experimental radiation sources for future medical applications requires effective methods for measuring three-dimensional (3D) dose distributions generated in a water-equivalent medium. Current dosimeters based on ionization chambers, diodes, thermoluminescence detectors, radiochromic films, or polymer gels exhibit various drawbacks: High quality 3D dose determination is either very sophisticated and expensive or requires high amounts of effort and time for the preparation or read out. New detectors based on scintillator blocks in combination with optical tomography are studied, since they have the potential to facilitate the desired cost-effective, transportable, and long-term stable dosimetry system that is able to determine 3D dose distributions with high spatial resolution in a short time.Methods: A portable detector prototype was set up based on a plastic scintillator block and four digital cameras. During irradiation the scintillator emits light, which is detected by the fixed cameras. The light distribution is then reconstructed by optical tomography, using maximum-likelihood expectation maximization. The result of the reconstruction approximates the 3D dose distribution. First performance tests of the prototype using laser light were carried out. Irradiation experiments were performed with ionizing radiation, i.e., bremsstrahlung (6 to 21 MV), electrons (6 to 21 MeV), and protons (68 MeV), provided by clinical and research accelerators.Results: Laser experiments show that the current imaging properties differ from the design specifications: The imaging scale of the optical systems is position dependent, ranging from 0.185 mm/pixel to 0.225 mm/pixel. Nevertheless, the developed dosimetry method is proven to be functional for electron and proton beams. Induced radiation doses of 50 mGy or more made 3D dose reconstructions possible. Taking the imaging properties into account, determined dose profiles are in agreement with reference measurements. An inherent drawback of the scintillator is the nonlinear light output for high stopping-power radiation due to the quenching effect. It impacts the depth dose curves measured with the dosimeter. For single Bragg peak distributions this leads to a peak to plateau ratio of 2.8 instead of 4.5 for the reference ionization chamber measurement. Furthermore, the transmission of the clinical bremsstrahlung beams through the scintillator leads to the saturation of one camera, making dose reconstructions in that case presently not feasible.Conclusions: It is shown that distributions of scintillation light generated by proton or electron beams can be reconstructed by the dosimetry system within minutes. The quenching apparent for proton irradiation, and the yet not precisely determined position dependency of the imaging scale, require further investigation and corrections. Upgrading the prototype with larger or inorganic scintillators would increase the detectable proton and electron energy range. The presented results show that the determination of 3D dose distributions using scintillator blocks and optical tomography is a promising dosimetry method.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-J-62: Estimating Plausible Treatment Course Dose Distributions by Accounting for Registration Uncertainty and Organ Motion

Purpose: Dose accumulation following deformable image registration (DIR) is challenging. In this study, we used a statistical sampling approach, which takes into account both DIR uncertainties and patient-specific organ motion, to study the distribution of possible true dose distributions. Methods: The study included ten patients (six CT scans/patient) treated with radiotherapy for prostate cancer. For each patient, the planned dose was re-calculated on the repeated geometries, following rigid registration based on fiducial markers. The dose re-calculated on the first CT served as our snapshot dose distribution (D1) and the average of the first five repeat scans as our treatment course reference dose distribution (Dref). Patient-specific motion and DIR-uncertainties, at each voxel in CT1, were assessed using a previously developed DIR performance measure, the distance discordance metric (DDM). To sample the distribution of possible true, predicted dose distributions (Dpred), we resampled D1 by perturbing the location of each voxel with the corresponding DDM. The three dose distribution approaches are compared for the rectum and the bladder. Results: The bladder generalized equivalent uniform dose (gEUD) from the averaged Dpred was closer to the gEUDref than to the gEUD1 (difference: 0.6 vs. 1.0 Gy). For both structures, the gEUDpred was higher than the gEUDref, and significantly higher (p≤0.05) for the rectum (average: 50.8 Gy vs. 48.0 Gy). Conclusion: We have shown that the bladder gEUD values resulting from our DIR-uncertainty inclusive dose sampling approach, Dpred, were closer to the gEUD from Dref than the gEUD values from D1. For the rectum, gEUDpred overestimated gEUDref. Theoretically however, gEUDpred values, sampled from DDM uncertainties are more representative of dose uncertainties.

60 APPLIED LIFE SCIENCES↗

SU-E-T-205: Improving Quality Assurance of HDR Brachytherapy: Verifying Agreement Between Planned and Delivered Dose Distributions Using DICOM RTDose and Advanced Film Dosimetry

Purpose: HDR brachytherapy is undergoing significant development, and quality assurance (QA) checks must keep pace. Current recommendations do not adequately verify delivered against planned dose distributions: This is particularly relevant for new treatment planning system (TPS) calculation algorithms (non TG-43 based), and an era of significant patient-specific plan optimisation. Full system checks are desirable in modern QA recommendations, complementary to device-centric individual tests. We present a QA system incorporating TPS calculation, dose distribution export, HDR unit performance, and dose distribution measurement. Such an approach, more common in external beam radiotherapy, has not previously been reported in the literature for brachytherapy. Methods: Our QA method was tested at 24 UK brachytherapy centres. As a novel approach, we used the TPS DICOM RTDose file export to compare planned dose distribution with that measured using Gafchromic EBT3 films placed around clinical brachytherapy treatment applicators. Gamma analysis was used to compare the dose distributions. Dose difference and distance to agreement were determined at prescription Point A. Accurate film dosimetry was achieved using a glass compression plate at scanning to ensure physically-flat films, simultaneous scanning of known dose films with measurement films, and triple-channel dosimetric analysis. Results: The mean gamma pass rate of RTDose compared to film-measured dose distributions was 98.1% at 3%(local), 2 mm criteria. The mean dose difference, measured to planned, at Point A was -0.5% for plastic treatment applicators and -2.4% for metal applicators, due to shielding not accounted for in TPS. The mean distance to agreement was 0.6 mm. Conclusion: It is recommended to develop brachytherapy QA to include full-system verification of agreement between planned and delivered dose distributions. This is a novel approach for HDR brachytherapy QA. A methodology using advanced film dosimetry and gamma comparison to DICOM RTDose files has been demonstrated as suitable to fulfil this need.

07 ISOTOPE AND RADIATION SOURCES↗