Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

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

Reactor Containment Passive Safety Analysis: Steam Condensation in Presence of Non-condensable Gas Scaled Experiment and Modeling

This study presents steam condensation scaled experiments and semi-empirical models in presence of nitrogen (N)—a noncondensable gas (NCG), simulating air in the reactor containment—to support water-cooled small modular reactors (SMRs) passive containment cooling system (PCCS) design and analysis. Previous experimental studies on PCCS are focused on fixed and smaller tube (mostly 2-in.) geometries and specific test condition variations, bringing challenges with geometric scaling and mismatching with SMR prototypic design. To address these challenges, this study presents steam condensation test dataset obtained from three scaled test sections of 1-, 2-, and 4-in.-diameter steam condensers with an annular/jacket cooling of 2-, 3-, and 6 in.-diameter tubes, respectively. Test data were collected for steam ranges from 58 to 63 kg/hr., and NCG flow of 4.4 to 13.3 kg/hr. Annular cooling water flow was varied to obtain required testing conditions of saturated steam inlet and fully condensed outlet. Axial temperature test data of bulk cooling water, steam and condensate were collected by thermocouples for three test sections and various steam-NCG mixing/testing conditions. A standard data reduction method was adopted—utilizing iterative and nodalized mass and heat transfer calculation—to estimate axial local heat fluxes, heat transfer coefficients (HTCs), condensation rates, film thickness, and Nusselt number. Based on the obtained dataset semi-empirical model results—a ratio of experimental and Nusselt’s theoretical HTC are presented. Such results and findings are supportive of developing scaled-up testing facility, to enable model validations and accelerate next generation of reactors development and deployment

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Presentation: Reactor Containment Passive Safety Analysis: Steam Condensation in Presence of Non-condensable Gas Scaled Experiment and Modeling

This study presents steam condensation scaled experiments and semi-empirical models in presence of nitrogen--a noncondensable gas (NCG), simulating air in the reactor containment--to support water-cooled small modular reactors (SMRs) passive containment cooling system (PCCS) design and analysis. Previous experimental studies on PCCS are focused on fixed and smaller tube (mostly 2-in.) geometries and specific test condition variations, bringing challenges with geometric scaling and mismatching with SMR prototypic design. To address these challenges, this study presents steam condensation test dataset obtained from three scaled test sections of 1-, 2-, and 4-in.-diameter steam condensers with an annular/jacket cooling of 2-, 3-, and 6 in.-diameter tubes, respectively. Test data were collected for steam ranges from 58 to 63 kg/hr., and NCG flow of 4.4 to 13.3 kg/hr. Annular cooling water flow was varied to obtain required testing conditions of saturated steam inlet and fully condensed outlet. Axial temperature test data of bulk cooling water, steam and condensate were collected by thermocouples for three test sections and various steam-NCG mixing/testing conditions. A standard data reduction method was adopted--utilizing iterative and nodalized mass and heat transfer calculation to estimate axial local heat fluxes, heat transfer coefficients (HTCs), condensation rates, film thickness, and Nusselt number. Based on the obtained dataset semi-empirical model results--a ratio of experimental and Nusselt's theoretical HTC are presented. Such results and findings are supportive of developing scaled-up testing facility, to enable model validations and accelerate next generation of reactors development and deployment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Strategies for automatic online treatment plan reoptimization using clinical treatment planning system: A planning parameters study

