Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model fitting”

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

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

At least 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-T-598: Parametric Equation for Quick and Reliable Estimate of Stray Neutron Doses in Proton Therapy and Application for Intracranial Tumor Treatments

Purpose: Development of a parametric equation suitable for a daily use in routine clinic to provide estimates of stray neutron doses in proton therapy. Methods: Monte Carlo (MC) calculations using the UF-NCI 1-year-old phantom were exercised to determine the variation of stray neutron doses as a function of irradiation parameters while performing intracranial treatments. This was done by individually changing the proton beam energy, modulation width, collimator aperture and thickness, compensator thickness and the air gap size while their impact on neutron doses were put into a single equation. The variation of neutron doses with distance from the target volume was also included in it. Then, a first step consisted in establishing the fitting coefficients by using 221 learning data which were neutron absorbed doses obtained with MC simulations while a second step consisted in validating the final equation. Results: The variation of stray neutron doses with irradiation parameters were fitted with linear, polynomial, etc. model while a power-law model was used to fit the variation of stray neutron doses with the distance from the target volume. The parametric equation fitted well MC simulations while establishing fitting coefficients as the discrepancies on the estimate of neutron absorbed doses were within 10%. The discrepancy can reach ∼25% for the bladder, the farthest organ from the target volume. Finally, the validation showed results in compliance with MC calculations since the discrepancies were also within 10% for head-and-neck and thoracic organs while they can reach ∼25%, again for pelvic organs. Conclusion: The parametric equation presents promising results and will be validated for other target sites as well as other facilities to go towards a universal method.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Harderian Gland Tumorigenesis: Low-Dose and LET Response

Increased cancer risk remains a primary concern for travel into deep space and may preclude manned missions to Mars due to large uncertainties that currently exist in estimating cancer risk from the spectrum of radiations found in space with the very limited available human epidemiological radiation-induced cancer data. Existing data on human risk of cancer from X-ray and gamma-ray exposure must be scaled to the many types and fluences of radiations found in space using radiation quality factors and dose-rate modification factors, and assuming linearity of response since the shapes of the dose responses at low doses below 100 mSv are unknown. The goal of this work was to reduce uncertainties in the relative biological effect (RBE) and linear energy transfer (LET) relationship for space-relevant doses of charged-particle radiation-induced carcinogenesis. The historical data from the studies of Fry et al. and Alpen et al. for Harderian gland (HG) tumors in the female CB6F1 strain of mouse represent the most complete set of experimental observations, including dose dependence, available on a specific radiation-induced tumor in an experimental animal using heavy ion beams that are found in the cosmic radiation spectrum. However, these data lack complete information on low-dose responses below 0.1 Gy, and for chronic low-dose-rate exposures, and there are gaps in the LET region between 25 and 190 keV/μm. In this study, we used the historical HG tumorigenesis data as reference, and obtained HG tumor data for 260 MeV/u silicon (LET ~70 keV/μm) and 1,000 MeV/u titanium (LET ~100 keV/μm) to fill existing gaps of data in this LET range to improve our understanding of the dose-response curve at low doses, to test for deviations from linearity and to provide RBE estimates. Animals were also exposed to five daily fractions of 0.026 or 0.052 Gy of 1,000 MeV/u titanium ions to simulate chronic exposure, and HG tumorigenesis from this fractionated study were compared to the results from single 0.13 or 0.26 Gy acute titanium exposures. Theoretical modeling of the data show that a nontargeted effect model provides a better fit than the targeted effect model, providing important information at space-relevant doses of heavy ions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2018

The Bright and the Slow—GRBs 100724B and 160509A with High-energy Cutoffs at ≲100 MeV

