Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level set method”

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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Implications of z ~ 6 Quasar Proximity Zones for the Epoch of Reionization and Quasar Lifetimes

In this paper, we study quasar proximity zones in the redshift range $5.77\leqslant z\leqslant 6.54$ by homogeneously analyzing 34 medium-resolution spectra, encompassing both archival and newly obtained data, and exploiting recently updated systemic redshift and magnitude measurements. Whereas previous studies found strong evolution of proximity zone sizes with redshift and argued that this provides evidence for a rapidly evolving intergalactic medium (IGM) neutral fraction during reionization, we measure a much shallower trend $\propto {(1+z)}^{-1.44}$. We compare our measured proximity zone sizes to predictions from hydrodynamical simulations post-processed with one-dimensional radiative transfer and find good agreement between observations and theory, irrespective of the ionization state of the ambient IGM. This insensitivity to IGM ionization state has been previously noted, and results from the fact that the definition of proximity zone size as the first drop of the smoothed quasar spectrum below the 10% flux transmission level probes locations where the ionizing radiation from the quasar is an order of magnitude larger than the expected ultraviolet ionizing background that sets the neutral fraction of the IGM. Our analysis also uncovered three objects with exceptionally small proximity zones (two have ${R}_{p}\lt 1$ proper Mpc), which constitute outliers from the observed distribution and are challenging to explain with our radiative transfer simulations. Finally, we consider various explanations for their origin, such as strong absorption line systems associated with the quasar or patchy reionization, but find that the most compelling scenario is that these quasars have been shining for ≲10 5 years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Electrical properties and electronic structure of Si-implanted hexagonal boron nitride films

Si ion implantation with a set of ion energies and ion doses was carried out to dope hexagonal boron nitride (hBN) thin films synthesized by radio-frequency magnetron sputtering. Hall effect measurements revealed n-type conduction with a low resistivity of 0.5 Ω cm at room temperature, corresponding to an electron concentration of 2.0 × 10{sup 19} cm{sup −3} and a mobility of 0.6 cm{sup 2}/V s. Temperature-dependent resistivity measurements in a wide temperature range from 50 to 800 K demonstrated two shallow donor levels in the hBN band gap induced by Si doping, which was in consistence with the theoretical calculation by density function theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Polarization contributions to intermolecular interactions revisited with fragment electric-field response functions

The polarization energy in intermolecular interactions treated by self-consistent field electronic structure theory is often evaluated using a constraint that the atomic orbital (AO) to molecular orbital transformation is blocked by fragments. This approach is tied to AO basis sets, overestimates polarization energies in the overlapping regime, particularly in large AO basis sets, and lacks a useful complete basis set limit. These problems are addressed by the construction of polarization subspaces based on the responses of isolated fragments to weak electric fields. These subspaces are spanned by fragment electric-field response functions, which can capture effects up to the dipole (D), or quadrupole (DQ) level, or beyond. Schemes are presented for the creation of both non-orthogonal and orthogonal fragment subspaces, and the basis set convergence of the polarization energies computed using these spaces is assessed. Numerical calculations for the water dimer, water–Na{sup +}, water–Mg{sup 2+}, water–F{sup −}, and water–Cl{sup −} show that the non-orthogonal DQ model is very satisfactory, with small differences relative to the orthogonalized model. Additionally, we prove a fundamental difference between the polarization degrees of freedom in the fragment-blocked approaches and in constrained density schemes. Only the former are capable of properly prohibiting charge delocalization during polarization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Adaptive experimental design for multi-fidelity surrogate modeling of multi-disciplinary systems

We present an adaptive algorithm for constructing surrogate models of multi-disciplinary systems composed of a set of coupled components. With this goal we introduce “coupling” variables with a priori unknown distributions that allow surrogates of each component to be built independently. Once built, the surrogates of the components are combined to form an integrated-surrogate that can be used to predict system-level quantities of interest at a fraction of the cost of the original model. The error in the integrated-surrogate is greedily minimized using an experimental design procedure that allocates the amount of training data, used to construct each component-surrogate, based on the contribution of those surrogates to the error of the integrated-surrogate. Here, the multi-fidelity procedure presented is a generalization of multi-index stochastic collocation that can leverage ensembles of models of varying cost and accuracy, for one or more components, to reduce the computational cost of constructing the integrated-surrogate. Extensive numerical results demonstrate that, for a fixed computational budget, our algorithm is able to produce surrogates that are orders of magnitude more accurate than methods that treat the integrated system as a black-box.