Purpose: Adaptive radiation therapy for prostate cancer using online reoptimization provides an improved control of interfractional anatomy variations. However, the clinical implementation of online reoptimization is currently limited by the low efficiency of current strategies and the difficulties associated with integration into the current treatment planning system. This study investigates the strategies for performing fast (∼2 min) automatic online reoptimization with a clinical fluence-map-based treatment planning system; and explores the performance with different input parameters settings: dose-volume histogram (DVH) objective settings, starting stage, and iteration number (in the context of real time planning).Methods: Simulated treatments of 10 patients were reoptimized daily for the first week of treatment (5 fractions) using 12 different combinations of optimization strategies. Options for objective settings included guideline-based RTOG objectives, patient-specific objectives based on anatomy on the planning CT, and daily-CBCT anatomy-based objectives adapted from planning CT objectives. Options for starting stages involved starting reoptimization with and without the original plan's fluence map. Options for iteration numbers were 50 and 100. The adapted plans were then analyzed by statistical modeling, and compared both in terms of dosimetry and delivery efficiency.Results: All online reoptimized plans were finished within ∼2 min with excellent coverage and conformity to the daily target. The three input parameters, i.e., DVH objectives, starting stage, and iteration number, contributed to the outcome of optimization nearly independently. Patient-specific objectives generally provided better OAR sparing compared to guideline-based objectives. The benefit in high-dose sparing from incorporating daily anatomy into objective settings was positively correlated with the relative change in OAR volumes from planning CT to daily CBCT. The use of the original plan fluence map as the starting stage reduced OAR dose at the mid-dose region, but increased the monitor units by 17%. Differences of only 2cc or less in OAR V50%/V70Gy/V76Gy were observed between 100 and 50 iterations.Conclusions: It is feasible to perform automatic online reoptimization in ∼2 min using a clinical treatment planning system. Selecting optimal sets of input parameters is the key to achieving high quality reoptimized plans, and should be based on the individual patient's daily anatomy, delivery efficiency, and time allowed for plan adaptation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Hierarchical Krylov and nested Krylov methods for extreme-scale computing

The solution of large, sparse linear systems is typically a dominant phase of computation for simulations based on partial differential equations, which are ubiquitous in scientific and engineering applications. While preconditioned Krylov methods are widely used and provide many advantages for solving sparse linear systems that do not have highly convergent, geometric multigrid solvers or specialized fast solvers, Krylov methods encounter well-known scaling difficulties for over 10,000 processor cores because each iteration requires at least one vector inner product, which in turn requires a global synchronization that scales poorly because of internode latency. To aid in overcoming these difficulties, we have developed hierarchical Krylov methods and nested Krylov methods in the PETSc library that reduce the number of global inner products required across the entire system (where they are expensive), though freely allow vector inner products across smaller subsets of the entire system (where they are inexpensive) or use inner iterations that do not invoke vector inner products at all. Nested Krylov methods are a generalization of inner-outer iterative methods with two or more layers. Hierarchical Krylov methods are a generalization of block Jacobi and overlapping additive Schwarz methods, where each block itself is solved by Krylov methods on smaller blocks. Conceptually, the hierarchy can continue recursively to an arbitrary number of levels of smaller and smaller blocks. As a specific case, we introduce the hierarchical FGMRES method, or h-FGMRES, and we demonstrate the impact of two-level h-FGMRES with a variable preconditioner on the PFLOTRAN subsurface flow application. We also demonstrate the impact of nested FGMRES, BiCGStab and Chebyshev methods. These hierarchical Krylov methods and nested Krylov methods significantly reduced overall PFLOTRAN simulation time on the Cray XK6 when using 10,000 through 224,000 cores through the combined effects of reduced global synchronization due to fewer global inner products and stronger inner hierarchical or nested preconditioners.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An investigation of Newton-Sketch and subsampled Newton methods

Sketching, a dimensionality reduction technique, has received much attention in the statistics community. In this paper, we study sketching in the context of Newton's method for solving finite-sum optimization problems in which the number of variables and data points are both large. In this work, we study two forms of sketching that perform dimensionality reduction in data space: Hessian subsampling and randomized Hadamard transformations. Each has its own advantages, and their relative tradeoffs have not been investigated in the optimization literature. Additionally, our study focuses on practical versions of the two methods in which the resulting linear systems of equations are solved approximately, at every iteration, using an iterative solver. The advantages of using the conjugate gradient method vs. a stochastic gradient iteration are revealed through a set of numerical experiments, and a complexity analysis of the Hessian subsampling method is presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Improving the efficiency of molecular replacement by utilizing a new iterative transform phasing algorithm

An iterative transform algorithm is proposed to improve the conventional molecular-replacement method for solving the phase problem in X-ray crystallography. Several examples of successful trial calculations carried out with real diffraction data are presented. An iterative transform method proposed previously for direct phasing of high-solvent-content protein crystals is employed for enhancing the molecular-replacement (MR) algorithm in protein crystallography. Target structures that are resistant to conventional MR due to insufficient similarity between the template and target structures might be tractable with this modified phasing method. Trial calculations involving three different structures are described to test and illustrate the methodology. The relationship of the approach to PHENIX Phaser-MR and MR-Rosetta is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

