Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 3”

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 379 records · Page 21

WE-AB-204-08: Uncertainties in Kinetic Modeling of Shortened Dynamic 18F-FMISO PET Due to Image Noise

Purpose: To assess uncertainties in voxel-wise kinetic modeling (KM) of dynamic 18F-FMISO PET (dPET) due to image noise for different lengths of dPET acquisitions. Methods: 12 tumor time activity curves (TACs) deduced from 6 head and neck cancer patient dPET datasets (45min dPET + 10min frames at ∼95min and ∼160min respectively) were modeled using an irreversible two-tissue compartment model to estimate kinetic rate constants (KRCs) k3, vB, K1, and K1/k2 (reference standard). For each modeled TAC, 1000 noisy TACs were simulated by adding Gaussian noise with standard deviation equivalent to that observed on a voxel level in GE DSTE PET-CT. KM was conducted for each set of noisy TACs using (i) full dataset (170min), and repeated using (ii) 105min and (iii) 45min shortened subsets, with an input function that was image-derived from the dPET dataset of matching length. Absolute value of percent difference between KRCs from noisy TACs as estimated from either of 3 datasets and reference standard KRCs was used to represent bias. Results: For all KRCs, bias was higher for shortened datasets. For k3, bias was lower when actual k3 was higher: from 7%±5% (170min), 9%±7% (105min), and 14%±11% (45min) (actual k3=0.0098), to 32%±25% (170min), 45%±34% (95min), and 79%±48% (45min) (actual k3=0.0008). Similar relationship was observed for vB, with bias being between ∼7% (actual vB=0.28) and ∼26–34% (actual vB=0.08). For K1 [K1/k2], bias was ∼6–10% [∼2–3%] (dependence on actual value of K1 [K1/k2] was not observed). Conclusion: Uncertainties in voxel-wise KM of shortened 18F-FMISO dPET due to image noisy only are larger for shorter acquisitions and, for k3 [vB], the bias was found to be inversely correlated with the actual values of k3 [vB]. Using 45min subset, k3, vB, K1, and K1/k2 could be estimated to within ∼15–80%, ∼10–35%, ∼10%, and ∼3%, respectively.

60 APPLIED LIFE SCIENCES↗

TU-EF-204-11: Impact of Using Multi-Slice Training Sets On the Performance of a Channelized Hotelling Observer in a Low-Contrast Detection Task in CT

Purpose: To investigate using multiple CT image slices from a single acquisition as independent training images for a channelized Hotelling observer (CHO) model to reduce the number of repeated scans for CHO-based CT image quality assessment. Methods: We applied a previously validated CHO model to detect low contrast disk objects formed from cross-sectional images of three epoxy-resin-based rods (diameters: 3, 5, and 9 mm; length: ∼5cm). The rods were submerged in a 35x 25 cm2 iodine-doped water filled phantom, yielding-15 HU object contrast. The phantom was scanned 100 times with and without the rods present. Scan and reconstruction parameters include: 5 mm slice thickness at 0.5 mm intervals, 120 kV, 480 Quality Reference mAs, and a 128-slice scanner. The CHO’s detectability index was evaluated as a function of factors related to incorporating multi-slice image data: object misalignment along the z-axis, inter-slice pixel correlation, and number of unique slice locations. In each case, the CHO training set was fixed to 100 images. Results: Artificially shifting the object’s center position by as much as 3 pixels in any direction relative to the Gabor channel filters had insignificant impact on object detectability. An inter-slice pixel correlation of >∼0.2 yielded positive bias in the model’s performance. Incorporating multi-slice image data yielded slight negative bias in detectability with increasing number of slices, likely due to physical variations in the objects. However, inclusion of image data from up to 5 slice locations yielded detectability indices within measurement error of the single slice value. Conclusion: For the investigated model and task, incorporating image data from 5 different slice locations of at least 5 mm intervals into the CHO model yielded detectability indices within measurement error of the single slice value. Consequently, this methodology would Result in a 5-fold reduction in number of image acquisitions. This project was supported by National Institutes of Health grants R01 EB017095 and U01 EB017185 from the National Institute of Biomedical Imaging and Bioengineering.

