Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistics”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Statistical Combination of ATLAS Run 2 Searches for Charginos and Neutralinos at the LHC

Statistical combinations of searches for charginos and neutralinos using various decay channels are performed using 139 fb −1 of 𝑝⁢𝑝 collision data at $\sqrt{𝑠}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. Searches targeting pure-wino chargino pair production, pure-wino chargino-neutralino production, or Higgsino production decaying via standard model 𝑊, 𝑍, or ℎ bosons are combined to extend the mass reach to the produced supersymmetric particles by 30–100 GeV. The depth of the sensitivity of the original searches is also improved by the combinations, lowering the 95% C.L. cross-section upper limits by 15%–40%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Statistically averaged molecular dynamics simulations of hydrogen diffusion in magnesium and magnesium hydrides

Magnesium has a different crystal structure from its dihydride with hydrogenation leading to a phase transition from the hexagonal closely packed Mg into a tetragonal α-MgH 2 rutile type structure. Such materials exhibit complex hydrogen uptake and release kinetics because hydrogen diffusivities significantly change when the crystal structure changes. To provide a foundational understanding of (de)hydrogenation kinetics that is applicable to all stages of the reaction, we performed statistically averaged molecular dynamics simulations to derive hydrogen diffusivities as a function of temperature and hydrogen content for both magnesium and magnesium hydride. Furthermore, our studies confirm that hydrogen diffusivities in magnesium hydride are much lower than in magnesium, in agreement with experimental data. Additionally, we observe that in either magnesium or magnesium hydride, higher hydrogen compositions result in reduced diffusivities. The latter was not revealed by prior experiments, which were conducted at fixed hydrogen composition. Finally, we discover a non-Arrhenius behavior in magnesium hydride. The physical origin of this behavior is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Statistics on magnetic properties of Co compounds: A database-driven method for discovering Co-based ferromagnets

The search for new ferromagnetic compounds is often targeted at known structure families, particularly those containing iron, cobalt, and manganese. Here, we propose a method to expand this search to lesser-known structure types, using a database of experimental Curie and Néel temperatures. This study, in particular, illustrates the case of compounds containing the element cobalt, and we demonstrate how the use of such a database can lead to the discovery of ferromagnetic materials that had previously been overlooked. We report statistics from a literature survey of the magnetic properties of Co-based compounds with more than 33 at. % Co. We classify more than 13 000 compounds by structure type, cobalt content, and magnetic ground state. From these data, compounds TaCo 2 Ga, La 6 Co 13 Bi, and Nd 2 Co 3 were identified as potential ferromagnets, and we confirm their ferromagnetic ordering theoretically via first-principles calculations and experimentally via synthesis and characterization measurements. In addition, the analysis is focused on the collection of data trends and discovery of ferromagnetic materials with easy-axis magnetic anisotropy. Both known ferromagnetic materials with unknown magnetic anisotropy, and unstudied compounds were considered. From the subset of known ferromagnets with unknown anisotropy, the compound Co 2 Mg was synthesized, characterized, and determined to have easy-axis magnetic anisotropy at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Statistical mechanics of kinks on a gliding screw dislocation

The ability of a body-centered-cubic metal to deform plastically is limited by the thermally activated glide motion of screw dislocations, which are line defects with a mobility exhibiting complex dependence on temperature, stress, and dislocation segment length. We derive an analytical expression for the velocity of dislocation glide, based on a statistical mechanics argument, and identify an apparent phase transition marked by a crossover temperature above which the activation energy for glide effectively halves, changing from the formation energy of a double kink to that of a single kink. The analysis is in quantitative agreement with direct kinetic Monte Carlo simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Employing Fracture Statistics to Track Cell Reliability Through Module Fabrication

