Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Stochastic joint inversion of hydrogeophysical data for salt tracer test monitoring and hydraulic conductivity imaging

The assessment of hydraulic conductivity of heterogeneous aquifers is a difficult task using traditional hydrogeological methods (e.g., steady state or transient pumping tests) due to their low spatial resolution. Geophysical measurements performed at the ground surface and in boreholes provide additional information for increasing the resolution and accuracy of the inverted hydraulic conductivity field. We used a stochastic joint inversion of Direct Current (DC) resistivity and Self-Potential (SP) data plus in situ measurement of the salinity in a downstream well during a synthetic salt tracer experiment to reconstruct the hydraulic conductivity field between two wells. The pilot point parameterization was used to avoid over-parameterization of the inverse problem. Bounds on the model parameters were used to promote a consistent Markov chain Monte Carlo sampling of the model parameters. To evaluate the effectiveness of the joint inversion process, we compared 8 cases in which the geophysical data are coupled or not to the in situ sampling of the salinity to map the hydraulic conductivity. We first tested the effectiveness of the inversion of each type of data alone (concentration sampling, self-potential, and DC resistivity), and then we combined the data two by two. We finally combined all the data together to show the value of each type of geophysical data in the joint inversion process because of their different sensitivity map. We also investigated a case in which the data were contaminated with noise and the variogram unknown and inverted stochastically. The results of the inversion revealed that incorporating the self-potential data improves the estimate of hydraulic conductivity field especially when the self-potential data were combined to the salt concentration measurement in the second well or to the time-lapse cross-well electrical resistivity data. Various tests were also performed to quantify the uncertainty in the inverted hydraulic conductivity field.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Evaluation of the Community Land Model simulated carbon and water fluxes against observations over ChinaFLUX sites

The Community Land Model (CLM) is an advanced process-based land surface model that simulates carbon, nitrogen, water vapor and energy exchanges between terrestrial ecosystems and the atmosphere at various spatial and temporal scales. We use observed carbon and water fluxes from five representative Chinese Terrestrial Ecosystem Flux Research Network (ChinaFLUX) eddy covariance tower sites to systematically evaluate the new version CLM4.5 and old version CLM4.0, and to generate insights that may inform future model developments. CLM4.5 underestimates the annual carbon sink at three forest sites and one alpine grassland site but overestimates the carbon sink of a semi-arid grassland site. The annual carbon sink underestimation for the deciduous-dominated forest site results from underestimated daytime carbon sequestration during summer and overestimated nighttime carbon emission during spring and autumn. Compared to CLM4.0, the bias of annual gross primary production (GPP) is reduced by 24% and 28% in CLM4.5 at two subtropical forest sites. However, CLM4.5 still presents a large positive bias in annual GPP. The improvement in net ecosystem exchange (NEE) is limited, although soil respiration bias decreases by 16%–43% at three forest sites. CLM4.5 simulates lower soil water content in the dry season than CLM4.0 at two grassland sites. Drier soils produce a significant drop in the leaf area index and in GPP and an increase in respiration for CLM4.5. The new fire parameterization approach in CLM4.5 causes excessive burning at the Changbaishan forest site, resulting in an unexpected underestimation of NEE, vegetation carbon, and soil organic carbon by 46%, 95%, and 87%, respectively. Altogether, our study reveals significant improvements achieved by CLM4.5 compared to CLM4.0, and suggests further developments on the parameterization of seasonal GPP and respiration, which will require a more effective representation of seasonal water conditions and the partitioning of net radiation between sensible and heat fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Photosynthesis phenology, as defined by solar-induced chlorophyll fluorescence, is overestimated by vegetation indices in the extratropical Northern Hemisphere

