Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EVAPORATION MODEL”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Theoretical analysis of mixing in liquid clouds – Part IV: DSD evolution and mixing diagrams

Evolution of droplet size distribution (DSD) due to mixing between cloudy and dry volumes is investigated for different values of the cloud fraction and for different initial DSD shapes. The analysis is performed using a diffusion–evaporation model which describes time-dependent processes of turbulent diffusion and droplet evaporation within a mixing volume. Time evolution of the DSD characteristics such as droplet concentration, LWC and mean volume radii is analyzed. The mixing diagrams are plotted for the final mixing stages. It is shown that the difference between the mixing diagrams for homogeneous and inhomogeneous mixing is insignificant and decreases with an increase in the DSD width. The dependencies of the normalized cube of the mean volume radius on the cloud fraction were compared with those on normalized droplet concentration and found to be quite different. If the normalized droplet concentration is used, mixing diagrams do not show any significant dependence on relative humidity in the dry volume. The main conclusion of the study is that traditional mixing diagrams cannot serve as a reliable tool for analysis of mixing type.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Trends and meteorological drivers of extreme daily reservoir evaporation events in the western United States

Extreme daily evaporation from reservoir surfaces can lead to significant short-term water losses, affecting water quality, water supply, and reservoir operation strategies. Historical trends in daily reservoir evaporation events have eluded the scientific and operational communities, largely due to a lack of long-term, consistent data record. This study quantifies trends in extreme daily reservoir evaporation events at 165 major reservoirs located in the western U.S. Here, we use the place-based energy balance and aerodynamic Daily Lake Evaporation Model (DLEM) driven by multiple meteorological data products (RTMA, gridMET, Daymet) to estimate daily evaporation rates at these reservoirs from 1981 to 2022. The results—while are based on different meteorological forcing datasets—consistently indicate that the California, Lower Colorado, and Rio Grande hydrologic regions are more prone to higher daily evaporation extremes. Compared to the rest of western U.S, these three regions also experience a more pronounced increasing trend in the annual maximum daily evaporation rate, at about 0.3 mm day -1 decade -1 during 1981-2022. The results show that heat and dryness are the main drivers to the increasing trend of extreme evaporation, while extreme wind speed is the primary contributor to exceptionally high daily evaporation events across all regions. This phenomenon is particularly prominent in the arid Lower Colorado region, but less significant in the cold and humid Pacific Northwest region. We also find that the correlation between extreme wind speed and extreme evaporation degrades as the time scale increases from daily, to monthly and seasonal. Our findings have strong implications for the pattern and distribution of extreme evaporation events at the western U.S. reservoirs, and illustrate how various drivers influence extreme evaporation across different time scales.

13 HYDRO ENERGY↗
DOE OSTI2022

Life Expectancy of Evaporating Capillary Bridges Predicted by Tertiary Creep Modeling

The evaporation of capillary bridges is experimentally investigated at the microscale through a three-grain capillary cluster. This setting provides the minimum viable description of Haines jumps during evaporation, that is, capillary instabilities stemming from air entry into a saturated granular material. The displacement profile of a meniscus is obtained via digital image correlation for different grain materials, geometries, and separations. While it is well known that Haines jumps are triggered at the pore throat, we find that these instabilities are of three types depending on the separation. We also provide a temporal characterization of Haines jumps; we find that they are accurately described, as tertiary creep instabilities, by Voight’s relation, similarly to landslides and volcanic eruptions. This finding extends the description of capillary instabilities beyond their onset predicted by Laplace equilibrium. Our contribution also paves the way for a microscopically-informed description of desiccation cracks, of which Haines jumps are the precursors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Data and scripts associated with a manuscript analyzing ELM-FATES parameter sensitivity under pre-fire and postfire scenarios using machine learning

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).

Aboveground biomass↗
DOE OSTI2025

Measurement of neutron production in atmospheric neutrino interactions at Super-Kamiokande

We present measurements of total neutron production from atmospheric neutrino interactions in water, analyzed as a function of electron-equivalent visible energy over a range of 30 MeV to 10 GeV. These results are based on 4,270 days of data collected by Super-Kamiokande, including 564 days with 0.011 wt% gadolinium added to enhance neutron detection. Neutron signal selection is based on a neural network trained on simulation, with its performance validated using an Am/Be neutron point source. The measurements are compared to predictions from neutrino event generators combined with various hadron-nucleus interaction models, which include an intranuclear cascade model and a nuclear deexcitation model. We observe significant variations in the predictions depending on the choice of hadron-nucleus interaction model. We discuss key factors that contribute to describing our data, such as in-medium effects in the intranuclear cascade and the accuracy of statistical evaporation modeling.