In this work we demonstrate how fracture statistics may be employed to track cell reliability through module fabrication steps. We were able to detect how cell cutting equipment could reduce cell fracture strain by 25 % as well as increase its variability. A smaller reduction was found to result from the ribbon soldering operation, and all cells measured were weaker when strained in an orientation perpendicular to the gridlines. Cell fracture measurements were found to be consistent with the number of cracked cells experienced in full sized modules exposed to mechanical loading. A structural mechanics model was created to elucidate how the requisite fracture strain is developed in a cell encapsulated within a PV module and demonstrates that this strain decreases more than 50 % from its center to edge when placed within the inner loading span of four-point flexure. The associated implications on the apparent cell strength is examined with an effective volume treatment.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2019

Estimation of Behind-the-Meter Solar Generation by Integrating Physical with Statistical Models

Accurate estimation of solar photovoltaic (PV) generation is crucial for distribution grid control and optimization. Unfortunately, most of the residential solar PV installations are behind-the-meter. Thus, utilities only have access to the net load readings. This paper presents an unsupervised framework for estimating solar PV generation by disaggregating the net load readings. The proposed framework synergistically combines a physical PV system performance model with a statistical model for load estimation. Specifically, our algorithm iteratively estimates solar PV generation with a physical model and electric load with the Hidden Markov model regression. The proposed algorithm is also capable of estimating the key technical parameters of the solar PV systems. Our proposed method is validated against net load and solar PV generation data gathered from residential customers located in Austin, Texas. The validation results show that our method reduces mean squared error by 44% compared to the state-of-the-art disaggregation algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Statistical and clustering analysis for disturbances: A case study of voltage dips in wind farms

This study proposes and evaluates an alternative statistical methodology to analyze a large number of voltage dips. For a given voltage dip, a set of lengths is first identified to characterize the root mean square (rms) voltage evolution along the disturbance, deduced from partial linearized time intervals and trajectories. Principal component analysis and K-means clustering processes are then applied to identify rms-voltage patterns and propose a reduced number of representative rms-voltage profiles from the linearized trajectories. This reduced group of averaged rms-voltage profiles enables the representation of a large amount of disturbances, which offers a visual and graphical representation of their evolution along the events, aspects that were not previously considered in other contributions. The complete process is evaluated on real voltage dips collected in intense field-measurement campaigns carried out in a wind farm in Spain among different years. The results are included in this paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Characterization and differentiation of aluminum powders used in improvised explosive devices. Part 2: Micromorphometric method refinement and preliminary statistical analysis

Abstract Aluminum (Al) powder is commonly encountered in improvised explosive devices (IEDs) as a metallic fuel due to its availability and low cost. Although available commercially in powder form, amateur bomb‐makers also produce their own Al powder via simple methods found online. In order to provide investigative leads and forensic intelligence, it is important to evaluate not only the composition of homemade devices, but also to distinguish between the various forms of Al powder they contain. To achieve this goal, a method using automated microscopy in combination with statistical techniques has been demonstrated to have the potential to provide source discrimination and investigative leads in source attribution of Al powders in IEDs. The present research refined this method and investigated 59 industrially and amateurly produced Al powder sources with seven subsamples per source using two traditional linear discriminant analyses (LDA), one with a standard data split for training and testing, and another using leave‐one‐out cross‐validation. Averaging the classification accuracies for the two LDA‐based analyses, LDA has the ability to correctly classify 59.26%, 83.35%, and 80.69% of the samples based on their powder source, type, and production method, respectively. This classification accuracy represents a 3407%, 317%, and 61.38% increase in accuracy from random class assignment, respectively. Further, in most instances of incorrect data attribution to a particular source, the subsample has been misidentified with another sample of the same powder type or production method.

Ommen, Danica M.↗
DOE OSTI2025

A statistical representation of bond coating oxidation under environmental barrier coatings

