Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FIT”

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,495 records · Page 83

Space Shuttle RCS Oxidizer Leak Repair for STS-26

Following propellant loading of the Space Shuttle's reaction control system (RCS) for mission STS 26, an oxidizer leak was detected in the left orbital maneuvering system (OMS) pod, where the RCS is located. Subsequent investigation determined that the leak was isolated at a mechanical Dynatube fitting near the RCS nitrogen tetroxide tank. An intense effort was initiated to design, fabricate, and qualify a sealing device to stop the oxidizer leak externally so that the Space Shuttle launch could proceed. It was discovered that sealing devices called clamshells were widely used throughout the petrochemical and power generation industries to stop leaks developed in large diameter pipes which carry steam or other hazardous fluids. These clamshells are available in different diameters and strengths and are placed around the pipe at the location of the leak. A sealing compound is then injected under high pressure into the clamshell to stop the leak. This technology was scaled down and applied to the problem of stopping the leak on the Orbiter, which was on a half-inch diameter line in a nearly inaccessible location. Many obstacles had to be overcome such as determining that the sealing material would be compatible with the nitrogen tetroxide and ensuring that the clamshell would actually fit around the Dynatube fitting without interfering with other lines which were in close proximity. The effort at the NASA Johnson Space Center included materials compatibility testing of several sealants, design of a clamshell to fit in the confined compartment, and manufacture and qualification of the flight hardware. A clamshell was successfully placed around the Dynatube fitting on the Orbiter and the oxidizer leak was terminated. Then it was decided to apply this technology further and design clamshells for other mechanical fittings onboard the Orbiter and develop sealing compounds which will be compatible with fuels such as monomethyl hydrazine (MMH). The potential exists for using this type of sealing device in numerous other applications throughout the aerospace industry.

Delventhal, R. A.↗

Studies of Vestibular Neurons in Normal, Hyper- and Hypogravity

During the past year, pre-, in- and postflight studies were conducted in association with the Axon project for Bion 10 (Cosmos 2229). Recordings were made during pre- and postflight studies, from 118 horizontal semicircular canal afferents and 27 vestibular nucleus neurons in 7 rhesus monkeys; 137 pulse rotation protocols alone were executed (548 acceleration and deceleration responses were curve fit). Usable data was obtained from 127 horizontal afferents concerning their spontaneous discharge. Curve fits and analysis was made of sinusoidal and sum of sinusoidal responses from 42 and 35 horizontal afferents, respectively. Also recordings were made from neurons inflight from the two flight animals. The mean spontaneous rate varied from 128 spikes/sec. during preflight to 92 spikes/sec during postflight (day 5) - a change of 28%. In direct contrast to the results of Cosmos 2044, the best fitted neural adaptation operator (k) and the gain of the pulse response were decreased during post flight when compared to preflight. Surprisingly, the best fitted gain and k values for the sum of sines were slightly elevated during post flight tests. The gain and phase of single sine responses were compared for pre- and post flight tests and compared to a larger population of afferents. In contrast to Cosmos 2044 results where on the first day of post flight testing the gains of the best fitted sine response were skewed toward the higher values of the Miles and Braitman distribution, the gain of the best fitted sine responses during the first day of post flight testing (day 2) during Cosmos 2229 were exactly on the mode of the Miles and Braitman distribution. Thus, at least for the periodic stimuli, (pulses and sine waves) we found no change in gain and neural adaptation during post flight testing following Cosmos 2229. This conclusion is different from the one derived following the Cosmos 2044 flight. Cosmos flight 2229 differed from Cosmos flight 2044 in several significant ways: For example, during preflight, (1) The animals preflight training was different (less well trained on the gaze task) and (2) the animals were exposed to more experimental manipulations (surgical and rotational). Inflight, (1) the animals were required to make a pointing gesture (motor response) in association with eye movements to obtain reward, (2) the inflight diet was different (more balanced), (3) the feeder for one of the animals clogged following 9 days of flight resulting in evident dehydration and probably less head motion exposure in that monkey and (4) there was limited video taping of the monkeys in space. During postflight, (1) we were unable to test the flight animals until 26 hours postflight as compared to 14.5 hours during Cosmos 2044, (2) the animals received significantly more exposure to motion stimuli during postflight testing than during Cosmos 2044. These differences in the vestibular environment will require analysis of several parameters other than just neural and eye movement responses. For example, computer programs will have to be written and used to recover and quantify the number of head movements made by each animal during flight. This activity is critical to the production of neural adaptation and increased gain.

