Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “volume estimation”

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 811 records · Page 45

Constraining the Thickness of Polar Ice Deposits on Mercury Using the Mercury Laser Altimeter and Small Craters in Permanently Shadowed Regions

Radar-bright deposits at the poles of Mercury are located in permanently shadowed regions, which pro- vide thermally stable environments for hosting and retaining water ice on the surface or in the near sub- surface for geologic timescales. While the areal distribution of these radar-bright deposits is well characterized, their thickness, and thus their total mass and volume, remain poorly constrained. Here we derive thickness estimates for selected water-ice deposits using small, simple craters visible within the permanently shadowed, radar-bright deposits. We examine two endmember scenarios: in Case I, these craters predate the emplacement of the ice, and in Case II, these craters postdate the emplacement of the ice. In Case I, we find the difference between estimated depths of the original unfilled craters and the measured depths of the craters to find the estimated infill of material. The average estimated infilled material for 9 craters assumed to be overlain with water ice is ∼41 +30 −14 m, where 1 −σ standard error of the mean is reported as uncertainty. Reported uncertainties are for statistical errors only. Additional systematic uncertainty may stem from georeferencing the images and topographic datasets, from the radial accuracy of the altimeter measurements, or from assumptions in our models including (1) ice is flat in the bowl- shaped crater and (2) there is negligible ice at the crater rims. In Case II, we derive crater excavation depths to investigate the thickness of the ice layer that may have been penetrated by the impact. While the absence of excavated regolith associated with the small craters observed suggests that impacts generally do not penetrate through the ice deposit, the spatial resolution and complex illumination geometry of images may limit the observations. Therefore, it is not possible to conclude whether the small craters in this study penetrate through the ice deposit, and thus Case II does not provide a constraint on the ice thickness. For Mercury’s polar deposits, we argue that Case I of the small craters predating the emplacement of the ice deposits is more likely, given other geologic evidence that suggests that these ice deposits are relatively young. Using the ice thickness estimates from Case I to calculate the total amount of water ice currently contained in Mercury’s polar deposits results in a value of ∼10 14 –10 15 kg. This is equivalent to ∼100–1000 km 3 ice in volume. This volume of water ice is consistent with delivery via micrometeorite bombardment, Jupiter-family comets, or potentially a single impactor.

Deutsch, Ariel N.↗

What's in WIPP? Packaging TRU Waste to Enhance WIPP's Capacity - 16515

Beginning in 1999, and continuing for 15 years, the US DOE Waste Isolation Pilot Plant (WIPP) repository successfully disposed of defense Transuranic (TRU) waste. Six of the originally planned ten disposal panels were filled before two unrelated accidents temporarily halted disposal operations in February 2014. Based on counting the gross volume of all containers emplaced in WIPP to date, a little more than half of the legislated volume capacity has been subscribed (91,000 cubic meters). This paper describes the many different packaging configurations that have been used to ship TRU waste to WIPP, and makes observations that may enhance the repository capacity when shipments resume. The contents of shipments to WIPP are limited by many requirements of the repository waste acceptance criteria. These criteria derive from three primary sources; 1) the Certificate of Compliance (CoC), issued by the NRC for Type B shipping packaging used by WIPP; 2) The Hazardous Waste Facility Permit, issued by the New Mexico Environment Department; and 3) the WIPP facility's own Documented Safety Analysis (DSA), developed under DOE's regulation 10CFR830. To meet these criteria, various TRU waste forms have been packaged in different ways, with varying packaging efficiency. Of course, there is an overall upper limit on WIPP's volume capacity that was legislated as part of the WIPP Land Withdrawal Act. The volume counted against this limit has been calculated by simply summing the volumes of individual payload containers shipped to and emplaced at WIPP. This method counts significant non-waste volumes inside the containers due to over-packing inner containers and partially filled containers. Without some change to how disposal volumes are calculated, or how waste is packaged, the amount of TRU waste yet to be shipped to WIPP could exceed the legislated capacity. Concepts to enhance the packaging efficiency of future waste shipments to WIPP are presented. Approximate estimates of the waste versus non-waste volumes that have been emplaced to date are discussed. This paper only considers contact handled TRU waste (i.e., contact dose rates < 2 mSv/hr). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparative Study of Hydrology and Icemelt in Three Nepal River Basins Using the Glacio-Hydrological Degree-Day Model (GDM) and Observations From the Advanced Scatterometer (ASCAT)