Environmental barrier coatings (EBCs) protect SiC-based ceramic matrix composites (CMCs) in turbine hot sections from high-temperature volatilization in combustion gases. The formation of a SiO 2 thermally grown oxide (TGO) is expected under the EBC after long-term operation. The oxidation resistance of the EBC is understood as a life-limiting factor for the CMC, and this work predicts long-term oxidation behavior under EBCs through a simple statistical approach. Specimens were exposed to 1350°C isothermal conditions for 100-h thermal cycles in flowing steam for up to 1000 h. The EBC morphology, SiO 2 thickness, and SiO 2 cracking behavior were assessed. Using thousands of SiO 2 thickness measurements across many millimeters of the interface, a realistic representation of the entire TGO was captured via a lognormal distribution. The lognormal fit parameters were extrapolated out to 25 000 h to assess the degree of SiO 2 growth, the spread of SiO 2 thicknesses related to the rough oxidizing interface, and percentages of the intermediate bond coating consumed. In conclusion, local interfacial defects from the coating deposition process are identified as local failure points for EBC—CMC systems.

SiO 2↗
DOE OSTI2021

Seeded Stage III glass dissolution behavior of a statistically designed glass matrix

The glass dissolution rate of some glasses accelerates after prolonged time spent at a slow, residual glass dissolution rate. This phenomenon is referred to as Stage III behavior. Here, the acceleration in glass dissolution rate linked to Stage III behavior is significant and may be the most impactful to long-term performance of glass in a repository. This work is aimed at understanding the effect of glass composition on Stage III behavior to add a level of technical defensibility to glass disposal. To this end, a set of twenty-four glass compositions were statistically designed, where eight glass components (SiO 2 , B 2 O 3 , Al 2 O 3 , CaO, Na 2 O, SnO 2 , ZrO 2 , and Others) have been independently varied in order to study the individual effects of each. These glasses have been subjected to static dissolution tests at 90 °C in deionized water and then seeded with zeolite Na-P2 28 days into the testing to induce Stage III behavior. The response of the glasses to the zeolite seeds fell into four primary types: 1) no response to seeds; 2) an immediate linear sustained acceleration in the rate; 3) an immediate linear acceleration in the rate followed by a decrease; and, 4) a progressive acceleration in the rate that is concurrent with the addition of the seeds. The main glass components observed to influence these behaviors were CaO, Al 2 O 3 , B 2 O 3 , and ZrO 2 , where: 1) CaO influenced which glasses showed a Stage III response to seeds (high CaO: Types 2, 3, and 4) or did not respond to seeds (low CaO: Type 1), 2) Al 2 O 3 and B 2 O 3 influenced which glasses showed a sustainable Stage III response (high Al 2 O 3 : Types 2 and 4) versus transitory response (low Al 2 O 3 and high B 2 O 3 : Type 3), and 3) ZrO 2 concentration influenced whether glasses showed a linear (high ZrO 2 : Type 2) versus progressive (low ZrO 2 : Type 4) response to seeds.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Investigation of statistical iterative reconstruction for dedicated breast CT

Purpose: Dedicated breast CT has great potential for improving the detection and diagnosis of breast cancer. Statistical iterative reconstruction (SIR) in dedicated breast CT is a promising alternative to traditional filtered backprojection (FBP). One of the difficulties in using SIR is the presence of free parameters in the algorithm that control the appearance of the resulting image. These parameters require tuning in order to achieve high quality reconstructions. In this study, the authors investigated the penalized maximum likelihood (PML) method with two commonly used types of roughness penalty functions: hyperbolic potential and anisotropic total variation (TV) norm. Reconstructed images were compared with images obtained using standard FBP. Optimal parameters for PML with the hyperbolic prior are reported for the task of detecting microcalcifications embedded in breast tissue.Methods: Computer simulations were used to acquire projections in a half-cone beam geometry. The modeled setup describes a realistic breast CT benchtop system, with an x-ray spectra produced by a point source and an a-Si, CsI:Tl flat-panel detector. A voxelized anthropomorphic breast phantom with 280 μm microcalcification spheres embedded in it was used to model attenuation properties of the uncompressed woman's breast in a pendant position. The reconstruction of 3D images was performed using the separable paraboloidal surrogates algorithm with ordered subsets. Task performance was assessed with the ideal observer detectability index to determine optimal PML parameters.Results: The authors' findings suggest that there is a preferred range of values of the roughness penalty weight and the edge preservation threshold in the penalized objective function with the hyperbolic potential, which resulted in low noise images with high contrast microcalcifications preserved. In terms of numerical observer detectability index, the PML method with optimal parameters yielded substantially improved performance (by a factor of greater than 10) compared to FBP. The hyperbolic prior was also observed to be superior to the TV norm. A few of the best-performing parameter pairs for the PML method also demonstrated superior performance for various radiation doses. In fact, using PML with certain parameter values results in better images, acquired using 2 mGy dose, than FBP-reconstructed images acquired using 6 mGy dose.Conclusions: A range of optimal free parameters for the PML algorithm with hyperbolic and TV norm-based potentials is presented for the microcalcification detection task, in dedicated breast CT. The reported values can be used as starting values of the free parameters, when SIR techniques are used for image reconstruction. Significant improvement in image quality can be achieved by using PML with optimal combination of parameters, as compared to FBP. Importantly, these results suggest improved detection of microcalcifications can be obtained by using PML with lower radiation dose to the patient, than using FBP with higher dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Dose impact in radiographic lung injury following lung SBRT: Statistical analysis and geometric interpretation