Correia, Manning J.↗

A Study of the Non-Thermal X-Ray Emission of Shell-Type Supernova Remnants

We present an analysis of the X-ray spectrum of the Galactic shell-type SNR G347.3-0.5 (RX 51713.7-3946). This SNR is a member of a growing class of SNRs which are dynamically young, shell-type sources that emit non-thermal X-rays from specific regions on their outer shells. By performing a joint spectral analysis of data from observations made of G347.3-0.5 using the ROSAT PSPC, the ASCA GIS and the RXTE PCA, we have fit the spectra of particular regions of this SNR (including the bright northwestern and southwestern rims, the northeast rim and the interior diffuse emission) over the approximate energy range of 0.5 through 30 keV. We find that fits to the spectra of this SNR over this energy range using the SRCUT model were superior to a simple power law model or the SRESC model. We find that the inclusion of a thermal model with the SRCUT model helps to improve the fit to the observed X-ray spectrum: this represents the first detection of thermal X-ray emission from G347.3-0.5. Thermal emission appears to be more clearly associated with the diffuse emission in the interior of the SNR than with the bright X-ray emitting rims. A weak emission feature seen near 6.4 keV in the RXTE PCA spectrum most likely originates from diffuse X-ray emission from the surrounding Galactic Ridge rather than from G347.3-0.5 itself. We have analyzed our RXTE PCA data to search for pulsations from a recently discovered radio pulsar (PSR 51713-3949) which may be associated with G347.3-0.5, and we do not detect any X-ray pulsations at the measured radio period of 392 ms. Using the best-fit parameters obtained from the SRCUT model, we estimate the maximum energy of cosmic-ray electrons accelerated by the rims of G347.3-0.5 to be 19-25 TeV (assuming a magnetic field strength of B = 10muG), consistent with the results of Ellison et al. We present a broadband (radio to gamma ray) photon energy-flux spectrum for the northwestern rim of G347.3-0.5, where we have fit the spectrum using a more sophisticated synchrotron-inverse Compton model with a variable magnetic field strength. Our fit derived from this model yields a maximum energy of only 8.8 (+4.1) (-3.4) TeV for the accelerated cosmic-ray electrons and a much greater magnetic field strength of 150 (+250) (-80). Finally, we compare the gross properties of G347.3-0.5 with other SNRs known to possess X-ray spectra dominated by non-thermal emission.

Allen, Glenn E.↗

Evolving Multi Rover Systems in Dynamic and Noisy Environments

