SEARCH · Engineering Papers

Results for “model weights”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Collapse analysis of a waffle plate strongback for Space Station Freedom

The structural integrity was studied of the Integrated Equipment Assembly (IEA) Strongback of the SSF for the launch environment. The strongback structure supports the electrical power system for SSF. To achieve minimum launch mass, it is essential that flight structures are designed as light as possible. A nonlinear structural analysis was conducted to determine the collapse load of the structure and the associated factor of safety against the service loads. A modeling technique is provided for simulating the load conditions and the buckling and post buckling (collapse) load is evaluated of the IEA Strongback structure, using the finite element computer code MARC. Two of four strongback panels were modeled and analyzed. The effects were examined of the following factors on the global behavior of the strongback panels: (1) load simplification and simulation; (2) type of support boundary conditions; and (3) the possibility of weight reduction of the original structure. For this purpose, several models of the two panels of the strongback were considered. The stress level and distribution in the panels for launch condition, the Eigenvalue critical buckling load and/or the collapse load were determined.

Monasa, Frank F.↗

NASA NTRS2018

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗

DOE OSTI2012

Predicting Chest Wall Pain From Lung Stereotactic Body Radiotherapy for Different Fractionation Schemes

Purpose: Recent studies with two fractionation schemes predicted that the volume of chest wall receiving >30 Gy (V30) correlated with chest wall pain after stereotactic body radiation therapy (SBRT) to the lung. This study developed a predictive model of chest wall pain incorporating radiobiologic effects, using clinical data from four distinct SBRT fractionation schemes. Methods and Materials: 102 SBRT patients were treated with four different fractionations: 60 Gy in three fractions, 50 Gy in five fractions, 48 Gy in four fractions, and 50 Gy in 10 fractions. To account for radiobiologic effects, a modified equivalent uniform dose (mEUD) model calculated the dose to the chest wall with volume weighting. For comparison, V30 and maximum point dose were also reported. Using univariable logistic regression, the association of radiation dose and clinical variables with chest wall pain was assessed by uncertainty coefficient (U) and C statistic (C) of receiver operator curve. The significant associations from the univariable model were verified with a multivariable model. Results: 106 lesions in 102 patients with a mean age of 72 were included, with a mean of 25.5 (range, 12-55) months of follow-up. Twenty patients reported chest wall pain at a mean time of 8.1 (95% confidence interval, 6.3-9.8) months after treatment. The mEUD models, V30, and maximum point dose were significant predictors of chest wall pain (p < 0.0005). mEUD improved prediction of chest wall pain compared with V30 (C = 0.79 vs. 0.77 and U = 0.16 vs. 0.11). The mEUD with moderate weighting (a = 5) better predicted chest wall pain than did mEUD without weighting (a = 1) (C = 0.79 vs. 0.77 and U = 0.16 vs. 0.14). Body mass index (BMI) was significantly associated with chest wall pain (p = 0.008). On multivariable analysis, mEUD and BMI remained significant predictors of chest wall pain (p = 0.0003 and 0.03, respectively). Conclusion: mEUD with moderate weighting better predicted chest wall pain than did V30, indicating that a small chest wall volume receiving a high radiation dose is responsible for chest wall pain. Independently of dose to the chest wall, BMI also correlated with chest wall pain.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2019

Frit Recommendation For Sludge Batch 9 Coupled Operation With The Salt Waste Processing Facility

Prior to the processing of Sludge Batch 9 (SB9) with Frit 803 at the Defense Waste Processing Facility (DWPF), assessments and experimental glass work were completed to support SB9 coupled processing with streams from the Actinide Removal Process (ARP) and Modular Caustic Side Solvent Extraction Unit (MCU). This previous work did not address coupled processing with streams from the Salt Waste Processing Facility (SWPF) because the Product Composition Control System (PCCS) models limited the TiO2 concentration in glass to 2 weight percent (wt.%). These high activity streams from SWPF include monosodium titanate (MST) and sludge solids from the Sludge Solids Receipt Tank (SSRT) as well as Cscontaining strip effluent (SE). The SWPF streams are expected to reach TiO2 concentrations in glass greater than 2 wt.% based on an MST strike of 0.4 g MST/L of salt solution. Therefore, revisions to the models contained in PCCS were completed to allow the evaluation of glasses containing greater than 2 wt.% TiO2, but less than 6 wt.%. Since the initial SB9 qualification for ARP-MCU coupled operation, additional Tank 40 analyses have been completed, which served as the basis for an updated projection of the SB9 blend composition in Tank 40. This projection will be used to support SB9 coupled operation with SWPF in 2019. Savannah River Remediation (SRR) has requested that the Savannah River National Laboratory (SRNL) select a second frit to optimize SB9 coupled operation with SWPF. The selected frit would provide operational flexibility by allowing an acceptable waste loading (WL) to be targeted while maintaining an operating window of at least ± 4 percentage points around this target. The objective of this study was to select a second frit based on a paper-study utilizing Measurement Acceptance Region (MAR) assessments with DWPF PCCS. To support these assessments, SRNL developed a new method for rigorous and directed optimization of frit compositions, which allowed for rapid screening and an assessment of a large array of candidate frits, and the ability to simultaneously satisfy the following single and double MST strike operating scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2013