An assessment of the water supply and its seasonal and annual changes over the century in the High Mountain Asia (HMA) region is of increasing interest due to its potential impact on one-sixth of the global population. In order to understand the changing hydrology and snow and ice melt, we used remotely sensed Advanced Scatterometer (ASCAT) observations of glacier melt (GM) and a distributed and gridded Glacio-hydrological Degree-day Model (GDM) in three river basins: Tamor, Trishuli and Marsyangdi. The GDM-estimated contribution of snowmelt, icemelt, rainfall and baseflow in river flows is found to be most accurate in the Trishuli River basin, with Nash-Sutcliffe efficiency (NSE) between the estimated and observed discharges of 0.81 and volume differences of -0.5%, and reasonably accurate in the Tamor River basin, with NSE of 0.69 and volume difference of -7.51%. Similarly, NSE of 0.81 and volume difference of 4.64% in Marsyangdi River basin. We find strong similarities in the timing of glacier melting using the GDM and from observations from the ASCAT GM, determining the seasonal start of glacier melting to within 6 days on average. In all basins ASCAT GM observes melting at higher elevations relative to GDM, average of 5,328 m a.s.l. Systematic differences in glacier melting area determined by modeling and satellite observations indicate ASCAT may have suboptimal resolution, view geometry and/or polarimetry for delineating glacier melting at the process-scale in complex topography, especially in the ablation zone. This is the first step in examining the remote sensing products that could potentially be incorporated into hydrologic models to increase the accuracy of the hydrologic flow as well as the ability to estimate river discharge in other basins with limited data.

54 ENVIRONMENTAL SCIENCES↗

Determining the IV fluids required for a ten day medical emergency on Space Station Freedom - Comparison of packaged vs. on-orbit produced solutions

To aid planning for the storage of supplies onboard Space Station Freedom, an estimate was made of the amount of intravenous (IV) fluid required to support a patient who has suffered a medical emergency for a period of up to 10 days. Six different medical scenarios were evaluated, and the volume of IV fluids required for each scenario was estimated. Up to 220 liters of fluid would be required to support a patient for all of the scenarios. When optimizing the volumes to support any single scenario, a total of 123 liters is required. Use of a water polishing system to produce sterile water for injection from potable supplies and on-station formulation of IV fluids results in a smaller mass and volume requirement for the Fluid Therapy Subsystem than carrying prepackaged bags of fluid.

Creager, Gerald J.↗

CHIRP-Like Signals: Estimation, Detection and Processing A Sequential Model-Based Approach

Chirp signals have evolved primarily from radar/sonar signal processing applications specifically attempting to estimate the location of a target in surveillance/tracking volume. The chirp, which is essentially a sinusoidal signal whose phase changes instantaneously at each time sample, has an interesting property in that its correlation approximates an impulse function. It is well-known that a matched-filter detector in radar/sonar estimates the target range by cross-correlating a replicant of the transmitted chirp with the measurement data reflected from the target back to the radar/sonar receiver yielding a maximum peak corresponding to the echo time and therefore enabling the desired range estimate. In this application, we perform the same operation as a radar or sonar system, that is, we transmit a “chirp-like pulse” into the target medium and attempt to first detect its presence and second estimate its location or range. Our problem is complicated by the presence of disturbance signals from surrounding broadcast stations as well as extraneous sources of interference in our frequency bands and of course the ever present random noise from instrumentation. First, we discuss the chirp signal itself and illustrate its inherent properties and then develop a model-based processing scheme enabling both the detection and estimation of the signal from noisy measurement data.

42 ENGINEERING↗

Technical challenges of sparing infrahyoid swallowing organs at risk in oropharynx squamous cell cancer treated with IMRT

