Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accountability”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Leveraging Extremal Dependence to Better Characterize the 2021 Pacific Northwest Heatwave

Abstract In late June, 2021, a devastating heatwave affected the US Pacific Northwest and western Canada, breaking numerous all-time temperature records by large margins and directly causing hundreds of fatalities. The observed 2021 daily maximum temperature across much of the U.S. Pacific Northwest exceeded upper bound estimates obtained from single-station temperature records even after accounting for anthropogenic climate change, meaning that the event could not have been predicted under standard univariate extreme value analysis assumptions. In this work, we utilize a flexible spatial extremes model that considers all stations across the Pacific Northwest domain and accounts for the fact that many stations simultaneously experience extreme temperatures. Our analysis incorporates the effects of anthropogenic forcing and natural climate variability in order to better characterize time-varying changes in the distribution of daily temperature extremes. We show that greenhouse gas forcing, drought conditions and large-scale atmospheric modes of variability all have significant impact on summertime maximum temperatures in this region. Our model represents a significant improvement over corresponding single-station analysis, and our posterior medians of the upper bounds are able to anticipate more than 96% of the observed 2021 high station temperatures after properly accounting for extremal dependence. Supplementary materials accompanying this paper appear online.

Zhang, Likun (ORCID:0000000154584556)↗
DOE OSTI2022

Aboveground and belowground contributions to ecosystem respiration in a temperate deciduous forest

In this study, we developed a three-way carbon dioxide (CO 2 ) flux-partitioning algorithm that separates net ecosystem exchange (NEE) into aboveground plant respiration (R above ), belowground root and soil respiration (R below ), and gross primary production (GPP). We applied this algorithm to a coupled dataset of continuous chamber-measured soil respiration and eddy covariance (EC)-measured NEE of CO 2 in an oak-hickory (Quercus-Carya) deciduous broadleaf forest from 2006 to 2015. We found that on annual time scale, R below dominated over R above with the former accounting for 66.9–86.4% and the latter 13.6–33.1%, of the total ecosystem respiration (R eco ). The ratio of R below to R above varied seasonally, ranging from 1.77 to 7.25 in growing season, and 1.02 to 4.57 in non-growing season. The temperature sensitivity (E 0 ) of R below was significantly higher than that of R above , and E 0 of R eco responded differently to air and soil temperature. Over the whole study period, annual mean R above , R below , and GPP were 243, 806, and 1170 g C m –2 , respectively, with annual R eco accounting for 89.6% of GPP, of which 68.8% was lost as R below and 20.8% lost as R above , and leaving only 10% of the carbon fixation in ecosystems. Furthermore, these estimates, however, did not consider potential light inhibition of leaf respiration. If we accept the presence of light inhibition, then the daytime three-way partitioning method would underestimate annual R above by 20.4% whereas the nighttime method would overestimate R above by 23.9% and GPP by 4.7%, compared with estimates accounting for light inhibition in leaves.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

An efficient parallel sampling technique for Multivariate Poisson-Lognormal model: Analysis with two crash count datasets

Our study investigates the Multivariate Poisson-lognormal (MVPLN) model that jointly models crash frequency and severity accounting for correlations. The ordinary univariate count models analyze crashes of different severity level separately ignoring the correlations among severity levels. The MVPLN model is capable to incorporate the general correlation structure and takes account of the over dispersion in the data that leads to a superior data fitting. But, the traditional estimation approach for MVPLN model is computationally expensive, which often limits the use of MVPLN model in practice. In this work, a parallel sampling scheme is introduced to improve the original Markov Chain Monte Carlo (MCMC) estimation approach of the MVPLN model, which significantly reduces the model estimation time. Two MVPLN models are developed using the pedestrian vehicle crash data collected in New York City from 2002 to 2006, and the highway-injury data from Washington State (5-year data from 1990 to 1994) The Deviance Information Criteria (DIC) is used to evaluate the model fitting. The estimation results show that the MVPLN models provide a superior fit over univariate Poisson-lognormal (PLN), univariate Poisson, and Negative Binomial models. Moreover, the correlations among the latent effects of different severity levels are found significant in both datasets that justifies the importance of jointly modeling crash frequency and severity accounting for correlations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Ontological models with epistemic constraints: Local reconstruction of a dual-rail qubit

