Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model based development”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

MODELING AND METHOD FOR OPTICAL PERFORMANCE OF A RECEIVER COATING AND ITS DEGRADATION DUE TO OPERATION

The optical properties of the coating used on a concentrated solar power (CSP)’s receiver directly impact the performance of the power plant. Receiver coatings are designed to absorb as much heat as possible and transfer it to the circulating fluid. Over time, the absorptivity of the receiver coating degrades due to environmental and operational conditions. This project aims to (1) create a test apparatus that allows coatings to undergo accelerated lifetime testing under conditions that mimic what the coating will experience on-tower; and (2) develop a mathematical model that based on the test data predicts absorptivity over time. To create the mathematical model, samples of BrightSource Energy’s Gen1 and Gen2 coatings were tested first looking at individual failure modes (such as high temperature cycling) and then for combined failure modes in the newly developed apparatus. In addition to taking optical measurements, the samples were destroyed to allow for metalogical evaluation to understand what was happening on a microscopic level. Pre-existing models from the literature that describe accelerated testing as well as individual phenomena observed in the coating were then combined to create the model. Additionally, the Gen2 and Gen3 coatings are being testing at the Plataforma Solar de Almería (PSA) facility in Almeria, Spain to confirm that the apparatus test results are similar to those received on-sun. (PSA’s facility closely mimics on-tower conditions). The goal is to enable industry players to test coatings using the new apparatus, gather absorptivity measurements (without needing to destroy samples), and input the information into the model to determine how the coating will perform over time, thus allowing the optical performance of coatings to be compared even before they are used on-tower. The model will be publicly available by the end of the project (currently projected as December 2024).

14 SOLAR ENERGY↗
DOE OSTI2018

SCC Aerosol Transport Model Summary Report

The purpose of the present work is to introduce a generic and reliable numerical model for prediction of aerosol transport, deposition, and plugging in leak paths similar to stress corrosion cracks while accounting for potential plugging formation. The model is dynamic (changing leak path geometry due to plugging) and it relies on the numerical solution of the aerosol transport equation in one dimension using finite differences. An extensive validation exercise (particle diameters: 0.01-10 μm and pressure difference up to 12 kPa) was conducted based on comparisons with experimental and theoretical data. The developed model is fairly general since it is based on a generic mechanistic description of the aerosol flow and particle deposition in the leak path, without requiring experimental fittings other than using the leak rate parameters commonly available in practical applications. Despite the approximate agreement of the model with the experimental data, the obtained predictions can be qualified as realistic. This is a major achievement of the model, which has been based on generic aerosol mechanics with no reference or adjustment to experimental data.

42 ENGINEERING↗
DOE OSTI2023

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Modeling thermochemistry of fuel and coupling to fuel performance codes

A thermochemical library, Thermochimica, utilizes new algorithms and solvers designed for calculating equilibria of multicomponent and multiphase systems and has been coupled to the nuclear fuel performance code BISON. Thermochimica utilizes a generalized thermochemical database for urania-based fuels that has been under development for a number of years and continues to be expanded. One of the key motivations for use of thermodynamic models is that they can better represent the physics of oxygen migration. The coupling of the BISON code with Thermochimica and examples of its use in representing fuel phenomena are described. In the current implementation in BISON, fluxes are driven by the elemental concentration gradient. However, equilibrium is characterized by the absence of spatial variations of chemical potentials of the system species/components. The driving force to equilibrium (species flux) is proportional to the deviation from equilibrium, i.e., the gradient in chemical potential. It is therefore necessary to integrate chemical potential driven diffusion in the BISON representation of oxygen transport. The diffusion kernels based on thermodynamic models have been developed to simulate multicomponent transport in light water reactor (LWR) materials. The Thermochimica library was used to calculate various thermodynamics properties needed for transport calculations. The library also can calculate properties that are specific to various transport mechanisms for different fuel materials. For example, it can calculate defect site fractions among the sublattices in oxide fuels as a function of temperature, burnup, and stoichiometry. The transport and thermodynamic models are being integrated with other simulation efforts such as calculation of mobility functions that couple fluxes and chemical gradients. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Benchmarks for Model Validation based on LiDAR Wake Measurements

