Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model checking”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

DETERMINATION OF THE POINT-SPREAD FUNCTION FOR THE FERMI LARGE AREA TELESCOPE FROM ON-ORBIT DATA AND LIMITS ON PAIR HALOS OF ACTIVE GALACTIC NUCLEI

The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope is a pair-conversion telescope designed to detect photons with energies from Almost-Equal-To 20 MeV to >300 GeV. The pre-launch response functions of the LAT were determined through extensive Monte Carlo simulations and beam tests. The point-spread function (PSF) characterizing the angular distribution of reconstructed photons as a function of energy and geometry in the detector is determined here from two years of on-orbit data by examining the distributions of {gamma} rays from pulsars and active galactic nuclei (AGNs). Above 3 GeV, the PSF is found to be broader than the pre-launch PSF. We checked for dependence of the PSF on the class of {gamma}-ray source and observation epoch and found none. We also investigated several possible spatial models for pair-halo emission around BL Lac AGNs. We found no evidence for a component with spatial extension larger than the PSF and set upper limits on the amplitude of halo emission in stacked images of low- and high-redshift BL Lac AGNs and the TeV blazars 1ES0229+200 and 1ES0347-121.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Gravitational Waves from Merging Intermediate-mass Black Holes. II. Event Rates at Ground-based Detectors

Based on a dynamical formation model of a supermassive black hole (SMBH), we estimate the expected observational profile of gravitational waves at ground-based detectors, such as KAGRA or advanced LIGO/VIRGO. Noting that the second generation of detectors have enough sensitivity from 10 Hz and up (especially with KAGRA owing to its location at less seismic noise), we are able to detect the ring-down gravitational wave of a BH with mass M<2×10{sup 3}M{sub ⊙}. This enables us to check the sequence of BH mergers to SMBHs via intermediate-mass BHs. We estimate the number density of galaxies from the halo formation model and estimate the number of BH mergers from the giant molecular cloud model assuming hierarchical growth of merged cores. At the designed KAGRA (and/or advanced LIGO/VIRGO), we find that the BH merger of its total mass M∼60M{sub ⊙} is at the peak of the expected mass distribution. With its signal-to-noise ratio ρ=10 (30), we estimate the event rate R∼200 (20) per year in the most optimistic case, and we also find that BH mergers in the range M<150M{sub ⊙} are R>1 per year for ρ=10. Thus, if we observe a BH with more than 100M{sub ⊙} in future gravitational-wave observations, our model naturally explains its source.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Genarris 2.0: A Random Structure Generator for Molecular Crystals

Genarris is an open source Python package for generating random molecular crystal structures with physical constraints for seeding crystal structure prediction algorithms and training machine learning models. Here we present a new version of the code, containing several major improvements. A MPI-based parallelization scheme has been implemented, which facilitates the seamless sequential execution of user-defined workflows. A new method for estimating the unit cell volume based on the single molecule structure has been developed using a machine-learned model trained on experimental structures. A new algorithm has been implemented for generating crystal structures with molecules occupying special Wyckoff positions. A new hierarchical structure check procedure has been developed to detect unphysical close contacts efficiently and accurately. New intermolecular distance settings have been implemented for strong hydrogen bonds. To demonstrate these new features, we study two specific cases: benzene and glycine. Genarris finds the experimental structures of the two polymorphs of benzene and the three polymorphs of glycine. Program summary Program Title: Genarris 2.0 Program Files doi: http://dx.doi.org/10.17632/grx6mz4pjn.1 Licensing provisions: BSD-3 Clause Programming language: Python, C External routines/libraries: Spglib, ASE, pymatgen, SciPy, mpi4py, scikit-learn, PyTorch, FHI-aims. Nature of problem: Molecular crystal structure prediction. Solution method: Genarris 2.0 generates molecular crystal structures over the 230 space groups, on general and special Wyckoff positions, using physical constraints. Down-sampling of the generated structures may be performed subsequently, based on molecular crystal packing descriptors and an unsupervised machine learning algorithm. Lastly, ab initio structure relaxation may be performed for the final pool. Depending on the user-defined workflow implemented, Genarris may be used to generate diverse molecular crystal datasets to seed evolutionary algorithms or to train machine learning algorithms or as a standalone crystal structure prediction method. Restrictions: For crystal structure generation, the molecule of interest must be semi-rigid with no bond rotational degrees of freedom. Unusual features: Genarris 2.0 is a highly distributed program, making use of MPI for Python parallelization. The user has the ability to design and implement workflows by executing a user-defined list of procedures. Genarris 2.0 offers new features including a machine learning model for estimating the molecular volume in the solid state from the single molecule structure, structure generation in special Wyckoff positions of space groups, hierarchical structure checks including rigorous treatment of non-orthogonal structures, and clustering and down-selection workflows combining first principles simulations with machine learning. (C) 2020 Elsevier B.V. All rights reserved.