Restrictions on gaining information may significantly affect perception of the system under study. In this paper we are concerned with mechanisms introducing epistemic constraints in the framework of ontological models and discuss a scheme leading to the minimal operational account which does not rely on the knowledge of the underlying ontology. As a case in point, we construct local ontological model which simulates behaviour of a single quantum particle in dual-rail interferometric circuits. We show that for agents analysing conceivable experimental circuits predictions of the model are indistinguishable from their quantum mechanical counterparts, that is operational account of the model is equivalent to a qubit. Note that it contrasts with the usual understanding the collapse of the wave function as non-local phenomenon, namely the model is local from the construction with any form of ‘non-locality’ being explained as an epistemic effect due to incomplete information about the system. This illustrates flexibility of the ontological model framework, and in particular the possibility of different interpretations behind the same operational account. - Highlights: • Effects of restricted access to information about a system are discussed. • Mechanisms of epistemic constraints in the ontological model framework. • Operational description from the perspective of agents with restricted resources. • Construction of a local ontological model of dual-rail interferometric circuits. • Predictions of the model are indistinguishable from the quantum behaviour of a qubit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

From chiral anomaly to two-fluid hydrodynamics for electronic vortices

Highlights: • Chiral transformations with account for chiral anomaly were applied to the CDW model. • Two-fluid hydrodynamics was constructed of the order parameter and the normal liquid. • Topologically nontrivial dynamics appears under applied fields or charge injection. • Spontaneous spacetime vorticity under the driving current was modeled numerically. • Periodic sequence of amplitude nodes with 2π phase increments was demonstrated. - Abstract: Many recent experiments addressed manifestations of electronic crystals, particularly the charge density waves, in nano-junctions, under electric field effect, at high magnetic fields, together with real space visualizations by STM and micro X-ray diffraction. This activity returns the interest to stationary or transient states with static and dynamic topologically nontrivial configurations: electronic vortices as dislocations, instantons as phase slip centers, and ensembles of microscopic solitons. Describing and modeling these states and processes calls for an efficient phenomenological theory which should take into account the degenerate order parameter, various kinds of normal carriers and the electric field. Here we notice that the commonly employed time-depend Ginzburg–Landau approach suffers with violation of the charge conservation law resulting in unphysical generation of particles which is particularly strong for nucleating or moving electronic vortices. We present a consistent theory which exploits the chiral transformations taking into account the principle contribution of the fermionic chiral anomaly to the effective action. The resulting equations clarify partitions of charges, currents and rigidity among subsystems of the condensate and normal carriers. On this basis we perform the numerical modeling of a spontaneously generated coherent sequence of phase slips – the spacetime vortices – serving for the conversion among the injected normal current and the collective one.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Closing the Carbon Balance for Fermentation by Clostridium thermocellum (ATCC 27405)

Our lab and most others have not been able to close a carbon balance for fermentation by the thermophilic, cellulolytic anaerobe, Clostridium thermocellum. We undertook a detailed accounting of product formation in C. thermocellum ATCC 27405. Elemental analysis revealed that for both cellulose (Avicel) and cellobiose, {>=}92% of the substrate carbon utilized could be accounted for in the pellet, supernatant and off-gas when including sampling. However, 11.1% of the original substrate carbon was found in the liquid phase and not in the form of commonly-measured fermentation products - ethanol, acetate, lactate, and formate. Further detailed analysis revealed all the products to be <720 da and have not usually been associated with C. thermocellum fermentation, including malate, pyruvate, uracil, soluble glucans, and extracellular free amino acids. By accounting for these products, 92.9% and 93.2% of the final product carbon was identified during growth on cellobiose and Avicel, respectively.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Effects of redox cycling of iron in nontronite on reduction of technetium

