Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LIMITING VALUES”

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 721 records · Page 40

A probabilistic graphical model based stochastic input model construction

Model reduction techniques have been widely used in modeling of high-dimensional stochastic input in uncertainty quantification tasks. However, the probabilistic modeling of random variables projected into reduced-order spaces presents a number of computational challenges. Due to the curse of dimensionality, the underlying dependence relationships between these random variables are difficult to capture. In this work, a probabilistic graphical model based approach is employed to learn the dependence by running a number of conditional independence tests using observation data. Thus a probabilistic model of the joint PDF is obtained and the PDF is factorized into a set of conditional distributions based on the dependence structure of the variables. The estimation of the joint PDF from data is then transformed to estimating conditional distributions under reduced dimensions. To improve the computational efficiency, a polynomial chaos expansion is further applied to represent the random field in terms of a set of standard random variables. This technique is combined with both linear and nonlinear model reduction methods. Numerical examples are presented to demonstrate the accuracy and efficiency of the probabilistic graphical model based stochastic input models. - Highlights: • Data-driven stochastic input models without the assumption of independence of the reduced random variables. • The problem is transformed to a Bayesian network structure learning problem. • Examples are given in flows in random media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A projection gradient method for computing ground state of spin-2 Bose–Einstein condensates

In this paper, a projection gradient method is presented for computing ground state of spin-2 Bose–Einstein condensates (BEC). We first propose the general projection gradient method for solving energy functional minimization problem under multiple constraints, in which the energy functional takes real functions as independent variables. We next extend the method to solve a similar problem, where the energy functional now takes complex functions as independent variables. We finally employ the method into finding the ground state of spin-2 BEC. The key of our method is: by constructing continuous gradient flows (CGFs), the ground state of spin-2 BEC can be computed as the steady state solution of such CGFs. We discretized the CGFs by a conservative finite difference method along with a proper way to deal with the nonlinear terms. We show that the numerical discretization is normalization and magnetization conservative and energy diminishing. Numerical results of the ground state and their energy of spin-2 BEC are reported to demonstrate the effectiveness of the numerical method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Augmented Lagrangian formulation of orbital-free density functional theory

We present an Augmented Lagrangian formulation and its real-space implementation for non-periodic Orbital-Free Density Functional Theory (OF-DFT) calculations. In particular, we rewrite the constrained minimization problem of OF-DFT as a sequence of minimization problems without any constraint, thereby making it amenable to powerful unconstrained optimization algorithms. Further, we develop a parallel implementation of this approach for the Thomas–Fermi–von Weizsacker (TFW) kinetic energy functional in the framework of higher-order finite-differences and the conjugate gradient method. With this implementation, we establish that the Augmented Lagrangian approach is highly competitive compared to the penalty and Lagrange multiplier methods. Additionally, we show that higher-order finite-differences represent a computationally efficient discretization for performing OF-DFT simulations. Overall, we demonstrate that the proposed formulation and implementation are both efficient and robust by studying selected examples, including systems consisting of thousands of atoms. We validate the accuracy of the computed energies and forces by comparing them with those obtained by existing plane-wave methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An iterative particle filter approach for coupled hydro-geophysical inversion of a controlled infiltration experiment

The modeling of unsaturated groundwater flow is affected by a high degree of uncertainty related to both measurement and model errors. Geophysical methods such as Electrical Resistivity Tomography (ERT) can provide useful indirect information on the hydrological processes occurring in the vadose zone. In this paper, we propose and test an iterated particle filter method to solve the coupled hydrogeophysical inverse problem. We focus on an infiltration test monitored by time-lapse ERT and modeled using Richards equation. The goal is to identify hydrological model parameters from ERT electrical potential measurements. Traditional uncoupled inversion relies on the solution of two sequential inverse problems, the first one applied to the ERT measurements, the second one to Richards equation. This approach does not ensure an accurate quantitative description of the physical state, typically violating mass balance. To avoid one of these two inversions and incorporate in the process more physical simulation constraints, we cast the problem within the framework of a SIR (Sequential Importance Resampling) data assimilation approach that uses a Richards equation solver to model the hydrological dynamics and a forward ERT simulator combined with Archie's law to serve as measurement model. ERT observations are then used to update the state of the system as well as to estimate the model parameters and their posterior distribution. The limitations of the traditional sequential Bayesian approach are investigated and an innovative iterative approach is proposed to estimate the model parameters with high accuracy. The numerical properties of the developed algorithm are verified on both homogeneous and heterogeneous synthetic test cases based on a real-world field experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An algebraic method for constructing stable and consistent autoregressive filters