MICADO: Parallel implementation of a 2D-1D iterative algorithm for the 3D neutron transport problem in prismatic geometries

The large increase in computing power over the past few years now makes it possible to consider developing 3D full-core heterogeneous deterministic neutron transport solvers for reference calculations. Among all approaches presented in the literature, the method first introduced in [1] seems very promising. It consists in iterating over resolutions of 2D and ID MOC problems by taking advantage of prismatic geometries without introducing approximations of a low order operator such as diffusion. However, before developing a solver with all industrial options at EDF, several points needed to be clarified. In this work, we first prove the convergence of this iterative process, under some assumptions. We then present our high-performance, parallel implementation of this algorithm in the MICADO solver. Benchmarking the solver against the Takeda case shows that the 2D-1D coupling algorithm does not seem to affect the spatial convergence order of the MOC solver. As for performance issues, our study shows that even though the data distribution is suited to the 2D solver part, the efficiency of the ID part is sufficient to ensure a good parallel efficiency of the global algorithm. After this study, the main remaining difficulty implementation-wise is about the memory requirement of a vector used for initialization. An efficient acceleration operator will also need to be developed. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Benchmark of numerical tools simulating beam propagation and secondary particles in ITER NBI

Injection of high energy beams of neutral particles is a method for plasma heating in fusion devices. The ITER injector, and its prototype MITICA (Megavolt ITER Injector and Concept Advancement), are large extrapolations from existing devices: therefore numerical modeling is needed to set thermo-mechanical requirements for all beam-facing components. As the power and charge deposition originates from several sources (primary beam, co-accelerated electrons, and secondary production by beam-gas, beam-surface, and electron-surface interaction), the beam propagation along the beam line is simulated by comprehensive 3D models. This paper presents a comparative study between two codes: BTR has been used for several years in the design of the ITER HNB/DNB components; SAMANTHA code was independently developed and includes additional phenomena, such as secondary particles generated by collision of beam particles with the background gas. The code comparison is valuable in the perspective of the upcoming experimental operations, in order to prepare a reliable numerical support to the interpretation of experimental measurements in the beam test facilities. The power density map calculated on the Electrostatic Residual Ion Dump (ERID) is the chosen benchmark, as it depends on the electric and magnetic fields as well as on the evolution of the beam species via interaction with the gas. Finally the paper shows additional results provided by SAMANTHA, like the secondary electrons produced by volume processes accelerated by the ERID fringe-field towards the Cryopumps.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Enabling Predictive Simulation and UQ of Complex Multiphysics PDE Systems by the Development of Goal-Oriented Variational Sensitivity Analysis and A Posteriori Error Estimation Methods

It was demonstrated that a posteriori analyses in general and in particular one that uses adjoint methods can accurately and efficiently compute numerical error estimates and sensitivity for critical Quantities of Interest (QoIs) that depend on a large number of parameters. Activities include: analysis and implementation of several time integration techniques for solving system of ODEs as typically obtained from spatial discretization of PDE systems; multirate integration methods for ordinary differential equations; formulation and analysis of an iterative multi-discretization Galerkin finite element method for multi-scale reaction-diffusion equations; investigation of an inexpensive postprocessing technique to estimate the error of finite element solution of the second-order quasi-linear elliptic problems measured in some global metrics; investigation of an application of the residual-based a posteriori error estimates to symmetric interior penalty discontinuous Galerkin method for solving a class of second order quasi-linear elliptic problems; a posteriori analysis of explicit time integrations for system of linear ordinary differential equations; derivation of accurate a posteriori goal oriented error estimates for a user-defined quantity of interest for two classes of first and second order IMEX schemes for advection-diffusion-reaction problems; Postprocessing finite element solution; and A Bayesian Framework for Uncertain Quantification of Porous Media Flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

An Iterative Optimization Strategy for Nuclear Reactor Design