In this chapter, we address how to evolve control strategies for a collective: a set of entities that collectively strives to maximize a global evaluation function that rates the performance of the full system. Addressing such problems by directly applying a global evolutionary algorithm to a population of collectives is unworkable because the search space is prohibitively large. Instead, we focus on evolving control policies for each member of the collective, where each member is trying to maximize the fitness of its own population. The main difficulty with this approach is creating fitness evaluation functions for the members of the collective that induce the collective to achieve high performance with respect to the system level goal. To overcome this difficulty, we derive member evaluation functions that are both aligned with the global evaluation function (ensuring that members trying to achieve high fitness results in a collective with high fitness) and sensitive to the fitness of each member (a member's fitness depends more on its own actions than on actions of other members). In a difficult rover coordination problem in dynamic and noisy environments, we show how to construct evaluation functions that lead to good collective behavior. The control policy evolved using aligned and member-sensitive evaluations outperforms global evaluation methods by up to a factor of four. in addition we show that the collective continues to perform well in the presence of high noise levels and when the environment is highly dynamic. More notably, in the presence of a larger number of rovers or rovers with noisy sensors, the improvements due to the proposed method become significantly more pronounced.

Tumer, Kagan↗

Spline-Based Smoothing of Airfoil Curvatures

Constrained fitting for airfoil curvature smoothing (CFACS) is a splinebased method of interpolating airfoil surface coordinates (and, concomitantly, airfoil thicknesses) between specified discrete design points so as to obtain smoothing of surface-curvature profiles in addition to basic smoothing of surfaces. CFACS was developed in recognition of the fact that the performance of a transonic airfoil is directly related to both the curvature profile and the smoothness of the airfoil surface. Older methods of interpolation of airfoil surfaces involve various compromises between smoothing of surfaces and exact fitting of surfaces to specified discrete design points. While some of the older methods take curvature profiles into account, they nevertheless sometimes yield unfavorable results, including curvature oscillations near end points and substantial deviations from desired leading-edge shapes. In CFACS as in most of the older methods, one seeks a compromise between smoothing and exact fitting. Unlike in the older methods, the airfoil surface is modified as little as possible from its original specified form and, instead, is smoothed in such a way that the curvature profile becomes a smooth fit of the curvature profile of the original airfoil specification. CFACS involves a combination of rigorous mathematical modeling and knowledge-based heuristics. Rigorous mathematical formulation provides assurance of removal of undesirable curvature oscillations with minimum modification of the airfoil geometry. Knowledge-based heuristics bridge the gap between theory and designers best practices. In CFACS, one of the measures of the deviation of an airfoil surface from smoothness is the sum of squares of the jumps in the third derivatives of a cubicspline interpolation of the airfoil data. This measure is incorporated into a formulation for minimizing an overall deviation- from-smoothness measure of the airfoil data within a specified fitting error tolerance. CFACS has been extensively tested on a number of supercritical airfoil data sets generated by inverse design and optimization computer programs. All of the smoothing results show that CFACS is able to generate unbiased smooth fits of curvature profiles, trading small modifications of geometry for increasing curvature smoothness by eliminating curvature oscillations and bumps (see figure).

Li, W.↗

Dust and Gas in the Magellanic Clouds from the Heritage Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm (Submillimeter) Excess Emission

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micromillimeters. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a powerlaw emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submillimeter excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micromillimeters. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micromillimeters submillimeter excesses of 27% and 43% for the Large and Small Magellanic Clouds, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our fitting results show that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submillimeter excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 plus or minus 1.7) x 10 (sup 5) and (8.3 plus or minus 2.1) x 10 (sup 4) solar masses for the Large and Small Magellanic Clouds, respectively. We find significant correlations between the submillimeter excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

Magellanic↗

Investigating the Differential Emission Measure and Energetics of Microflares with Combined SDO/AIA and RHESSI Observations

An important question in solar physics is whether solar microflares, the smallest currently observable flare events in X-rays, possess the same energetic properties as large flares. Recent surveys have suggested that microflares may be less efficient particle accelerators than large flares, and hence contribute less non-thermal energy, which may have implications for coronal heating mechanisms. We therefore explore the energetic properties of microflares by combining EUV and X-ray measurements. We present forward-fitting differential emission measure (DEM) analysis of 10 microflares. The fitting is constrained by combining, for the first time, high-temperature Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations and flux data from the Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA). Two fitting models are tested for the DEM; a Gaussian distribution and a uniform DEM profile. A Gaussian fit proved unable to explain the observations for any of the studied microflares. However, 8 of 10 events studied were reasonably fit by a uniform DEM profile. Hence microflare plasma can be considered to be significantly multi-thermal, and may not be significantly peaked or contain resolvable fine structure, within the uncertainties of the observational instruments. The thermal and non-thermal energy is estimated for each microflare, comparing the energy budget with an isothermal plasma assumption. From the multi-thermal fits the minimum non-thermal energy content was found to average approximately 30% of the estimated thermal energy. By comparison, under an isothermal model the non-thermal and thermal energy estimates were generally comparable. Hence, multi-thermal plasma is an important consideration for solar microflares that substantially alters their thermal and non-thermal energy content.

Sun: flares↗

Light-Curve Modelling Constraints on the Obliquities and Aspect Angles of the Young Fermi Pulsars