Eulerian adaptive finite-difference method for high-velocity impact and penetration problems

Owing to the complex processes involved, faithful prediction of high-velocity impact events demands a simulation method delivering efficient calculations based on comprehensively formulated constitutive models. Such an approach is presented herein, employing a weighted essentially non-oscillatory (WENO) method within an adaptive mesh refinement (AMR) framework for the numerical solution of hyperbolic partial differential equations. Applied widely in computational fluid dynamics, these methods are well suited to the involved locally non-smooth finite deformations, circumventing any requirement for artificial viscosity functions for shock capturing. Application of the methods is facilitated through using a model of solid dynamics based upon hyper-elastic theory comprising kinematic evolution equations for the elastic distortion tensor. The model for finite inelastic deformations is phenomenologically equivalent to Maxwell s model of tangential stress relaxation. Closure relations tailored to the expected high-pressure states are proposed and calibrated for the materials of interest. Sharp interface resolution is achieved by employing level-set functions to track boundary motion, along with a ghost material method to capture the necessary internal boundary conditions for material interactions and stress-free surfaces. The approach is demonstrated for the simulation of high velocity impacts of steel projectiles on aluminium target plates in two and three dimensions.

Barton, Philip T.↗

DOE OSTI2026

ReMU: regional minimal updating for model-based derivative-free optimization

Derivative-free optimization (DFO) problems are optimization problems where derivative information is unavailable or extremely difficult to obtain. Model-based DFO solvers have been applied extensively in scientific computing. Powell's NEWUOA (2004) [Powell, The NEWUOA software for unconstrained optimization without derivatives, in Large-Scale Nonlinear Optimization, Nonconvex Optimization and its Applications Vol. 83, G. Di Pillo and M. Roma, eds., Springer, 2006, pp. 255–297] and Wild's POUNDerS (2014) [Wild, Solving derivative-free nonlinear least squares problems with POUNDERS, in Advances and Trends in Optimization with Engineering Applications, T. Terlaky, M.F. Anjos, and S. Ahmed, eds., SIAM, 2017, pp. 529–540] explore the numerical power of the minimal norm Hessian (MNH) model for DFO and contributed to the open discussion on building better models with fewer data to achieve faster numerical convergence. Another decade later, we propose the regional minimal updating (ReMU) models, and extend the previous models into a broader class, including the H 2 norm models [Xie and Yuan, Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms, IMA J. Numer. Anal. 46 (2025), pp. 21–50]. This paper shows motivation behind ReMU models, computational details, theoretical and numerical results on particular extreme points and the barycentre of ReMU's weight coefficient region, and the associated KKT matrix error and distance. Novel metrics, such as the truncated Newton step error, are proposed to numerically understand the new models' properties. A new algorithmic strategy, based on iteratively adjusting the ReMU model type, is also proposed, and shows numerical advantages by combining and switching between the barycentric model and the classic least Frobenius norm model in an online fashion.

derivative-free trust-region methods↗

NASA NTRS1993

Cost model relationships between textile manufacturing processes and design details for transport fuselage elements