Purpose: To demonstrate a new method of evaluating dose response of treatment-induced lung radiographic injury post-SBRT (stereotactic body radiotherapy) treatment and the discovery of bimodal dose behavior within clinically identified injury volumes. Methods: Follow-up CT scans at 3, 6, and 12 months were acquired from 24 patients treated with SBRT for stage-1 primary lung cancers or oligometastic lesions. Injury regions in these scans were propagated to the planning CT coordinates by performing deformable registration of the follow-ups to the planning CTs. A bimodal behavior was repeatedly observed from the probability distribution for dose values within the deformed injury regions. Based on a mixture-Gaussian assumption, an Expectation-Maximization (EM) algorithm was used to obtain characteristic parameters for such distribution. Geometric analysis was performed to interpret such parameters and infer the critical dose level that is potentially inductive of post-SBRT lung injury. Results: The Gaussian mixture obtained from the EM algorithm closely approximates the empirical dose histogram within the injury volume with good consistency. The average Kullback-Leibler divergence values between the empirical differential dose volume histogram and the EM-obtained Gaussian mixture distribution were calculated to be 0.069, 0.063, and 0.092 for the 3, 6, and 12 month follow-up groups, respectively. The lower Gaussian component was located at approximately 70% prescription dose (35 Gy) for all three follow-up time points. The higher Gaussian component, contributed by the dose received by planning target volume, was located at around 107% of the prescription dose. Geometrical analysis suggests the mean of the lower Gaussian component, located at 35 Gy, as a possible indicator for a critical dose that induces lung injury after SBRT. Conclusions: An innovative and improved method for analyzing the correspondence between lung radiographic injury and SBRT treatment dose has been demonstrated. Bimodal behavior was observed in the dose distribution of lung injury after SBRT. Novel statistical and geometrical analysis has shown that the systematically quantified low-dose peak at approximately 35 Gy, or 70% prescription dose, is a good indication of a critical dose for injury. The determined critical dose of 35 Gy resembles the critical dose volume limit of 30 Gy for ipsilateral bronchus in RTOG 0618 and results from previous studies. The authors seek to further extend this improved analysis method to a larger cohort to better understand the interpatient variation in radiographic lung injury dose response post-SBRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

A statistical method for retrospective cardiac and respiratory motion gating of interventional cardiac x-ray images