This study reports clinical performance in the sparing of infrahyoid swallowing organs at risk (SWOARs) in oropharynx cancer intensity-modulated radiation therapy (IMRT) plans. Rates of meeting dose-volume planning goals are reported and compared with geometry-based estimates of what is achievable. This study also develops 3 measures of target-SWOAR geometry and tests their usefulness in providing geometry-based dose-volume planning goals. A total of 50 oropharynx cancer IMRT plans were reviewed. Success rates in meeting institutional dose-volume goals were determined for the glottic larynx (G), postcricoid pharynx (P), and esophagus (E). The following 3 measures of target-SWOAR geometry were investigated as methods of identifying geometry-based planning goals: presence of gross disease in neck levels 3 to 4, target-SWOAR overlap, and a 3-dimensional (3D) measure of target-SWOAR geometry. Locally advanced disease was predominant in this patient population with target volumes overlapping SWOARs in 68% to 98% of cases. Clinical rates of success in meeting dose-volume goals varied by SWOAR (16% to 82%) but compared well with estimated potentially achievable rates in most cases (14% average difference between clinical and potential). Cases grouped by the presence of levels 3 to 4 neck nodes or target-SWOAR overlap did not have significantly different SWOAR doses. Cases grouped using a 3D measure of target-SWOAR geometry differed significantly, providing useful geometry-based planning goals (e.g., mean Glottis dose <45 Gy was achieved 19%, 44%, or 81% of the time in each of 3 groups). This study describes the technical challenge of sparing SWOARs and investigates several potential methods for grouping cases to assist with treatment plan evaluation. Quantifying the 3-D relationship between the targets and SWOARs is a promising way of approaching this complex problem. Data presented in this paper may be useful to evaluate treatment plans using objective geometry-based goals.

61 RADIATION PROTECTION AND DOSIMETRY↗

Tris (isodecyl) guanidine degradation in the MCU system

The current solvent blend consists of four components; an extractant, the modifier, a suppressor, and the diluent. Of the four components, only the suppressor – tris(isodecyl)guanidine (TiDG) has exhibited an appreciable depletion rate during facility operations.∇ Using data derived from MCU process samples, SRNL derived a method to predict the TiDG depletion based upon time and volume of feed processed. With the current data set, the following formula can estimate the TiDG concentrations after processing: A t = [A 0 ×exp(-0.000900×t)] – [5.39E-06×volume] where A t is the TiDG concentration (mM) at time t (days), A 0 is the TiDG concentration (mM) at time = 0, and volume is the amount of salt solution processed in gallons from time 0 to time t. The ability to use this formula as a predictive tool is limited due to the number of data points obtained for this scope of work. As such, this formula should not be used to precisely predict future TiDG concentrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TU-H-207A-08: Estimating Radiation Dose From Low-Dose Lung Cancer Screening CT Exams Using Tube Current Modulation

Purpose: The purpose of this work is to estimate effective and lung doses from a low-dose lung cancer screening CT protocol using Tube Current Modulation (TCM) across patient models of different sizes. Methods: Monte Carlo simulation methods were used to estimate effective and lung doses from a low-dose lung cancer screening protocol for a 64-slice CT (Sensation 64, Siemens Healthcare) that used TCM. Scanning parameters were from the AAPM protocols. Ten GSF voxelized patient models were used and had all radiosensitive organs identified to facilitate estimating both organ and effective doses. Predicted TCM schemes for each patient model were generated using a validated method wherein tissue attenuation characteristics and scanner limitations were used to determine the TCM output as a function of table position and source angle. The water equivalent diameter (WED) was determined by estimating the attenuation at the center of the scan volume for each patient model. Monte Carlo simulations were performed using the unique TCM scheme for each patient model. Lung doses were tallied and effective doses were estimated using ICRP 103 tissue weighting factors. Effective and lung dose values were normalized by scanspecific 32 cm CTDIvol values based upon the average tube current across the entire simulated scan. Absolute and normalized doses were reported as a function of WED for each patient. Results: For all ten patients modeled, the effective dose using TCM protocols was below 1.5 mSv. Smaller sized patient models experienced lower absolute doses compared to larger sized patients. Normalized effective and lung doses showed some dependence on patient size (R2 = 0.77 and 0.78, respectively). Conclusion: Effective doses for a low-dose lung screening protocol using TCM were below 1.5 mSv for all patient models used in this study. Institutional research agreement, Siemens Healthcare; Past recipient, research grant support, Siemens Healthcare; Consultant, Toshiba America Medical Systems; Consultant, Samsung Electronics.

60 APPLIED LIFE SCIENCES↗

Comparison between target margins derived from 4DCT scans and real-time tumor motion tracking: Insights from lung tumor patients treated with robotic radiosurgery