Vegetation phenology is highly sensitive to climate change, although the data and methods used to estimate key phenological states can influence this sensitivity. Because of its direct relation to leaf photosynthetic carbon uptake, remotely sensed solar-induced chlorophyll fluorescence (SIF) can provide new insight assessing changes in vegetation phenology. In this work, we investigated the potential of using a SIF time series product named contiguous SIF (CSIF) to estimate spring, summer, and autumn phenology in the extratropical Northern Hemisphere (>30°N) and compared the results with those based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) for the period 2001–2018. Overall, we found similar spatial patterns in phenological states. However, specific dates of key phenological events differed when using CSIF vs. MODIS NDVI data. NDVI data indicated that the growing season started earlier (by 10.1 days on average) and ended later (11.5 days on average) relative to CSIF data. This implies that actual periods for photosynthetic activity are shorter (by 21.6 days on average) than those estimated from vegetation indices more directly related to changes in canopy structure. These large differences between results from NDVI and that from CSIF suggest that vegetation indices such as NDVI seem to overestimate the period for active photosynthesis over the extratropical Northern Hemisphere. Furthermore, while phenology of the early growing season is dominated by temperature for both NDVI and CSIF data, phenology of the late growing season is mainly controlled by temperature for NDVI but by precipitation for CSIF. Our findings were further confirmed by other SIF (GOME-2 SIF) and vegetation index (MODIS EVI) datasets. Phenology modes in Earth system modelling are often parameterized using leaf unfolding and senescence from either station or satellite observations. Our results imply that canopy structure-based parameterization schemes may have overestimated photosynthesis active period, and thus productivity responses. We conclude that SIF data offers a novel and unique approach for assessing phenological change - one that is more directly tied to the carbon cycle and how it is being influenced by climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

The double ITCZ syndrome in GCMs: A coupled feedback problem among convection, clouds, atmospheric and ocean circulations

The appearance of a spurious double ITCZ south of the equator in coupled general circulation models has been a stubborn problem ever since the beginning of coupled model development. This article reviews the past research in this area, with a focus on three possible major contributors to the double ITCZ biases: 1) the southeastern Pacific marine stratus clouds and associated warm sea surface temperature (SST) biases; 2) the extratropical shortwave absorption biases over Southern Ocean; and 3) convective parameterization. The negative biases in marine boundary layer clouds in the southeastern Pacific lead to warm SST biases. Through coupled atmosphere-ocean interactions, it contributes to the double ITCZ bias. Positive shortwave absorption biases over Southern Ocean is believed to be an extratropical contributor to the double ITCZ biases in models. Since the heat from excess shortwave absorption must be transported to the northern hemisphere by the Hadley circulation, the position of the ITCZ must shift southward. However, later research finds that reducing the extratropical shortwave absorption has little effect on the ITCZ position. One possibility is that some feedback mechanism involving subtropical low-level clouds as intermediary is missing. For convective parameterization, changes in various elements in convection schemes can have large impacts on the ITCZ simulation. It involves a complex chain of interactions among convection, large-scale atmospheric circulation, SST, and upper ocean circulation. In one scheme, changes in the scheme lead to the elimination of double ITCZ in all seasons, giving hope that we can finally take on the double ITCZ problem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Atmospheric ice nuclei concentration measurements over a high altitude-station in the Western Ghats, India

Mixed-phase and ice clouds play a significant role in modulating the Indian summer monsoon rainfall, and modeling studies suggest that inaccurate representation of ice nucleating particles (INP) is one of the factors, which have resulted in a deficiency in the prediction of cloud properties and hydrological cycle. We do not understand the ice nucleation mechanisms to represent the ice formation in the model, and therefore more INP measurements at various supercooled temperatures are needed. This study reports for the first time real-time INP concentration at a high altitude station in the Western Ghats region of India, and these results will be used to constrain the ice nucleation parameterizations in the model. The variation of INP for any particular cloud condition can have strong microphysical responses of mixed-phase precipitation formation processes. Also, to improve the representation of heterogeneous ice nucleation, INP measurements of various types of aerosols commonly observed over the Indian region are needed. This study presents, the INP concentrations and daily concentration of non-refractory chemical composition (organics, sulfate, nitrate, ammonium, and chloride), cloud condensation nuclei (CCN), and aerosol size distribution. Spectrometer for Ice Nuclei (SPIN) was operated to measure INP concentration at the mountain site from August to December 2018. INP measurements were performed at three different temperatures (-25, -30, and - 34 °C) in the immersion freezing mode relevant for mixed-phase cloud conditions. The average INP concentrations approximately varied from 0.18 to 12.4 L-1, 0.39 to 24 L-1, 1.1 to 40.2 L-1 at -25, -30, and - 34 °C, respectively. These concentrations are within the concentrations reported in the literature from different regions across the globe. The air mass back trajectory analysis showed that continental air mass contains more INP compared to maritime air mass. The information on INP concentration and chemical characteristics of aerosol will help to improve heterogeneous ice nucleation parameterization in numerical models.