Cherenkov detectors↗
DOE OSTI2019

Evaluating Integrated Surface/Subsurface Permafrost Thermal Hydrology Models in ATS (v0.88) Against Observations from a Polygonal Tundra Site: Modeling Archive

This Modeling Archive supports the NGEE Arctic publication "Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site" available in the Geoscientific Model Development journal. We evaluated the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations in polygonal tundra at the Barrow Environmental Observatory. We simulated thermal hydrology of three-dimensional ice-wedge polygons with generic but broadly representative surface microtopography. We drove the simulations with meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of the soil properties and evaporation model parameters, the simulations were found to be consistent with observed water tables, active layer thickness, snowpack depth, evaporation, and observed soil temperatures at several depths in trough, rim, and center. The study builds confidence in the emerging class of integrated surface/subsurface permafrost simulators and provides an optimized set of model parameters for use in projections of permafrost dynamics in a warming climate. This dataset includes python, html, pdf, EXO, H5, xlsx, and csv files. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Volatile propellant droplet evaporation measurements in metered dose inhaler sprays

Many aerosol products rely on the rapid vaporization of volatile propellants to produce a fine spray. In the simplest case, these are binary mixtures of propellant and a delivered product which undergo a flash-evaporation process leaving only the less volatile product in the resultant droplet. In more complex applications, such as pressurized metered-dose inhalers, the non-propellant component may contain dissolved or suspended drug which precipitates or dries to form a matured particle. The size and morphology of the particles depend strongly on the time-history of the droplet as the propellant evaporates. However, measuring the dynamic evaporation processes that occur in dense sprays containing millions of droplets is challenging. In this paper, we demonstrate a novel application of Ultra Small Angle X-ray Scattering to measure the bulk composition of volatile HFC134a–ethanol sprays and compare the obtained results with simple evaporation models in a dry nitrogen environment. The data reveal that diffusion-limiting processes inside the droplet are equally important as external convection and mixing-limited factors in determining evaporative timescales.

42 ENGINEERING↗
DOE OSTI2021

Inclusive J /$\psi$ production at midrapidity in pp collisions at $\sqrt{s}$ = 13 TeV

We report on the inclusive J/ψ production cross section measured at the CERN Large Hadron Collider in proton–proton collisions at a center-of-mass energy $\sqrt{s}$ = 13 TeV. The J/ψ mesons are reconstructed in the e + e - decay channel and the measurements are performed at midrapidity (|y|<0.9) in the transverse-momentum interval 0 < p T < 40 GeV/c, using a minimum-bias data sample corresponding to an integrated luminosity L int =32.2 nb -1 and an Electromagnetic Calorimeter triggered data sample with L int =8.3 pb -1 . The p T -integrated J/ψ production cross section at midrapidity, computed using the minimum-bias data sample, is dσ/dy| y = 0 = 8.97±0.24 (stat)±0.48 (syst)±0.15 (lumi) μb. An approximate logarithmic dependence with the collision energy is suggested by these results and available world data, in agreement with model predictions. The integrated and p T -differential measurements are compared with measurements in pp collisions at lower energies and with several recent phenomenological calculations based on the non-relativistic QCD and Color Evaporation models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

MODEX ─ A Program for Calculation of the Energy Spectra of Particle Emitted in the Reactions of Pre-Equilibrium and Equilibrium Statistical Decays

The program is described intended for calculation of the energy distribution of particles, emitted by excited nucleus, using the Modified Exciton Model (MEM). The complete text of the program is presented in FORTRAN66 language. We reproduce here this 1974 JINR Communication as it was never published in a journal but it represents a historical interest as the first (in the world) Monte Carlo realization of the preequilibrium plus evaporation model, and was and still is used in dozens of countries all over the world. E.g., the current CEM03.03 and LAQGSM03.03 event-generators of the latest LANL Monte Carlo transport code MCNP6 have today modified parts from MODEX, the first Monte Carlo preequilibrium code in the world. [Abstract is in English. The 1974 JINR Communication is in Russian.]

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

CHROMOSPHERIC EVAPORATION IN AN X1.0 FLARE ON 2014 MARCH 29 OBSERVED WITH IRIS AND EIS