42 ENGINEERING↗
DOE OSTI2023

Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT

Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.

58 GEOSCIENCES↗
DOE OSTI2016

Radiation Dose Escalation in Esophageal Cancer Revisited: A Contemporary Analysis of the National Cancer Data Base, 2004 to 2012

Purpose: To evaluate the effect of radiation dose escalation on overall survival (OS) for patients with nonmetastatic esophageal cancer treated with concurrent radiation and chemotherapy. Methods and Materials: Patients diagnosed with stage I to III esophageal cancer treated from 2004 to 2012 were identified from the National Cancer Data Base. Patients who received concurrent radiation and chemotherapy with radiation doses of ≥50 Gy and did not undergo surgery were included. OS was compared using Cox proportional hazards regression and propensity score matching. Results: A total of 6854 patients were included; 3821 (55.7%) received 50 to 50.4 Gy and 3033 (44.3%) received doses >50.4 Gy. Univariate analysis revealed no significant difference in OS between patients receiving 50 to 50.4 Gy and those receiving >50.4 Gy (P=.53). The dose analysis, binned as 50 to 50.4, 51 to 54, 55 to 60, and >60 Gy, revealed no appreciable difference in OS within any group compared with 50 to 50.4 Gy. Subgroup analyses investigating the effect of dose escalation by histologic type and in the setting of intensity modulated radiation therapy also failed to reveal a benefit. Propensity score matching confirmed the absence of a statistically significant difference in OS among the dose levels. The factors associated with improved OS on multivariable analysis included female sex, lower Charlson-Deyo comorbidity score, private insurance, cervical/upper esophagus location, squamous cell histologic type, lower T stage, and node-negative status (P<.01 for all analyses). Conclusions: In this large national cohort, dose escalation >50.4 Gy did not result in improved OS among patients with stage I to III esophageal cancer treated with definitive concurrent radiation and chemotherapy. These data suggest that despite advanced contemporary treatment techniques, OS for patients with esophageal cancer remains unaltered by escalation of radiation dose >50.4 Gy, consistent with the results of the INT-0123 trial. Furthermore, these data highlight that many radiation oncologists have not embraced the concept that dose escalation does not improve OS. Although local control, not investigated in the present study, might benefit from dose escalation, novel therapies are needed to improve the OS of patients with esophageal cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Identifying Sociomarkers of Pediatric Asthma Patients at Risk of Hospital Revisiting