In this paper, we introduce an algebraic method to construct stable and consistent univariate autoregressive (AR) models of low order for filtering and predicting nonlinear turbulent signals with memory depth. By stable, we refer to the classical stability condition for the AR model. By consistent, we refer to the classical consistency constraints of Adams–Bashforth methods of order-two. One attractive feature of this algebraic method is that the model parameters can be obtained without directly knowing any training data set as opposed to many standard, regression-based parameterization methods. It takes only long-time average statistics as inputs. The proposed method provides a discretization time step interval which guarantees the existence of stable and consistent AR model and simultaneously produces the parameters for the AR models. In our numerical examples with two chaotic time series with different characteristics of decaying time scales, we find that the proposed AR models produce significantly more accurate short-term predictive skill and comparable filtering skill relative to the linear regression-based AR models. These encouraging results are robust across wide ranges of discretization times, observation times, and observation noise variances. Finally, we also find that the proposed model produces an improved short-time prediction relative to the linear regression-based AR-models in forecasting a data set that characterizes the variability of the Madden–Julian Oscillation, a dominant tropical atmospheric wave pattern.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Variational Bayesian strategies for high-dimensional, stochastic design problems

This paper is concerned with a lesser-studied problem in the context of model-based, uncertainty quantification (UQ), that of optimization/design/control under uncertainty. The solution of such problems is hindered not only by the usual difficulties encountered in UQ tasks (e.g. the high computational cost of each forward simulation, the large number of random variables) but also by the need to solve a nonlinear optimization problem involving large numbers of design variables and potentially constraints. We propose a framework that is suitable for a class of such problems and is based on the idea of recasting them as probabilistic inference tasks. To that end, we propose a Variational Bayesian (VB) formulation and an iterative VB–Expectation-Maximization scheme that is capable of identifying a local maximum as well as a low-dimensional set of directions in the design space, along which, the objective exhibits the largest sensitivity. We demonstrate the validity of the proposed approach in the context of two numerical examples involving thousands of random and design variables. In all cases considered the cost of the computations in terms of calls to the forward model was of the order of 100 or less. The accuracy of the approximations provided is assessed by information-theoretic metrics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying and reducing model-form uncertainties in Reynolds-averaged Navier–Stokes simulations: A data-driven, physics-informed Bayesian approach

Despite their well-known limitations, Reynolds-Averaged Navier–Stokes (RANS) models are still the workhorse tools for turbulent flow simulations in today's engineering analysis, design and optimization. While the predictive capability of RANS models depends on many factors, for many practical flows the turbulence models are by far the largest source of uncertainty. As RANS models are used in the design and safety evaluation of many mission-critical systems such as airplanes and nuclear power plants, quantifying their model-form uncertainties has significant implications in enabling risk-informed decision-making. In this work we develop a data-driven, physics-informed Bayesian framework for quantifying model-form uncertainties in RANS simulations. Uncertainties are introduced directly to the Reynolds stresses and are represented with compact parameterization accounting for empirical prior knowledge and physical constraints (e.g., realizability, smoothness, and symmetry). An iterative ensemble Kalman method is used to assimilate the prior knowledge and observation data in a Bayesian framework, and to propagate them to posterior distributions of velocities and other Quantities of Interest (QoIs). We use two representative cases, the flow over periodic hills and the flow in a square duct, to evaluate the performance of the proposed framework. Both cases are challenging for standard RANS turbulence models. Simulation results suggest that, even with very sparse observations, the obtained posterior mean velocities and other QoIs have significantly better agreement with the benchmark data compared to the baseline results. At most locations the posterior distribution adequately captures the true model error within the developed model form uncertainty bounds. The framework is a major improvement over existing black-box, physics-neutral methods for model-form uncertainty quantification, where prior knowledge and details of the models are not exploited. This approach has potential implications in many fields in which the governing equations are well understood but the model uncertainty comes from unresolved physical processes. - Highlights: • Proposed a physics–informed framework to quantify uncertainty in RANS simulations. • Framework incorporates physical prior knowledge and observation data. • Based on a rigorous Bayesian framework yet fully utilizes physical model. • Applicable for many complex physical systems beyond turbulent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Virtual Environment for Reactor Applications (VERA): Design and architecture