60 APPLIED LIFE SCIENCES↗

CTF (Subchannel) Calculations and Validation L3:VVI.H2L.P15.01

The goal of the Verification and Validation Implementation (VVI) High to Low (Hi2Lo) process is utilizing a validated model in a high resolution code to generate synthetic data for improvement of the same model in a lower resolution code. This process is useful in circumstances where experimental data does not exist or it is not sufficient in quantity or resolution. Data from the high-fidelity code is treated as calibration data (with appropriate uncertainties and error bounds) which can be used to train parameters that affect solution accuracy in the lower-fidelity code model, thereby reducing uncertainty. This milestone presents a demonstration of the Hi2Lo process derived in the VVI focus area. The majority of the work performed herein describes the steps of the low-fidelity code used in the process with references to the work detailed in the companion high-fidelity code milestone (Reference 1). The CASL low-fidelity code used to perform this work was Cobra Thermal Fluid (CTF) and the high-fidelity code was STAR-CCM+ (STAR). The master branch version of CTF (pulled May 5, 2017 – Reference 2) was utilized for all CTF analyses performed as part of this milestone. The statistical and VVUQ components of the Hi2Lo framework were performed using Dakota version 6.6 (release date May 15, 2017 – Reference 3). Experimental data from Westinghouse Electric Company (WEC – Reference 4) was used throughout the demonstrated process to compare with the high-fidelity STAR results. A CTF parameter called Beta was chosen as the calibration parameter for this work. By default, Beta is defined as a constant mixing coefficient in CTF and is essentially a tuning parameter for mixing between subchannels. Since CTF does not have turbulence models like STAR, Beta is the parameter that performs the most similar function to the turbulence models in STAR. The purpose of the work performed in this milestone is to tune Beta to an optimal value that brings the CTF results closer to those measured in the WEC experiments.

97 MATHEMATICS AND COMPUTING↗

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗

The inner zone electron model AE-5

A description is given of the work performed in the development of the inner radiation zone electron model, AE-5. A complete description of the omnidirectional flux model is given for energy thresholds E sub T in the range 4.0 E sub T/(MeV) 0.04 and for L values in the range 2.8 L 1.2 for an epoch of October 1967. Confidence codes for certain regions of B-L space and certain energies are given based on data coverage and the assumptions made in the analysis. The electron model programs that can be supplied to a user are referred to. One of these, a program for accessing the model flux at arbitrary points in B-L space and arbitrary energies, includes the latest outer zone electron model and proton model. The model AE-5, is based on data from five satellites, OGO 1, OGO 3, 1963-38C, OV3-3, and Explorer 26, spanning the period December 1964 to December 1967.

Teague, M. J.↗

Magnetospheric substorms

The generation of magnetospheric substorms as a magnetospheric response to a rectangular wave of a component of the interplanetary magnetic field is discussed. The development and decay of auroral substorms (the only visible manifestations of magnetic substorms) are described with reference to auroral particle precipitation, joule heat dissipation, and ring current injection. Various models of substorm phenomena are reviewed, including: (1) the conversion of magnetotail magnetic energy, (2) hot plasma injection from the plasma sheet into the Van Allen belt and ring current formation, (3) field-aligned currents and the auroral electrojet, and (4) the nature of interplanetary magnetic field fluctuations.

Akasofu, S.-I.↗

3-D inelastic analysis methods for hot section components (base program)

A 3-D inelastic analysis methods program consists of a series of computer codes embodying a progression of mathematical models (mechanics of materials, special finite element, boundary element) for streamlined analysis of combustor liners, turbine blades, and turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (dynamics, buckling) structural behavior of the three selected components. These models are used to solve 3-D inelastic problems using linear approximations in the sense that stresses/strains and temperatures in generic modeling regions are linear functions of the spatial coordinates, and solution increments for load, temperature and/or time are extrapolated linearly from previous information. Three linear formulation computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (MARC-Hot Section Technology), and BEST (Boundary Element Stress Technology), were developed and are described.