Crystal structure prediction↗
DOE OSTI2019

A framework for assessment of predictive capability maturity and its application in nuclear thermal hydraulics

Here this work presents a formalized and computerized framework for the assessment of decision regarding the adequacy of a simulation tool for a nuclear reactor application. The decision regarding a code’s adequacy for an application is dependent on the assessment of different attributes that govern verification, validation, and uncertainty quantification of the code. In this work, the focus is on code validation. Therefore, the framework is developed and illustrated from the perspective of decision regarding the validation assessment of a code. Code validation assessment is performed based on the validation test results, data applicability, and process quality assurance factors. The process quality assurance factors warrant the trustworthiness of the evidence and help in checking people and process compliance with respect to the standard requirements. The proposed framework is developed using an argument modeling technique called Goal Structuring Notation (GSN). Goal structuring notation facilitates structural knowledge representation, information abstraction, evidence incorporation, and provides a skeletal structure for quantitative maturity assessment. The decision schema for the development of the decision model is based on the Predictive Capability Maturity Model (PCMM) and Analytic Hierarchy Process (AHP) and formalized using Goal structuring notation. Each decision attribute is formulated as a claim, where the degree of validity of the claim (attribute’s assessment) is expressed using different maturity levels. The GSN representation of the decision model is transformed into a confidence network to provide evidence-based quantitative maturity assessment using the Bayesian network. A metric based on the expected utility of maturity levels, called expected distance metric, is proposed to measure the distance between target maturity and achieved maturity on a scale of zero to one. Expected distance metric helps in comparing the assessment of different attributes and identification of major areas of concern in terms of modeling capability, data needs, and quality of assessment process. The practical application of the framework is demonstrated by a case study on validation assessment of a thermal-hydraulic code for a challenge problem called Departure from Nucleate Boiling (DNB).

42 ENGINEERING↗
DOE OSTI2020

Enhancement of risk informed validation framework for external hazard scenario

In recent years, the U.S. Nuclear Regulatory Commission (USNRC) and the International Atomic Energy Agency (IAEA) have developed methodologies to assess the vulnerabilities of nuclear plants against site specific extreme hazards. In many cases, advanced simulation tools are being considered to simulate multi-physics, multi-scale phenomena and to evaluate vulnerability of nuclear facilities. The credibility of advanced simulation tools is assessed based on a formal verification, validation, and uncertainty quantification procedure. One of the key limitations in validation is the lack of relevant experimental data at system-level. This limitation leads to a decrease in the confidence of system-level risk predictions. Therefore, a robust validation framework is needed to formalize the confidence in predictive capability of advanced simulation results. Additionally, this study enhances the existing risk informed validation methodology, originally proposed by Kwag et al. [1] and Bodda et al. [2], by developing additional attributes and a new set of validation indicies for a complete and wider applicability of the framework. In this manuscript, the methodology to identify the critical path that leads to the system-level failure is illustrated. The overall validation is checked for completeness and consistency by comparing the critical path for both the system-level simulation and experimental models. The applicability of the code for an intended application is represented in terms of various maturity levels and helps in the process of decision making.

42 ENGINEERING↗
DOE OSTI2015

Excitation correlation photoluminescence in the presence of Shockley-Read-Hall recombination