Textile manufacturing processes offer potential cost and weight advantages over traditional composite materials and processes for transport fuselage elements. In the current study, design cost modeling relationships between textile processes and element design details were developed. Such relationships are expected to help future aircraft designers to make timely decisions on the effect of design details and overall configurations on textile fabrication costs. The fundamental advantage of a design cost model is to insure that the element design is cost effective for the intended process. Trade studies on the effects of processing parameters also help to optimize the manufacturing steps for a particular structural element. Two methods of analyzing design detail/process cost relationships developed for the design cost model were pursued in the current study. The first makes use of existing databases and alternative cost modeling methods (e.g. detailed estimating). The second compares design cost model predictions with data collected during the fabrication of seven foot circumferential frames for ATCAS crown test panels. The process used in this case involves 2D dry braiding and resin transfer molding of curved 'J' cross section frame members having design details characteristic of the baseline ATCAS crown design.

Metschan, Stephen L.↗

NASA NTRS2008

Oxidation of Carbon/Carbon through Coating Cracks

Reinforced carbon/carbon (RCC) is used to protect the wing leading edge and nose cap of the Space Shuttle Orbiter on re-entry. It is composed of a lay-up of carbon/carbon fabric protected by a SiC conversion coating. Due to the thermal expansion mismatch of the carbon/carbon and the SiC, the SiC cracks on cool-down from the processing temperature. The cracks act as pathways for oxidation of the carbon/carbon. A model for the diffusion controlled oxidation of carbon/carbon through machined slots and cracks is developed and compared to laboratory experiments. A symmetric cylindrical oxidation cavity develops under the slots, confirming diffusion control. Comparison of cross sectional dimensions as a function of oxidation time shows good agreement with the model. A second set of oxidation experiments was done with samples with only the natural craze cracks, using weight loss as an index of oxidation. The agreement of these rates with the model is quite reasonab

Jacobson, N. S.↗

DOE OSTI2021

Dark Energy Survey year 3 results: covariance modelling and its impact on parameter estimation and quality of fit

ABSTRACT We describe and test the fiducial covariance matrix model for the combined two-point function analysis of the Dark Energy Survey Year 3 (DES-Y3) data set. Using a variety of new ansatzes for covariance modelling and testing, we validate the assumptions and approximations of this model. These include the assumption of Gaussian likelihood, the trispectrum contribution to the covariance, the impact of evaluating the model at a wrong set of parameters, the impact of masking and survey geometry, deviations from Poissonian shot noise, galaxy weighting schemes, and other sub-dominant effects. We find that our covariance model is robust and that its approximations have little impact on goodness of fit and parameter estimation. The largest impact on best-fitting figure-of-merit arises from the so-called fsky approximation for dealing with finite survey area, which on average increases the χ2 between maximum posterior model and measurement by $3.7{{\ \rm per\ cent}}$ (Δχ2 ≈ 18.9). Standard methods to go beyond this approximation fail for DES-Y3, but we derive an approximate scheme to deal with these features. For parameter estimation, our ignorance of the exact parameters at which to evaluate our covariance model causes the dominant effect. We find that it increases the scatter of maximum posterior values for Ωm and σ8 by about $3{{\ \rm per\ cent}}$ and for the dark energy equation-of-state parameter by about $5{{\ \rm per\ cent}}$.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1979

Effect of lip and centerbody geometry on aerodynamic performance of inlets for tilting-nacelle VTOL aircraft

Inlets for tilt-nacelle VTOL aircraft must operate over a wide range of incidence angles and engine weight flows without internal flow separation. Wind tunnel tests of scale model inlets were conducted to evaluate the effectiveness of three geometric variables to provide this capability. Increasing the lip contraction ratio increased the separation angle at all engine weight flows. The optimum axial location of the centerbody occurred when its leading edge was located just downstream of the inlet lip. Compared with a short centerbody, the optimum location of the centerbody resulted in an increase in separation angle at all engine weight flows. Decreasing the lip major-to-minor-axis ratio increased the separation angle at the lower engine weight flows.

Burley, R. R.↗

NASA NTRS1979

Effect of lip and centerbody geometry on aerodynamic performance of inlets for tilting-nacelle VTOL aircraft

Burley, R. R.↗

DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2017

Aeroservoelastic Modeling of Body Freedom Flutter for Control System Design