Technology development and design decisions in wind energy are often based on results from simulations performed for individual wind turbines or entire wind plants. It is therefore critical to ensure that the models being used for research and industry applications in wind energy be thoroughly validated against measurements. A full-system validation of wind plant simulations must consider the atmospheric inflow, the response of the wind turbines, and their wakes. This task is complicated by the lack of freely available, quality-controlled, high-quality measurements. Here, such measurements are used to offer a validation exercise that can be used to assess the accuracy of models of any fidelity level. When it comes to real-world measurements, the dataset considered herein is simple in terms of terrain but exhibits pronounced diurnal cycles. Instead of a full-scale wind plant, we consider an individual research-scale, utility wind turbine instrumented for power and loads measurements. Three benchmarks are defined, with increasing levels of complexity: near neutral, slightly unstable, and very stable atmospheric stratification. Through comparisons between observations and simulations, the benchmarks provide complementary information about the model performance and its ability to reproduce mean and dynamic wake characteristics. This article describes the measurements and methodology used to define these benchmarks and provides the information required to perform simulations and conduct the model-measurement comparison. The objective is to provide a robust wake model validation exercise open to anyone, which will serve to minimize uncertainty in model validation practices related to varying methodologies across simulation tools and users.

17 WIND ENERGY↗
DOE OSTI2012

Toward a mechanistic modeling of nitrogen limitation on vegetation dynamics

Nitrogen is a dominant regulator of vegetation dynamics, net primary production, and terrestrial carbon cycles; however, most ecosystem models use a rather simplistic relationship between leaf nitrogen content and photosynthetic capacity. Such an approach does not consider how patterns of nitrogen allocation may change with differences in light intensity, growing-season temperature and CO{sub 2} concentration. To account for this known variability in nitrogen-photosynthesis relationships, we develop a mechanistic nitrogen allocation model based on a trade-off of nitrogen allocated between growth and storage, and an optimization of nitrogen allocated among light capture, electron transport, carboxylation, and respiration. The developed model is able to predict the acclimation of photosynthetic capacity to changes in CO{sub 2} concentration, temperature, and radiation when evaluated against published data of V{sub c,max} (maximum carboxylation rate) and J{sub max} (maximum electron transport rate). A sensitivity analysis of the model for herbaceous plants, deciduous and evergreen trees implies that elevated CO{sub 2} concentrations lead to lower allocation of nitrogen to carboxylation but higher allocation to storage. Higher growing-season temperatures cause lower allocation of nitrogen to carboxylation, due to higher nitrogen requirements for light capture pigments and for storage. Lower levels of radiation have a much stronger effect on allocation of nitrogen to carboxylation for herbaceous plants than for trees, resulting from higher nitrogen requirements for light capture for herbaceous plants. As far as we know, this is the first model of complete nitrogen allocation that simultaneously considers nitrogen allocation to light capture, electron transport, carboxylation, respiration and storage, and the responses of each to altered environmental conditions. We expect this model could potentially improve our confidence in simulations of carbon-nitrogen interactions and the vegetation feedbacks to climate in Earth system models.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Determination of equilibrium electron temperature and times using an electron swarm model with BOLSIG+ calculated collision frequencies and rate coefficients

Electromagnetic pulse (EMP) events produce low-energy conduction electrons from Compton electron or photoelectron ionizations with air. It is important to understand how conduction electrons interact with air in order to accurately predict EMP evolution and propagation. An electron swarm model can be used to monitor the time evolution of conduction electrons in an environment characterized by electric field and pressure. Here a swarm model is developed that is based on the coupled ordinary differential equations (ODEs) described by Higgins et al. (1973), hereinafter HLO. The ODEs characterize the swarm electric field, electron temperature, electron number density, and drift velocity. Important swarm parameters, the momentum transfer collision frequency, energy transfer collision frequency, and ionization rate, are calculated and compared to the previously reported fitted functions given in HLO. These swarm parameters are found using BOLSIG+, a two term Boltzmann solver developed by Hagelaar and Pitchford (2005), which utilizes updated cross sections from the LXcat website created by Pancheshnyi et al. (2012). We validate the swarm model by comparing to experimental effective ionization coefficient data in Dutton (1975) and drift velocity data in Ruiz-Vargas et al. (2010). In addition, we report on electron equilibrium temperatures and times for a uniform electric field of 1 StatV/cm for atmospheric heights from 0 to 40 km. We show that the equilibrium temperature and time are sensitive to the modifications in the collision frequencies and ionization rate based on the updated electron interaction cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A combined experimental and computational analysis of failure mechanisms in open-hole cross-ply laminates under flexural loading