In more than four years of observation the Large Area Telescope on board the Fermi satellite has identified pulsed gamma-ray emission from more than 80 young or middle-aged pulsars, in most cases providing light curves with high statistics. Fitting the observed profiles with geometrical models can provide estimates of the magnetic obliquity alpha and of the line of sight angle zeta, yielding estimates of the radiation beaming factor and radiated luminosity. Using different gamma-ray emission geometries (Polar Cap, Slot Gap, Outer Gap, One Pole Caustic) and core plus cone geometries for the radio emission, we fit gamma-ray light curves for 76 young or middle-aged pulsars and we jointly fit their gamma-ray plus radio light curves when possible. We find that a joint radio plus gamma-ray fit strategy is important to obtain (alpha, zeta) estimates that can explain simultaneously detectable radio and gamma-ray emission: when the radio emission is available, the inclusion of the radio light curve in the fit leads to important changes in the (alpha, gamma) solutions. The most pronounced changes are observed for Outer Gap and One Pole Caustic models for which the gamma-ray only fit leads to underestimated alpha or zeta when the solution is found to the left or to the right of the main alpha-zeta plane diagonal respectively. The intermediate-to-high altitude magnetosphere models, Slot Gap, Outer Gap, and One pole Caustic, are favored in explaining the observations. We find no apparent evolution of α on a time scale of 106 years. For all emission geometries our derived gamma-ray beaming factors are generally less than one and do not significantly evolve with the spin-down power. A more pronounced beaming factor vs. spin-down power correlation is observed for Slot Gap model and radio-quiet pulsars and for the Outer Gap model and radio-loud pulsars. The beaming factor distributions exhibit a large dispersion that is less pronounced for the Slot Gap case and that decreases from radio-quiet to radio-loud solutions. For all models, the correlation between gamma-ray luminosity and spin-down power is consistent with a square root dependence. The gamma-ray luminosities obtained by using the beaming factors estimated in the framework of each model do not exceed the spin-down power. This suggests that assuming a beaming factor of one for all objects, as done in other studies, likely overestimates the real values. The data show a relation between the pulsar spectral characteristics and the width of the accelerator gap. The relation obtained in the case of the Slot Gap model is consistent with the theoretical prediction.

neutron aEuro" pulsars:↗

Wind Tunnel Strain-Gage Balance Calibration Data Analysis Using a Weighted Least Squares Approach

A new approach is presented that uses a weighted least squares fit to analyze wind tunnel strain-gage balance calibration data. The weighted least squares fit is specifically designed to increase the influence of single-component loadings during the regression analysis. The weighted least squares fit also reduces the impact of calibration load schedule asymmetries on the predicted primary sensitivities of the balance gages. A weighting factor between zero and one is assigned to each calibration data point that depends on a simple count of its intentionally loaded load components or gages. The greater the number of a data point's intentionally loaded load components or gages is, the smaller its weighting factor becomes. The proposed approach is applicable to both the Iterative and Non-Iterative Methods that are used for the analysis of strain-gage balance calibration data in the aerospace testing community. The Iterative Method uses a reasonable estimate of the tare corrected load set as input for the determination of the weighting factors. The Non-Iterative Method, on the other hand, uses gage output differences relative to the natural zeros as input for the determination of the weighting factors. Machine calibration data of a six-component force balance is used to illustrate benefits of the proposed weighted least squares fit. In addition, a detailed derivation of the PRESS residuals associated with a weighted least squares fit is given in the appendices of the paper as this information could not be found in the literature. These PRESS residuals may be needed to evaluate the predictive capabilities of the final regression models that result from a weighted least squares fit of the balance calibration data.

calibration analysis↗

Program to Optimize Simulated Trajectories II (POST2) Surrogate Models for Mars Ascent Vehicle (MAV) Performance Assessment