In situ technetium-99 ( 99 Tc) immobilization by Fe(II) associated with clay minerals is a potential cost-effective method for Tc remediation at contaminated sites. The oxidation state of Fe is constantly cycled in sedimentary environments, however the effect of this cycle on Tc reduction and immobilization has not yet been investigated. The objective of this project was therefore to study how multiple cycles of reduction-reoxidation of Fe-rich clay mineral, nontronite, affected its reactivity toward Tc (VII) reduction. Iron-rich nontronite NAu-2 was used as a model clay mineral. NAu-2 suspension was first bioreduced by Shewanella putrefaciens CN32, which was subsequently re-oxidized by air. Three cycles of reduction-oxidation were conducted and bioreduced NAu-2 samples from all three cycles were collected and used for Tc(VII) reduction experiments. Each redox cycle resulted in a small fraction of dissolution of small size and/or poorly crystalline NAu-2. The released Fe(II) from the dissolution was likely adsorbed onto NAu-2 surface/edge sites with a high reactivity. Upon exposure to O2, this reactive Fe(II) fraction was oxidized more rapidly than structural Fe(II) and may have accounted for a two-step reoxidation kinetics of NAu-2 associated Fe(II): rapid oxidation over first few hours followed by slow oxidation. Progressive increase of this reactive fraction of Fe(II), from increased dissolution, accounted for the successively higher rate of bioreduction and reoxidation with increased redox cycles. The same Fe redistribution accounted for two-step Tc(VII) reduction kinetics as well. Rapid Tc(VII) reduction in the first few hours may be attributed to a small fraction of highly reactive Fe(II) at the NAu-2 surface/edge sites, and more steady Tc(VII) reduction over longer time may be carried out by structural Fe(II). Similar to the increased rates of Fe(III) reduction and Fe(II) oxidation, the Tc(VII) reduction rate also increased with redox cycles and could be explained by progressive increase of the reactive Fe(II) on NAu-2 surface/edges. Iron-rich clay minerals undergo important changes after redox cycles, but eventually reach a steady state with continued reactivity toward heavy metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A MPET 2 -mPBPK model for subcutaneous injection of biotherapeutics with different molecular weights: From local scale to whole-body scale

Subcutaneous injection of biotherapeutics has attracted considerable attention in the pharmaceutical industry. However, there is limited understanding of the mechanisms underlying the absorption of drugs with different molecular weights and the delivery of drugs from the injection site to the targeted tissue. Here, we propose the MPET 2 -mPBPK model to address this issue. This multiscale model couples the MPET 2 model, which describes subcutaneous injection at the local tissue scale from a biomechanical view, with a post-injection absorption model at injection site and a minimal physiologically-based pharmacokinetic (mPBPK) model at whole-body scale. Utilizing the principles of tissue biomechanics and fluid dynamics, the local MPET 2 model provides solutions that account for tissue deformation and drug absorption in local blood vessels and initial lymphatic vessels during injection. Additionally, we introduce a model accounting for the molecular weight effect on the absorption by blood vessels, and a nonlinear model accounting for the absorption in lymphatic vessels. The post-injection model predicts drug absorption in local blood vessels and initial lymphatic vessels, which are integrated into the whole-body mPBPK model to describe the pharmacokinetic behaviors of the absorbed drug in the circulatory and lymphatic system. We establish a numerical model which links the biomechanical process of subcutaneous injection at local tissue scale and the pharmacokinetic behaviors of injected biotherapeutics at whole-body scale. With the help of the model, we propose an explicit relationship between the reflection coefficient and the molecular weight and predict the bioavalibility of biotherapeutics with varying molecular weights via subcutaneous injection. The considered drug absorption mechanisms enable us to study the differences in local drug absorption and whole-body drug distribution with varying molecular weights. This model enhances the understanding of drug absorption mechanisms and transport routes in the circulatory system for drugs of different molecular weights, and holds the potential to facilitate the application of computational modeling to drug formulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

The role of correlations in uncertainty quantification of transportation relevant fuel models