Purpose: Image-guided cardiac interventions involve the use of fluoroscopic images to guide the insertion and movement of interventional devices. Cardiorespiratory gating can be useful for 3D reconstruction from multiple x-ray views and for reducing misalignments between 3D anatomical models overlaid onto fluoroscopy. Methods: The authors propose a novel and potentially clinically useful retrospective cardiorespiratory gating technique. The principal component analysis (PCA) statistical method is used in combination with other image processing operations to make our proposed masked-PCA technique suitable for cardiorespiratory gating. Unlike many previously proposed techniques, our technique is robust to varying image-content, thus it does not require specific catheters or any other optically opaque structures to be visible. Therefore, it works without any knowledge of catheter geometry. The authors demonstrate the application of our technique for the purposes of retrospective cardiorespiratory gating of normal and very low dose x-ray fluoroscopy images. Results: For normal dose x-ray images, the algorithm was validated using 28 clinical electrophysiology x-ray fluoroscopy sequences (2168 frames), from patients who underwent radiofrequency ablation (RFA) procedures for the treatment of atrial fibrillation and cardiac resynchronization therapy procedures for heart failure. The authors established end-systole, end-expiration, and end-inspiration success rates of 97.0%, 97.9%, and 97.0%, respectively. For very low dose applications, the technique was tested on ten x-ray sequences from the RFA procedures with added noise at signal to noise ratio (SNR) values of√(5)0, √(1)0, √(8), √(6), √(5), √(2), and √(1) to simulate the image quality of increasingly lower dose x-ray images. Even at the low SNR value of √(2), representing a dose reduction of more than 25 times, gating success rates of 89.1%, 88.8%, and 86.8% were established. Conclusions: The proposed technique can therefore extract useful information from interventional x-ray images while minimizing exposure to ionizing radiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-CAMPUS-T-04: Statistical Process Control for Patient-Specific QA in Proton Beams

Purpose: To evaluate and improve the reliability of proton QA process, to provide an optimal customized level using the statistical process control (SPC) methodology. The aim is then to suggest the suitable guidelines for patient-specific QA process. Methods: We investigated the constancy of the dose output and range to see whether it was within the tolerance level of daily QA process. This study analyzed the difference between the measured and calculated ranges along the central axis to suggest the suitable guidelines for patient-specific QA in proton beam by using process capability indices. In this study, patient QA plans were classified into 6 treatment sites: head and neck (41 cases), spinal cord (29 cases), lung (28 cases), liver (30 cases), pancreas (26 cases), and prostate (24 cases). Results: The deviations for the dose output and range of daily QA process were ±0.84% and ±019%, respectively. Our results show that the patient-specific range measurements are capable at a specification limit of ±2% in all treatment sites except spinal cord cases. In spinal cord cases, comparison of process capability indices (Cp, Cpm, Cpk ≥1, but Cpmk ≤1) indicated that the process is capable, but not centered, the process mean deviates from its target value. The UCL (upper control limit), CL (center line) and LCL (lower control limit) for spinal cord cases were 1.37%, −0.27% and −1.89%, respectively. On the other hands, the range differences in prostate cases were good agreement between calculated and measured values. The UCL, CL and LCL for prostate cases were 0.57%, −0.11% and −0.78%, respectively. Conclusion: SPC methodology has potential as a useful tool to customize an optimal tolerance levels and to suggest the suitable guidelines for patient-specific QA in clinical proton beam.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Feasibility study of using statistical process control to customized quality assurance in proton therapy

Purpose: To evaluate and improve the reliability of proton quality assurance (QA) processes and, to provide an optimal customized tolerance level using the statistical process control (SPC) methodology. Methods: The authors investigated the consistency check of dose per monitor unit (D/MU) and range in proton beams to see whether it was within the tolerance level of the daily QA process. This study analyzed the difference between the measured and calculated ranges along the central axis to improve the patient-specific QA process in proton beams by using process capability indices. Results: The authors established a customized tolerance level of ±2% for D/MU and ±0.5 mm for beam range in the daily proton QA process. In the authors’ analysis of the process capability indices, the patient-specific range measurements were capable of a specification limit of ±2% in clinical plans. Conclusions: SPC methodology is a useful tool for customizing the optimal QA tolerance levels and improving the quality of proton machine maintenance, treatment delivery, and ultimately patient safety.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Multileaf collimator performance monitoring and improvement using semiautomated quality control testing and statistical process control