The primary purpose of the multiPOST tool is to enable the execution of much larger sets of vehicle cases to allow for broader trade space exploration. However, this exploration is not achieved solely with the increased case throughput. The multiPOST tool is applied to carry out a Design of Experiments (DOE), which is a set of cases that have been structured to capture a maximum amount of information about the design space with minimal computational effort. The results of the DOE are then used to fit a surrogate model, ultimately enabling parametric design space exploration. The approach used for the MAV study includes both DOE and surrogate modeling. First, the primary design considerations for the vehicle were used to develop the variables and ranges for the multiPOST DOE. The final set of DOE variables were carefully selected in order to capture the desired vehicle trades and take into account any special considerations for surrogate modeling. Next, the DOE sets were executed through multiPOST. Following successful completion of the DOE cases, a manual verification trial was performed. The trial involved randomly selecting cases from the DOE set and running them by hand. The results from the human analyst's run and multiPOST were then compared to ensure that the automated runs were being executed properly. Completion of the verification trials was then followed by surrogate model fitting. After fits to the multiPOST data were successfully created, the surrogate models were used as a stand-in for POST2 to carry out the desired MAV trades. Using the surrogate models in lieu of POST2 allowed for visualization of vehicle sensitivities to the input variables as well as rapid evaluation of vehicle performance. Although the models introduce some error into the output of the trade study, they were very effective at identifying areas of interest within the trade space for further refinement by human analysts. The next section will cover all of the ground rules and assumptions associated with DOE setup and multiPOST execution. Section 3.1 gives the final DOE variables and ranges, while section 3.2 addresses the POST2 specific assumptions. The results of the verification trials are given in section 4. Section 5 gives the surrogate model fitting results, including the goodness-of-fit metrics for each fit. Finally, the MAV specific results are discussed in section 6.

Zwack, M. R.↗

The Atacama Cosmology Telescope: delensed power spectra and parameters

We present ΛCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (T T), 30% (EE), 26% (T E) and 20% (BB); this results in detections of the delensing effect at 8.7σ (T T), 5.1σ (EE), 2.6σ (T E), and 2.4σ (BB) significance. The combination of 150 and 98 GHz T T, EE, and T E delensed spectra is well fit by a standard ΛCDM model. We also measure the shift in best-fit parameters when fitting delensed versus lensed spectra; while this shift does not inform our ability to measure cosmological parameters, it does provide a three-way consistency check among the lensing inferred from the best-fit parameters, the lensing in the CMB power spectrum, and the reconstructed lensing map. This shift is predicted to be zero when fitting with the correct model since both lensed and delensed spectra originate from the same region of sky. Fitting with a ΛCDM model and marginalizing over foregrounds, we find that the shift in cosmological parameters is consistent with zero. Our results show that gravitational lensing of the microwave background is internally consistent within the framework of the standard cosmological model.

Dongwon Han↗

Growth Curve Parameterization of Metabolic Activity of Yeast Cells for BioSentinel

The goal of the BioSentinel small satellite payload is to measure the effect of deep space radiation on the growth and metabolic activity of yeast cells. Raw test data is generated by fluidics cards containing yeast cells rehydrated at different periods, with metabolic activity measured by the reduction of alamarBlue. Each card well has a sensor array that measures the amount of red, green, and infrared light transmitted through the yeast culture. This illumination data is then converted to absorbance values, which are further converted into concentrations. The ultimate objective is to convert these concentrations into biologically-relevant metrics that can be compared against one another to determine changes due to differential radiation exposure. Beginning with IR absorbance data (corresponding to cell density) from ground studies, three parameters from a sigmoidal growth curve were extracted and analyzed: 𝜆 (lag phase), 𝜇 (max growth rate), and A (max cell growth). The data was fit to the Gompertz model of microbial growth using non-linear regression (Minitab), as the fit error was reduced compared to the simpler logistic growth curve. Graphs showed that the data contained a discrepancy (drift) in the lag phase that is attributable to a slow, constant loss of moisture. Correcting this discrepancy by fitting the first 25 hours of the data to a power function and subtracting these values from the absorbance readings obtained a better statistical fit to the growth curve in the lag phase. A power fit was selected over a linear fit because it reflected the effects of constant volume loss. This correction to the BioSentinel data analysis pipeline will enable quantitative statistical analysis of the effect of different levels of deep space radiation on yeast cells. Future work includes automation of drift correction and curve modeling to extract these parameters directly from data.