In this work, integrated experimental tests and computational modeling are proposed to investigate the failure mechanisms of open-hole cross-ply carbon fiber reinforced polymer (CFRP) laminated composites. In particular, we propose two effective methods, which include width-tapered double cantilever beam (WTDCB) and fixed-ratio mixed-mode end load split (FRMMELS) tests, to obtain the experimental data more reliably. We then calibrate the traction-separation laws of cohesive zone model (CZM) used among laminas of the composites by leveraging these two methods. The experimental results of fracture energy, i.e. G Ic and G Tc , obtained from WTDCB and FRMMELS tests are generally insensitive to the crack length thus requiring no effort to accurately measure the crack tip. Moreover, FRMMELS sample contains a fixed mixed-mode ratio of G IIc /G Tc depending on the width taper ratio. Examining comparisons between experimental results of FRMMELS tests and failure surface of B–K failure criterion predicted from a curve fitting, good agreement between the predictions and experimental data has been found, indicating that FRMMELS tests are an effective method to determine mixed-mode fracture criterion. In addition, a coupled experimental-computational modeling of WTDCB, edge notched flexure, and FRMMELS tests are adopted to calibrate and validate the interfacial strengths. Finally, failure mechanisms of open-hole cross-ply CFRP laminates under flexural loading have been studied systematically using experimental and multi-scale computational analyses based on the developed CZM model. The initiation and propagation of delamination, the failure of laminated layers as well as load-displacement curves predicted from computational analyses are in good agreement with what we have observed experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling suggests that virion production cycles within individual cells is key to understanding acute hepatitis B virus infection kinetics

Hepatitis B virus (HBV) infection kinetics in immunodeficient mice reconstituted with humanized livers from inoculation to steady state is highly dynamic despite the absence of an adaptive immune response. To recapitulate the multiphasic viral kinetic patterns, we developed an agent-based model that includes intracellular virion production cycles reflecting the cyclic nature of each individual virus lifecycle. The model fits the data well predicting an increase in production cycles initially starting with a long production cycle of 1 virion per 20 hours that gradually reaches 1 virion per hour after approximately 3–4 days before virion production increases dramatically to reach to a steady state rate of 4 virions per hour per cell. Together, modeling suggests that it is the cyclic nature of the virus lifecycle combined with an initial slow but increasing rate of HBV production from each cell that plays a role in generating the observed multiphasic HBV kinetic patterns in humanized mice.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Pharmacokinetic modeling: Prediction and evaluation of route dependent dosimetry of bisphenol A in monkeys with extrapolation to humans

A physiologically based pharmacokinetic (PBPK) model was developed for bisphenol A (BPA) in adult rhesus monkeys using intravenous (iv) and oral bolus doses of 100 {mu}g d6-BPA/kg (). This calibrated PBPK adult monkey model for BPA was then evaluated against published monkey kinetic studies with BPA. Using two versions of the adult monkey model based on monkey BPA kinetic data from and , the aglycone BPA pharmacokinetics were simulated for human oral ingestion of 5 mg d16-BPA per person (Voelkel et al., 2002). Voelkel et al. were unable to detect the aglycone BPA in plasma, but were able to detect BPA metabolites. These human model predictions of the aglycone BPA in plasma were then compared to previously published PBPK model predictions obtained by simulating the Voelkel et al. kinetic study. Our BPA human model, using two parameter sets reflecting two adult monkey studies, both predicted lower aglycone levels in human serum than the previous human BPA PBPK model predictions. BPA was metabolized at all ages of monkey (PND 5 to adult) by the gut wall and liver. However, the hepatic metabolism of BPA and systemic clearance of its phase II metabolites appear to be slower in younger monkeys than adults. The use of the current non-human primate BPA model parameters provides more confidence in predicting the aglycone BPA in serum levels in humans after oral ingestion of BPA. -- Highlights: Black-Right-Pointing-Pointer A bisphenol A (BPA) PBPK model for the infant and adult monkey was constructed. Black-Right-Pointing-Pointer The hepatic metabolic rate of BPA increased with age of the monkey. Black-Right-Pointing-Pointer The systemic clearance rate of metabolites increased with age of the monkey. Black-Right-Pointing-Pointer Gut wall metabolism of orally administered BPA was substantial across all ages of monkeys. Black-Right-Pointing-Pointer Aglycone BPA plasma concentrations were predicted in humans orally given oral doses of deuterated BPA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Advanced Modeling of Beam Physics and Performance Optimization for Nuclear Physics Colliders