Chromospheric evaporation refers to dynamic mass motions in flare loops as a result of rapid energy deposition in the chromosphere. These motions have been observed as blueshifts in X-ray and extreme-ultraviolet (EUV) spectral lines corresponding to upward motions at a few tens to a few hundreds of km s{sup −1}. Past spectroscopic observations have also revealed a dominant stationary component, in addition to the blueshifted component, in emission lines formed at high temperatures (∼10 MK). This is contradictory to evaporation models predicting predominant blueshifts in hot lines. The recently launched Interface Region Imaging Spectrograph (IRIS) provides high-resolution imaging and spectroscopic observations that focus on the chromosphere and transition region in the UV passband. Using the new IRIS observations, combined with coordinated observations from the EUV Imaging Spectrometer, we study the chromospheric evaporation process from the upper chromosphere to the corona during an X1.0 flare on 2014 March 29. We find evident evaporation signatures, characterized by Doppler shifts and line broadening, at two flare ribbons that are separating from each other, suggesting that chromospheric evaporation takes place in successively formed flaring loops throughout the flare. More importantly, we detect dominant blueshifts in the high-temperature Fe xxi line (∼10 MK), in agreement with theoretical predictions. We also find that, in this flare, gentle evaporation occurs at some locations in the rise phase of the flare, while explosive evaporation is detected at some other locations near the peak of the flare. There is a conversion from gentle to explosive evaporation as the flare evolves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Theoretical analysis of mixing in liquid clouds – Part 3: Inhomogeneous mixing

An idealized diffusion–evaporation model of time-dependent mixing between a cloud volume and a droplet-free volume is analyzed. The initial droplet size distribution (DSD) in the cloud volume is assumed to be monodisperse. It is shown that evolution of the microphysical variables and the final equilibrium state are unambiguously determined by two non-dimensional parameters. The first one is the potential evaporation parameter R, proportional to the ratio of the saturation deficit to the liquid water content in the cloud volume, that determines whether the equilibrium state is reached at 100 % relative humidity, or is characterized by a complete evaporation of cloud droplets. The second parameter Da is the Damkölher number equal to the ratio of the characteristic mixing time to the phase relaxation time. Parameters R and Da determine the type of mixing. The results are analyzed within a wide range of values of R and Da. It is shown that there is no pure homogeneous mixing, since the first mixing stage is always inhomogeneous. The mixing type can change during the mixing process. Any mixing type leads to formation of a tail of small droplets in DSD and, therefore, to DSD broadening that depends on Da. At large Da, the final DSD dispersion can be as large as 0.2. The total duration of mixing varies from several to 100 phase relaxation time periods, depending on R and Da. The definitions of homogeneous and inhomogeneous types of mixing are reconsidered and clarified, enabling a more precise delimitation between them. The paper also compares the results obtained with those based on the classic mixing concepts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Formation trajectories of solution-processed perovskite thin films from mixed solvents

The engineering of mixed-solvent formulations and their evaporation conditions are key to reproducible perovskite coatings for high-performance photovoltaics. Here, we report a lumped-parameter evaporation model to predict the evolution of a perovskite ink liquid film over time (solvent ratio, solute concentration, and film thickness). The drying-rate model is validated via in situ film-thickness measurements, and the predicted transient liquid film state is mapped as a process path. These methods allow for the prediction of process sensitivity to local environmental factors and the understanding and visualization of a broader processing parameter space enabled through the coupling of process and ink engineering. Process maps are applied to create a new framework for scalable perovskite coating development with a goal of improving the reproducibility and transferability of perovskite fabrication. This approach is demonstrated with blade-coated FA 0.83 Cs 0.17 PbI 3 photovoltaic devices, improving the photovoltaic conversion efficiency from 17.5% ± 1.7% to 20.3% ± 0.6%.

14 SOLAR ENERGY↗
DOE OSTI2022

On the hot surface ignition of a wall-stagnating spray flame

