Search NASA⌕ Search

SEARCH · Search NASA

Results for “error fields”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Motion tracking in the liver: Validation of a method based on 4D ultrasound using a nonrigid registration technique

Purpose: Treatments like radiotherapy and focused ultrasound in the abdomen require accurate motion tracking, in order to optimize dosage delivery to the target and minimize damage to critical structures and healthy tissues around the target. 4D ultrasound is a promising modality for motion tracking during such treatments. In this study, the authors evaluate the accuracy of motion tracking in the liver based on deformable registration of 4D ultrasound images. Methods: The offline analysis was performed using a nonrigid registration algorithm that was specifically designed for motion estimation from dynamic imaging data. The method registers the entire 4D image data sequence in a groupwise optimization fashion, thus avoiding a bias toward a specifically chosen reference time point. Three healthy volunteers were scanned over several breathing cycles (12 s) from three different positions and angles on the abdomen; a total of nine 4D scans for the three volunteers. Well-defined anatomic landmarks were manually annotated in all 96 time frames for assessment of the automatic algorithm. The error of the automatic motion estimation method was compared with interobserver variability. The authors also performed experiments to investigate the influence of parameters defining the deformation field flexibility and evaluated how well the method performed with a lower temporal resolution in order to establish the minimum frame rate required for accurate motion estimation. Results: The registration method estimated liver motion with an error of 1 mm (75% percentile over all datasets), which was lower than the interobserver variability of 1.4 mm. The results were only slightly dependent on the degrees of freedom of the deformation model. The registration error increased to 2.8 mm with an eight times lower temporal resolution. Conclusions: The authors conclude that the methodology was able to accurately track the motion of the liver in the 4D ultrasound data. The authors believe that the method has potential in interventions on moving abdominal organs such as MR or ultrasound guided focused ultrasound therapy and radiotherapy, pending the method is enabled to run in real-time. The data and the annotations used for this study are made publicly available for those who would like to test other methods on 4D liver ultrasound data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Measurement of the electron antineutrino mass in tritium beta decay in the Troitsk nu-mass experiment

The results obtained in the Troitsk nu-mass experiment by measuring the electron-antineutrino mass in tritium beta decay are presented. The facility used consists of a gaseous windowless tritium source and an electrostatic electron spectrometer involving an adiabatic magnetic collimation. Runs in which measurement conditions were reliably established were thoroughly selected in analyzing data obtained from 1994 to 2004. All known systematic effects were taken into account. For the square of the electron-antineutrino mass, the treatment of measured spectra yielded the following result: m{sub {nu}}{sup 2} = -0.67 {+-} 1.89{sub stat.} {+-} 1.68{sub syst.} eV{sup 2}. The use of the Bayesian method and the Feldman-Cousins unified approach made it possible to obtain the following upper limits on the mass: m{sub {nu}} < 2.12 eV (at a 95% C.L.; Bayesian method) and m{sub {nu}} < 2.05 eV (at a 95% C.L., Feldman-Cousins method). At the same time, an estimation of the sensitivity limit without allowance for negative values of the square of the mass leads to m{sub {nu}} < 2.2 eV (at a 95% C.L.). Measured spectra were analyzed for the possible existence of an additional structure (step) in the electron spectrum near the boundary energy. The conclusion drawn from this analysis was that, within the existing statistical errors, there are no reasons for introducing such a feature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Evaluation of three lidar scanning strategies for turbulence measurements