Wilson, R. B.↗

The NASA MSFC Earth Global Reference Atmospheric Model-2007 Version

Reference or standard atmospheric models have long been used for design and mission planning of various aerospace systems. The NASA/Marshall Space Flight Center (MSFC) Global Reference Atmospheric Model (GRAM) was developed in response to the need for a design reference atmosphere that provides complete global geographical variability, and complete altitude coverage (surface to orbital altitudes) as well as complete seasonal and monthly variability of the thermodynamic variables and wind components. A unique feature of GRAM is that, addition to providing the geographical, height, and monthly variation of the mean atmospheric state, it includes the ability to simulate spatial and temporal perturbations in these atmospheric parameters (e.g. fluctuations due to turbulence and other atmospheric perturbation phenomena). A summary comparing GRAM features to characteristics and features of other reference or standard atmospheric models, can be found Guide to Reference and Standard Atmosphere Models. The original GRAM has undergone a series of improvements over the years with recent additions and changes. The software program is called Earth-GRAM2007 to distinguish it from similar programs for other bodies (e.g. Mars, Venus, Neptune, and Titan). However, in order to make this Technical Memorandum (TM) more readable, the software will be referred to simply as GRAM07 or GRAM unless additional clarity is needed. Section 1 provides an overview of the basic features of GRAM07 including the newly added features. Section 2 provides a more detailed description of GRAM07 and how the model output generated. Section 3 presents sample results. Appendices A and B describe the Global Upper Air Climatic Atlas (GUACA) data and the Global Gridded Air Statistics (GGUAS) database. Appendix C provides instructions for compiling and running GRAM07. Appendix D gives a description of the required NAMELIST format input. Appendix E gives sample output. Appendix F provides a list of available parameters to enable the user to generate special output. Appendix G gives an example and guidance on incorporating GRAM07 as a subroutine in other programs such as trajectory codes or orbital propagation routines.

Leslie, F.W.↗

Phonons, Diffusons, and the Boson Peak in Two-Dimensional Lattices with Random Bonds

Within the model of stable random matrices possessing translational invariance, a two-dimensional (on a square lattice) disordered oscillatory system with random strongly fluctuating bonds is considered. By a numerical analysis of the dynamic structure factor S(q, ω), it is shown that vibrations with frequencies below the Ioffe-Regel frequency ω{sub IR} are ordinary phonons with a linear dispersion law ω(q) ∝ q and a reciprocal lifetime δ~ q{sup 3}. Vibrations with frequencies above ω{sub IR}, although being delocalized, cannot be described by plane waves with a definite dispersion law ω(q). They are characterized by a diffusion structure factor with a reciprocal lifetime δ ~ q{sup 2}, which is typical of a diffusion process. In the literature, they are often referred to as diffusons. It is shown that, as in the three-dimensional model, the boson peak at the frequency ωb in the reduced density of vibrational states g(ω)/ω is on the order of the frequency ω{sub IR}. It is located in the transition region between phonons and diffusons and is proportional to the Young’s modulus of the lattice, ω{sub b}≃E.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

End-Use Savings Shapes Upgrade Package Documentation: Wall and Roof Insulation and New Windows

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock and ComStock models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) upgrade measures, or upgrades. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility upgrade applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on an upgrade package of three end-use savings shapes upgrades - Window Replacement, Exterior Wall Insulation, and Roof Insulation, which we will refer to collectively as the "High Efficiency Envelope" package. More details on the individual upgrades can be found on the ComStock Measures Documentation page. An upgrade package applies two or more EUSS upgrades to a single building model simulation. Since ComStock is a bottom-up physics-based model, an upgrade package will go beyond aggregating or summing the individual upgrade results and produce novel results by simulating interactions between the upgrades. For example, pairing an envelope upgrade with an electrification upgrade would likely result in higher savings results than the sum of these upgrades individually, and the size of the heating, ventilating, and air conditioning (HVAC) equipment may be reduced if the envelope upgrade reduces the loads significantly.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