High energy colliders provide a critical tool in nuclear physics study by probing the fundamental structure and dynamics of matter. To maximize the potential of scientific discovery in nuclear physics study, it is important to optimize the parameters of these colliders to attain the best performance. The performance of a collider is typically measured by its integrated luminosity of colliding beams since the probability of a new event is proportional to the integrated luminosity. However, the achievable luminosity is limited by the electromagnetic interactions (beam-beam effects) of two colliding beams at higher energy, and the interplay between the space-charge effects and the beam-beam effects at lower energy. To achieve the best performance of a collider means to attain the highest luminosity of the collider with optimized collider parameters. Optimizing the collider’s machine parameters is both computationally and experimentally expensive. A fast and robust computational framework including beam-beam and space-charge effects will be critical to attaining the best performance of the collider. In this project, we will study the beam dynamics challenges, specifically the interplay of the space-charge and the beam-beam effects, and the machine tuning models for maximizing the performance of RHIC experiments. We will develop an advanced modeling framework based on first-principles physical simulations, lattice models and the state-of-the-art machine learning methods and apply this framework to performance improvement of the RHIC in operation. We will build data manipulation packages to connect the simulation data and the experimental data with the framework, develop a self-consistent hybrid model of space-charge and beam-beam effects, study underlying physics mechanisms, build surrogate models using the labeled data, integrate the models into the advanced modeling framework, and apply the framework to RHIC luminosity (STAR and sPHENIX) optimization. The success of this project would substantially improve the performance of existing and future colliders and increase the opportunity for scientific discovery.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Multiscale Modeling of Accident Tolerant Fuels under the NEAMS ATF Program

The U.S. Department of Energy's Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has been developing computational analysis tools for nuclear fuel simulation. These codes include BISON and Marmot for investigation of fuel behavior at the continuum and microstructural scales, respectively. In 2014, NEAMS introduced what it calls High Impact Problems (HIPs) into its development plan. These HIPs are intended to make significant advancements in a particular area of nuclear power research in a short period of time (3 years or less). NEAMS has chosen an accident tolerant fuel (ATF) project, which emphasizes utilizing BISON and Marmot to model proposed materials, as its first HIP. Since the events that occurred at the Fukushima-Daiichi nuclear power plant in 2011, the U.S. Department of Energy's Office of Nuclear Energy accelerated research in nuclear fuel concepts with enhanced accident tolerance as part of its Fuel Cycle Research and Development (FCRD) Advanced Fuels Campaign (AFC). Accident tolerant fuels are expected to perform at least as well as conventional fuel during normal conditions and provide significantly increased response time in the event of an accident. Given the many unknowns associated with proposed ATF concepts, the research and development being undertaken by AFC in this area is a challenge, particularly considering the desire to have a lead test rod or assembly in a commercial plant by 2022. With the AFC as its primary stakeholder, the NEAMS ATF HIP is developing material models (using a multi-scale approach) for multiple candidate materials. Through a variety of analyses employing these models, the ATF HIP is providing insight and support to the AFC as it strives to understand these materials and prioritize further development. Use of computational tools such as BISON and Marmot are a vital component of the AFC's strategy since experimental data will not be available broadly and quickly enough to understand the proposed materials fully by the time the AFC must make further development decisions. This summary reviews the overall approach of the NEAMS ATF HIP. First, we give an overview of the current materials chosen for investigation. We also review the overall multi-scale approach to model development. We briefly review the work underway at multiple national laboratories to develop mechanistic engineering scale models based on lower length scale insight. Note that coordinated work is also underway at multiple universities, but those efforts are not detailed here. The NEAMS ATF HIP is assisting the AFC in assessing two accident tolerant materials, uranium silicide fuel and FeCrAl cladding. A multi-scale material model development approach is being undertaken to formulate the models necessary to investigate the materials over a range of normal and off-normal scenarios. This approach includes molecular dynamics, density functional theory, mesoscale modeling, and engineering scale modeling. Sensitivity studies are underway to understand the most critical scenarios and material model input parameters. The information generated is being given to the AFC for their use in selecting promising ATF concepts. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗
DOE OSTI2018