Here, we analyze the prompt emission of GRB 100724B and GRB 160509A, two of the brightest gamma-ray bursts (GRBs) observed by Fermi at ≲ MeV energies but surprisingly faint at ≳100 MeV energies. Time-resolved spectroscopy reveals a sharp high-energy cutoff at energies E c ~ 20–60 MeV for GRB 100724B and E c ~ 80–150 MeV for GRB 160509A. We first characterize phenomenologically the cutoff and its time evolution. We then fit the data to two models where the high-energy cutoff arises from intrinsic opacity to pair production within the source (τ γγ ): (i) a Band spectrum with τ γγ from the internal-shocks-motivated model of Granot et al. (2008) and (ii) the photospheric model of Gill & Thompson (2014). Alternative explanations for the cutoff, such as an intrinsic cutoff in the emitting electron energy distribution, appear to be less natural. Both models provide a good fit to the data with very reasonable physical parameters, providing an estimate of bulk Lorentz factors in the range Γ ~ 100–400, on the lower end of what is generally observed in Fermi GRBs. Surprisingly, their lower cutoff energies E c compared to other Fermi/LAT GRBs arise not only predominantly from the lower Lorentz factors, but also at a comparable level from differences in variability time, luminosity, and high-energy photon index. Lastly, particularly low E c values may prevent detection by Fermi/LAT, thus introducing a bias in the Fermi/LAT GRB sample against GRBs with low Lorentz factors or variability times.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Unveiling acoustic physics of the CMB using nonparametric estimation of the temperature angular power spectrum for Planck

Estimation of the angular power spectrum is one of the important steps in Cosmic Microwave Background (CMB) data analysis. Here, we present a nonparametric estimate of the temperature angular power spectrum for the Planck 2013 CMB data. The method implemented in this work is model-independent, and allows the data, rather than the model, to dictate the fit. Since one of the main targets of our analysis is to test the consistency of the ΛCDM model with Planck 2013 data, we use the nuisance parameters associated with the best-fit ΛCDM angular power spectrum to remove foreground contributions from the data at multipoles ℓ ≥50. We thus obtain a combined angular power spectrum data set together with the full covariance matrix, appropriately weighted over frequency channels. Our subsequent nonparametric analysis resolves six peaks (and five dips) up to ℓ ∼1850 in the temperature angular power spectrum. We present uncertainties in the peak/dip locations and heights at the 95% confidence level. We further show how these reflect the harmonicity of acoustic peaks, and can be used for acoustic scale estimation. Based on this nonparametric formalism, we found the best-fit ΛCDM model to be at 36% confidence distance from the center of the nonparametric confidence set—this is considerably larger than the confidence distance (9%) derived earlier from a similar analysis of the WMAP 7-year data. Another interesting result of our analysis is that at low multipoles, the Planck data do not suggest any upturn, contrary to the expectation based on the integrated Sachs-Wolfe contribution in the best-fit ΛCDM cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Environmental controls on the light use efficiency of terrestrial gross primary production

Abstract Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf‐level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy‐level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C 3 species; however, additional assumptions are required to “scale up” from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed‐effects model. The same model design was fitted to site‐based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple “goodness of fit” comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

An improved method to increase the predictive accuracy of the ECR technique

The technique of electrical conductivity relaxation (ECR) is a widely utilized method for obtaining oxygen chemical surface exchange (k) and bulk diffusion (D) coefficients for mixed ionic electronic conductors. Although the technique is conceptually simple, its practical application does not frequently produce reliably accurate values of k and D, and many examples of poor agreement can be found in the published literature. In this present paper, methods for enhancing the accuracy of the ECR technique relative to data processing are described in detail. Physically accurate determination of the actual starting time (t{sub 0}) and improved methods for data fitting of ECR measurements are discussed. To test the methods described, generated relaxation data sets are analyzed by using non-linear-least-square (NLLS) methods. The results showed that fixing t{sub 0} based on an ideally stirred reactor model improves the fitted results compared to treating t{sub 0} as a fitted parameter. Additionally, evaluation of the global minimum in the fitted error is made by analyzing multiple samples of varying thickness, which improves parameter evaluation compared to analysis of a single relaxation data set. The improved accuracy of the computed values of k and D is described in comparison to conventional methods.

Bulk diffusion coefficient↗
DOE OSTI2017

Pati-Salam SUSY GUT with Yukawa unification