Excitation correlation photoluminescence (ECPL) measurements are often analyzed in the approximation of a cross correlation of charge carrier populations generated by the two delayed pulses. In semiconductors, this approach is valid for a linear non-radiative recombination path, but not for a non-linear recombination rate as in the general Shockley-Read-Hall recombination scenario. Here, the evolution of the ECPL signal was studied for deep trap recombination following Shockley-Read-Hall statistics. Analytic solutions can be obtained for a fast minority trapping regime and steady state recombination. For the steady state case, our results show that the quadratic radiative term plays only a minor role, and that the shape of the measured signal is mostly determined by the non-linearity of the recombination itself. We find that measurements with unbalanced intense pump and probe pulses can directly provide information about the dominant non-radiative recombination mechanism. The signal traces follow the charge carrier concentrations, despite the complex origins of the signal, thus showing that ECPL can be applied to study charge carrier dynamics in semiconductors without requiring elaborate calculations. The model is compared with measurements on a reference sample with alternating layers of InGaAs/InAlAs that were additionally cross-checked with time resolved optical pump terahertz probe measurements and found to be in excellent agreement.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Probing the PeV region in the astrophysical neutrino spectrum using 𝜈 𝜇 from the Southern sky

IceCube has observed a diffuse astrophysical neutrino flux over the energy region from a few TeV to a few PeV. At PeV energies, the spectral shape is not yet well measured due to the low statistics of the data. This analysis probes the gap between 1 and 10 PeV by using high-energy downgoing muon neutrinos. Here, to reject the large atmospheric muon background, two complementary techniques are combined. The first technique selects events with high stochasticity to reject atmospheric muon bundles whose stochastic energy losses are smoothed due to high muon multiplicity. The second technique vetoes atmospheric muons with the IceTop surface array. Using 9 yrs of data, we found two neutrino candidate events in the signal region, consistent with expectation from background, each with relatively high signal probabilities. A joint maximum likelihood estimation is performed using this sample and an independent 9.5-yr sample of tracks to measure the neutrino spectrum. A likelihood ratio test is done to compare the single power-law (SPL) vs SPL+cutoff hypothesis; the SPL+cutoff model is not significantly better than the SPL. High-energy astrophysical objects from four source catalogs are also checked around the direction of the two events. No significant coincidence was found.