The ignition of fuel sprays due to interaction with hot surfaces is an important phenomenon in the safety analysis of many engineering systems. We perform a parametric study of the hot surface ignition (HSI) of a fuel spray approaching a heated surface caused by the accidental leakage of a fuel line. Here, to this end, we employ a one-dimensional Eulerian-Eulerian formulation with a non-equilibrium evaporation model and a realistic chemical mechanism to describe $n$-dodecane fuel chemistry. We first describe and analyze the phenomenology of the unsteady processes leading to ignition using non-dimensionalized quantities. Through consideration of the temporal development of the most reactive mixture, we demonstrate that ignition occurs at a fuel-lean composition in a premixed region near the hot surface. Using non-dimensional parameters identified from the governing equations, we perform a parametric study of the time, location and local mixture composition at ignition and determine the ignition limits. We then identify the most important parametric sensitivities for physical analysis using a data-driven classification method. Our analysis demonstrates a contraction of the ignition limits with increased Stokes number and a regime of parametric insensitivity of igniting mixture composition. We also show that at high Damköhler numbers, the ignition location conforms to the parametric behavior of the thermal boundary layer, whereas at low Damköhler numbers approaching the ignition limit it reaches a near-unity value of the quenching Peclet number. We then compare the demonstrated parametric dependencies to the results of the quasi-steady asymptotic ignition literature, showing that our results are consistent with those obtained analytically within the limitations imposed by the simplified formulation of the latter.

42 ENGINEERING↗
DOE OSTI2022

Accessing subnuclear fluctuations and saturation with multiplicity dependent J/ψ production in p+p and p+Pb collisions

We study the production of J/ψ vector mesons as a function of charged hadron multiplicity in p+p and p+Pb collisions at LHC energies. We employ the color glass condensate framework, using running coupling Balitsky-Kovchegov evolved dipole amplitudes, to compute gluon and cc¯-pair production. We use fragmentation functions to obtain charged hadrons, and explore two different hadronization schemes for the J/ψ: non-relativistic quantum chromodynamics and the improved color evaporation model. In our framework, event-by-event multiplicity fluctuations of both hadrons and J/ψ are driven by geometric and saturation scale normalization fluctuations. Studying the correlation between J/ψ and hadron multiplicity, we show that the characteristic difference between forward and backward rapidity in p+Pb collisions is a result of different degrees of saturation probed at different rapidities. We demonstrate that experimental data on heavy-flavor production as a function of event activity provide stringent constraints on the fluctuating proton structure

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Volatility characterization of nanoparticles from single and dual-fuel low temperature combustion in compression ignition engines

Our work explores the volatility of particles produced from two diesel low temperature combustion (LTC) modes proposed for high-efficiency compression ignition engines. It also explores mechanisms of particulate formation and growth upon dilution in the near-tailpipe environment. Moreover, the number distribution of exhaust particles from low- and mid-load dual-fuel reactivity controlled compression ignition (RCCI) and single-fuel premixed charge compression ignition (PPCI) modes were experimentally studied over a gradient of dilution temperature. Particle volatility of select particle diameters was investigated using volatility tandem differential mobility analysis (V-TDMA). Evaporation rates for exhaust particles were compared with V-TDMA results for candidate pure n-alkanes to identify species with similar volatility characteristics. The results show that LTC particles are mostly comprised of material with volatility similar to engine oil alkanes. V-TDMA results were used as inputs to an aerosol condensation and evaporation model to support the finding that smaller particles in the distribution are comprised of lower volatility material than large particles under primary dilution conditions. Although the results show that saturation levels are high enough to drive condensation of alkanes onto existing particles under the dilution conditions investigated, they are not high We conclude that observed particles from LTC operation must grow from low concentrations of highly non-volatile compounds present in the exhaust.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2012

THE STRUCTURE AND SPECTRAL FEATURES OF A THIN DISK AND EVAPORATION-FED CORONA IN HIGH-LUMINOSITY ACTIVE GALACTIC NUCLEI

We investigate the accretion process in high-luminosity active galactic nuclei (HLAGNs) in the scenario of the disk evaporation model. Based on this model, the thin disk can extend down to the innermost stable circular orbit (ISCO) at accretion rates higher than 0.02 M-dot{sub Edd} while the corona is weak since part of the coronal gas is cooled by strong inverse Compton scattering of the disk photons. This implies that the corona cannot produce as strong X-ray radiation as observed in HLAGNs with large Eddington ratio. In addition to the viscous heating, other heating to the corona is necessary to interpret HLAGN. In this paper, we assume that a part of accretion energy released in the disk is transported into the corona, heating up the electrons, and is thereby radiated away. For the first time, we compute the corona structure with additional heating, fully taking into account the mass supply to the corona, and find that the corona could indeed survive at higher accretion rates and that its radiation power increases. The spectra composed of bremsstrahlung and Compton radiation are also calculated. Our calculations show that the Compton-dominated spectrum becomes harder with the increase of energy fraction (f) liberating in the corona, and the photon index for hard X-ray (2-10 keV) is 2.2 < {Gamma} < 2.7. We discuss possible heating mechanisms for the corona. Combining the energy fraction transported to the corona with the accretion rate by magnetic heating, we find that the hard X-ray spectrum becomes steeper at a larger accretion rate and the bolometric correction factor (L{sub bol}/L{sub 2-10keV}) increases with increasing accretion rate for f < 8/35, which is roughly consistent with the observational results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Impacts of Pasture Conversion to Sugarcane on Water Fluxes and Water Use Efficiency in the Southeastern US