Purpose: A unique capability of the CyberKnife system is dynamic target tracking. However, not all patients are eligible for this approach. Rather, their tumors are tracked statically using the vertebral column for alignment. When using static tracking, the internal target volume (ITV) is delineated on the four-dimensional (4D) CT scan and an additional margin is added to account for setup uncertainty [planning target volume (PTV)]. Treatment margins are difficult to estimate due to unpredictable variations in tumor motion and respiratory pattern during the course of treatment. The inability to track the target and detect changes in respiratory characteristics might result in geographic misses and local tumor recurrences. The purpose of this study is to develop a method to evaluate the adequacy of ITV-to-PTV margins for patients treated in this manner. Methods: Data from 24 patients with lesions in the upper lobe (n = 12), middle lobe (n = 3), and lower lobe (n = 9) were included in this study. Each patient was treated with dynamic tracking and underwent 4DCT scanning at the time of simulation. Data including the 3D coordinates of the target over the course of treatment were extracted from the treatment log files and used to determine actual target motion in the superior–inferior (S–I), anterior–posterior (A–P), and left–right (L–R) directions. Different approaches were used to calculate anisotropic and isotropic margins, assuming that the tumor moves as a rigid body. Anisotropic margins were calculated by separating target motion in the three anatomical directions, and a uniform margin was calculated by shifting the gross tumor volume contours in the 3D space and by computing the percentage of overlap with the PTV. The analysis was validated by means of a theoretical formulation. Results: The three methods provided consistent results. A uniform margin of 4.5 mm around the ITV was necessary to assure 95% target coverage for 95% of the fractions included in the analysis. In the case of anisotropic margins, the expansion required in the S–I direction was larger (8.1 mm) than those in the L–R (4.9 mm) and A–P (4.5 mm) directions. This margin accounts for variations of target position within the same treatment fraction. Conclusions: The use of bony alignment for CyberKnife lung stereotactic body radiation therapy requires careful considerations, in terms of the potential for increased toxicity or local miss. Our method could be used by other centers to determine the adequacy of ITV-to-PTV margins for their patients.

60 APPLIED LIFE SCIENCES↗

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗

Dynamics of the Bingham Canyon rock avalanches (Utah, USA) resolved from topographic, seismic, and infrasound data: Bingham Canyon Rock Avalanches

The 2013 Bingham Canyon Mine rock avalanches represent one of the largest cumulative landslide events in recorded U.S. history and provide a unique opportunity to test remote analysis techniques for landslide characterization. We combine aerial photogrammetry surveying, topographic reconstruction, numerical runout modeling, and analysis of broadband seismic and infrasound data to extract salient details of the dynamics and evolution of the multiphase landslide event. Our results reveal a cumulative intact rock source volume of 52 Mm 3 , which mobilized in two main rock avalanche phases separated by 1.5 h. We estimate that the first rock avalanche had 1.5–2 times greater volume than the second. Each failure initiated by sliding along a gently dipping (21°), highly persistent basal fault before transitioning to a rock avalanche and spilling into the inner pit. The trajectory and duration of the two rock avalanches were reconstructed using runout modeling and independent force history inversion of intermediate-period (10–50 s) seismic data. Intermediate- and shorter-period (1–50 s) seismic data were sensitive to intervals of mass redirection and constrained finer details of the individual slide dynamics. Back projecting short-period (0.2–1 s) seismic energy, we located the two rock avalanches within 2 and 4 km of the mine. Further analysis of infrasound and seismic data revealed that the cumulative event included an additional 11 smaller landslides (volumes ~10 4 –10 5 m 3 ) and that a trailing signal following the second rock avalanche may result from an air-coupled Rayleigh wave. These results demonstrate new and refined techniques for detailed remote characterization of the dynamics and evolution of large landslides.

58 GEOSCIENCES↗

Isovector and isoscalar tensor charges of the nucleon from lattice QCD

We present results for the iso-vector and flavor diagonal tensor charges $g$$^{u−d}_{T}$, $g$$^{u}_{T}$, $g$$^{d}_{T}$, and $g$$^{s}_{T}$ needed to probe novel tensor interactions at the TeV scale in neutron and nuclear β-decays and the contribution of the quark electric dipole moment (EDM) to the neutron EDM. The lattice QCD calculations were done using nine ensembles of gauge configurations generated by the MILC collaboration using the HISQ action with 2+1+1 dynamical flavors. These ensembles span three lattice spacings a ≈ 0.06, 0.09 and 0.12 fm and three quark masses corresponding to the pion masses M π ≈ 130, 220 and 310 MeV. Using estimates from these ensembles, we quantify all systematic uncertainties and perform a simultaneous extrapolation in the lattice spacing, volume and light quark masses for the connected contributions. The final estimates of the connected nucleon (proton) tensor charge for the iso-vector combination is $g$$^{u−d}_{T}$ = 1.020(76) in the MS scheme at 2 GeV. The additional disconnected quark loop contributions needed for the flavor-diagonal matrix elements are calculated using a stochastic estimator employing the truncated solver method with the all-mode-averaging technique. We find that the size of the disconnected contribution is smaller than the statistical error in the connected contribution. This allows us to bound the disconnected contribution and include it as an additional uncertainty in the flavor-diagonal charges. After a continuum extrapolation, we find $g$$^{u}_{T}$ = 0.774(66), $g$$^{d}_{T}$ = −0.233(28) and $g$$^{u+d}_{T}$ = 0.541(67). The strangeness tensor charge, that can make a significant contribution to the neutron EDM due to the large ratio ms/mu,d, is $g$$^{s}_{T}$ = 0.008(9) in the continuum limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Financial Analysis of the Smallmouth Bass Flows implemented at the Glen Canyon Dam during Water Year 2024