Large reaction mechanisms are often used to describe the combustion behavior of transportation-relevant fuels like gasoline, where these are typically represented by surrogate blends, e.g., n-heptane/iso-octane/toluene. We describe efforts to quantify the uncertainty in the predictions of such mechanisms at realistic engine conditions, seeking to better understand the robustness of the model as well as the important reaction pathways and their impacts on combustion behavior. In this work, we examine the importance of taking into account correlations among reactions that utilize the same rate rules and those with multiple product channels on forward propagation of uncertainty by Monte Carlo simulations. Automated means are developed to generate the uncertainty factor assignment for a detailed chemical kinetic mechanism, by first uniquely identifying each reacting species, then sorting each of the reactions based on the rate rule utilized. Simulation results reveal that in the low temperature combustion regime for iso-octane, the majority of the uncertainty in the model predictions can be attributed to low temperature reactions of the fuel sub-mechanism. The foundational, or small-molecule chemistry (C 0 -C 4 ) only contributes significantly to uncertainties in the predictions at the highest temperatures (Tc=900 K). Accounting for correlations between important reactions is shown to produce non-negligible differences in the estimates of uncertainty. Including correlations among reactions that use the same rate rules increases uncertainty in the model predictions, while accounting for correlations among reactions with multiple branches decreases uncertainty in some cases. Significant non-linear response is observed in the model predictions depending on how the probability distributions of the uncertain rate constants are defined.Finally, we concluded that care must be exercised in defining these probability distributions in order to reduce bias, and physically unrealistic estimates in the forward propagation of uncertainty for a range of UQ activities.

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

Multimodal Bayesian registration of noisy functions using Hamiltonian Monte Carlo

Functional data registration is a necessary processing step for many applications. The observed data can be inherently noisy, often due to measurement error or natural process uncertainty; which most functional alignment methods cannot handle. A pair of functions can also have multiple optimal alignment solutions, which is not addressed in current literature. In this paper, a flexible Bayesian approach to functional alignment is presented, which appropriately accounts for noise in the data without any pre-smoothing required. Additionally, by running parallel MCMC chains, the method can account for multiple optimal alignments via the multi-modal posterior distribution of the warping functions. To most efficiently sample the warping functions, the approach relies on a modification of the standard Hamiltonian Monte Carlo to be well-defined on the infinite-dimensional Hilbert space. In this work, this flexible Bayesian alignment method is applied to both simulated data and real data sets to show its efficiency in handling noisy functions and successfully accounting for multiple optimal alignments in the posterior; characterizing the uncertainty surrounding the warping functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Statistical change detection of building energy consumption: Applications to savings estimation

The surge in interval meter data availability and associated activity in energy data analytics has inspired new interest in advanced methods for building efficiency savings estimation. Statistical and machine learning approaches are being explored to improve the energy baseline models used to measure and verify savings. One outstanding challenge is the ability to identify and account for operational changes that may confound savings estimates. In the measurement and verification (M&V) context, ‘non-routine events’ (NREs) cause changes in building energy use that are not attributable to installed efficiency measures, and not accounted for in the baseline model’s independent variables. In the M&V process NREs must be accounted for as ‘adjustments’ to appropriately attribute the estimated energy savings to the specific efficiency interventions that were implemented. Currently this is a manual and custom process, conducted using professional judgment and engineering expertise. As such it remains a barrier in scaling and standardizing meter-based savings estimation. Here in this work, a data driven methodology was developed to (partially) automate, and therefore streamline the process of detecting NREs in the post-retrofit period and making associated savings adjustments. The proposed NRE detection algorithm is based on a statistical change point detection method and a dissimilarity metric. The dissimilarity metric measures the proximity between the actual time series of the post-retrofit energy consumption and the projected baseline, which is generated using a statistical baseline model. The suggested approach for NRE adjustment involves the NRE detection algorithm, the M&V practitioner, and a regression modeling algorithm. The performance of the detection and adjustment algorithm was evaluated using a simulation-generated test data set, and two benchmark algorithms. Results show a high true positive detection rate (75%-100% across the test cases), higher than ideal false positive detection rates (20%-70%), and low errors in energy adjustment (<0.7%). These results indicate that the algorithm holds for helping M&V practitioners to streamline the process of handling NREs. Moreover, the change point algorithm and underlying statistical principles could prove valuable for other building analytics applications such as anomaly detection and fault diagnostics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

The impact of behavioral and geographic heterogeneity on residential-sector carbon abatement costs