Abbasi, R. [Loyola University Chicago] (ORCID:0000↗
DOE OSTI2022

TRIPOLI-4 neutronics calculations for IAEA-CRP benchmark of CEFR start-up tests using new libraries JEFF-3.3 and ENDF/B-VIII

Sodium-cooled Fast Reactors (SFR) are one of the promising Generation IV fast reactors. The China Experimental Fast Reactor (CEFR) is a 65 MWth pool-type SFR with a high neutron leakage core using high enriched uranium oxide fuel. The CEFR start-up tests in 2010 consist of series of neutronics experiments. Essential experimental data are available from the 2018 IAEA-CEFR-CRP benchmark document and useful for the validation of neutron transport codes and nuclear data libraries. The TRIPOLI-4 Monte Carlo transport code is a general-purpose neutronics code using continuous-energy nuclear data libraries. It has a rich validation database covering different computational and experimental benchmark data sets to assure the accuracy and credibility of numerical studies. Previous TRIPOLI-4 SFR core physics calculations used mainly MOX fuels. SFR control rod worth studies with TRIPOLI-4 were mainly on big-size cores. To test the modeling capability of TRIPOLI-4 on fuel loading patterns and on variable control rods positions for SFRs and to check recent nuclear data libraries, it is interesting to investigate the CEFR start-up tests in this work, including core states from subcritical to supercritical, control rods and rod group worth, point kinetics parameters, and radial foil activation measurements. TRIPOLI-4 calculation results using new data libraries JEFF-3.3, ENDF/B-VIII, and those of ENDF/B-VII.1 for different CEFR core states are reported here. By means of different options of TRIPOLI-4, control-rod SAs reactivity worth, point kinetics parameters, and radial {sup 237}Np(n,f) fission rate distribution were successfully evaluated. Both ENDF/B-VIII and JEFF-3.3 nuclear data libraries provided reliable results with a difference of 170 +/- 11 pcm in k{sub eff}.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Design and Analysis of Components of a Pressure/Vacuum Vessel System [Slides]

Performed design and analysis tasks to support design and testing of a prototype Experimental Physics Package for Shock Physics Experiments at Proton Radiography. Testing of low-fidelity vessel system prototype is in progress. General finite element analysis (FEA) process: identify critical components, simplify model geometry, take advantage of symmetry, apply proper boundary conditions, create a mesh that will solve quickly but accurately, check for solution convergence.

42 ENGINEERING↗
DOE OSTI2024

Initial PIP-II Beam Current Monitor Fault Case Analyses & Beam Position Monitor Linearity Studies in CST Studio Suite

The use of non-invasive sensors & systems to measure particle beam characteristics is a crucial part of modern accelerator control systems due to their ability to return real time beam data while minimizing negative effects on beam quality. To ensure that one can be reasonably confident these sensors will behave as desired upon be-ing implemented within the beamline, simulations pre-dicting the performance of these sensors under beamline conditions can be used as a valuable tool for checking sensor functionality without a physical test bench. This paper details the design, testing, and results of two sensor models developed using CST Studio Suite soft-ware designed to mimic two sensors to be implemented within the PIP-III beamline: an elliptical, large-aperture beam position monitor (BPM) for which vertical & hori-zontal position signal linearity was analyzed, and an AC current transformer (ACCT) beam current monitor (BCM) used to search for potential fault cases within the BCM and beam pipe flange gaps. Special focus is given to the discovery of linearity variations within the BPM and the use of frequency domain techniques in the BCM fault case analyses.

Rouzky, A. R.↗
DOE OSTI2026

Diagnosing the representation of surface and layered soil moisture in Earth system models

Surface soil moisture (mrsos) and vertically integrated soil moisture (mrsol) over the top 10 cm should, by definition, be physically consistent in Earth System Models (ESMs). However, an evaluation of nine CMIP6 models reveals substantial inconsistencies: in some models, mrsos and integrated mrsol agree globally; in others, they align only in specific regions; and in a few, they diverge across all grid cells. These discrepancies arise from a combination of factors, including metadata errors, inconsistent variable definitions, or diagnostic sequencing within the model. We demonstrate how such issues can lead to significant biases, even when both variables are present and seemingly well-defined. As model complexity increases and multi-model comparisons become more common, assumptions about variable equivalence may lead to flawed conclusions. This study highlights the need for routine consistency checks, improved metadata standards, and community-wide practices that ensure reliability of derived variables across ESM outputs, particularly in preparation for CMIP7.

Earth system models↗
DOE OSTI2021

Distributed deep reinforcement learning for simulation control

Abstract Several applications in the scientific simulation of physical systems can be formulated as control/optimization problems. The computational models for such systems generally contain hyperparameters, which control solution fidelity and computational expense. The tuning of these parameters is non-trivial and the general approach is to manually ‘spot-check’ for good combinations. This is because optimal hyperparameter configuration search becomes intractable when the parameter space is large and when they may vary dynamically. To address this issue, we present a framework based on deep reinforcement learning (RL) to train a deep neural network agent that controls a model solve by varying parameters dynamically. First, we validate our RL framework for the problem of controlling chaos in chaotic systems by dynamically changing the parameters of the system. Subsequently, we illustrate the capabilities of our framework for accelerating the convergence of a steady-state computational fluid dynamics solver by automatically adjusting the relaxation factors of the discretized Navier–Stokes equations during run-time. The results indicate that the run-time control of the relaxation factors by the learned policy leads to a significant reduction in the number of iterations for convergence compared to the random selection of the relaxation factors. Our results point to potential benefits for learning adaptive hyperparameter learning strategies across different geometries and boundary conditions with implications for reduced computational campaign expenses 4 4 Data and codes available at https://github.com/Romit-Maulik/PAR-RL . .

42 ENGINEERING↗
DOE OSTI2022

The HUNTER Dynamic Human Reliability Analysis Tool: Coupling an External Plant Code

The Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) is a framework to support dynamic human reliability analysis (HRA) with the aim to develop standalone software to perform the dynamic HRA calculations. Within the HRA, human actions in nuclear power plants (NPPs) are predicated by plant states, and human actions influence the plant. In other words, plant operations are necessarily recursive, and it becomes challenging to model complex human-plant interactions. Consequently, we have linked two software simulations that complement those shortcomings. RELAP5-3D—the Reactor Excursion and Leak Analysis Program (RELAP; Aumiller, Tomlinson, and Bauer 2001) is the foundational thermal-hydraulic software used to model nuclear systems. Using RELAP5-3D, we have simulated the plant operations proceeding according to procedures developed to address emergent situations in NPPs. Plant operations include various actions such as the operator checking plant parameters, as well as actions that are continuously performed over time until a specific parameter reaches certain criteria. That means that HUNTER and RELAP5-3D exchange information with each other and should be carried out simultaneously over time. To simulate plant operations, which represent the actual operator checks of plant parameters and corresponding manual control actions, changes in plant status are identified through simulation and performed according to the criteria and order of the procedure. Thus, the goal of coupling HUNTER with RELAP5-3D is to facilitate synchronous coupling, where human and plant models provide iterative feedback loops that drive the course of actions. The advantage of coupling with RELAP5-3D to serve as the external environment module in HUNTER is the ability to customize the plant model and streamline for particular model applications. In this paper, we will address the key features of the coupling and the coupling structure built to perform the feedback loops.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Global River Topology (GRIT): A Bifurcating River Hydrography