The Glen Canyon Dam (GCD) is a Colorado River Storage Project (CRSP) power resource that is a component of the Salt Lake City Area Integrated Projects (SLCA/IP). The 2016 record of decision (ROD) for the GCD long-term experimental and management plan (LTEMP) final Environmental Impact Statement (EIS) specifies criteria for GCD monthly water releases, daily and hourly operating limits, and experimental releases. This report presents a financial analysis of the Smallmouth bass (Micropterus dolomieu) (SMB) flows implemented at GCD during Water Year (WY) 2024. These bypass flows were introduced by the U.S. Bureau of Reclamation (USBR) as an emergency response to the growing threat posed by invasive SMB in the Colorado River ecosystem downstream of the dam. SMB are a non-native predatory species that pose a significant threat to native fish populations, including the endangered humpback chub (Gila cypha). The thermal regime below GCD, typically cold due to hypolimnetic releases from Lake Powell, has historically served as a thermal barrier limiting SMB establishment. However, persistently low reservoir levels in recent years have reduced stratification in Lake Powell, allowing warmer water to be released downstream. This has enabled SMB to spawn successfully below the dam, prompting urgent ecological concerns. To mitigate the risk of SMB proliferation, the USBR implemented a series of bypass flows in WY 2024. Drawn from a lower elevation than the penstocks, the bypass structures released cooler water downstream. These short-duration bypass flows aimed to keep temperatures cool enough to prevent SMB from spawning, thereby reducing the ecological threat posed by this invasive species. Although motivated by ecological objectives, these bypass flows came with financial tradeoffs. Releasing water through the bypass structures instead of the turbines at GCD reduced hydropower generation, resulting in a significantly lower financial position for Western Area Power Administration (WAPA), which is responsible for marketing the electricity produced by the GCD Powerplant. This report analyzes the financial impact of the SMB flows implemented from July to November 2024. These experimental releases led to an estimated financial cost of approximately $18.9 million, primarily driven by the substantial volume of water diverted through the bypass structures. This study applies an integrated set of tools to estimate WAPA financial impacts by simulating GCD under two types of cases; namely, (1) a “With Experiment” case that mimics the water operations that actually occurred, including the SMB bypass flows, and (2) a “Without Experiment” case that simulates operations under the assumption that the SMB flows did not occur. Both cases comply with LTEMP hourly and daily operating criteria, and the monthly water release volumes are assumed to be identical under both cases. The Colorado River Storage Project Python-based model (CRiSPPy) model was the main modeling tool used to simulate the dispatch of the GCD hydropower plant and associated water releases from Lake Powell. In the modeling process, the research team used extensive data sets and historical information on SLCA/IP power plant characteristics, hydrologic conditions, and WAPA’s power purchases and sales prices.

13 HYDRO ENERGY↗

Changes in cardiac output and tibial artery flow during and after progressive LBNP

A 3.0 MHz Pulsed Doppler velocity meter (PD) was used to determine blood velocities in the ascending aorta from the suprasternal notch before, during and after progressive 5 min stages of lower body negative pressure (LBNP) in 7 subjects. Changes in stroke volume were calculated from the systolic velocity integrals. A unique 20 MHz PD was used to estimate bloodflow in the posterior tibial artery. With -20 torr mean stroke volume fell 11% and then continued to decline by 48% before LBNP was terminated. Mean tibial flow fell progressively with LBNP stress, due to an increase in reverse flow component and a reduction in peak forward flow and diameter. Stroke volume increased and heart rate fell dramatically during the first 15 sec of recovery. The LBNP was terminated early in 2 subjects because of vasovagal symptons (V). During V the stroke volume rose 86% which more than compensated for the drop in heart rate. This implies that V is accompanied by a paradoxical increase in venous return and that the reduction in HR is the primary cardiovascular event. During the first 15 sec of recovery these 2 subjects had a distinctive marked rise to heart rate reminiscent of the Bainbridge reflex.