Development of Application Function Blocks for Power-Hardware-in-the-Loop Testing of Grid-Connected Inverters: Preprint

Stability and accuracy of the PHIL platform has been a common problem and challenging issue for a long time, thus continuous efforts and research have been devoted to developing interface algorithms to improve upon this. Some of existing works exhibit promising solutions to compensate the errors in the PHIL platform. However, there is a need to harmonize various solutions into a standard framework that provides a general solution for various PHIL test beds and that can be reused with minimum modifications to avoid extra effort. A good way to do this is by using application function blocks (AFBs), which are model-based and can be developed as a standard library component, thus achieving satisfactory reusability and reconfigurability for various platforms. Since PHIL simulation has the flexibility of the software simulation, proper functions can be easily implemented in the digital real-time simulator (DRTS) to preprocess the reference signal before sending it to the power amplifier. Therefore, this paper aims to develop standard AFBs in DRTS for the PHIL interface to test grid-connected inverters. The PHIL interface algorithm is separated into five AFBs that will be presented in detail in this paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

A Causal Inference Model Based on Random Forests to Identify the Effect of Soil Moisture on Precipitation

Soil moisture influences precipitation mainly through its impact on land–atmosphere interactions. Understanding and correctly modeling soil moisture–precipitation (SM–P) coupling is crucial for improving weather forecasting and subseasonal to seasonal climate predictions, especially when predicting the persistence and magnitude of drought. However, the sign and spatial structure of SM–P feedback are still being debated in the climate research community, mainly due to the difficulty in establishing causal relationships and the high degree of nonlinearity in land–atmosphere processes. To this end, we developed a causal inference model based on the Granger causality analysis and a nonlinear machine learning model. This model includes three steps: nonlinear anomaly decomposition, nonlinear Granger causality analysis, and evaluation of the quality of SM–P feedback, which eliminates the nonlinear response of interannual and seasonal variability and the memory effects of climatic factors and isolates the causal relationship of local SM–P feedback. We applied this model by using National Climate Assessment–Land Data Assimilation System (NCA-LDAS) datasets over the United States. Here, the results highlight the importance of nonlinear atmosphere responses in land–atmosphere interactions. In addition, the strong feedback over the southwestern United States and the Great Plains both highlight the impacts of topographic factors rather than only the sensitivity of evapotranspiration to soil moisture. Furthermore, the SM–P index defined by our framework is used to benchmark Earth system models (ESMs), which provides a new metric for efficiently identifying potential model biases in modeling local land–atmosphere interactions and may help the development of ESMs in improving simulations of water cycle variability and extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling the effects of carbon payments and forest owner cooperatives on carbon storage and revenue in Pacific Northwest forestlands

The forestlands of the Pacific Northwest store a significant amount of carbon, provide wildlife habitat, generate revenue for forestland owners, and provide for sustainable generation of forest products. Promoting forest carbon storage is increasingly recognized as a key management objective to mitigate climate change. We developed an agent-based model to examine the effectiveness of current land management practices in sustaining forest carbon storage and revenue flows, and the utility of market-based carbon payments in improving carbon storage while also improving revenue generation in the Pacific Northwest. We find that while market based incentives can improve both forest carbon storage and revenue, they are limited in their benefit to small forest owners. Our results highlight the crucial need for additional policy interventions to reduce barriers to market access for small private owners and ensure sustainable forest management in this region. We further show that a forest owners’ cooperative, which aggregates small owners to reduce the costs of participating in a carbon market, can help overcome these barriers and result in both greater carbon storage and a more equitable benefit from carbon payments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