Prior studies of a three family Yukawa unified supersymmetric grand unified theory (SUSY GUT) with SO(10) or Pati-Salam (PS) gauge symmetry proposed by Raby and students show that this model is able to fit low energy and inflation observables. Yet, the fit to low energy observables is not great especially for $\sin2\beta$, and up and down quark masses. In this paper, we show that by choosing PS as the gauge group and modifying the Yukawa sector, the low energy fit improves significantly while other qualities of the model are maintained. In particular, the lightest SUSY particle is the neutralino with mass of order 300 - 500 GeV, the lightest stop and sbottom have mass of order 3 - 5 TeV and the CP odd Higgs mass is of order 5 - 6 TeV, so we are in the decoupling limit for the light Higgs. In addition, we reinterpret gluino simplified model analyses by the ATLAS and CMS collaborations and find that the most stringent gluino mass bound for our model is $M_{\tilde{g}}\sim1.9\,\text{TeV}$. The current best fit point, consistent with this bound, with gluino mass $M_{\tilde{g}}=1.9\,\text{TeV}$ has $\chi^2/\text{dof}\approx1.12$, compared to the best fit point of the previous model with $\chi^2/\text{dof}=1.90$. We find that predictions for the electric dipole moment of the electron, the branching ratio $BR(\mu \rightarrow e \gamma)$ and the CP violating angle in the lepton sector, $\sin\delta$, are affected significantly as compared to previous results. In conclusion, we are unable to rule out this model or place an upper bound on gluino mass accessible by this run of the LHC because the $\chi^2/\text{dof}$ of this model is well below $2\sigma$ even for a gluino mass as high as 2.7 TeV. On the bright side, this means that this model is still viable and we might find low energy SUSY particles in the near future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

THE PHYSICAL CHARACTERIZATION OF THE POTENTIALLY HAZARDOUS ASTEROID 2004 BL86: A FRAGMENT OF A DIFFERENTIATED ASTEROID

The physical characterization of potentially hazardous asteroids (PHAs) is important for impact hazard assessment and evaluating mitigation options. Close flybys of PHAs provide an opportunity to study their surface photometric and spectral properties that enable the identification of their source regions in the main asteroid belt. We observed PHA (357439) 2004 BL86 during a close flyby of the Earth at a distance of 1.2 million km (0.0080 AU) on 2015 January 26, with an array of ground-based telescopes to constrain its photometric and spectral properties. Lightcurve observations showed that the asteroid was a binary and subsequent radar observations confirmed the binary nature and gave a primary diameter of 300 m and a secondary diameter of 50–100 m. Our photometric observations were used to derive the phase curve of 2004 BL86 in the V-band. Two different photometric functions were fitted to this phase curve, the IAU H–G model and the Shevchenko model. From the fit of the H–G function we obtained an absolute magnitude of H = 19.51 ± 0.02 and a slope parameter of G = 0.34 ± 0.02. The Shevchenko function yielded an absolute magnitude of H = 19.03 ± 0.07 and a phase coefficient b = 0.0225 ± 0.0006. The phase coefficient was used to calculate the geometric albedo (Ag) using the relationship found by Belskaya and Schevchenko, obtaining a value of Ag = 40% ± 8% in the V-band. With the geometric albedo and the absolute magnitudes derived from the H–G and the Shevchenko functions we calculated the diameter (D) of 2004 BL86, obtaining D = 263 ± 26 and D = 328 ± 35 m, respectively. 2004 BL86 spectral band parameters and pyroxene chemistry are consistent with non-cumulate eucrite meteorites. A majority of these meteorites are derived from Vesta and are analogous with surface lava flows on a differentiated parent body. A non-diagnostic spectral curve match using the Modeling for Asteroids tool yielded a best-match with non-cumulate eucrite Bereba. Three other near-Earth asteroids (1993 VW, 1998 KK17, and 2000 XH44) that were observed by Burbine et al. also have spectral properties similar to 2004 BL86. The presence of eucrites with anomalous oxygen isotope ratios compared to the howardites, eucrites, and diogenites meteorites from Vesta suggests the possible presence of multiple differentiated bodies in the inner main belt or the contamination of Vesta’s surface with exogenic material. The spectral properties of both anomalous and Vestan eucrites are degenerate, making it difficult to identify the parent bodies of anomalous eucrites in the main belt and the NEO population using remote sensing. This makes it difficult to link 2004 BL86 directly to Vesta, although the Vesta family is the largest contributor of V-types to near-Earth space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A generalized linear-quadratic model incorporating reciprocal time pattern of radiation damage repair