Analyses of monetary and emissions savings from residential efficiency upgrades usually neglect behavioral differences between consumers. Variations in behavior are frequently large: 13% of the U.S. population sets thermostats for cooling at 20 °C (68 °F) or lower and 15% at 25.5 °C (78 °F) or higher. Efficiency analyses should account for behavioral heterogeneity as well as variations in climate and housing characteristics. We model energy and economic savings of efficiency upgrades for U.S. single-family detached houses, accounting for differences in thermostat settings, climate zone, fuel prices, and home characteristics. We consider five efficiency interventions: wall insulation, attic insulation, air sealing, high-efficiency furnaces, and high-efficiency air conditioners. Energy and economic savings vary widely by consumer; the average net present value for wall insulation is $3,020 United States Dollar (USD), but considering heterogeneity, it has a 10th percentile value of -$231 and 90th percentile value of $5,000. Neglecting heterogeneity, technology upgrades can be prioritized from lower to higher carbon abatement costs: air conditioners, wall insulation, furnace, air sealing, and attic insulation. Accounting for heterogeneity reorders this prioritization, and consumer-specific conditions affect mitigation costs such that technologies are mixed in with one another and are no longer arranged one after another. These results indicate that interventions to improve efficiency should consider differences between consumers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2012

Lifestyle Factors in U.S. Residential Electricity Consumption

A multivariate statistical approach to lifestyle analysis of residential electricity consumption is described and illustrated. Factor analysis of selected variables from the 2005 U.S. Residential Energy Consumption Survey (RECS) identified five lifestyle factors reflecting social and behavioral choices associated with air conditioning, laundry usage, personal computer usage, climate zone of residence, and TV use. These factors were also estimated for 2001 RECS data. Multiple regression analysis using the lifestyle factors yields solutions accounting for approximately 40% of the variance in electricity consumption for both years. By adding the associated household and market characteristics of income, local electricity price and access to natural gas, variance accounted for is increased to approximately 54%. Income contributed only {approx}1% unique variance to the 2005 and 2001 models, indicating that lifestyle factors reflecting social and behavioral choices better account for consumption differences than income. This was not surprising given the 4-fold range of energy use at differing income levels. Geographic segmentation of factor scores is illustrated, and shows distinct clusters of consumption and lifestyle factors, particularly in suburban locations. The implications for tailored policy and planning interventions are discussed in relation to lifestyle issues.

03 NATURAL GAS↗
DOE OSTI2017

Evaluation of accessible mineral surface areas for improved prediction of mineral reaction rates in porous media

The rates of mineral dissolution reactions in porous media are difficult to predict, in part because of a lack of understanding of mineral reactive surface area in natural porous media. Common estimates of mineral reactive surface area used in reactive transport models for porous media are typically ad hoc and often based on average grain size, increased to account for surface roughness or decreased by several orders of magnitude to account for reduced surface reactivity of field as opposed to laboratory samples. In this paper, accessible mineral surface areas are determined for a sample from the reservoir formation at the Nagaoka pilot CO 2 injection site (Japan) using a multi-scale image analysis based on synchrotron X-ray microCT, SEM QEMSCAN, XRD, SANS, and FIB-SEM. This analysis not only accounts for accessibility of mineral surfaces to macro-pores, but also accessibility through connected micro-pores in smectite, the most abundant clay mineral in this sample. While the imaging analysis reveals that most of the micro- and macro-pores are well connected, some pore regions are unconnected and thus inaccessible to fluid flow and diffusion. To evaluate whether mineral accessible surface area accurately reflects reactive surface area a flow-through core experiment is performed and modeled at the continuum scale. The core experiment is performed under conditions replicating the pilot site and the evolution of effluent solutes in the aqueous phase is tracked. Various reactive surface area models are evaluated for their ability to capture the observed effluent chemistry, beginning with parameter values determined as a best fit to a disaggregated sediment experiment (Beckingham et al., 2016) described previously. Simulations that assume that all mineral surfaces are accessible (as in the disaggregated sediment experiment) over-predict the observed mineral reaction rates, suggesting that a reduction of RSA by a factor of 10–20 is required to match the core flood experimental data. While the fit of the effluent chemistry (and inferred mineral dissolution rates) greatly improve when the pore-accessible mineral surface areas are used, it was also necessary to include highly reactive glass phases to match the experimental observations, in agreement with conclusions from the disaggregated sediment experiment. Finally, it is hypothesized here that the 10–20 reduction in reactive surface areas based on the limited pore accessibility of reactive phases in core flood experiment may be reasonable for poorly sorted and cemented sediments like those at the Nagaoka site, although this reflects pore rather than larger scale heterogeneity.