SU-E-T-578: MCEBRT, A Monte Carlo Code for External Beam Treatment Plan Verifications

Purpose: Present a new Monte Carlo code (MCEBRT) for patient-specific dose calculations in external beam radiotherapy. The code MLC model is benchmarked and real patient plans are re-calculated using MCEBRT and compared with commercial TPS. Methods: MCEBRT is based on the GEPTS system (Med. Phys. 29 (2002) 835–846). Phase space data generated for Varian linac photon beams (6 – 15 MV) are used as source term. MCEBRT uses a realistic MLC model (tongue and groove, rounded ends). Patient CT and DICOM RT files are used to generate a 3D patient phantom and simulate the treatment configuration (gantry, collimator and couch angles; jaw positions; MLC sequences; MUs). MCEBRT dose distributions and DVHs are compared with those from TPS in absolute way (Gy). Results: Calculations based on the developed MLC model closely matches transmission measurements (pin-point ionization chamber at selected positions and film for lateral dose profile). See Fig.1. Dose calculations for two clinical cases (whole brain irradiation with opposed beams and lung case with eight fields) are carried out and outcomes are compared with the Eclipse AAA algorithm. Good agreement is observed for the brain case (Figs 2-3) except at the surface where MCEBRT dose can be higher by 20%. This is due to better modeling of electron contamination by MCEBRT. For the lung case an overall good agreement (91% gamma index passing rate with 3%/3mm DTA criterion) is observed (Fig.4) but dose in lung can be over-estimated by up to 10% by AAA (Fig.5). CTV and PTV DVHs from TPS and MCEBRT are nevertheless close (Fig.6). Conclusion: A new Monte Carlo code is developed for plan verification. Contrary to phantombased QA measurements, MCEBRT simulate the exact patient geometry and tissue composition. MCEBRT can be used as extra verification layer for plans where surface dose and tissue heterogeneity are an issue.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-EF-207-03: Design and Optimization of a CBCT Head Scanner for Detection of Acute Intracranial Hemorrhage

Purpose: To design a dedicated x-ray cone-beam CT (CBCT) system suitable to deployment at the point-of-care and offering reliable detection of acute intracranial hemorrhage (ICH), traumatic brain injury (TBI), stroke, and other head and neck injuries. Methods: A comprehensive task-based image quality model was developed to guide system design and optimization of a prototype head scanner suitable to imaging of acute TBI and ICH. Previously reported models were expanded to include the effects of x-ray scatter correction necessary for detection of low contrast ICH and the contribution of bit depth (digitization noise) to imaging performance. Task-based detectablity index provided the objective function for optimization of system geometry, x-ray source, detector type, anti-scatter grid, and technique at 10–25 mGy dose. Optimal characteristics were experimentally validated using a custom head phantom with 50 HU contrast ICH inserts imaged on a CBCT imaging bench allowing variation of system geometry, focal spot size, detector, grid selection, and x-ray technique. Results: The model guided selection of system geometry with a nominal source-detector distance 1100 mm and optimal magnification of 1.50. Focal spot size ∼0.6 mm was sufficient for spatial resolution requirements in ICH detection. Imaging at 90 kVp yielded the best tradeoff between noise and contrast. The model provided quantitation of tradeoffs between flat-panel and CMOS detectors with respect to electronic noise, field of view, and readout speed required for imaging of ICH. An anti-scatter grid was shown to provide modest benefit in conjunction with post-acquisition scatter correction. Images of the head phantom demonstrate visualization of millimeter-scale simulated ICH. Conclusions: Performance consistent with acute TBI and ICH detection is feasible with model-based system design and robust artifact correction in a dedicated head CBCT system. Further improvements can be achieved with incorporation of model-based iterative reconstruction techniques also within the scope of the task-based optimization framework. David Foos and Xiaohui Wang are employees of Carestream Health.

60 APPLIED LIFE SCIENCES↗