Purpose: It has been conventionally assumed that the repair rate for sublethal damage (SLD) remains constant during the entire radiation course. However, increasing evidence from animal studies suggest that this may not the case. Rather, it appears that the repair rate for radiation-induced SLD slows down with increasing time. Such a slowdown in repair would suggest that the exponential repair pattern would not necessarily accurately predict repair process. As a result, the purpose of this study was to investigate a new generalized linear-quadratic (LQ) model incorporating a repair pattern with reciprocal time. The new formulas were tested with published experimental data. Methods: The LQ model has been widely used in radiation therapy, and the parameter G in the surviving fraction represents the repair process of sublethal damage with T{sub r} as the repair half-time. When a reciprocal pattern of repair process was adopted, a closed form of G was derived analytically for arbitrary radiation schemes. The published animal data adopted to test the reciprocal formulas. Results: A generalized LQ model to describe the repair process in a reciprocal pattern was obtained. Subsequently, formulas for special cases were derived from this general form. The reciprocal model showed a better fit to the animal data than the exponential model, particularly for the ED50 data (reduced {chi}{sup 2}{sub min} of 2.0 vs 4.3, p = 0.11 vs 0.006), with the following gLQ parameters: {alpha}/{beta} = 2.6-4.8 Gy, T{sub r} = 3.2-3.9 h for rat feet skin, and {alpha}/{beta} = 0.9 Gy, T{sub r} = 1.1 h for rat spinal cord. Conclusions: These results of repair process following a reciprocal time suggest that the generalized LQ model incorporating the reciprocal time of sublethal damage repair shows a better fit than the exponential repair model. These formulas can be used to analyze the experimental and clinical data, where a slowing-down repair process appears during the course of radiation therapy.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2015

Galactic center γ-ray excess in hidden sector DM models with dark gauge symmetries: local Z{sub 3} symmetry as an example

We show that hidden sector dark matter (DM) models with local dark gauge symmetries make a natural playground for the possible γ-ray excess from the galactic center (GC). We first discuss in detail the GC γ-ray excess in a scalar dark matter (DM) model with local Z{sub 3} symmetry which was recently proposed by the present authors. Within this model, scalar DM with mass 30–70 GeV is allowed due to the newly-opened (semi-)annihilation channels of a DM pair into dark Higgs ϕ and/or dark photon Z′ pair, and the γ-ray spectrum from the GC can be fit within this model. Then we argue that the GC gamma ray excess can be easily accommodated within hidden sector dark matter models where DM is stabilized by local gauge symmetries, due to the presence of dark Higgs (and also dark photon for Abelian dark gauge symmetry)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Buzzard to Cardinal: Improved Mock Catalogs for Large Galaxy Surveys

Abstract We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including the Dark Energy Survey (DES), DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of z = 2.35 to a depth of m r = 27. Compared to the Buzzard mocks, the Cardinal mocks include an updated subhalo abundance matching model that considers orphan galaxies and includes mass-dependent scatter between galaxy luminosity and halo properties. This model can simultaneously fit galaxy clustering and group–galaxy cross-correlations measured in three different luminosity threshold samples. The Cardinal mocks also feature a new color assignment model that can simultaneously fit color-dependent galaxy clustering in three different luminosity bins. We have developed an algorithm that uses photometric data to further improve the color assignment model and have also developed a novel method to improve small-scale lensing below the ray-tracing resolution. These improvements enable the Cardinal mocks to accurately reproduce the abundance of galaxy clusters and the properties of lens galaxies in the DES data. As such, these simulations will be a valuable tool for future cosmological analyses based on large sky surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Identification of broadly neutralizing antibody epitopes in the HIV-1 envelope glycoprotein using evolutionary models