Growth Curve↗

Solar Temperature Variations Computed from SORCE SIM Irradiances Observed During 2003-2020

NASA’s Solar Radiation and Climate Experiment (SORCE) Spectral Irradiance Monitor (SIM) instrument produced about 17 years of daily average Spectral Solar Irradiance ( SSI ) data for wavelengths 240 nm – 2416 nm. We choose a day of minimal solar activity, 2008-08-24, during the 2008 − 2009 minimum between cycles 23 and 24, and compute the brightness temperature (𝑇 o ) from that day’s solar spectral irradiance (𝑆𝑆𝐼 o ). We consider small variations of T and SSI about these reference values, and derive linear and quadratic analytic approximations by Taylor expansion about the reference day values. To determine approximation accuracy, we compare to exact brightness temperatures T computed from the Planck spectrum, by solving analytically for T , or equivalent root-finding in Wolfram Mathematica. We find that the linear analytic approximation overestimates, while the quadratic underestimates the exact result. This motivates search for statistical “fit” models “in between” the two analytic models, with minimum root-mean-square-error RMSE. We make this search using open-source statistical R software, determine coefficients for linear and quadratic fit models, and compare statistical with analytic RMSE’s. When only linear analytic and fit models are compared, the fit model is superior at ultraviolet, visible, and near infrared wavelengths. This again holds true when comparing only quadratic models. Quadratic is superior to linear for both analytic and statistical models, and statistical fits give smallest RMSE’s. Lastly, we use linear analytic and fit models to find an interpolating function in wavelength, useful in case the SIM results need adjustment to another choices of wavelengths, to compare or extend to any other instrument.

SORCE↗

Plasma Physics Approximations in Ares

Lee & More derived analytic forms for the transport properties of a plasma. Many hydro-codes use their formulae for electrical and thermal conductivity. The coefficients are complex functions of Fermi-Dirac integrals, F n ( μ/θ ), the chemical potential, μ or ζ = ln(1+e μ/θ ), and the temperature, θ = kT. Since these formulae are expensive to compute, rational function approximations were fit to them. Approximations are also used to find the chemical potential, either μ or ζ . The fits use ζ as the independent variable instead of μ/θ . New fits are provided for A α (ζ ),A β (ζ ), ζ, f(ζ ) = (1 + e -μ/θ )F 1/2 (μ/θ), F 1/2 '/F 1/2 , F c α , and F c β . In each case the relative error of the fit is minimized since the functions can vary by many orders of magnitude. The new fits are designed to exactly preserve the limiting values in the non-degenerate and highly degenerate limits or as ζ→ 0 or ∞. The original fits due to Lee & More and George Zimmerman are presented for comparison.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A High-performance Track Fitter for Use in Ultra-fast Electronics

This article describes a new charged-particle track fitting algorithm designed for use in high-speed electronics applications such as hardware-based triggers in high-energy physics experiments. Following a novel technique designed for fast electronics, the positions of the hits on the detector are transformed before being passed to a linearized track parameter fit. This transformation results in fitted track parameters with a very linear dependence on the hit positions. The approach is demonstrated in a representative detector geometry based on the CMS detector at the Large Hadron Collider. The fit is implemented in FPGA chips and optimized for track fitting throughput and obtains excellent track parameter performance. Such an algorithm is potentially useful in any high-speed track-fitting application.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of Schedule Duration on Head and Neck Radiotherapy: Accelerated Tumor Repopulation Versus Compensatory Mucosal Proliferation