One of the most severe forms of coupling between aeroelasticity and flight dynamics is an instability called freedom flutter. The existing tools often assume relatively weak coupling, and are therefore unable to accurately model body freedom flutter. Because the existing tools were developed from traditional flutter analysis models, inconsistencies in the final models are not compatible with control system design tools. To resolve these issues, a number of small, but significant changes have been made to the existing approaches. A frequency domain transformation is used with the unsteady aerodynamics to ensure a more physically consistent stability axis rational function approximation of the unsteady aerodynamic model. The aerodynamic model is augmented with additional terms to account for limitations of the baseline unsteady aerodynamic model and to account for the gravity forces. An assumed modes method is used for the structural model to ensure a consistent definition of the aircraft states across the flight envelope. The X-56A stiff wing flight-test data were used to validate the current modeling approach. The flight-test data does not show body-freedom flutter, but does show coupling between the flight dynamics and the aeroelastic dynamics and the effects of the fuel weight.

aeroservoelastic↗

DOE OSTI2024

Mitigation of distortion of Al/steel part under simulated paint baking condition: Experiment and numerical model studies

Multi-material joining of lightweight structures is essential to reduce vehicle weight for more energy savings and less greenhouse gas emission. However, mismatch of thermal expansion coefficient for dissimilar materials during the paint baking process can induce part distortion and joint failure for adhesive bonding. Here, in the present work, a thermomechanical model based on contact mechanics and large deformation theory was developed for dissimilar high-strength Al alloy and steel components to study the distortion mechanism and influential factors of the residual gap. The established model was used to optimize joint conditions, such as pitch distance and part geometry. When a weld pitch is shorter than 100 mm, the maximum gap between Al and steel part can be greatly reduced to 0.1 mm, and the local stress and plastic strain around the joint during the oven heating and cooling cycle are also substantially reduced compared with the long pitch case (900 mm). The numerical modeling results revealed that a comparable bending stiffness ratio between the steel and Al cross sections is critical to the minimization of gap and distortion under paint baking condition. Digital image correlation technique was used to measure the overall part distortion and local strain distribution that were used to validate the model prediction. Weld bonding (adhesive bonding with friction bit joining) process was successfully employed to join Al to steel component without gap opening in adhesive after the paint baking and cooling.

36 MATERIALS SCIENCE↗

DOE OSTI2023

The star-formation history in the last 10 billion years from CIB cross-correlations

ABSTRACT The cosmic infrared background (CIB) traces the emission of star-forming galaxies throughout all cosmic epochs. Breaking down the contribution from galaxies at different redshifts to the observed CIB maps would allow us to probe the history of star formation. In this paper, we cross-correlate maps of the CIB with galaxy samples covering the range z ≲ 2 to measure the bias-weighted star-formation rate (SFR) density 〈bρSFR〉 as a function of time in a model independent way. This quantity is complementary to direct measurements of the SFR density ρSFR, giving a higher weight to more massive haloes, and thus provides additional information to constrain the physical properties of star formation. Using cross-correlations of the CIB with galaxies from the DESI Legacy Survey and the extended Baryon Oscillation Spectroscopic Survey, we obtain high signal-to-noise ratio measurements of 〈bρSFR〉, which we then use to place constraints on halo-based models of the star-formation history. We fit halo-based SFR models to our data and compare the recovered ρSFR with direct measurements of this quantity. We find a qualitatively good agreement between both independent data sets, although the details depend on the specific halo model assumed. This constitutes a useful robustness test for the physical interpretation of the CIB, and reinforces the role of CIB maps as valuable astrophysical probes of the large-scale structure. We report our measurements of 〈bρSFR〉 as well as a thorough account of their statistical uncertainties, which can be used to constrain star-formation models in combination with other data.

Astronomy & Astrophysics↗

DOE OSTI2012

SEARCH FOR DARK MATTER ANNIHILATION SIGNALS FROM THE FORNAX GALAXY CLUSTER WITH H.E.S.S

The Fornax galaxy cluster was observed with the High Energy Stereoscopic System for a total live time of 14.5 hr, searching for very high energy (VHE; E > 100GeV) {gamma}-rays from dark matter (DM) annihilation. No significant signal was found in searches for point-like and extended emissions. Using several models of the DM density distribution, upper limits on the DM velocity-weighted annihilation cross-section ({sigma}v) as a function of the DM particle mass are derived. Constraints are derived for different DM particle models, such as those arising from Kaluza-Klein and supersymmetric models. Various annihilation final states are considered. Possible enhancements of the DM annihilation {gamma}-ray flux, due to DM substructures of the DM host halo, or from the Sommerfeld effect, are studied. Additional {gamma}-ray contributions from internal bremsstrahlung and inverse Compton radiation are also discussed. For a DM particle mass of 1 TeV, the exclusion limits at 95% of confidence level reach values of ({sigma}v){sup 95%C.L.} {approx} 10{sup -23} cm{sup 3} s{sup -1}, depending on the DM particle model and halo properties. Additional contribution from DM substructures can improve the upper limits on ({sigma}v) by more than two orders of magnitude. At masses around 4.5 TeV, the enhancement by substructures and the Sommerfeld resonance effect results in a velocity-weighted annihilation cross-section upper limit at the level of ({sigma}v){sup 95%C.L.} {approx}10{sup -26} cm{sup 3} s{sup -1}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