Kumar, Anil V.↗
DOE OSTI2020

An efficient method to integrate polynomials over polytopes and curved solids

Here in this paper, we present an efficient approach to compute the integral of monomials and polynomials over polyhedra and regions defined by parametric curved boundary surfaces. We use Euler's theorem for homogeneous functions in combination with Stokes's theorem to reduce the integration of a monomial over a three-dimensional solid to its boundary. If the solid is a polytope, through a recursive application of these theorems, the integral is further reduced to just the evaluation of the monomial and its derivatives at the vertices of the polytope. The present approach is simpler than existing techniques that rely on repeated use of the divergence theorem, which require the antiderivative of the monomials and the projection of these functions onto hyperplanes. For convex and nonconvex polytopes, our approach does not introduce any approximation for the integration of monomials. For curved solid regions bounded by surfaces that admit a parameterization, the same approach yields simplified formulas to compute the integral of any homogeneous function, including monomials. For surfaces parameterized by polynomial surfaces (such as Bezier surface triangles and B-spline patches), the method yields machine-precision accuracy for the volumetric integration of monomials with an appropriate quadrature rule. Numerical examples over regions bounded by polynomial surfaces and rational surfaces are presented to establish the accuracy and efficiency of the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Optimal toolpath design of additive manufactured composite cylindrical structures

Structural design optimization for additive manufacturing is primarily focused on planar layer-by-layer processes and design of composite cylindrical structures does not often accommodate manufacturing constraints. In this study, we propose to optimize the toolpath trajectory of additively manufactured composite cylinders comprised of multiple thin cylindrical annuli. Additionally, our printing process is based on direct ink writing in which short carbon fibers in a thermoset resin are extruded through a moving nozzle to build up a structure. Since the fibers are primarily aligned in the flow direction of the extrudate, the printing trajectory influences the material properties of the composite structure. To accommodate this, we define the toolpaths as contours of parameterized level-set functions. This parameterization allows us to optimize the material properties and impose manufacturing constraints such as no-overlap, no-sag, minimum radius of curvature, and continuity of the toolpaths. Several examples illustrate the optimization procedure.

42 ENGINEERING↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multifidelity uncertainty quantification with models based on dissimilar parameters

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. Here we present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

42 ENGINEERING↗
DOE OSTI2023

Learning subgrid-scale models with neural ordinary differential equations

We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Thermo4PFM: Facilitating Phase-field simulations of alloys with thermodynamic driving forces

Phase-field modeling is a popular front-tracking approach used to model solidification. Its time-evolution equations are often coupled to alloy composition and/or thermal diffusion in high-resolution multiphysics approaches. Materials thermodynamic properties tabulated in CALPHAD databases can be used for phase-field modeling to parameterize bulk energies of alloys. In addition, they can be naturally integrated into models such as the Kim-Kim-Suzuki (KKS) model where driving forces depend on the differences between chemical potentials of co-existing phases. In that case, a small system of coupled nonlinear equations needs to be solved at every point in space where the phase-field order parameter is to be updated and evolved in time. Here we present Thermo4PFM, a solver for the KKS equations for binary and ternary alloys, with two or three phases, and parameterized with CALPHAD models. Thermo4PFM is open source, written in C++, and can take advantage of Graphics Processing Units (GPU) accelerators. Using OpenMP offload capabilities for C++ classes, an excellent performance is demonstrated on GPU using the LLVM compiler. CALPHAD data is read from simple JSON files using an open source parser from the boost library.