End-Use Savings Shapes Upgrade Package Documentation: LED Lighting, HP-RTU and ASHP-Boiler

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock and ComStock models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) upgrade measures, or upgrades. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility upgrade applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on an upgrade package of three end-use savings shapes upgrades - Light Emitting Diode (LED) Lighting, Heat Pump Rooftop Unit (RTU) (HP-RTU), and Air-Source Heat Pump (ASHP) Boiler, which we will refer to collectively as the "Interior Lighting and Heat Pump" package. More details on the individual upgrades can be found on the ComStock Measures Documentation page. An upgrade package applies two or more EUSS upgrades to a single building model simulation. Since ComStock is a bottom-up physics-based model, an upgrade package will go beyond aggregating or summing the individual upgrade results and produce novel results by simulating interactions between the upgrades. For example, pairing an envelope upgrade with an electrification upgrade would likely result in higher savings results than the sum of these upgrades individually, and the size of the heating, ventilating, and air conditioning (HVAC) equipment may be reduced if the envelope upgrade reduces the loads significantly.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Central and Eastern North America Ground-Motion Characterization (NGA-East)

This document is the final project report of the Next Generation Attenuation for Central and Eastern North America (CENA) project (NGA-East). The NGA-East objective was to develop a new ground-motion characterization (GMC) model for the CENA region. The GMC model consists of a set of new ground-motion models (GMMs) for median and standard deviation of ground motions and their associated weights to be used with logic-trees in probabilistic seismic hazard analyses (PSHA). NGA-East is a large multidisciplinary project coordinated by the Pacific Earthquake Engineering Research Center (PEER), at the University of California. The project has two components: (1) a set of scientific research tasks, and (2) a model-building component following the framework of the “Seismic Senior Hazard Analysis Committee (SSHAC) Level 3” (Budnitz et al. 1997; NRC 2012). Component (2) is built on the scientific results of component (1) of the NGA-East project. This report documents the tasks under component (2) of the project. Under component (1) of NGA-East, several scientific issues were addressed, including: (a) development of a new database of ground motion data recorded in CENA; (b) development of a regionalized ground-motion map for CENA, (c) definition of the reference site condition; (d) simulations of ground motions based on different methodologies; and (e) development of numerous GMMs for CENA. The scientific tasks of NGA-East were all documented as a series of PEER reports. The scope of component (2) of NGA-East was to develop the complete GMC. This component was designed as a SSHAC Level 3 study with the goal of capturing the ground motions’ center, body, and range of the technically defensible interpretations in light of the available data and models. The SSHAC process involves four key tasks: evaluation, integration, formal review by the Participatory Peer Review Panel (PPRP), and documentation (this report). Key tasks documented in this report include review and evaluation of the empirical groundmotion database, the regionalization of ground motions, and screening sets of candidate GMMs. These are followed by the development of new median and standard deviation GMMs, the development of new analyses tools for quantifying the epistemic uncertainty in ground motions, and the documentation of implementation guidelines of the complete GMC for PSHA computations. Appendices include further documentation of the relevant SSHAC process and additional supporting technical documentation of numerous sensitivity analyses results. The PEER reports documenting component (1) of NGA-East are also considered “attachments” to the current report and are all available online on the PEER website (https://peer.berkeley.edu/). The final NGA-East GMC model includes a set of 17 GMMs defined for 24 ground-motion intensity measures, applicable to CENA in the moment magnitude range of 4.0 to 8.2 and covering distances up to 1500 km. Standard deviation models are also provided for site-specific analysis (single-station standard deviation) and for general PSHA applications (ergodic standard deviation). Adjustment factors are provided for consideration of source-depth effects and hanging-wall effects, as well as for hazard computations at sites in the Gulf Coast region.

58 GEOSCIENCES↗

Modelling the Impact of Fractionation on Late Urinary Toxicity After Postprostatectomy Radiation Therapy

Purpose: To fit urinary toxicity data of patients treated with postprostatectomy radiation therapy with the linear quadratic (LQ) model with/without introducing a time factor. Methods and Materials: Between 1993 and 2010, 1176 patients were treated with conventional fractionation (1.8 Gy per fraction, median 70.2 Gy, n=929) or hypofractionation (2.35-2.90 Gy per fraction, n=247). Data referred to 2004-2010 (when all schemes were in use, n=563; conventional fractionation: 316; hypofractionation: 247) were fitted as a logit function of biological equivalent dose (BED), according to the LQ model with/without including a time factor γ (fixing α/β = 5 Gy). The 3-year risks of severe urethral stenosis, incontinence, and hematuria were considered as endpoints. Best-fit parameters were derived, and the resulting BEDs were taken in multivariable backward logistic models, including relevant clinical variables, considering the whole population. Results: The 3-year incidences of severe stenosis, incontinence, and hematuria were, respectively, 6.6%, 4.8%, and 3.3% in the group treated in 2004-2010. The best-fitted α/β values were 0.81 Gy and 0.74 Gy for incontinence and hematuria, respectively, with the classic LQ formula. When fixing α/β = 5 Gy, best-fit values for γ were, respectively, 0.66 Gy/d and 0.85 Gy/d. Sensitivity analyses showed reasonable values for γ (0.6-1.0 Gy/d), with comparable goodness of fit for α/β values between 3.5 and 6.5 Gy. Likelihood ratio tests showed that the fits with/without including γ were equivalent. The resulting multivariable backward logistic models in the whole population included BED, pT4, and use of antihypertensives (area under the curve [AUC] = 0.72) for incontinence and BED, pT4, and year of surgery (AUC = 0.80) for hematuria. Stenosis data could not be fitted: a 4-variable model including only clinical factors (acute urinary toxicity, pT4, year of surgery, and use of antihypertensives) was suggested (AUC = 0.73). Conclusions: The unexpected impact of moderate hypofractionation on severe incontinence and hematuria after postprostatectomy radiation therapy may be explained by a bladder α/β value <1 Gy or, radiobiologically more plausible, by introducing a time factor likely to represent a previously hypothesized consequential component of late effect.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

{sup 18}F-Choline Positron Emission Tomography/Computed Tomography and Multiparametric Magnetic Resonance Imaging for the Detection of Early Local Recurrence of Prostate Cancer Initially Treated by Radiation Therapy: Comparison With Systematic 3-Dimensional Transperineal Mapping Biopsy

Purpose: To compare the diagnostic performance of {sup 18}F-fluorocholine positron emission tomography/computed tomography (FCH-PET/CT), multiparametric prostate magnetic resonance imaging (mpMRI), and a combination of both techniques for the detection of local recurrence of prostate cancer initially treated by radiation therapy. Methods and Materials: This was a retrospective, single-institution study of 32 patients with suspected prostate cancer recurrence who underwent both FCH-PET/CT and 3T mpMRI within 3 months of one another for the detection of recurrence. All included patients had to be cleared for metastatic recurrence. The reference procedure was systematic 3-dimensional (3D)-transperineal prostate biopsy for the final assessment of local recurrence. Both imaging modalities were analyzed by 2 experienced readers blinded to clinical data. The analysis was made per-patient and per-segment using a 4-segment model. Results: The median prostate-specific antigen value at the time of imaging was 2.92 ng/mL. The mean prostate-specific antigen doubling time was 14 months. Of the 32 patients, 31 had a positive 3D-transperineal mapping biopsy for a local relapse. On a patient-based analysis, the detection rate was 71% (22 of 31) for mpMRI and 74% (23 of 31) for FCH-PET/CT. On a segment-based analysis, the sensitivity and specificity were, respectively, 32% and 87% for mpMRI, 34% and 87% for FCH-PET/CT, and 43% and 83% for the combined analysis of both techniques. Accuracy was 64%, 65%, and 66%, respectively. The interobserver agreement was κ = 0.92 for FCH-PET/CT and κ = 0.74 for mpMRI. Conclusions: Both mpMRI and FCH-PET/CT show limited sensitivity but good specificity for the detection of local cancer recurrence after radiation therapy, when compared with 3D-transperineal mapping biopsy. Prostate biopsy still seems to be mandatory to diagnose local relapse and select patients who could benefit from local salvage therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-F-T-31: Shape and Isodose Distributions in Co60 HDR Brachytherapy for Different Utero-Vaginal Time Ratios

Purpose: To optimize the dose in bladder and rectum and show the different shapes of the isodose volumes in Co60-HDR brachytherapy, considering different utero and vaginal sources dwell ratio times (TU:TV). Methods: Besides Ir192-HDR, new Co60-HDR sources are being incorporated. We considered different TU:TV times and computed the dosis in bladder, rectum and at the reference points of the Manchester system. Also, we calculated the isodose volume and shape in each case. We used a EZAG-BEBIG Co0.A86 model with TPS HDRplus3.0.4. and LCT42-7, LCT42-2(R,L) applicators. A reference dose RA= 1.00 Gy was given to the A-right point. We considered the TU:TV dwell time ratios 1:0.25, 1:0.33, 1:0.5, 1:1, 1:2, 1:3, and 1:4. Given TU:TV, the stop time at each dwell position is fixed for each applicator. Results: Increasing TU:TV systematically results in a decreasing of the dose in bladder and rectum, e.g. 9% and 27% reduction were found in 1:0.25 with respect to 1:1, while 12% and 34% increase were found in 1:4 with respect to 1:1. Also, the isodose volume parameters height (h), width (w), thickness (t) and volume (hwt) increased from the 1:0.25 case to the 1:4 value: hwt is 25% lower and 31% higher than the 1:1 reference volume in these cases. Also w decreased for higher TU:TV and may compromise the tumoral volume coverage, decreasing 17% in the 1:0.25 case compared to the 1:1 case. The shape of the isodose volume was obtained for the different TU:TV considered. Conclusion: We obtained the shape of isodose volumes for different TU:TV values in gynecological Co60-HDR. We studied the dose reduction in bladder and rectum for different TU:TV ratios. The volume parameters and hwt are strongly dependent on this ratio. This information is useful for a quantitative check of the TPS and as a starting point towards optimization.

60 APPLIED LIFE SCIENCES↗

Nonlinear stability and control study of highly maneuverable high performance aircraft

This project is intended to research and develop new nonlinear methodologies for the control and stability analysis of high-performance, high angle-of-attack aircraft such as HARV (F18). Past research (reported in our Phase 1, 2, and 3 progress reports) is summarized and more details of final Phase 3 research is provided. While research emphasis is on nonlinear control, other tasks such as associated model development, system identification, stability analysis, and simulation are performed in some detail as well. An overview of various models that were investigated for different purposes such as an approximate model reference for control adaptation, as well as another model for accurate rigid-body longitudinal motion is provided. Only a very cursory analysis was made relative to type 8 (flexible body dynamics). Standard nonlinear longitudinal airframe dynamics (type 7) with the available modified F18 stability derivatives, thrust vectoring, actuator dynamics, and control constraints are utilized for simulated flight evaluation of derived controller performance in all cases studied.

Mohler, R. R.↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 3: Illustrative test problems

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 3 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 3 explains the kinetics and kinetics-plus-sensitivity analysis problems supplied with LSENS and presents sample results. These problems illustrate the various capabilities of, and reaction models that can be solved by, the code and may provide a convenient starting point for the user to construct the problem data file required to execute LSENS. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Bittker, David A.↗

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