VERA, the Virtual Environment for Reactor Applications, is the system of physics capabilities being developed and deployed by the Consortium for Advanced Simulation of Light Water Reactors (CASL). CASL was established for the modeling and simulation of commercial nuclear reactors. VERA consists of integrating and interfacing software together with a suite of physics components adapted and/or refactored to simulate relevant physical phenomena in a coupled manner. VERA also includes the software development environment and computational infrastructure needed for these components to be effectively used. We describe the architecture of VERA from both software and numerical perspectives, along with the goals and constraints that drove major design decisions, and their implications. We explain why VERA is an environment rather than a framework or toolkit, why these distinctions are relevant (particularly for coupled physics applications), and provide an overview of results that demonstrate the use of VERA tools for a variety of challenging applications within the nuclear industry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Vector tomography for reconstructing electric fields with non-zero divergence in bounded domains

In vector tomography (VT), the aim is to reconstruct an unknown multi-dimensional vector field using line integral data. In the case of a 2-dimensional VT, two types of line integral data are usually required. These data correspond to integration of the parallel and perpendicular projection of the vector field along the integration lines and are called the longitudinal and transverse measurements, respectively. In most cases, however, the transverse measurements cannot be physically acquired. Therefore, the VT methods are typically used to reconstruct divergence-free (or source-free) velocity and flow fields that can be reconstructed solely from the longitudinal measurements. In this paper, we show how vector fields with non-zero divergence in a bounded domain can also be reconstructed from the longitudinal measurements without the need of explicitly evaluating the transverse measurements. To the best of our knowledge, VT has not previously been used for this purpose. In particular, we study low-frequency, time-harmonic electric fields generated by dipole sources in convex bounded domains which arise, for example, in electroencephalography (EEG) source imaging. We explain in detail the theoretical background, the derivation of the electric field inverse problem and the numerical approximation of the line integrals. We show that fields with non-zero divergence can be reconstructed from the longitudinal measurements with the help of two sparsity constraints that are constructed from the transverse measurements and the vector Laplace operator. As a comparison to EEG source imaging, we note that VT does not require mathematical modeling of the sources. By numerical simulations, we show that the pattern of the electric field can be correctly estimated using VT and the location of the source activity can be determined accurately from the reconstructed magnitudes of the field. - Highlights: • Vector tomography is used to reconstruct electric fields generated by dipole sources. • Inverse solutions are based on longitudinal and transverse line integral measurements. • Transverse line integral measurements are used as a sparsity constraint. • Numerical procedure to approximate the line integrals is described in detail. • Patterns of the studied electric fields are correctly estimated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A fast mollified impulse method for biomolecular atomistic simulations

Classical integration methods for molecular dynamics are inherently limited due to resonance phenomena occurring at certain time-step sizes. The mollified impulse method can partially avoid this problem by using appropriate filters based on averaging or projection techniques. However, existing filters are computationally expensive and tedious in implementation since they require either analytical Hessians or they need to solve nonlinear systems from constraints. In this work we follow a different approach based on corotation for the construction of a new filter for (flexible) biomolecular simulations. The main advantages of the proposed filter are its excellent stability properties and ease of implementation in standard softwares without Hessians or solving constraint systems. By simulating multiple realistic examples such as peptide, protein, ice equilibrium and ice–ice friction, the new filter is shown to speed up the computations of long-range interactions by approximately 20%. The proposed filtered integrators allow step sizes as large as 10 fs while keeping the energy drift less than 1% on a 50 ps simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exactly energy conserving semi-implicit particle in cell formulation