36 MATERIALS SCIENCE↗
DOE OSTI2022

QUIC-URB and QUIC-fire extension to complex terrain: Development of a terrain-following coordinate system

Ensemble-based approaches to prescribed fire planning cannot be supported by CFD-based models like FIRETEC and WFDS because they are too computationally expensive and cannot leverage LES approaches like CAWFE and WRF-SFIRE because too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, typical for instance of the Western United States. In this paper, we describe the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. In particular, the paper presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captured the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. We report the model agreed well with measurements on the upwind side of the peak, but overestimated speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field was obtained in 6 s, making the solver appealing for coupled fire–atmosphere simulations. The Socorro Mountain was highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect was addressed in this work). Results highlight that future work should focus on the implementation of parameterizations of wake-eddies, similar to QUIC-URB’s building parameterizations, and on thermodynamic-driven flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Molecular thermodynamic modeling of aqueous Na + -K + -Mg 2+ -Ca 2+ -SO 4 2- quinary system

As part of the effort to develop a comprehensive molecular thermodynamic model for high salinity produced water based on the electrolyte nonrandom two-liquid theory, this work presents the model results for the aqueous Na + -K + -Mg 2+ -Ca 2+ -SO 4 2- quinary system. With the binary interaction parameters fully parameterized for the water-electrolyte pairs and the electrolyte-electrolyte pairs, the model satisfactorily describes the thermodynamic properties for the aqueous quinary system and its binary, ternary, and quaternary subsystems for the temperature range 273.15–473.15 K and the salt concentrations up to saturation. The model results are compared with the salt solubility and thermodynamic data available for the subsystems with special attention to systems involving (Ca 2+ SO 4 2- ).

Chemistry↗
DOE OSTI2018

The dissolution behavior of borosilicate glasses in far-from equilibrium conditions

In far-from-equilibrium conditions, the dissolution of borosilicate glasses used to immobilize nuclear waste is known to be a function of both temperature and pH. The aim of this paper is to study effects of these variables on three model waste glasses (SON68, ISG, AFCI). To do this, experiments were conducted at temperatures of 23, 40, 70, and 90 °C and pH(RT) values of 9, 10, 11, and 12 with the single-pass flow-through (SPFT) test method. The results from these tests were then used to parameterize a kinetic rate model based on transition state theory. Both the absolute dissolution rates and the rate model parameters are compared with previous results. Discrepancies in the absolute dissolution rates as compared to those obtained using other test methods are discussed. Rate model parameters for the three glasses studied here are nearly equivalent within error and in relative agreement with previous studies. The results were analyzed with a linear multivariate regression (LMR) and a nonlinear multivariate regression performed with the use of the Glass Corrosion Modeling Tool (GCMT), which is capable of providing a robust uncertainty analysis. This robust analysis highlights the high degree of correlation of various parameters in the kinetic rate model. As more data are obtained on borosilicate glasses with varying compositions, the effect of glass composition on the rate parameter values could possibly be obtained. This would allow for the possibility of predicting the forward dissolution rate of glass based solely on composition

Neeway, James J.↗
DOE OSTI2023

Uncertain characterization of reservoir fluids due to brittleness of equation of state regression