Purpose: High-quality radiation therapy using highly conformal dose distributions and image-guided techniques requires optimum machine delivery performance. In this work, a monitoring system for multileaf collimator (MLC) performance, integrating semiautomated MLC quality control (QC) tests and statistical process control tools, was developed. The MLC performance monitoring system was used for almost a year on two commercially available MLC models. Control charts were used to establish MLC performance and assess test frequency required to achieve a given level of performance. MLC-related interlocks and servicing events were recorded during the monitoring period and were investigated as indicators of MLC performance variations. Methods: The QC test developed as part of the MLC performance monitoring system uses 2D megavoltage images (acquired using an electronic portal imaging device) of 23 fields to determine the location of the leaves with respect to the radiation isocenter. The precision of the MLC performance monitoring QC test and the MLC itself was assessed by detecting the MLC leaf positions on 127 megavoltage images of a static field. After initial calibration, the MLC performance monitoring QC test was performed 3–4 times/week over a period of 10–11 months to monitor positional accuracy of individual leaves for two different MLC models. Analysis of test results was performed using individuals control charts per leaf with control limits computed based on the measurements as well as two sets of specifications of ±0.5 and ±1 mm. Out-of-specification and out-of-control leaves were automatically flagged by the monitoring system and reviewed monthly by physicists. MLC-related interlocks reported by the linear accelerator and servicing events were recorded to help identify potential causes of nonrandom MLC leaf positioning variations. Results: The precision of the MLC performance monitoring QC test and the MLC itself was within ±0.22 mm for most MLC leaves and the majority of the apparent leaf motion was attributed to beam spot displacements between irradiations. The MLC QC test was performed 193 and 162 times over the monitoring period for the studied units and recalibration had to be repeated up to three times on one of these units. For both units, rate of MLC interlocks was moderately associated with MLC servicing events. The strongest association with the MLC performance was observed between the MLC servicing events and the total number of out-of-control leaves. The average elapsed time for which the number of out-of-specification or out-of-control leaves was within a given performance threshold was computed and used to assess adequacy of MLC test frequency. Conclusions: A MLC performance monitoring system has been developed and implemented to acquire high-quality QC data at high frequency. This is enabled by the relatively short acquisition time for the images and automatic image analysis. The monitoring system was also used to record and track the rate of MLC-related interlocks and servicing events. MLC performances for two commercially available MLC models have been assessed and the results support monthly test frequency for widely accepted ±1 mm specifications. Higher QC test frequency is however required to maintain tighter specification and in-control behavior.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-E-J-177: A Computational Approach for Determination of Anisotropic PTV Margins Based On Statistical Shape Analysis for Prostate Cancer Radiotherapy

Purpose: Our aim of this study was to propose a computational approach for determination of anisotropic planning target volume (PTV) margins based on statistical shape analysis with taking into account time variations of clinical target volume (CTV) shapes for the prostate cancer radiation treatment planning (RTP). Methods: Systematic and random setup errors were measured using orthogonal projection and cone beam computed tomography (CBCT) images for data of 20 patients, who underwent the intensity modulated radiation therapy for prostate cancer. The low-risk, intermediate-risk, and high-risk CTVs were defined as only a prostate, a prostate plus proximal 1-cm seminal vesicles, and a prostate plus proximal 2-cm seminal vesicles, respectively. All CTV regions were registered with a reference CTV region with a median volume to remove the effect of the setup errors, and converted to a point distribution models. The systematic and random errors for translations of CTV regions were automatically evaluated by analyzing the movements of centroids of CTV regions. The random and systematic errors for shape variations of CTV regions were obtained from covariance matrices based on point distributions for the CTV contours on CBCT images of 72 fractions of 10 patients. Anisotropic PTV margins for 6 directions (right, left, anterior, posterior, superior and inferior) were derived by using Yoda’s PTV margin model. Results: PTV margins with and without shape variations were 5.75 to 8.03 mm and 5.23 to 7.67 mm for low-risk group, 5.87 to 8.33 mm and 5.23 to 7.67 mm for intermediate-risk group, and 5.88 to 8.25 mm and 5.29 to 7.82 mm for highrisk group, respectively. Conclusion: The proposed computational approach could be feasible for determination of the anisotropic PTV margins with taking into account CTV shape variations for the RTP.

60 APPLIED LIFE SCIENCES↗