Source record↗

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SU-F-T-450: The Investigation of Radiotherapy Quality Assurance and Automatic Treatment Planning Based On the Kernel Density Estimation Method

Purpose: To develop a fast automatic algorithm based on the two dimensional kernel density estimation (2D KDE) to predict the dose-volume histogram (DVH) which can be employed for the investigation of radiotherapy quality assurance and automatic treatment planning. Methods: We propose a machine learning method that uses previous treatment plans to predict the DVH. The key to the approach is the framing of DVH in a probabilistic setting. The training consists of estimating, from the patients in the training set, the joint probability distribution of the dose and the predictive features. The joint distribution provides an estimation of the conditional probability of the dose given the values of the predictive features. For the new patient, the prediction consists of estimating the distribution of the predictive features and marginalizing the conditional probability from the training over this. Integrating the resulting probability distribution for the dose yields an estimation of the DVH. The 2D KDE is implemented to predict the joint probability distribution of the training set and the distribution of the predictive features for the new patient. Two variables, including the signed minimal distance from each OAR (organs at risk) voxel to the target boundary and its opening angle with respect to the origin of voxel coordinate, are considered as the predictive features to represent the OAR-target spatial relationship. The feasibility of our method has been demonstrated with the rectum, breast and head-and-neck cancer cases by comparing the predicted DVHs with the planned ones. Results: The consistent result has been found between these two DVHs for each cancer and the average of relative point-wise differences is about 5% within the clinical acceptable extent. Conclusion: According to the result of this study, our method can be used to predict the clinical acceptable DVH and has ability to evaluate the quality and consistency of the treatment planning.

60 APPLIED LIFE SCIENCES↗

Estimation of Rain Intensity Spectra over the Continental US Using Ground Radar-Gauge Measurements

A high-resolution surface rainfall product is used to estimate rain characteristics over the continental US as a function of rain intensity. By defining each data at 4-km horizontal resolutions and 1-h temporal resolutions as an individual precipitating/nonprecipitating sample, statistics of rain occurrence and rain volume including their geographical and seasonal variations are documented. Quantitative estimations are also conducted to evaluate the impact of missing light rain events due to satellite sensors' detection capabilities. It is found that statistics of rain characteristics have large seasonal and geographical variations across the continental US. Although heavy rain events (> 10 mm/hr.) only occupy 2.6% of total rain occurrence, they may contribute to 27% of total rain volume. Light rain events (< 1.0 mm/hr.), occurring much more frequently (65%) than heavy rain events, can also make important contributions (15%) to the total rain volume. For minimum detectable rain rates setting at 0.5 and 0.2 mm/hr which are close to sensitivities of the current and future space-borne precipitation radars, there are about 43% and 11% of total rain occurrence below these thresholds, and they respectively represent 7% and 0.8% of total rain volume. For passive microwave sensors with their rain pixel sizes ranging from 14 to 16 km and the minimum detectable rain rates around 1 mm/hr., the missed light rain events may account for 70% of train occurrence and 16% of rain volume. Statistics of rain characteristics are also examined on domains with different temporal and spatial resolutions. Current issues in estimates of rain characteristics from satellite measurements and model outputs are discussed.

Lin, Xin↗

Design of the protoDUNE raw data management infrastructure

The Deep Underground Neutrino Experiment (DUNE) will employ a set of Liquid Argon Time Projection Chambers (LArTPC) with a total mass of 40 kt as the main components of its Far Detector. In order to validate this technology and characterize the detector performance at full scale, an ambitious experimental program (called “protoDUNE”) has been initiated which includes a test of the large-scale prototypes for the single-phase and dual-phase LArTPC technologies, which will run in a beam at CERN. The total raw data volume that is slated to be collected during the scheduled 3-month beam run is estimated to be in excess of 2.5 PB for each detector. This data volume will require that the protoDUNE experiment carefully design the DAQ, data handling and data quality monitoring systems to be capable of dealing with challenges inherent with peta-scale data management while simultaneously fulfilling the requirements of disseminating the data to a worldwide collaboration and DUNE associated computing sites. Here in this paper, we present our approach to solving these problems by leveraging the design, expertise and components created for the LHC and Intensity Frontier experiments into a unified architecture that is capable of meeting the needs of protoDUNE.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