Background: Identification of the epitopes targeted by antibodies that can neutralize diverse HIV-1 strains can provide important clues for the design of a preventative vaccine. Methods: We have developed a computational approach that can identify key amino acids within the HIV-1 envelope glycoprotein that influence sensitivity to broadly cross-neutralizing antibodies. Given a sequence alignment and neutralization titers for a panel of viruses, the method works by fitting a phylogenetic model that allows the amino acid frequencies at each site to depend on neutralization sensitivities. Sites at which viral evolution influences neutralization sensitivity were identified using Bayes factors (BFs) to compare the fit of this model to that of a null model in which sequences evolved independently of antibody sensitivity. Conformational epitopes were identified with a Metropolis algorithm that searched for a cluster of sites with large Bayes factors on the tertiary structure of the viral envelope. Results: We applied our method to ID 50 neutralization data generated from seven HIV-1 subtype C serum samples with neutralization breadth that had been tested against a multi-clade panel of 225 pseudoviruses for which envelope sequences were also available. For each sample, between two and four sites were identified that were strongly associated with neutralization sensitivity (2ln(BF) > 6), a subset of which were experimentally confirmed using site-directed mutagenesis. Conclusions: Our results provide strong support for the use of evolutionary models applied to cross-sectional viral neutralization data to identify the epitopes of serum antibodies that confer neutralization breadth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Geometric compatibility of IceCube TeV-PeV neutrino excess and its galactic dark matter origin

Here, we perform a geometric analysis for the sky map of the IceCube TeV-PeV neutrino excess and test its compatibility with the sky map of decaying dark matter signals in our galaxy. Furthermore, we have found that a galactic decaying dark matter component in general improve the goodness of the fit of our model, although the pure isotropic hypothesis has a better fit than the pure dark matter one. Finally, we also consider several representative decaying dark matter, which can provide a good t to the observed spectrum at IceCube with a dark matter lifetime of around 12 orders of magnitude longer than the age of the universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Prediction of Cooling Curves for Controlled Unidirectional Solidification Under the Influence of Shrinkage: A Semi-analytical Approach

Quality of casting products can be significantly improved by controlling the crystal growth rate. Both dendritic and equiaxed grain structures depend strongly on crystal growth rate along the crystallographic direction. Uncontrolled solidification processes may lead to the formation of porous and columnar material with an extremely non-homogeneous composition distribution. These undesirable micro-structures can be avoided by adopting an optimal cooling rate (K/s) during the crystal growth. The present work focuses on determining suitable cooling curves to obtain desired unidirectional crystal growth rates, by implementing a semi-analytical model. Here, we defined cooling curve as the time history of temperature evolution at the cold boundary. The proposed semi-analytical heat transfer model is diffusion driven and accounts for shrinkage effect during solidification. The proposed model is found capable of predicting cooling curves for prescribed interface velocities. The model is self-sufficient in validating itself by means of solving inverse and forward problems. The cooling curve prediction by the model is validated by using it as the transient boundary condition for an existing enthalpy updating scheme-based numerical model. The proposed model is also applied to empirical data reported in literature, to validate the model against experimental results. All the validations showcased excellent fit. The proposed model will facilitate controlled unidirectional crystal growth rate by providing cooling rate (K/s) as an input for the experiments involving directional solidification.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Amplitude analysis of resonant production in three pions

We present some results on the analysis of three pion resonances. The analyses are motivated by the recent release of the largest data set on diffractively produced three pions by the COMPASS collaboration. We construct reaction amplitudes that satisfy fundamental $S$-matrix principles, which allows the use of models that have physical constraints to be used in fitting data. The models are motivated by the isobar model that satisfy unitarity constraints. The model consist of a Deck production amplitude with which final state interactions are constrained by unitarity. We employ the isobar model where two of the pions form a quasi-stable particle. The analysis is performed in the high-energy, single Regge limit. We specifically discuss the examples of the three pion $J^{PC}=2^{-+}$ resonance in the $\rho\pi$ and $f_2\pi$ channels.

Jackura, Andrew↗
DOE OSTI2012

Soil physical and hydrological properties under three biofuel crops in Ohio