Purpose: To determine how modelled maximum tumor control rates, achievable without exceeding mucositis tolerance (tcp{sub max-early}) vary with schedule duration for head and neck squamous cell carcinoma (HNSCC). Methods and materials: Using maximum-likelihood techniques, we have fitted a range of tcp models to two HNSCC datasets (Withers' and British Institute of Radiology [BIR]), characterizing the dependence of tcp on duration and equivalent dose in 2 Gy fractions (EQD{sub 2}). Models likely to best describe future data have been selected using the Akaike information criterion (AIC) and its quasi-AIC extension to overdispersed data. Setting EQD{sub 2}s in the selected tcp models to levels just tolerable for mucositis, we have plotted tcp{sub max-early} against schedule duration. Results: While BIR dataset tcp fits describe dose levels isoeffective for tumor control as rising significantly with schedule protraction, indicative of accelerated tumor repopulation, repopulation terms in fits to Withers' dataset do not reach significance after accounting for overdispersion of the data. The tcp{sub max-early} curves calculated from tcp fits to the overall Withers' and BIR datasets rise by 8% and 0-4%, respectively, between 20 and 50 days duration; likewise, tcp{sub max-early} curves calculated for stage-specific cohorts also generally rise slowly with increasing duration. However none of the increases in tcp{sub max-early} calculated from the overall or stage-specific fits reach significance. Conclusions: Local control rates modeled for treatments which lie just within mucosal tolerance rise slowly but insignificantly with increasing schedule length. This finding suggests that whereas useful gains may be made by accelerating unnecessarily slow schedules until they approach early reaction tolerance, little is achieved by shortening schedules further while reducing doses to remain within mucosal tolerance, an approach that may slightly worsen outcomes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Assessment of Lung Tumour Motion and Volume Size Dependencies Using Various Evaluation Measures

The purpose of the study was to assess internal target volume changes through the breathing cycle and associated tumour motion for lung patients and to establish possible correlations between different parameters. Respiration-induced volume changes with breathing cycle and the associated tumour motion were analyzed for 11 patients. Selected phases were the maximum and average intensity projections and the 10 phases of equal duration and separation obtained through the respiratory cycle. Tumour centre of mass (COM) motion planes were generated using least square fitting, and their angles and orientations were then compared between the cases studied. Trajectories that are composed by the points of COM location in different phases were identified, and their interrelation was assessed through different similarity measures. The results were used to determine if there is any correlation between parameters chosen and if the margins conventionally used for the planning target volume creation successfully encompassed lung tumour motion and volume change. The results show that the extent of tumour motion was related to its volume and location. The tumour displacement was predominantly left and inferior. Planar fitting to COM motion data through respiratory phases demonstrated some correlation in best fit motion plane positions between different data sets. In the plane fit comparison, for each patient, the lower root mean square error values showed that a good planar fit can be achieved for the COM motion path. The evaluation of the inhale and exhale trajectories may allow, for certain tumour locations and size, contouring on only inhalation or exhalation phases, knowing that tumour motion will be adequately covered on the other phases. Taking all the data into account and knowing the tumour size and location, a good estimate can be made of the motion plane position in the three-dimensional space and the required dosimetric margins.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Five ab initio potential energy and dipole moment surfaces for hydrated NaCl and NaF. I. Two-body interactions

We report full-dimensional, ab initio-based potentials and dipole moment surfaces for NaCl, NaF, Na{sup +}H{sub 2}O, F{sup −}H{sub 2}O, and Cl{sup −}H{sub 2}O. The NaCl and NaF potentials are diabatic ones that dissociate to ions. These are obtained using spline fits to CCSD(T)/aug-cc-pV5Z energies. In addition, non-linear least square fits using the Born-Mayer-Huggins potential are presented, providing accurate parameters based strictly on the current ab initio energies. The long-range behavior of the NaCl and NaF potentials is shown to go, as expected, accurately to the point-charge Coulomb interaction. The three ion-H{sub 2}O potentials are permutationally invariant fits to roughly 20 000 coupled cluster CCSD(T) energies (awCVTZ basis for Na{sup +} and aVTZ basis for Cl{sup −} and F{sup −}), over a large range of distances and H{sub 2}O intramolecular configurations. These potentials are switched accurately in the long range to the analytical ion-dipole interactions, to improve computational efficiency. Dipole moment surfaces are fits to MP2 data; for the ion-ion cases, these are well described in the intermediate- and long-range by the simple point-charge expression. The performance of these new fits is examined by direct comparison to additional ab initio energies and dipole moments along various cuts. Equilibrium structures, harmonic frequencies, and electronic dissociation energies are also reported and compared to direct ab initio results. These indicate the high fidelity of the new PESs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