The expansion of sugarcane (cane), a high-yielding perennial crop, will likely reshape the bioenergy landscape in the Southeastern US. However, its ecohydrological implications, particularly following conversion from grazed pastures, a dominant land use in the region, remain highly uncertain. We investigated the impact of cane expansion on evapotranspiration (ET) and its partitioning, and the mechanisms influencing both ET components and water use efficiency (WUE) across multiple scales and growth cycles in subtropical Florida. We combined eddy covariance, biometric measurements, and process-based stomatal conductance (g s ) models. ET was 1.7% lower in cane than in improved pasture (IMP) but exceeded that in semi-native pasture (SN) by 21%. Transpiration (T) followed a similar pattern, consistent with lower g s in cane relative to IMP. Cane had more conservative water use and greater sensitivity of g s to vapor pressure deficit (VPD) compared to IMP pasture, suggesting cane may be more tolerant of increasing atmospheric water demand. In contrast, SN showed lower g s and weaker stomatal sensitivity to VPD compared to cane, resulting in lower T. In cane, stomatal regulation and T varied across growth cycles, with stomata becoming less water conservative as stands matured, highlighting the importance of incorporating stand age-dependent stomatal regulation into hydrological models. Evaporation (E) was higher in cane than pastures (19%–26%), partially offsetting WUE gains. Cane exhibited higher intrinsic WUE (GPP/g s ; Gross Primary Productivity), ecosystem WUE (GPP/ET), and harvest WUE (harvest/ET) than both pasture types. Large-scale pasture-to-cane conversion could produce widely contrasting hydrological outcomes. The net regional impact will depend on the proportion of each pasture type converted and on cane's high g s sensitivity to VPD, which triggers tight stomatal regulation and conservative water use, both of which will become increasingly consequential under intensifying atmospheric water demand.

bioenergy↗
DOE OSTI2016

TU-AB-BRC-09: Fast Dose-Averaged LET and Biological Dose Calculations for Proton Therapy Using Graphics Cards

Purpose: To demonstrate fast and accurate Monte Carlo (MC) calculations of proton dose-averaged linear energy transfer (LETd) and biological dose (BD) on a Graphics Processing Unit (GPU) card. Methods: A previously validated GPU-based MC simulation of proton transport was used to rapidly generate LETd distributions for proton treatment plans. Since this MC handles proton-nuclei interactions on an event-by-event using a Bertini intranuclear cascade-evaporation model, secondary protons were taken into account. The smaller contributions of secondary neutrons and recoil nuclei were ignored. Recent work has shown that LETd values are sensitive to the scoring method. The GPU-based LETd calculations were verified by comparing with a TOPAS custom scorer that uses tabulated stopping powers, following recommendations by other authors. Comparisons were made for prostate and head-and-neck patients. A python script is used to convert the MC-generated LETd distributions to BD using a variety of published linear quadratic models, and to export the BD in DICOM format for subsequent evaluation. Results: Very good agreement is obtained between TOPAS and our GPU MC. Given a complex head-and-neck plan with 1 mm voxel spacing, the physical dose, LETd and BD calculations for 10{sup 8} proton histories can be completed in ∼5 minutes using a NVIDIA Titan X card. The rapid turnover means that MC feedback can be obtained on dosimetric plan accuracy as well as BD hotspot locations, particularly in regards to their proximity to critical structures. In our institution the GPU MC-generated dose, LETd and BD maps are used to assess plan quality for all patients undergoing treatment. Conclusion: Fast and accurate MC-based LETd calculations can be performed on the GPU. The resulting BD maps provide valuable feedback during treatment plan review. Partially funded by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