While biofuel crops are widely studied and compared for their energy and carbon footprints, less is known about their effects on other soil properties, particularly hydrologic characteristics. Soils under three biofuel crops, corn (Zea mays), switchgrass (Panicum virgatum), and willow (Salix spp.), were analyzed seven years after establishment to assess the effects on soil bulk density ({rho}{sub b}), penetration resistance (PR), water-holding capacity, and infiltration characteristics. The PR was the highest under corn, along with the lowest associated water content, while PR was 50-60% lower under switchgrass. In accordance with PR data, surface (0-10 cm) bulk density also tended to be lower under switchgrass. Both water infiltration rates and cumulative infiltration amounts varied widely among and within the three crops. Because the Philip model did not fit the data, results were analyzed using the Kostiakov model instead. Switchgrass plots had an average cumulative infiltration of 69 cm over 3 hours with a constant infiltration rate of 0.28 cm min{sup -1}, compared with 37 cm and 0.11 cm min{sup -1} for corn, and 26 cm and 0.06 cm min{sup -1} for willow, respectively. Results suggest that significant changes in soil physical and hydrologic properties may require more time to develop. Soils under switchgrass may have lower surface bulk density, higher field water capacity, and a more rapid water infiltration rate than those under corn or willow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Deuterium Permeation Through a Self-Supported Palladium-Silver Membrane in Helium Gas Mixtures

Palladium membranes and membrane reactors can separate and purify tritium from impurities in the plasma exhaust processing system for the fusion energy fuel cycle. Membranes can also act as a continuous separation method to remove tritium from helium streams in the breeder blanket tritium extraction system, such as from the purge gas of solid breeders. To investigate the potential of these membranes for these applications, we performed a deuterium permeation campaign with a self-supported palladium-silver (Pd-25Ag wt%) tube of 0.15 m length, 76 µ m wall thickness, and 3.0 mm inner diameter. A gas mixture of 3.95% D 2 and a balance of He was delivered to the inside of the tube and D permeated radially outwards through the membrane into a vacuum chamber. Permeation experiments were conducted at 300 ° C, 350 ° C, 400° C, and 450°C, from 100 to 1000 sccm total flow rate, and with total pressures of 90, 150, 190, and 250 kPa. Further, deuterium permeation was consistently lower than predicted from diffusion-limited permeation models, thus we developed a transport model that included gas-phase mass transfer and surface reactions to model experimental results. The dissociation constant was optimized to fit the developed model to experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

TU-AB-303-11: Predict Parotids Deformation Applying SIS Epidemiological Model in H&N Adaptive RT

Purpose: The aim is to investigate the use of epidemiological models to predict morphological variations in patients undergoing radiation therapy (RT). The susceptible-infected-susceptible (SIS) deterministic model was applied to simulate warping within a focused region of interest (ROI). Hypothesis is to consider each voxel like a single subject of the whole sample and to treat displacement vector fields like an infection. Methods: Using Raystation hybrid deformation algorithms and automatic re-contouring based on mesh grid, we post-processed 360 MVCT images of 12 H&N patients treated with Tomotherapy. Study focused on parotid glands, identified by literature and previous analysis, as ROI more susceptible to warping in H&N region. Susceptible (S) and infectious (I) cases were identified in voxels with inter-fraction movement respectively under and over a set threshold. IronPython scripting allowed to export positions and displacement data of surface voxels for every fraction. A MATLAB homemade toolbox was developed to model the SIS. Results: SIS model was validated simulating organ motion on QUASAR phantom. Applying model in patients, within a [0–1cm] range, a single voxel movement of 0.4cm was selected as displacement threshold. SIS indexes were evaluated by MATLAB simulations. Dynamic time warping algorithm was used to assess matching between model and parotids behavior days of treatments. The best fit of the model was obtained with contact rate of 7.89±0.94 and recovery rate of 2.36±0.21. Conclusion: SIS model can follow daily structures evolutions, making possible to compare warping conditions and highlighting challenges due to abnormal variation and set-up errors. By epidemiology approach, organ motion could be assessed and predicted not in terms of average of the whole ROI, but in a voxel-by-voxel deterministic trend. Identifying anatomical region subjected to variations, would be possible to focus clinic controls within a cohort of pre-selected patients eligible for adaptive RT. The research is partially co-funded by the Italian Research Grant: Dose warping methods for IGRT and Adaptive RT: dose accumulation based on organ motion and anatomical variations of the patients during radiation therapy treatments,MoH (GR-2010-2318757) and Tecnologie Avanzate S.r.l.(Italy)

60 APPLIED LIFE SCIENCES↗