ObjectiveAsthma is one of the most common chronic childhood diseases in the United States [2, 3]. In addition to its pervasiveness, pediatric asthma shows high sensitivity to the environment. Combining medical-social dataset with machine learning methods we demonstrate how socio-markers play an important role in identifying patients at risk of hospital revisits due to pediatric asthma within a year.IntroductionA socio-marker is a measurable indicator of social conditions where a patient is embedded in and exposed to, being analogous with a biomarker indicating the severity or presence of some disease state. Social factors are one of the most clinical health determinants [1], which play a critical role in explaining health outcomes. Socio-markers can help medical practitioners and researchers to reliably identify high-risk individuals in a timely manner.MethodsWe collected data from three different sources: pediatric asthma encounter records from Jan 1st, 2016 to Dec 31st, 2016 at a children’s hospital, the 2010 U.S census data and neighborhood quality survey data by Memphis Property Hub. After merging these datasets we examine the effect of social features in identifying the patients who visited the hospital more than once during the observation period. We only use the first time visit (3,678 cases) to avoid over-counting of the same patients. In addition to demographic features (age, gender, insurance type, and race (African American and White)), we incorporate the social features such as the proportion of individuals living below the federal poverty level, blight prevalence, neighborhood quality, neighborhood quality inequality, trash dumping presence, the broken window pervasiveness within the zip code area of patients’ residence are included.We then implemented a Support Vector Machine (SVM) based classification model using abovementioned 11 social features. The classification outcome is either patient visits the hospital only one-time (class 0) or revisits the hospital within a year (class 1). Among 3,678 unique patients in the dataset, only 823 pediatric patients revisited hospital with asthma. So, to overcome the class imbalance issue, we have used 823 patients’ data (randomly selected in 1,000 iterations) from each class. Further, to avoid overfitting and ensure generalizability, we divided the dataset as training, test, and validation with a proportion of 60%, 20%, and 20%, respectively. The reported test (5-folds cross-validation using training and testing data) and validation accuracy of the SVM method are averaged over 1,000 iterations to avoid sampling error and bias.ResultsThe proposed socio-marker model resulted in an average classification accuracy of 63.70% for the test set and 63.67 % for the validation set. Further, the average specificity (the total true negative cases divided by the sum of true negative and false positive) and sensitivity (the total number of true positive cases divided by the sum of positive predicted cases) is found to be 62.79% and 64.77%, respectively for the test set and 62.79% and 64.83%, respectively for the validation set. Results of this study suggest that socio-marker features that are not directly related to a patient’s medical conditions can still predict whether the patient will come back to the hospital within a year or not with approximately 64% accuracy.ConclusionsBringing the socio-marker features in the surveillance system may ease the burden of detecting the patients at risk of revisiting the hospital. The results should be interpreted with caution because we only used 12-month period of observation and the visit beyond the observation window is not considered. Also the patients may have visited different hospitals which are not captured in the data.References1. Booske BC, Athens JK, Kindig DA, Park H, Remington PL: Different perspectives for assigning weights to determinants of health. University of Wisconsin: Population Health Institute 2010.2. Subbarao P, Mandhane PJ, Sears MR: Asthma: epidemiology, etiology and risk factors. Canadian Medical Association Journal 2009, 181(9):E181-E190.3. Gold DR, Wright R: Population disparities in asthma. Annu Rev Public Health 2005, 26:89-113.

Shin, Eun Kyong↗
DOE OSTI2023

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes

Here, the thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller–Peterson–Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Interval-Valued Rank in Finite Ordered Sets

We consider the concept of rank as a measure of the vertical levels and positions of elements of partially ordered sets (posets). We are motivated by the need for algorithmic measures on large, real-world hierarchically-structured data objects like the semantic hierarchies of ontolog- ical databases. These rarely satisfy the strong property of gradedness, which is required for traditional rank functions to exist. Representing such semantic hierarchies as finite, bounded posets, we recognize the duality of ordered structures to motivate rank functions which respect verticality both from the bottom and from the top. Our rank functions are thus interval-valued, and always exist, even for non-graded posets, providing order homomorphisms to an interval order on the interval-valued ranks. The concept of rank width arises naturally, allowing us to identify the poset region with point-valued width as its longest graded portion (which we call the “spindle”). A standard interval rank function is naturally motivated both in terms of its extremality and on pragmatic grounds. Its properties are examined, including the relation- ship to traditional grading and rank functions, and methods to assess comparisons of standard interval-valued ranks.

Joslyn, Cliff↗
DOE OSTI2014

Probing ground and low-lying excited states for HIO{sub 2} isomers

We present a computational study on HIO{sub 2} molecules. Ground state properties such as equilibrium structures, relative energetics, vibrational frequencies, and infrared intensities were obtained for all the isomers at the coupled-cluster with single and double excitations as well as perturbative inclusion of triples (CCSD(T)) level of theory with the aug-cc-pVTZ-PP basis set and ECP-28-PP effective core potential for iodine and the aug-cc-pVTZ basis set for hydrogen and oxygen atoms. The HOIO structure is confirmed as the lowest energy isomer. The relative energies are shown to be HOIO < HOOI < HI(O)O. The HO(O)I isomer is only stable at the density functional theory (DFT) level of theory. The transition states determined show interconversion of the isomers is possible. In order to facilitate future experimental identification, vibrational frequencies are also determined for all corresponding deuterated species. Vertical excitation energies for the three lowest-lying singlet and triplet excited states were determined using the configuration interaction singles, time-dependent density functional theory (TD-DFT)/B3LYP, TD-DFT/G96PW91, and equation of motion-CCSD approaches with the LANL2DZ basis set plus effective core potential for iodine and the aug-cc-pVTZ basis set for hydrogen and oxygen atoms. It is shown that HOIO and HOOI isomers have excited states accessible at solar wavelengths (<4.0 eV) but these states have very small oscillator strengths (<2 × 10{sup −3})

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