Optimization of a complex system involves the determination of optimum values for a set of design parameters to meet a specific set of objectives on the quantities of interest (QOI) wherein, the design parameters are a subset of the input parameters, and the QOIs are determined from the output parameters. Particularly, when the parameter space is large, optimization necessitates a significant number of executions of the simulator to obtain a desired solution in tolerance limits. When the simulations are expensive in terms of computation time, an emulator, based on regression methods is useful for predictions. This work presents a novel methodology that uses an iterative hybrid global optimization method (GOM) using genetic algorithms (GA), and simulated annealing (SA) model coupled (HYBGASA) with a Gaussian process regression method based emulator (GPMEM) to optimize a set of input parameters based on a set of defined objectives in a nuclear reactor power system. Hereafter this iterative hybrid method comprising of HYBGASA, and GPMEM would be called as 'IHGOM'. In addition to optimization, IHGOM iteratively updates the trial data, obtained from the neighborhood of the near optimal solution, used to train the GPMEM in order to reduce regression errors. The objective is to develop, model and analyze IHGOM and apply it to an optimization problem in the design of a nuclear reactor. Development and analysis of IHGOM and its implementation in a nuclear reactor power system problem is a significant contribution to the optimization and the nuclear engineering communities. The objective of this work is to develop IHGOM and apply it to an optimization problem in the design of a nuclear reactor. The research involves the development and analysis of a novel optimization strategy, IHGOM. The feasibility analysis involves the optimization of a set of defined parameters in the design and development of a complex system. A GCFBR reactor power system with a set of defined objectives and parameters based on criticality and depletion studies have been characterized as the complex system. Development and analysis of IHGOM and its implementation in an optimization problem in the design of a GCFBR is a significant contribution to the optimization and the nuclear engineering communities. In particular based on the literature review, a coupled GOM and GPMEM method for optimization and regression error reduction in nuclear engineering domain has never been explored. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

MODELING THE Ly α FOREST IN COLLISIONLESS SIMULATIONS

Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining the high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present “Iteratively Matched Statistics” (IMS), a novel method to accurately model the Ly α forest with collisionless N -body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Ly α forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N -body simulation, which we construct from the matter density. We demonstrate that our method can reproduce the PDF, line of sight and 3D power spectra of the Ly α forest with good accuracy (7%, 4%, and 7% respectively). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20%–80%), especially for N -body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic “mock” skies for Ly α forest surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Development of iteration strategies for a practical implementation of a higher order transverse leakage approximation

Transverse integrated nodal diffusion methods currently represent the standard in full core neutronic simulation. The primary shortcoming in this approach is the utilization of the quadratic transverse leakage approximation. This approach, although proven to work well for typical LWR problems, is not consistent with the formulation of nodal methods and can cause accuracy and convergence problems. In previous work, an improved, consistent quadratic leakage approximation was formulated, which derived from the class of higher order nodal methods developed some years ago. In this paper a number of iteration schemes are developed around this consistent quadratic leakage approximation which yield accurate node average results in much improved calculational times. The developed consistent leakage approximation is extended in this work via a number of numerical schemes, the most promising of which results from utilizing the consistent leakage approximation as a correction method to the standard quadratic leakage approximation. Numerical results are demonstrated on a set of benchmark problems, such as the 3D IAEA LWR and MOX C5 problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Communication: Application of state-specific multireference coupled cluster methods to core-level excitations

The concept of the model space underlying multireference coupled-cluster (MRCC) formulations is a powerful tool to deal with complex correlation effects for various electronic states. Here, we demonstrate that iterative state-specific MRCC methods (SS-MRCC) based on properly defined model spaces can be used to describe core-level excited states even when canonical Hartree-Fock orbitals are utilized. We show that the SS-MRCC models with single and double excitations (SSMRCCSD) are comparable in accuracy to high-level single reference equation-of-motion coupled cluster (EOMCC) formalism. We also demonstrate that the SS-MRCC methods are capable of providing high accuracy results without experiencing numerical problems of the EOMCC methods.

coupled cluster↗
DOE OSTI2014

Field of view extension and truncation correction for MR-based human attenuation correction in simultaneous MR/PET imaging