Abstract. Several errors occur when a traditional Doppler beam swinging (DBS) or velocity–azimuth display (VAD) strategy is used to measure turbulence with a lidar. To mitigate some of these errors, a scanning strategy was recently developed which employs six beam positions to independently estimate the u, v, and w velocity variances and covariances. In order to assess the ability of these different scanning techniques to measure turbulence, a Halo scanning lidar, WindCube v2 pulsed lidar, and ZephIR continuous wave lidar were deployed at field sites in Oklahoma and Colorado with collocated sonic anemometers.Results indicate that the six-beam strategy mitigates some of the errors caused by VAD and DBS scans, but the strategy is strongly affected by errors in the variance measured at the different beam positions. The ZephIR and WindCube lidars overestimated horizontal variance values by over 60 % under unstable conditions as a result of variance contamination, where additional variance components contaminate the true value of the variance. A correction method was developed for the WindCube lidar that uses variance calculated from the vertical beam position to reduce variance contamination in the u and v variance components. The correction method reduced WindCube variance estimates by over 20 % at both the Oklahoma and Colorado sites under unstable conditions, when variance contamination is largest. This correction method can be easily applied to other lidars that contain a vertical beam position and is a promising method for accurately estimating turbulence with commercially available lidars.

Newman, Jennifer F.↗
DOE OSTI2012

A Beam-Specific Planning Target Volume (PTV) Design for Proton Therapy to Account for Setup and Range Uncertainties

Purpose: To report a method for explicitly designing a planning target volume (PTV) for treatment planning and evaluation in heterogeneous media for passively scattered proton therapy and scanning beam proton therapy using single-field optimization (SFO). Methods and Materials: A beam-specific PTV (bsPTV) for proton beams was derived by ray-tracing and shifting ray lines to account for tissue misalignment in the presence of setup error or organ motion. Range uncertainties resulting from inaccuracies in computed tomography-based range estimation were calculated for proximal and distal surfaces of the target in the beam direction. The bsPTV was then constructed based on local heterogeneity. The bsPTV thus can be used directly as a planning target as if it were in photon therapy. To test the robustness of the bsPTV, we generated a single-field proton plan in a virtual phantom. Intentional setup and range errors were introduced. Dose coverage to the clinical target volume (CTV) under various simulation conditions was compared between plans designed based on the bsPTV and a conventional PTV. Results: The simulated treatment using the bsPTV design performed significantly better than the plan using the conventional PTV in maintaining dose coverage to the CTV. With conventional PTV plans, the minimum coverage to the CTV dropped from 99% to 67% in the presence of setup error, internal motion, and range uncertainty. However, plans using the bsPTV showed minimal drop of target coverage from 99% to 94%. Conclusions: The conventional geometry-based PTV concept used in photon therapy does not work well for proton therapy. We investigated and validated a beam-specific PTV method for designing and evaluating proton plans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗
DOE OSTI2011

Constraints on the four-generation quark mixing matrix from a fit to flavor-physics data