SU-E-T-433: Calibration Accuracy in Mailed High-Resolution 3D Dosimetry Service for SRS/SBRT QA

Purpose: SRS/SBRT combines hypofractionation with excellent dose distributions. However, extremely steep gradients across the target along with dose escalation, if not administered accurately, may lead to serious complications, recurrences, or even fatalities. Existing commercial QA products either lack adequate spatial resolution or the 3D aspect. By contrast, the new CrystalBall™ mailed high-resolution 3D dosimetry service removes the above limitations while reducing the overall workload on medical physics staff. The exposed dosimeters, which change optical density in proportion to local dose, are sent back to the manufacturer (MGS Research Inc., Madison, CT) for sub-millimeter-resolution laser-CT scanning and QA data analysis. QA report is returned electronically within 24 hours. The purpose of this study was to evaluate the dose calibration accuracy in this system. Methods: Two spherical CrystalBall™ polymer gel dosimeters from the same batch, 166 mm diameter, with embedded 3D image registration markers, were mounted in a special phantom designed for reproducible positioning. For full end to end testing, the optical guidance array was mounted onto the phantom and a CT was taken. Two separate Rapid Arc SRS plans were designed. Varian Medical Systems optical guidance system was used to position the phantom and the SRS treatment plans were delivered to the two spheres on Varian's Trilogy Accelerator. Exposed dosimeters were mailed back to the manufacturer for laser CT scanning and analysis. Results: For each plan, 3D gamma passing rate was 100% for 2%/2mm distance-to-agreement criteria above 50% isodose level. The two calibration curves, generated using volumetric dose and optical density data, showed excellent mutual agreement (max difference 2.2%, median difference 0.75%). Conclusion: The clinical utility of new CrystalBall™ mailed QA service for SRS/SBRT and high accuracy of dose calibration have been validated. The workflow associated with the use of the CrystalBall™ in clinical setting was found to be minimal. The presenting author is the founder of and has an ownership in MGS Research Inc which manufactures the CrystalBall system for 3D dosimetry.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Mortandad Canyon seismic modeling (Progress Report)

The goal of this project is to assess site amplification for LANL through modeling of site and topography effects. Such assessment is difficult due to the lack of data documenting site effects for large normal faulting earthquakes such as the Pajarito Fault System that was found to be governing the seismic hazard at LANL. LANL site is also characterized by unique geologic settings, a series of mesa at the edge of the huge volcanic complex of the Jemez caldera. The latter has produced a series of successive and intertwined deposits of volcanic flow and sedimentary layer. Such layering of material with different elastic properties is expected to create large amplification effects magnified further by the presence at the surface of tuff materials. In 2014, a seismic campaign was performed in the Mortandad Canyon with a network of sensors on and near a local mesa. That campaign was designed to provide experimental data to estimation and prediction of amplification and site effects at LANL. An analysis on the seismic noise records of that campaign shows evidence of amplification of the seismic amplitude for the stations on the top of the mesa. Our goal is to reproduce the effect observed (both retrieving the resonance frequencies and amplification levels). We also plan to test if amplification effects observed for seismic noise can be extended to shaking produced by earthquakes. We are using a 3D modeling tool based on the Spectral Element Method (SEM) that is a particular type of high-order finite element modeling. Two types of seismic sources will be considered: (1) earthquakes for purpose of seismic hazard evaluation, (2) ambient noise. We show 2 series of numerical computations: (1) seismic waves are generated by random sources to simulate ambient noise. (2) seismic waves are planar waves to simulate a distant earthquake. We plan in the future to model rupture scenario on the Pajarito Fault System.

58 GEOSCIENCES↗
DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

PHISICS/RELAP5-3D RESULTS FOR EXERCISES II-1 AND II-2 OF THE OECD/NEA MHTGR-350 BENCHMARK