We report a new particle in cell (PIC) method based on the semi-implicit approach. The novelty of the new method is that unlike any of its semi-implicit predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. Recent research has presented fully implicit methods where energy conservation is obtained as part of a non-linear iteration procedure. The new method (referred to as Energy Conserving Semi-Implicit Method, ECSIM), instead, does not require any non-linear iteration and its computational cycle is similar to that of explicit PIC. The properties of the new method are: i) it conserves energy exactly to round-off for any time step or grid spacing; ii) it is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency and allowing the user to select any desired time step; iii) it eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length; iv) the particle mover has a computational complexity identical to that of the explicit PIC, only the field solver has an increased computational cost. The new ECSIM is tested in a number of benchmarks where accuracy and computational performance are tested. - Highlights: • We present a new fully energy conserving semi-implicit particle in cell (PIC) method based on the implicit moment method (IMM). The new method is called Energy Conserving Implicit Moment Method (ECIMM). • The novelty of the new method is that unlike any of its predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. • The new method is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency. • The new method eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length. • These features are achieved at a reduced cost compared with either previous IMM or fully implicit implementation of PIC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The influence of alloy composition on residual stresses in heat treated aluminium alloys

The as quenched properties of eight different heat treatable aluminium alloys are related to residual stress magnitudes with the objective being to establish if there is a relationship between the residual stress and the as quenched alloy hardness and strength. Near surface residual stresses were assessed with X-ray diffraction using both the established sin{sup 2}ψ method and the more recent cos α technique. Through thickness residual stresses were also characterised using neutron diffraction. The alloys were chosen to encompass a wide range of strengths. The low to medium strength alloys were 6060 and 6082, medium to high strength 2618A, 2014A, 7075, 7010 and two variants of 7449, while the very high strength alloy was the powder metallurgy alloy N707. To assess the as quenched strength, dynamic hardness and tensile properties were determined from samples tested immediately after quenching to minimise the influence of precipitation hardening by natural aging. In addition, hot hardness measurements were made in situ on samples cooled to simulate quench paths. Within the experimental constraints of the investigation, the distribution of residual stress through the thickness was found to follow the same pattern for all the alloys investigated, varying from tensile in the interior to surface compression. The influence of alloy strength was manifested as a change in the observed residual stress magnitudes, and surface residual stresses were found to vary linearly with as quenched hardness and strength. - Highlights: • As quenched aluminium alloys contain high magnitude residual stresses. • Surface is compressive balance by a tensile core. • As quenched surface residual stress is linear function of alloy strength. • In situ hot hardness demonstrates rapid change in intrinsic hardness during rapid cooling.

36 MATERIALS SCIENCE↗

Image-guided intensity-modulated radiotherapy for prostate cancer: Dose constraints for the anterior rectal wall to minimize rectal toxicity

Rectal adverse events (AEs) are a major concern with definitive radiotherapy (RT) treatment for prostate cancer. The anterior rectal wall is at the greatest risk of injury as it lies closest to the target volume and receives the highest dose of RT. This study evaluated the absolute volume of anterior rectal wall receiving a high dose to identify potential ideal dose constraints that can minimize rectal AEs. A total of 111 consecutive patients with Stage T1c to T3a N0 M0 prostate cancer who underwent image-guided intensity-modulated RT at our institution were included. AEs were graded according to the Common Terminology Criteria for Adverse Events, version 4.0. The volume of anterior rectal wall receiving 5 to 80 Gy in 2.5-Gy increments was determined. Multivariable Cox regression models were used to identify cut points in these volumes that led to an increased risk of early and late rectal AEs. Early AEs occurred in most patients (88%); however, relatively few of them (13%) were grade ≥2. At 5 years, the cumulative incidence of late rectal AEs was 37%, with only 5% being grade ≥2. For almost all RT doses, we identified a threshold of irradiated absolute volume of anterior rectal wall above which there was at least a trend toward a significantly higher rate of AEs. Most strikingly, patients with more than 1.29, 0.73, or 0.45 cm{sup 3} of anterior rectal wall exposed to radiation doses of 67.5, 70, or 72.5 Gy, respectively, had a significantly increased risk of late AEs (relative risks [RR]: 2.18 to 2.72; p ≤ 0.041) and of grade ≥ 2 early AEs (RR: 6.36 to 6.48; p = 0.004). Our study provides evidence that definitive image-guided intensity-modulated radiotherapy (IG-IMRT) for prostate cancer is well tolerated and also identifies dose thresholds for the absolute volume of anterior rectal wall above which patients are at greater risk of early and late complications.

61 RADIATION PROTECTION AND DOSIMETRY↗

Volumetric tumor burden and its effect on brachial plexus dosimetry in head and neck intensity-modulated radiotherapy