A Grid Intensity Map-Based Source Model for Monte Carlo Modeling of IMRT Field Irradiation - Paper 5

An accurate source model is important for TPS (Treatment Planning System) dose calculation and verification, especially for IMRT (Intensity Modulated Radiotherapy) irregular and inhomogeneous fields. Monte Carlo simulation is the most accurate way to calculate dose distributions in the heterogeneous media like human bodies. A full radiotherapy modeling generally includes Linac modeling and phantom dose calculation. Full Linac modeling is very time-consuming and requires highly detailed knowledge of the accelerator composition and structure, however, which is very difficult to obtain for commercial interests. Often all relevant parameters are not known with sufficient accuracy, and thus a time-consuming tuning of e.g. the electron beam spot size and energy must be performed. A multiple source model (MSM) regards particles coming from different parts of an accelerator as if they are coming from different virtual sources. MSM uses a set of virtual sources representing the treatment beam (e.g. primary photons, head scatter, electron contamination) by a set of analytic functions or histograms with parameters derived from dosimetric measurements or a full MC simulation. Tillikainen et al built a three sub-source model by matching with measurement data using an automatic process to optimize the free model parameters for a Varian 2100 C/D accelerator. Lin et al have built a convenient five-photon-sub-source model for Varian 2300C 6 MV-X beam based on Monte Carlo simulated Linac for ARTS (Accurate Radiotherapy System). The particle generation of MSM is very fast compared to full Linac simulation, and its sampling efficiency is not limited by the size of the phase space file (PSF). However, generally MSM is generated based on regular field PSF and needs an irregular component to form the beam shape meanwhile apt to bringing error Linac modeling and MSM producing from PSF analysis actually both are multi-freedom optimization process by adjusting several parameters (e.g. electron energy and angular divergence). Linac modeling adjusts its parameters from the beginning of electron hitting target through all components to fit the simulation to the measurement. MSM tunes its PSF-derived parameters by optimizing to minimize the deviations between the simulation and the measurement. Since the full Linac modeling and the MSM model both need to fine tune their parameters to match the simulation to the measurement, and sometimes maybe offsets the real accelerator for its multi-parameter optimization process. So, why not fine tune directly all source model parameters just behind all Linac components to fit the simulation to the measurement? A grid intensity-based surface source model was suggested here, which assumed the particle emitting from the bottom surface of Linac last beam trimmer MLC (Multi Leaf Collimator). The particle emitting position is sampled on the grid intensity map optimized by TPS, and the particle weight is adjusted by its original grid intensity. This source model not only takes advantage of the Monte Carlo method in the heterogeneous media, but also avoids effectively the transportation speed shortcoming of Monte Carlo simulation through Linac components. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

NASA NTRS1987

Robust time and frequency domain estimation methods in adaptive control

A robust identification method was developed for use in an adaptive control system. The type of estimator is called the robust estimator, since it is robust to the effects of both unmodeled dynamics and an unmeasurable disturbance. The development of the robust estimator was motivated by a need to provide guarantees in the identification part of an adaptive controller. To enable the design of a robust control system, a nominal model as well as a frequency-domain bounding function on the modeling uncertainty associated with this nominal model must be provided. Two estimation methods are presented for finding parameter estimates, and, hence, a nominal model. One of these methods is based on the well developed field of time-domain parameter estimation. In a second method of finding parameter estimates, a type of weighted least-squares fitting to a frequency-domain estimated model is used. The frequency-domain estimator is shown to perform better, in general, than the time-domain parameter estimator. In addition, a methodology for finding a frequency-domain bounding function on the disturbance is used to compute a frequency-domain bounding function on the additive modeling error due to the effects of the disturbance and the use of finite-length data. The performance of the robust estimator in both open-loop and closed-loop situations is examined through the use of simulations.

Lamaire, Richard Orville↗