The Idaho National Laboratory (INL) Advanced Reactor Technologies (ART) High-Temperature Gas-Cooled Reactor (HTGR) Methods group currently leads the Modular High-Temperature Gas-Cooled Reactor (MHTGR) 350 benchmark. The benchmark consists of a set of lattice-depletion, steady-state, and transient problems that can be used by HTGR simulation groups to assess the performance of their code suites. The paper summarizes the results obtained for the first two transient exercises defined for Phase II of the benchmark. The Parallel and Highly Innovative Simulation for INL Code System (PHISICS), coupled with the INL system code RELAP5-3D, was used to generate the results for the Depressurized Conduction Cooldown (DCC) (exercise II-1a) and Pressurized Conduction Cooldown (PCC) (exercise II-2) transients. These exercises require the time-dependent simulation of coupled neutronics and thermal-hydraulics phenomena, and utilize the steady-state solution previously obtained for exercise I-3 of Phase I. This paper also includes a comparison of the benchmark results obtained with a traditional system code “ring” model against a more detailed “block” model that include kinetics feedback on an individual block level and thermal feedbacks on a triangular sub-mesh. The higher spatial fidelity that can be obtained by the block model is illustrated with comparisons of the maximum fuel temperatures, especially in the case of natural convection conditions that dominate the DCC and PCC events. Differences up to 125 K (or 10%) were observed between the ring and block model predictions of the DCC transient, mostly due to the block model’s capability of tracking individual block decay powers and more detailed helium flow distributions. In general, the block model only required DCC and PCC calculation times twice as long as the ring models, and it therefore seems that the additional development and calculation time required for the block model could be worth the gain that can be obtained in the spatial resolution

Advanced Reactor Technologies↗
DOE OSTI2014

Improved detection of rough defects for ultrasonic NDE inspections based on finite element modeling of elastic wave scattering

Defects which posses rough surfaces greatly affect ultrasonic wave scattering behaviour, often reducing the magnitude of reflected signals. Ultrasonic inspections rely upon this response for detecting and sizing flaws. For safety critical components reliable characterisation is crucial. Therefore, providing an accurate means to predict reductions in signal amplitude is essential. An extension of Kirchhoff theory has formed the basis for the UK power industry inspection justifications. However, it is widely recognised that these predictions are pessimistic owing to analytical approximations. A numerical full field modelling approach does not fall victim to such limitations. Here, a Finite Element model is used to aid in setting a non-conservative reporting threshold during the inspection of a large pressure vessel forging that might contain embedded rough defects. The ultrasonic response from multiple rough surfaces defined by the same statistical class is calculated for normal incident compression waves. The approach is validated by comparing coherent scattering with predictions made by Kirchhoff theory. At lower levels of roughness excellent agreement is observed, whilst higher values confirm the pessimism of Kirchhoff theory. Furthermore, the mean amplitude in the specular direction is calculated. This represents the information obtained during an inspection, indicating that reductions due to increasing roughness are significantly less than the coherent component currently being used.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

SU-E-T-492: Implementing a Method for Brain Irradiation in Rats Utilizing a Commercially Available Radiosurgery Irradiator

Purpose: The purpose of the study was to implement a method for accurate rat brain irradiation using the Gamma Knife Perfexion unit. The system needed to be repeatable, efficient, and dosimetrically and spatially accurate. Methods: A platform (“rat holder”) was made such that it is attachable to the Leskell Gamma Knife G Frame. The rat holder utilizes two ear bars contacting bony anatomy and a front tooth bar to secure the rat. The rat holder fits inside of the Leskell localizer box, which utilizes fiducial markers to register with the GammaPlan planning system. This method allows for accurate, repeatable setup.A cylindrical phantom was made so that film can be placed axially in the phantom. We then acquired CT image sets of the rat holder and localizer box with both a rat and the phantom. Three treatment plans were created: a plan on the rat CT dataset, a phantom plan with the same prescription dose as the rat plan, and a phantom plan with the same delivery time as the rat plan. Results: Film analysis from the phantom showed that our setup is spatially accurate and repeatable. It is also dosimetrically accurate, with an difference between predicted and measured dose of 2.9%. Film analysis with prescription dose equal between rat and phantom plans showed a difference of 3.8%, showing that our phantom is a good representation of the rat for dosimetry purposes, allowing for +/- 3mm diameter variation. Film analysis with treatment time equal showed an error of 2.6%, which means we can deliver a prescription dose within 3% accuracy. Conclusion: Our method for irradiation of rat brain has been shown to be repeatable, efficient, and accurate, both dosimetrically and spatially. We can treat a large number of rats efficiently while delivering prescription doses within 3% at millimeter level accuracy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Pollutant Exposures from Natural Gas Cooking Burners: A Simulation-Based Assessment for Southern California