58 GEOSCIENCES↗
DOE OSTI2015

A sub-pore model for multi-scale reaction–diffusion problems in porous media

Applications of reaction–diffusion systems in porous media pose a challenging problem for computational modeling approaches due to their multi-physics and multi-scale nature. The length scales usually span 3–4 orders of magnitude while physical phenomena involved include heat and mass transfer processes, and surface reactions. In this paper, a novel methodology that accounts for all the length scales and physical phenomena involved in a single framework is described. A length scale based dual approach is proposed – the larger pore channels (macro-pores) are resolved using conventional numerical techniques and a novel ‘sub-pore’ model is used to account for the unresolved pore channels (sub-pores) and the important physics therein. The porous network in the sub-pore system is composed of a fractal-like hierarchical system of straight cylindrical pores. Simplified governing equations for mass and energy transport are solved within the sub-pore system along with a reaction kinetics model to account for surface adsorption. An implicit coupling strategy is used to couple the macro-pore and the sub-pore systems so as to ensure conservation. The developed methodology is then applied to a few test cases and it is established that the proposed framework is necessary for problems where the adsorption time scale is much smaller than (diffusion-limited) or comparable to the diffusion time scale. Furthermore, it is also demonstrated that the framework can be potentially used to model the network of porous channels in its entirety thus significantly reducing computational costs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A method for predicting failure statistics for steady state elevated temperature structural components

This paper presents the initial development of a high temperature life prediction method that accounts for the variability in the material properties of Grade 91 steel. The method accounts for material variability by fitting a variable 3-parameter Weibull distribution to experimental rupture data and accounts for the variability of creep deformation on the steady-state stresses via a Monte Carlo approach. To ensure reasonable computational times, the model represents the material as an extremely viscous Stokes fluid with a non-Newtonian viscosity, therefore solving the stress relaxation problem with a steady, static, instead of transient, analysis. Furthermore, the complete statistical analysis combines this model for creep deformation with a probabilistic model for creep rupture to evaluate the probability of premature failure for a set of sample problems, comparing the predicted failure statistics to the design life predicted by the ASME Boiler and Pressure Vessel Code rules.

42 ENGINEERING↗
DOE OSTI2020

Determination of permeability for hydrocarbon release due to excavation-induced stress redistribution in rock salt

Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. Specifically, a multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure buildups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.

58 GEOSCIENCES↗
DOE OSTI2011

Patient-Specific Monte Carlo-Based Dose-Kernel Approach for Inverse Planning in Afterloading Brachytherapy

Purpose: Brachytherapy planning software relies on the Task Group report 43 dosimetry formalism. This formalism, based on a water approximation, neglects various heterogeneous materials present during treatment. Various studies have suggested that these heterogeneities should be taken into account to improve the treatment quality. The present study sought to demonstrate the feasibility of incorporating Monte Carlo (MC) dosimetry within an inverse planning algorithm to improve the dose conformity and increase the treatment quality. Methods and Materials: The method was based on precalculated dose kernels in full patient geometries, representing the dose distribution of a brachytherapy source at a single dwell position using MC simulations and the Geant4 toolkit. These dose kernels are used by the inverse planning by simulated annealing tool to produce a fast MC-based plan. A test was performed for an interstitial brachytherapy breast treatment using two different high-dose-rate brachytherapy sources: the microSelectron iridium-192 source and the electronic brachytherapy source Axxent operating at 50 kVp. Results: A research version of the inverse planning by simulated annealing algorithm was combined with MC to provide a method to fully account for the heterogeneities in dose optimization, using the MC method. The effect of the water approximation was found to depend on photon energy, with greater dose attenuation for the lower energies of the Axxent source compared with iridium-192. For the latter, an underdosage of 5.1% for the dose received by 90% of the clinical target volume was found. Conclusion: A new method to optimize afterloading brachytherapy plans that uses MC dosimetric information was developed. Including computed tomography-based information in MC dosimetry in the inverse planning process was shown to take into account the full range of scatter and heterogeneity conditions. This led to significant dose differences compared with the Task Group report 43 approach for the Axxent source.

62 RADIOLOGY AND NUCLEAR MEDICINE↗