Equations of state (EoS) play a central role in modeling the phase equilibrium of fluid mixtures. Their parameterization involves fitting a model to experimental data, i.e., solving a nonlinear, non-convex, multivariate optimization problem. The latter requires one to select design variables, domains of definition for each variable, and weights assigned to individual measurements. We demonstrate that subjective choices of an optimization algorithm and an initial guess also impact the regression process. Consequently, EoS predictions are fundamentally uncertain even after the EoS tuning to a limited set of experimental data points. We demonstrate this observation for two hydrocarbon reservoir fluids, in which five properties of the heaviest carbon fraction are treated as design variables. While all the optimization algorithms and initial guesses match experimental data for the gas and liquid properties, the resulting EoS parameterizations lead to dramatically different predictions of the fluid’s thermophysical behavior in the unsampled pressure and temperature regions. In conclusion, we propose the probabilistic treatment of design variables to quantify the predictive uncertainty of the resulting fluid models.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Comparative verification of hydro-mechanical fracture behavior: Task G of international research project DECOVALEX–2023

Numerical simulations become a necessity when experimental approaches cannot cover the required physical and time scale of interest. One of such area is a simulation of long-term host rock behaviors for nuclear waste disposal and simulation tools involved in the assessment must go through rigorous validation tests. The DECOVALEX project (Development of COupled models and their VAlidation against EXperiments) is dedicated to this purpose by international participants. This work is part of the ongoing phase DECOVALEX–2023 (D–2023, Task G) particularly aiming to simulate fracture behaviors under various conditions. Here, in this study, we cross-verified a variety of numerical methods including continuous and discontinuous approaches against four benchmark exercises with emphasis on numerical accuracy and parameterization of the various numerical approaches. The systematic inter-comparisons of test cases highlight advantages and disadvantages of the different numerical models. Numerical details on discretization effects (e.g. mesh density and orientation) and domain size were investigated in detail for practical applications. It became evident that meticulous attention to mesh resolution and domain size is imperative for achieving accurate numerical simulations, even for static cracks. Moreover, when comparing numerical methods to closed-form solutions for static cracks, all models successfully reproduced the maximum crack opening but encountered challenges near the crack tips. Finally, the paper discusses how to convert between and therefore compare parameters of various numerical approaches. Our benchmark studies reveal that each model necessitates a distinct number of parameters, even in simple scenarios like static crack aperture benchmarks. It is generally more practical to employ fewer parameters to mitigate model over-parameterization and enhance experimental feasibility.

58 GEOSCIENCES↗
DOE OSTI2016

Investigating the Implications of a Variable RBE on Proton Dose Fractionation Across a Clinical Pencil Beam Scanned Spread-Out Bragg Peak

Purpose: To investigate the clinical implications of a variable relative biological effectiveness (RBE) on proton dose fractionation. Using acute exposures, the current clinical adoption of a generic, constant cell killing RBE has been shown to underestimate the effect of the sharp increase in linear energy transfer (LET) in the distal regions of the spread-out Bragg peak (SOBP). However, experimental data for the impact of dose fractionation in such scenarios are still limited. Methods and Materials: Human fibroblasts (AG01522) at 4 key depth positions on a clinical SOBP of maximum energy 219.65 MeV were subjected to various fractionation regimens with an interfraction period of 24 hours at Proton Therapy Center in Prague, Czech Republic. Cell killing RBE variations were measured using standard clonogenic assays and were further validated using Monte Carlo simulations and parameterized using a linear quadratic formalism. Results: Significant variations in the cell killing RBE for fractionated exposures along the proton dose profile were observed. RBE increased sharply toward the distal position, corresponding to a reduction in cell sparing effectiveness of fractionated proton exposures at higher LET. The effect was more pronounced at smaller doses per fraction. Experimental survival fractions were adequately predicted using a linear quadratic formalism assuming full repair between fractions. Data were also used to validate a parameterized variable RBE model based on linear α parameter response with LET that showed considerable deviations from clinically predicted isoeffective fractionation regimens. Conclusions: The RBE-weighted absorbed dose calculated using the clinically adopted generic RBE of 1.1 significantly underestimates the biological effective dose from variable RBE, particularly in fractionation regimens with low doses per fraction. Coupled with an increase in effective range in fractionated exposures, our study provides an RBE dataset that can be used by the modeling community for the optimization of fractionated proton therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