To determine the effect of gross tumor volume of the primary (GTV-P) and nodal (GTV-N) disease on planned radiation dose to the brachial plexus (BP) in head and neck intensity-modulated radiotherapy (IMRT). Overall, 75 patients underwent definitive IMRT to a median total dose of 69.96 Gy in 33 fractions. The right BP and left BP were prospectively contoured as separate organs at risk. The GTV was related to BP dose using the unpaired t-test. Receiver operating characteristics curves were constructed to determine optimized volumetric thresholds of GTV-P and GTV-N corresponding to a maximum BP dose cutoff of > 66 Gy. Multivariate analyses were performed to account for factors associated with a higher maximal BP dose. A higher maximum BP dose (> 66 vs ≤ 66 Gy) correlated with a greater mean GTV-P (79.5 vs 30.8 cc; p = 0.001) and ipsilateral GTV-N (60.6 vs 19.8 cc; p = 0.014). When dichotomized by the optimized nodal volume, patients with an ipsilateral GTV-N ≥ 4.9 vs < 4.9 cc had a significant difference in maximum BP dose (64.2 vs 59.4 Gy; p = 0.001). Multivariate analysis confirmed that an ipsilateral GTV-N ≥ 4.9 cc was an independent predictor for the BP to receive a maximal dose of > 66 Gy when adjusted individually for BP volume, GTV-P, the use of a low anterior neck field technique, total planned radiation dose, and tumor category. Although both the primary and the nodal tumor volumes affected the BP maximal dose, the ipsilateral nodal tumor volume (GTV-N ≥ 4.9 cc) was an independent predictor for high maximal BP dose constraints in head and neck IMRT.

61 RADIATION PROTECTION AND DOSIMETRY↗

Dosimetric evaluation of simultaneous integrated boost during stereotactic body radiation therapy for pancreatic cancer

Stereotactic body radiation therapy (SBRT) provides a promising way to treat locally advanced pancreatic cancer and borderline resectable pancreatic cancer. A simultaneous integrated boost (SIB) to the region of vessel abutment or encasement during SBRT has the potential to downstage otherwise likely positive surgical margins. Despite the potential benefit of using SIB-SBRT, the ability to boost is limited by the local geometry of the organs at risk (OARs), such as stomach, duodenum, and bowel (SDB), relative to tumor. In this study, we have retrospectively replanned 20 patients with 25 Gy prescribed to the planning target volume (PTV) and 33~80 Gy to the boost target volume (BTV) using an SIB technique for all patients. The number of plans and patients able to satisfy a set of clinically established constraints is analyzed. The ability to boost vessels (within the gross target volume [GTV]) is shown to correlate with the overlap volume (OLV), defined to be the overlap between the GTV + a 1(OLV1)- or 2(OLV2)-cm margin with the union of SDB. Integral dose, boost dose contrast (BDC), biologically effective BDC, tumor control probability for BTV, and normal tissue complication probabilities are used to analyze the dosimetric results. More than 65% of the cases can deliver a boost to 40 Gy while satisfying all OAR constraints. An OLV2 of 100 cm{sup 3} is identified as the cutoff volume: for cases with OLV2 larger than 100 cm{sup 3}, it is very unlikely the case could achieve 25 Gy to the PTV while successfully meeting all the OAR constraints.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Spine stereotactic body radiation therapy plans: Achieving dose coverage, conformity, and dose falloff

We report our experience of establishing planning objectives to achieve dose coverage, conformity, and dose falloff for spine stereotactic body radiation therapy (SBRT) plans. Patients with spine lesions were treated using SBRT in our institution since September 2009. Since September 2011, we established the following planning objectives for our SBRT spine plans in addition to the cord dose constraints: (1) dose coverage—prescription dose (PD) to cover at least 95% planning target volume (PTV) and 90% PD to cover at least 99% PTV; (2) conformity index (CI)—ratio of prescription isodose volume (PIV) to the PTV < 1.2; (3) dose falloff—ratio of 50% PIV to the PTV (R{sub 50%}); (4) and maximum dose in percentage of PD at 2 cm from PTV in any direction (D{sub 2cm}) to follow Radiation Therapy Oncology Group (RTOG) 0915. We have retrospectively reviewed 66 separate spine lesions treated between September 2009 and December 2012 (31 treated before September 2011 [group 1] and 35 treated after [group 2]). The χ{sup 2} test was used to examine the difference in parameters between groups. The PTV V{sub 100%} {sub PD} ≥ 95% objective was met in 29.0% of group 1 vs 91.4% of group 2 (p < 0.01) plans. The PTV V{sub 90%} {sub PD} ≥ 99% objective was met in 38.7% of group 1 vs 88.6% of group 2 (p < 0.01) plans. Overall, 4 plans in group 1 had CI > 1.2 vs none in group 2 (p = 0.04). For D{sub 2cm}, 48.3% plans yielded a minor violation of the objectives and 16.1% a major violation for group 1, whereas 17.1% exhibited a minor violation and 2.9% a major violation for group 2 (p < 0.01). Spine SBRT plans can be improved on dose coverage, conformity, and dose falloff employing a combination of RTOG spine and lung SBRT protocol planning objectives.