Purpose: In quantitative PET imaging, it is critical to accurately measure and compensate for the attenuation of the photons absorbed in the tissue. While in PET/CT the linear attenuation coefficients can be easily determined from a low-dose CT-based transmission scan, in whole-body MR/PET the computation of the linear attenuation coefficients is based on the MR data. However, a constraint of the MR-based attenuation correction (AC) is the MR-inherent field-of-view (FoV) limitation due to static magnetic field (B{sub 0}) inhomogeneities and gradient nonlinearities. Therefore, the MR-based human AC map may be truncated or geometrically distorted toward the edges of the FoV and, consequently, the PET reconstruction with MR-based AC may be biased. This is especially of impact laterally where the patient arms rest beside the body and are not fully considered. Methods: A method is proposed to extend the MR FoV by determining an optimal readout gradient field which locally compensates B{sub 0} inhomogeneities and gradient nonlinearities. This technique was used to reduce truncation in AC maps of 12 patients, and the impact on the PET quantification was analyzed and compared to truncated data without applying the FoV extension and additionally to an established approach of PET-based FoV extension. Results: The truncation artifacts in the MR-based AC maps were successfully reduced in all patients, and the mean body volume was thereby increased by 5.4%. In some cases large patient-dependent changes in SUV of up to 30% were observed in individual lesions when compared to the standard truncated attenuation map. Conclusions: The proposed technique successfully extends the MR FoV in MR-based attenuation correction and shows an improvement of PET quantification in whole-body MR/PET hybrid imaging. In comparison to the PET-based completion of the truncated body contour, the proposed method is also applicable to specialized PET tracers with little uptake in the arms and might reduce the computation time by obviating the need for iterative calculations of the PET emission data beyond those required for reconstructing images.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-C-201-06: Utility of Quantitative 3D SPECT/CT Imaging in Patient Specific Internal Dosimetry of 153-Samarium with GATE Monte Carlo Package

Purpose: Patient-specific 3-dimensional (3D) internal dosimetry in targeted radionuclide therapy is essential for efficient treatment. Two major steps to achieve reliable results are: 1) generating quantitative 3D images of radionuclide distribution and attenuation coefficients and 2) using a reliable method for dose calculation based on activity and attenuation map. In this research, internal dosimetry for 153-Samarium (153-Sm) was done by SPECT-CT images coupled GATE Monte Carlo package for internal dosimetry. Methods: A 50 years old woman with bone metastases from breast cancer was prescribed 153-Sm treatment (Gamma: 103keV and beta: 0.81MeV). A SPECT/CT scan was performed with the Siemens Simbia-T scanner. SPECT and CT images were registered using default registration software. SPECT quantification was achieved by compensating for all image degrading factors including body attenuation, Compton scattering and collimator-detector response (CDR). Triple energy window method was used to estimate and eliminate the scattered photons. Iterative ordered-subsets expectation maximization (OSEM) with correction for attenuation and distance-dependent CDR was used for image reconstruction. Bilinear energy mapping is used to convert Hounsfield units in CT image to attenuation map. Organ borders were defined by the itk-SNAP toolkit segmentation on CT image. GATE was then used for internal dose calculation. The Specific Absorbed Fractions (SAFs) and S-values were reported as MIRD schema. Results: The results showed that the largest SAFs and S-values are in osseous organs as expected. S-value for lung is the highest after spine that can be important in 153-Sm therapy. Conclusion: We presented the utility of SPECT-CT images and Monte Carlo for patient-specific dosimetry as a reliable and accurate method. It has several advantages over template-based methods or simplified dose estimation methods. With advent of high speed computers, Monte Carlo can be used for treatment planning on a day to day basis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

High Precision In-Pipe Robot Localization with Reciprocal Sensor Fusion - 19516

Manual measurement of U-235 deposits in uranium enrichment piping is a costly, time consuming, and labor-intensive process that is a well-known cost and schedule driver. Autonomous, robotic innovation enabled by this research is revolutionizing the measurement speed, quality and safety of this important operation. The huge advantage is the robot's ability to measure from inside the pipes. The upside is sensing the geometry, appearance and radiometry directly. The downside is the inability to know precise, absolute position of the measurements in very long pipe runs. This paper develops the unprecedented localization required for this purpose. This paper presents the precise localization method designed for in-pipe radiation measurement robots. Unlike other in-pipe localization methods, this approach achieves millimeter-level accuracy in hundred-foot runs and does everything on one robot without relying on tether cable or communicating devices out of pipe, which is an important feature for nuclear in-pipe robot. It overcomes challenges usually encountered by in-pipe localization such as long travelling distance and operation time, narrow view angle, and featureless environment. The robot is equipped with encoders embedded in tracks which can record travelled distance, but encoder odometry drifts unacceptably due to slip and non-linearity. The robot also incorporates a laser range finder measuring the absolute distance, but range finder in long pipes has many misleading false measurements that are difficult to filter. The innovation here is to use one to filter/calibrate the other so that data used from the two sensors improve each other iteratively and reciprocally. The survey measurements accuracy of this localization method is proven by experiments comparing to ground truth and 'Zippering error' experiments that compare measurements to fixed features in pipes. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Six Degrees-of-Freedom Prostate and Lung Tumor Motion Measurements Using Kilovoltage Intrafraction Monitoring