Background: Residential natural gas cooking burners (NGCBs) can emit substantial quantities of pollutants, and they are typically used without venting range hoods. Objective: We quantified pollutant concentrations and occupant exposures resulting from NGCB use in California homes. Methods: A mass-balance model was applied to estimate time-dependent pollutant concentrations throughout homes in Southern California and the exposure concentrations experienced by individual occupants. We estimated nitrogen dioxide (NO2), carbon monoxide (CO), and formaldehyde (HCHO) concentrations for 1 week each in summer and winter for a representative sample of Southern California homes. The model simulated pollutant emissions from NGCBs as well as NO2 and CO entry from outdoors, dilution throughout the home, and removal by ventilation and deposition. Residence characteristics and outdoor concentrations of NO2 and CO were obtained from available databases. We inferred ventilation rates, occupancy patterns, and burner use from household characteristics. We also explored proximity to the burner(s) and the benefits of using venting range hoods. Replicate model executions using independently generated sets of stochastic variable values yielded estimated pollutant concentration distributions with geometric means varying by < 10%. Results: The simulation model estimated that—in homes using NGCBs without coincident use of venting range hoods—62%, 9%, and 53% of occupants are routinely exposed to NO2, CO, and HCHO levels that exceed acute health-based standards and guidelines. NGCB use increased the sample median of the highest simulated 1-hr indoor concentrations by 100, 3,000, and 20 ppb for NO2, CO, and HCHO, respectively. Conclusions: Reducing pollutant exposures from NGCBs should be a public health priority. Simulation results suggest that regular use of even moderately effective venting range hoods would dramatically reduce the percentage of homes in which concentrations exceed health-based standards.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Exploring ice sheet model sensitivity to ocean thermal forcing and basal sliding using the Community Ice Sheet Model (CISM)

Abstract. Multi-meter sea level rise (SLR) is thought to be possible within the next few centuries, with most of the uncertainty originating from the Antarctic land ice contribution. One source of uncertainty relates to the ice sheet model initialization. Since ice sheets have a long response time (compared to other Earth system components such as the atmosphere), ice sheet model initialization methods can have significant impacts on how the ice sheet responds to future forcings. To assess this, we generated 25 different ice sheet spin-ups, using the Community Ice Sheet Model (CISM) at a 4 km resolution. During each spin-up, we varied two key parameters known to impact the sensitivity of the ice sheet to future forcing: one related to the sensitivity of the ice shelf melt rate to ocean thermal forcing (TF) and the other related to the basal friction. The spin-ups all nudge toward observed thickness and enforce a no-advance calving criterion, such that all final spin-up states resemble observations but differ in their melt and friction parameter settings. Each spin-up was then forced with future ocean thermal forcings from 13 different CMIP6 models under the Shared Socioeconomic Pathway (SSP)5-8.5 emissions scenario and modern climatological surface mass balance data. Our results show that the effects of the ice sheet and ocean parameter settings used during the spin-up are capable of impacting multi-century future SLR predictions by as much as 2 m. By the end of this century, the effects of these choices are more modest, but still significant, with differences of up to 0.2 m of SLR. We have identified a combined ocean and ice parameter space that leads to widespread mass loss within 500 years (low friction and high melt rate sensitivity). To explore temperature thresholds, we also ran a synthetically forced CISM ensemble that is focused on the Amundsen region only. Given certain ocean and ice parameter choices, Amundsen mass loss can be triggered with thermal forcing anomalies between 1.5 and 2 ∘C relative to the spin-up. Our results emphasize the critical importance of considering ice sheet and ocean parameter choices during spin-up for SLR predictions and suggest the importance of including glacial isostatic adjustment in ice sheet simulations.

54 ENVIRONMENTAL SCIENCES↗