In the scenario with four quark generations, we perform a fit using flavor-physics data and determine the allowed values - preferred central values and errors - of all of the elements of the 4x4 quark-mixing matrix. In addition to the direct measurements of some of the elements, we include in the fit the present measurements of several flavor-changing observables in the K and B systems that have small hadronic uncertainties, and also consider the constraints from the vertex corrections to Z{yields}bb. The values taken for the masses of the fourth-generation quarks are consistent with the measurements of the oblique parameters and perturbativity of the Yukawa couplings. We find that |V-tilde{sub tb}|{>=}0.98 at 3{sigma}, so that a fourth generation cannot account for any large deviation of |V{sub tb}| from unity. The fit also indicates that all the new-physics parameters are consistent with zero, and the mixing of the fourth generation with the other three is constrained to be small: we obtain |V-tilde{sub ub}{sup '}|<0.06, |V-tilde{sub cb}{sup '}|<0.027, and |V-tilde{sub tb}{sup '}|<0.31 at 3{sigma}. Still, this does allow for the possibility of new-physics signals in B{sub d}, B{sub s} and rare K decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Intelligent Process Visualization through Nuclear Operation Process Modeling, Reasoning, and Object Detection from Field Videos (Final Report)

This report is a deliverable for the “Final Report” task of DOE NEET Project 19-16790, "Context-Aware Safety Information Display for Nuclear Field Workers." This project's overall goal is to test the hypothesis that integrating computer vision and process reasoning methods will enable proactive visualization of the safe operation and maintenance processes of Nuclear Power Plants (NPP) for field workers. Augmented Reality (AR) glasses adopting such proactive safety information visualization techniques can significantly increase personnel safety and reduce the NPP’s operating costs. The current practice of monitoring NPPs requires workers to switch between digital models, data, and physical workspaces in identifying relevant but potentially occluded objects and in assessing the risks of operation and maintenance processes. On the other hand, frequently changed field conditions require field workers to report to supervisors for real-time guidance. Such guidance is essential to ensure that changing conditions will not invalidate or endanger the work order and other ongoing processes that may jeopardize NPP operations. Additionally, incorrect recognition of equipment objects can result in communication errors and safety problems. AR techniques can assist engineers in viewing the physical workspaces with objects labeled with detailed operation procedures and safety reminders during field operations. The project team developed an “Intelligent Context-Aware Safety Information Display” (ICAD) for supporting Nuclear Power Plant (NPP) field workers in achieving safe and efficient execution of a series of operational tasks in uncertain and changing workspaces of an NPP. Before designing the ICAD-AR prototype, the project team synthesized NPP operational knowledge models through literature review studies, surveys, interviews with domain experts, and knowledge modeling. The project team conducted an extensive study of the operational procedures of various NPPs, and digital technologies that can support the safe and efficient execution of those procedures in different NPP operational contexts. This literature review helped the project team conduct surveys and interviews with nuclear engineers and field workers to identify three categories of information. The NPP knowledge modeling efforts reveal that the three categories of information identified have different levels of importance in a typical procedure of carrying out a series of tasks to achieve a specific NPP operation goal (e.g., shutdown, mode changes). These three categories of information include 1) Workspace dynamics – the changing spatial arrangements of workspaces, tools, protection equipment, and supporting materials, 2) Workflow prognostics – the dynamic dependencies between different parts of an NPP that functionally support and influence each other in terms of safety and efficiency, and 3) Hazards – objects and spaces that contain hazardous materials or physical conditions that can pose risks to workers or mechanical systems. The project team has profiled the importance levels of these categories of information into a knowledge model. This knowledge model specifies what types of information are more critical for a given task in a given workspace so that computers can automatically identify critical objects and sensors in a scene for delivering context-ware safety information to field workers through AR devices. Significant research development of this project results in technical research outcomes and a prototyping system that illustrates the technical feasibility of establishing an ICAD-AR system supporting the proactive safety information display for nuclear field workers. This final report summarizes the project team’s technological achievements in the past three years. Overall, the project team completed the development and integration of five techniques into a prototype ICAD Augmented Reality (ICAD-AR) system and demonstrated the developed system’s real-time execution in a mechanical room. The project team completed the analysis of using this prototype in other types of workspaces based on 3D image data and digital design models collected from two additional workspaces (a water treatment plant and a flow loop training facility). The integrated techniques include 1) Natural Language Processing (NLP) algorithms supporting the generation and updates of nuclear fieldwork process models based on text analysis of work packages and operation manuals; 2) sensor log analysis for predicting control actions in given sensor reading contexts; 3) computer vision algorithms for automatic localization and navigation of workers; 4) object detection algorithms for identifying task-related objects and correlated sensors for safety checking; 5) AR technique as a platform for supporting the integration. The testing results of these five techniques have shown that 1) the sensor log analysis model can predict the next control action with an accuracy of 0.883; 2) the trained natural language processing model can extract more than 80% of the critical information from paper-based procedures (PBPs); 3) the navigation algorithm with the integration of Visual Inertial Odometry (VIO) and Non-Recursive Bayesian Filter methods make operator’s trajectory estimation resilient to drift error; 4) the computer vision algorithm can detect task-specific and safety-critical objects with an average accuracy of 95.3%. The project team used work procedures collected from a flow loop training facility and two datasets collected from two mechanical rooms simulating the workspaces of NPPs to demonstrate the technical capabilities of the developed ICAD-AR prototype. The demonstration validated the technical feasibility of establishing the ICAD-AR system for nuclear field workers and identified the challenges in 1) automatic text analysis of work packages; 2) use of limited samples of sensor logs for predicting the proper timings of control actions; 3) reliably tracking workers and their task progress in mechanical rooms with many similar objects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Code-verification techniques for the method-of-moments implementation of the magnetic-field integral equation

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementation of the magnetic-field integral equation, are frequently used to solve Maxwell's equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches for cases with and without coding errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Superadditivity at large charge

The weak gravity conjecture has been invoked to conjecture that the dimensions of charged operators in a CFT should obey a superadditivity relation (sometimes referred to as convexity). In this paper, we study superadditivity of the operator spectrum in theories expanded about the semi-classical saddle point that dominates correlators of large charge operators. We explore this in two contexts. The first is a model with two scalar fields that carry different charges, at a non-trivial Wilson-Fisher fixed point. A careful analysis of the semi-classics for this two field model demonstrates that ‘quantum’ violations of superadditivity (those not forbidden by the conjecture) persist in the large charge regime. We then turn to study the general properties of CFTs at large charge as bottom-up EFTs. By a trial and error procedure we come up with a seemingly consistent family of examples violating the conjecture. In so doing the presence of a genuine dilaton field appears necessary. On the one hand our result demonstrates that the superadditivity conjecture cannot be proven purely on the basis of a bottom-up analysis. On the other hand, the need for a dilaton, with the corresponding infinite fine tuning, indicates the conjecture-violating EFTs are unlikely to be UV completable.

effective field theories↗
DOE OSTI2011

Learning statistical correlation for fast prostate registration in image-guided radiotherapy

Purpose: In adaptive radiation therapy of prostate cancer, fast and accurate registration between the planning image and treatment images of the patient is of essential importance. With the authors' recently developed deformable surface model, prostate boundaries in each treatment image can be rapidly segmented and their correspondences (or relative deformations) to the prostate boundaries in the planning image are also established automatically. However, the dense correspondences on the nonboundary regions, which are important especially for transforming the treatment plan designed in the planning image space to each treatment image space, are remained unresolved. This paper presents a novel approach to learn the statistical correlation between deformations of prostate boundary and nonboundary regions, for rapidly estimating deformations of the nonboundary regions when given the deformations of the prostate boundary at a new treatment image. Methods: The main contributions of the proposed method lie in the following aspects. First, the statistical deformation correlation will be learned from both current patient and other training patients, and further updated adaptively during the radiotherapy. Specifically, in the initial treatment stage when the number of treatment images collected from the current patient is small, the statistical deformation correlation is mainly learned from other training patients. As more treatment images are collected from the current patient, the patient-specific information will play a more important role in learning patient-specific statistical deformation correlation to effectively reflect prostate deformation of the current patient during the treatment. Eventually, only the patient-specific statistical deformation correlation is used to estimate dense correspondences when a sufficient number of treatment images have been acquired from the current patient. Second, the statistical deformation correlation will be learned by using a multiple linear regression (MLR) model, i.e., ridge regression (RR) model, which has the best prediction accuracy than other MLR models such as canonical correlation analysis (CCA) and principal component regression (PCR). Results: To demonstrate the performance of the proposed method, we first evaluate its registration accuracy by comparing the deformation field predicted by our method with the deformation field estimated by the thin plate spline (TPS) based correspondence interpolation method on 306 serial prostate CT images of 24 patients. The average predictive error on the voxels around 5 mm of prostate boundary is 0.38 mm for our method of RR-based correlation model. Also, the corresponding maximum error is 2.89 mm. We then compare the speed for deformation interpolation by different methods. When considering the larger region of interest (ROI) with the size of 512 x 512 x 61, our method takes 24.41 seconds to interpolate the dense deformation field while TPS method needs 6.7 minutes; when considering a small ROI (surrounding prostate) with size of 112 x 110 x 93, our method takes 1.80 seconds, while TPS method needs 25 seconds. Conclusions: Experimental results show that the proposed method can achieve much faster registration speed yet with comparable registration accuracy, compared to the TPS-based correspondence (or deformation) interpolation approach.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Determination of the beam quality index of high-energy photon beams under nonstandard reference conditions

Purpose: At some modern radiotherapy machines it is not possible to achieve reference conditions for the measurement of beam quality indices used in dosimetry codes of practice, such as IAEA TRS-398 and AAPM TG-51. This work aims at providing self-consistent and simpler expressions and more accurate fits for a limited range of beams of interest than have been proposed previously for deriving these beam quality indices from measurements. Methods: The starting point is a formula proposed by Sauer [Med. Phys. 36, 4168-4172 (2009)] for deriving the beam quality index used in IAEA TRS-398, TPR{sub 20,10}, from a measurement of the tissue phantom ratio at depths of 20 cm and 10 cm in water for an s Multiplication-Sign s cm{sup 2} (equivalent) square field, TPR{sub 20,10}(s). First, a self-consistent version of this formula is established followed by a simpler version by making a linear approximation. A similar approach is proposed to derive the beam quality index used in AAPM TG-51, %dd(10){sub X}, from a measurement of PDD{sub 10}(s), the percentage depth dose at 10 cm for a square field with size s. All models were fitted to subsets of relevant data from BJR supplement 25. Results: The linear models for TPR{sub 20,10}(s) and exponential models for PDD{sub 10}(s) as a function of the (equivalent) square field size can reproduce the beam quality within 0.3% and beam quality correction factors within 0.05% for square field sizes ranging from 4 cm to 12 cm and nominal photon energies from 4 MV to 12 MV. For higher energy beams the errors are only slightly worse but for %dd(10){sub X}, an additional uncertainty component has to be considered for the electron contamination correction. Conclusions: The models proposed here can be used in practical recommendations for the dosimetry of small and nonstandard fields.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗
DOE OSTI2016

SU-F-T-313: Clinical Results of a New Customer Acceptance Test for Elekta VMAT

Purpose: To report the results of a customer acceptance test (CAT) for VMAT treatments for two matched Elekta linear accelerators. Methods: The CAT tests were performed on two clinically matched Elekta linear accelerators equipped with a 160-leaf MLC. Functional tests included performance checks of the control system during dynamic movements of the diaphragms, MLC, and gantry. Dosimetric tests included MLC picket fence tests at static and variable dose rates and a diaphragm alignment test, all performed using the on-board EPID. Additionally, beam symmetry during arc delivery was measured at the four cardinal angles for high and low dose rate modes using a 2D detector array. Results of the dosimetric tests were analyzed using the VMAT CAT analysis tool. Results: Linear accelerator 1 (LN1) met all stated CAT tolerances. Linear accelerator 2 (LN2) passed the geometric, beam symmetry, and MLC position error tests but failed the relative dose average test for the diaphragm abutment and all three picket fence fields. Though peak doses in the abutment regions were consistent, the average dose was below the stated tolerance corresponding to a leaf junction that was too narrow. Despite this, no significant differences in patient specific VMAT quality assurance measured were observed between the accelerators and both passed monthly MLC quality assurance performed with the Hancock test. Conclusion: Results from the CAT showed LN2 with relative dose averages in the abutment regions of the diaphragm and MLC tests outside the tolerances resulting from differences in leaf gap distances. Tolerances of the dose average tests from the CAT may be small enough to detect MLC errors which do not significantly affect patient QA or the routine MLC tests.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Rapid multiphase-field model development using a modular free energy based approach with automatic differentiation in MOOSE/MARMOT

In this paper, we present a novel rapid phase-field model development capability using a modular free energy based approach in the finite element based MOOSE framework. The phase-field equations are implemented in a general form as separate code form the model specific data such as the thermodynamic free energy density and mobility functions. This isolates researchers from complicated code and allows them to focus on implementing new material models. In addition, using a native symbolic differentiation algorithm, all required derivatives of the free energy density are automatically determined, removing sources for potential human error, and guaranteeing optimally preconditioned equation systems through fully populated Jacobian matrices. Free energy terms contributing to a phase-field model are abstracted into self contained objects that can be combined at will at simulation run time. Combining multiple chemical and mechanical free energy contributions allows the construction of coupled phase-field, mechanics, and multiphase models. Through just-in-time compilation, we greatly reduce the computational time required to evaluate the differentiated expressions. Lastly, we demonstrate the new capability on a variety of example simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
DOE OSTI2015

WE-D-BRA-05: Pseudo In Vivo Patient Dosimetry Using a 3D-Printed Patient-Specific Phantom

Purpose: To test the feasibility of using 3D-printed patient-specific phantoms for intensity-modulated radiation therapy (IMRT) quality assurance (QA). Methods: We created a patient-specific whole-head phantom using a 3D printer. The printer data file was created from high-resolution DICOM computed tomography (CT) images of 3-year old child treated at our institution for medulloblastoma. A custom-modified extruder system was used to create tissue-equivalent materials. For the printing process, the Hounsfield Units from the CT images were converted to proportional volumetric densities. A 5-field IMRT plan was created from the patient CT and delivered to the 3D- phantom. Dose was measured by an ion chamber placed through the eye. The ion chamber was placed at the posterior edge of the planning target volume in a high dose gradient region. CT scans of the patient and 3D-phantom were fused by using commercial treatment planning software (TPS). The patient’s plan was calculated on the phantom CT images. The ion chamber’s active volume was delineated in the TPS; dose per field and total dose were obtained. Measured and calculated doses were compared. Results: The 3D-phantom dimensions and tissue densities were in good agreement with the patient. However, because of a printing error, there was a large discrepancy in the density in the frontal cortex. The calculated and measured treatment plan doses were 1.74 Gy and 1.72 Gy, respectively. For individual fields, the absolute dose difference between measured and calculated values was on average 3.50%. Conclusion: This study demonstrated the feasibility of using 3D-printed patient-specific phantoms for IMRT QA. Such phantoms would be particularly advantageous for complex IMRT treatment plans featuring high dose gradients and/or for anatomical sites with high variation in tissue densities. Our preliminary findings are promising. We anticipate that, once the printing process is further refined, the agreement between measured and calculated doses will improve.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2010

BOOTSTRAPPING THE CORONAL MAGNETIC FIELD WITH STEREO: UNIPOLAR POTENTIAL FIELD MODELING

We investigate the recently quantified misalignment of {alpha}{sub mis} {approx} 20{sup 0}-40{sup 0} between the three-dimensional geometry of stereoscopically triangulated coronal loops observed with STEREO/EUVI (in four active regions (ARs)) and theoretical (potential or nonlinear force-free) magnetic field models extrapolated from photospheric magnetograms. We develop an efficient method of bootstrapping the coronal magnetic field by forward fitting a parameterized potential field model to the STEREO-observed loops. The potential field model consists of a number of unipolar magnetic charges that are parameterized by decomposing a photospheric magnetogram from the Michelson Doppler Imager. The forward-fitting method yields a best-fit magnetic field model with a reduced misalignment of {alpha}{sub PF} {approx} 13{sup 0}-20{sup 0}. We also evaluate stereoscopic measurement errors and find a contribution of {alpha}{sub SE} {approx} 7{sup 0}-12{sup 0}, which constrains the residual misalignment to {alpha}{sub NP} {approx} 11{sup 0}-17{sup 0}, which is likely due to the nonpotentiality of the ARs. The residual misalignment angle, {alpha}{sub NP}, of the potential field due to nonpotentiality is found to correlate with the soft X-ray flux of the AR, which implies a relationship between electric currents and plasma heating.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