Purpose: Tumor positional uncertainty has been identified as a major issue that deteriorates the efficacy of radiation therapy. Tumor rotational movement, which is not well understood, can result in significant geometric and dosimetric inaccuracies. The objective of this study was to measure 6 degrees-of-freedom (6 DoF) prostate and lung tumor motion, focusing on the more novel rotation, using kilovoltage intrafraction monitoring (KIM). Methods and Materials: Continuous kilovoltage (kV) projections of tumors with gold fiducial markers were acquired during radiation therapy for 267 fractions from 10 prostate cancer patients and immediately before or after radiation therapy for 50 fractions from 3 lung cancer patients. The 6 DoF motion measurements were determined from the individual 3-dimensional (3D) marker positions, after using methods to reject spurious and smooth noisy data, using an iterative closest point algorithm. Results: There were large variations in the magnitude of the tumor rotation among different fractions and patients. Various rotational patterns were observed. The average prostate rotation angles around the left-right (LR), superior-inferior (SI), and anterior-posterior (AP) axes were 1.0 ± 5.0°, 0.6 ± 3.3°, and 0.3 ± 2.0°, respectively. For 35% of the time, the prostate rotated more than 5° about the LR axis, indicating the need for intrafractional adaptation during radiation delivery. For lung patients, the average LR, SI, and AP rotation angles were 0.8 ± 4.2°, −0.8 ± 4.5°, and 1.7 ± 3.1°, respectively. For about 30% of the time, the lung tumors rotated more than 5° around the SI axis. Respiration-induced rotation was detected in 2 of the 3 lung patients. Conclusions: The prostate and lung tumors were found to undergo rotations of more than 5° for about a third of the time. The lung tumor data represent the first 6 DoF tumor motion measured by kV images. The 6 DoF KIM method can enable rotational and translational adaptive radiation therapy and potentially reduce treatment margins.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Neutronics and Depletion Methods for Parametric Studies of Fluoride Salt Cooled High Temperature Reactors with Slab Fuel Geometry and Multi-Batch Fuel Management Schemes

The Advanced High Temperature Reactor (AHTR) is a 3400 MWth fluoride salt cooled high temperature reactor (FHR) that uses TRISO particle fuel compacted into slabs rather than spherical fuel pebbles or cylindrical fuel compacts. Simplified methods are required for parametric design studies such that analyzing the entire feasible design space for an AHTR is tractable. These simplifications include fuel homogenization techniques to increase the speed of neutron transport calculations in depletion analysis and equilibrium depletion analysis methods to analyze systems with multi-batch fuel management schemes. This paper presents three elements of significant novelty. First, the reactivity-equivalent physical transformation (RPT) methodology usually applied in systems with coated particle fuel in cylindrical and spherical geometries was extended to slab geometries. Secondly, based on this newly developed RPT method for slab geometries, a methodology that uses Monte Carlo depletion approaches was further developed to search for the maximum discharge burnup in a multi-batch system by iteratively estimating the beginning of equilibrium cycle composition and sampling different discharge burnups. This iterative equilibrium depletion search (IEDS) method fully defines an equilibrium fuel cycle (keff, power, flux and composition evolutions across space and time), but is computationally demanding, although feasible on single-processor workstations. Finally, an analytical method, the non-linear reactivity model, was developed by expanding the linear reactivity model to include an arbitrary number of higher order terms to extrapolate single-batch depletion results to estimate the maximum discharge burnup and BOEC keff in systems with multi-batch fuel management schemes. Results from this method were benchmarked against equilibrium depletion analysis results using the IEDS.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