61 RADIATION PROTECTION AND DOSIMETRY↗

Effect of stereotactic dosimetric end points on overall survival for Stage I non–small cell lung cancer: A critical review

Stereotactic body radiation therapy (SBRT) delivers a high biologically effective dose while minimizing toxicities to surrounding tissues. Within the scope of clinical trials and local practice, there are inconsistencies in dosimetrics used to evaluate plan quality. The purpose of this critical review was to determine if dosimetric parameters used in SBRT plans have an effect on local control (LC), overall survival (OS), and toxicities. A database of relevant trials investigating SBRT for patients with early-stage non–small cell lung cancer was compiled, and a table of dosimetric variables used was created. These parameters were compared and contrasted for LC, OS, and toxicities. Dosimetric end points appear to have no effect on OS or LC. Incidences of rib fractures correlate with a lack of dose-volume constraints (DVCs) reported. This review highlights the great disparity present in clinical trials reporting dosimetrics, DVCs, and toxicities for lung SBRT. Further evidence is required before standard DVCs guidelines can be introduced. Dosimetric end points specific to stereotactic treatment planning have been proposed but require further investigation before clinical implementation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Dosimetric quality, accuracy, and deliverability of modulated radiotherapy treatments for spinal metastases

Cancer often metastasizes to the vertebra, and such metastases can be treated successfully using simple, static posterior or opposed-pair radiation fields. However, in some cases, including when re-irradiation is required, spinal cord avoidance becomes necessary and more complex treatment plans must be used. This study evaluated 16 sample intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT) treatment plans designed to treat 6 typical vertebral and paraspinal volumes using a standard prescription, with the aim of investigating the advantages and limitations of these treatment techniques and providing recommendations for their optimal use in vertebral treatments. Treatment plan quality and beam complexity metrics were evaluated using the Treatment And Dose Assessor (TADA) code. A portal-imaging–based quality assurance (QA) system was used to evaluate treatment delivery accuracy, and radiochromic film measurements were used to provide high-resolution verification of treatment plan dose accuracy, especially in the steep dose gradient regions between each vertebral target and spinal cord. All treatment modalities delivered approximately the same doses and the same levels of dose heterogeneity to each planning target volume (PTV), although the minimum PTV doses in the vertebral plans were substantially lower than the prescription, because of the requirement that the plans meet a strict constraint on the dose to the spinal cord and cord planning risk volume (PRV). All plans met required dose constraints on all organs at risk, and all measured PTV-cord dose gradients were steeper than planned. Beam complexity analysis suggested that the IMRT treatment plans were more deliverable (less complex, leading to greater QA success) than the VMAT treatment plans, although the IMRT plans also took more time to deliver. The accuracy and deliverability of VMAT treatment plans were found to be substantially increased by limiting the number of monitor units (MU) per beam at the optimization stage, and thereby limiting beam modulation complexity. The VMAT arcs that were optimized with MU limitation had higher QA pass rates as well as higher modulation complexity scores (less complexity), lower modulation indices (less modulation), lower MU per beam, larger beam segments, and fewer small apertures than the VMAT arcs that were optimized without MU limitation. It is recommended that VMAT treatments for vertebral volumes, where the PTV abuts or surrounds the spinal cord, should be optimized with MU limitation. IMRT treatments may be preferable to the VMAT treatments, for dosimetry and deliverability reasons, but may be inappropriate for some patients because of their increased treatment delivery time.

61 RADIATION PROTECTION AND DOSIMETRY↗