Existing global river networks underpin a wide range of hydrological applications but do not represent channels with divergent river flows (bifurcations, multi‐threaded channels, canals), as these features defy the convergent flow assumption that elevation‐derived networks (e.g., HydroSHEDS, MERIT Hydro) are based on. Yet, bifurcations are important features of the global river drainage system, especially on large floodplains and river deltas, and are also often found in densely populated regions. Here we developed the first raster and vector‐based Global RIver Topology that not only represents the tributaries of the global drainage network but also the distributaries, including multi‐threaded rivers, canals and deltas. We achieve this by merging a 30 m Landsat‐based river mask with elevation‐generated streams to ensure a homogeneous drainage density outside of the river mask for rivers narrower than approximately 30 m. Crucially, we employ the new 30 m digital terrain model, FABDEM, based on TanDEM‐X, which shows greater accuracy over the traditionally used SRTM derivatives. After vectorization and pruning, directionality is assigned by a series of elevation, flow angle and continuity approaches. The new global network and its attributes are validated using gauging stations, comparison with existing networks, and randomized manual checks. The new network represents 19.6 million km of streams and rivers with drainage areas greater than 50 km 2 and includes 67,495 bifurcations. With the advent of hyper‐resolution modeling and artificial intelligence, GRIT is expected to greatly improve the accuracy of many river‐based applications such as flood forecasting, water availability and quality simulations, or riverine habitat mapping.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Testing the Kerr Black Hole Hypothesis Using X-Ray Reflection Spectroscopy

We present the first X-ray reflection model for testing the assumption that the metric of astrophysical black holes is described by the Kerr solution. We employ the formalism of the transfer function proposed by Cunningham. The calculations of the reflection spectrum of a thin accretion disk are split into two parts: the calculation of the transfer function and the calculation of the local spectrum at any emission point in the disk. The transfer function only depends on the background metric and takes into account all the relativistic effects (gravitational redshift, Doppler boosting, and light bending). Our code computes the transfer function for a spacetime described by the Johannsen metric and can easily be extended to any stationary, axisymmetric, and asymptotically flat spacetime. Transfer functions and single line shapes in the Kerr metric are compared to those calculated from existing codes to check that we reach the necessary accuracy. We also simulate some observations with NuSTAR and LAD/eXTP and fit the data with our new model to show the potential capabilities of current and future observations to constrain possible deviations from the Kerr metric.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Development of a Discrepancy Checker for the Digital Twin in a Supervisory Control System for a Thermal Energy Delivery System

Defined as a virtual representation of a physical object, process, or service, and used to support real-world decision-making, a digital twin (DT) can be utilized to combine classical and novel frameworks in sensors, state predictions, and multi-input/multi-output systems, and to enable optimal autonomous operations. However, a DT’s usefulness largely depends on its ability to adequately mirror the state of its physical counterpart, and this adequacy should be reflected by the level of uncertainty in the underlying simulation models when estimating and predicting quantities of interest (QOIs). Moreover, simulation models in a DT may involve multiple fidelities of representations—ranging from physics-based models to data-driven ones—but classical uncertainty quantification (UQ) methods struggle to handle numerous uncertainty sources, nor are they designed for real-time applications. This work presents a UQ-based discrepancy checking and diagnosis tool for a DT-based supervisory control system applied to a thermal energy delivery system (TEDS) at Idaho National Laboratory. The discrepancy checker was developed using metadata from an automated DT development process, and these metadata included different combinations of physical model forms and model parameters, training data and hyperparameters for surrogate models, and design parameters for supervisory control systems. Next, correlations between the uncertainty results and the metadata were established and then applied to the DT operations. The discrepancy checker evaluates the discrepancies between model predictions from virtual and sensor measurements and backtraces them to the corresponding major sources of uncertainty. The discrepancy checker showed reasonable performance in detecting discrepancies and diagnosing sources of uncertainty in testing scenarios.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Solid core dipoles and switching power supplies: Lower cost light sources?

As a result of improvements in power semiconductors, moderate frequency switching supplies can now provide the hundreds of amps typically required by accelerators with zero-to-peak noise in the kHz region ~ 0.06% in current or voltage mode. Modeling was undertaken using a finite electromagnetic program to determine if eddy currents induced in the solid steel of CEBAF magnets and small supplemental additions would bring the error fields down to the 5ppm level needed for beam quality. The expected maximum field of the magnet under consideration is 0.85 T and the DC current required to produce that field is used in the calculations. An additional 0.1% current ripple is added to the DC current at discrete frequencies 360 Hz, 720 Hz or 7200 Hz. Over the region of the pole within 0.5% of the central integrated BdL the resulting AC field changes can be reduced to less than 1% of the 0.1% input ripple for all frequencies, and a sixth of that at 7200 Hz. Doubling the current, providing 1.5 T central field, yielded the same fractional reduction in ripple at the beam for the cases checked. A small dipole was measured at 60, 120, 360 and 720 Hz in two conditions and the results compared to the larger model for the latter two frequencies with surprisingly good agreement. Thus, for light sources with aluminum vacuum vessels and full energy linac injection, the combination of solid core dipoles and switching power supplies may result in significant cost savings.

acceleration cavities and magnets superconducting ↗
DOE OSTI2016

Predictors of Trigeminal Neuropathy After Radiosurgery for Vestibular Schwannomas

Purpose: To analyze the relationship between dosimetric characteristics and symptoms related to trigeminal neuropathy (TN) observed after radiosurgery (RS) for vestibular schwannomas (VS); to propose guidelines to optimize planification in VS RS regarding TN preservation; and to detail the mechanism of TN impairment after VS RS. Methods and Materials: One hundred seventy-nine patients treated between 2011 and 2013 for VS RS and without trigeminal impairment before RS were included in a retrospective study. Univariate and multivariate analyses were performed to determine predictors of TN among characteristics of the patients, the dosimetry, and the VS. Results: There were 20 Koos grade 1, 99 grade 2, 57 grade 3, and 3 grade 4. Fourteen patients (7.8%) presented a transitory or permanent TN. Between the patients with and without TN after VS RS, there was no significant difference regarding dosimetry or VS volume itself. Significant differences (univariate analysis P<.05, Mann-Whitney test) were found for parameters related to the cisternal portion of the trigeminal nerve: total integrated dose, maximum dose, mean dose, volume of the Vth nerve (Vol{sub v}), and volume of the Vth nerve receiving at least 11 Gy (Vol{sub Vcist>11Gy}), but also for maximal dose to the Vth nerve nucleus and intra-axial portion (Dose max{sub Vax}). After multivariate analysis, the best model predicting TN included Vol{sub Vcist>11Gy} (P=.0045), Dose max{sub Vax} (P=.0006), and Vol{sub v} (P=.0058). The negative predictive value of this model was 97%. Conclusions: The parameters Vol{sub Vcist>11Gy}, Dose max{sub Vax}, and Vol{sub v} should be checked when designing dosimetry for VS RS.

62 RADIOLOGY AND NUCLEAR MEDICINE↗