Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large deviations”

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 55 records · Page 3

SU-F-J-21: Clinical Evaluation of Surface Scanning Systems in Different Treatment Locations

Purpose: To reduce imaging dose in fractionated IGRT, the ability of optical surface imaging systems (OSIS) to detect setup errors was tested. Therefore, clinical studies to evaluate for different treatment locations setup corrections derived by OSIS in comparison with x-ray image guidance in fractionated radiation therapy was performed. Methods: The setup correction accuracy of an OSIS system (AlignRT, VisionRT, London, UK) will be analysed for the 4 tumour locations Pelvis, Upper Abdomen, Thorax and Breast, 20 patients for each location in comparison to a different system (Sentinel, C-RAD, SE). For each patient, the setup corrections of the cone-beam computed tomography (CBCT) of an Elekta Versa HD linear accelerator (Elekta, Crawley, UK) is considered as gold-standard and then compared with those of the OSIS for the first ten fractions retrospectively. There were no clinical decisions made based on the surrogate system. For the OSIS, the reference surface is highly important as it represents the actual ground truth. It can be obtained either with the system itself or the surface structure delineated in the planning CT can be imported via DICOM interface. In this paper, the first results for the treatment region thorax are presented. The reference image modalities were compared. Results: Table 1 displays the difference between the setup corrections obtained with OSIS and CBCT in lateral (LAT), longitudinal (LNG) and vertical (VRT) direction for the DICOM reference image. While the median deviations are within a few millimeters, some outliers showed large deviations. Generally, the mean deviation as well as the spread was smallest in lateral and largest in vertical direction. Conclusion: Although the system allows fast, simple and non-invasive determination of setup corrections, it should be evaluated treatment region dependant. Therefore, the study is ongoing. The application of OSIS may help to reduce the imaging dose for the patient. We gratefully acknowledge the support by our colleagues from the Radiological University Clinic Heidelberg, where the study was performed. This work was funded by the Federal Ministry of Education and Research (BMBF) Germany, grant number 01IB13001B.

60 APPLIED LIFE SCIENCES↗

Incorporating Active Runway Crossings in Airport Departure Scheduling

A mixed integer linear program is presented for deterministically scheduling departure and ar rival aircraft at airport runways. This method addresses different schemes of managing the departure queuing area by treating it as first-in-first-out queues or as a simple par king area where any available aircraft can take-off ir respective of its relative sequence with others. In addition, this method explicitly considers separation criteria between successive aircraft and also incorporates an optional prioritization scheme using time windows. Multiple objectives pertaining to throughput and system delay are used independently. Results indicate improvement over a basic first-come-first-serve rule in both system delay and throughput. Minimizing system delay results in small deviations from optimal throughput, whereas minimizing throughput results in large deviations in system delay. Enhancements for computational efficiency are also presented in the form of reformulating certain constraints and defining additional inequalities for better bounds.

Gupta, Gautam↗

The stochastic variability of asteroidal regolith depths

Modeling the depth of regolith on asteroids is approached from a statistical point of view. It is demonstrated that average values are not good descriptors of regolith depth on asteroids. Large deviations from the average can be expected to occur due to both large variations in depth over the surface of a body and to the fact that each asteroid has a unique regolith. The utility of the average depth is not significantly increased by excluding the parts of a surface which are occupied by large craters; a procedure adopted in existing regolith models. Although an asteroid's surface may be 'smoothed out' by movement of debris into gravitationally low spots, the regolith depth retains its variability because of variations in topography at the bottom of the regolith layer. The large variability associated with regolith depth severely limits the power of regolith models in predicting parent body size for the brecciated meteorites.

Housen, K. R.↗

Speed of Sound Measurements of R-1130(E) and an Azeotropic Blend of R-1336mzz(Z)/1130(E)

Sound speed data measured using a dual-path pulse-echo instrument are reported for pure trans-1,2-dichloroethene (R-1130(E)) and an azeotropic blend of cis-1,1,1,4,4,4-hexafluorobutene (R-1336mzz(Z)) and R-1130(E) with a composition of 74.8 mass % R-1336mzz(Z) with the balance being R-1130(E). The azeotropic blend of R-1336mzz(Z)/1130(E) is classified as R-514A in ANSI/ASHRAE standard 34. Liquid phase speed of sound data are reported from just above the saturation pressure of pure R-1130(E) or the bubble point pressure of R-514A to a maximum pressure of 26.7 MPa. The relative combined expanded uncertainty in the speed of sound varies from 0.032 % to 0.148 % with the greatest deviations occurring at the lowest sound speeds. At present, no reference Helmholtz-energy-explicit equation of state (EOS) is available for R-1130(E). Therefore, the reported data for pure R-1130(E) are compared to an extended corresponding states (ECS) model. Deviations between the pure R-1130(E) sound speed data and the ECS model were found to be consistently negative ranging between − 4.1 % and − 3.5 %. The R-514A data are compared to a multifluid model inclusive of the established reference Helmholtz-energy-explicit EOS for R-1336mzz(Z) and ECS model for R-1130(E) with estimated binary interaction parameters. Deviations between the experimental speed of sound data and the multifluid model were also found to be consistently negative. However, deviations from the multifluid model were found to be as great as − 17.1 %. The large deviations from the ECS model and multifluid model underscore the need for a robust Helmholtz-energy-explicit EOS for R-1130(E).

42 ENGINEERING↗

Gravitational Test beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

In this study, the 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-643: Pure Alanine Dosimeter for Verification Dosimetry in IMRT

Purpose: The objective of this study was evaluation of accuracy of pure alanine dosimeters measuring intensity-modulated radiation therapy (IMRT) dose distributions in a thorax phantom. Methods: Alanine dosimeters were prepared in the form of 110 mg pure L-α-alanine powder filled into clear tissue-equivalent polymethylmethacrylate (PMMA) plastic tubes with the dimensions 25 mm length, 3 mm inner diameter, and 1 mm wall thickness. A dose-response calibration curve was established for the alanine by placing the dosimeters at 1.5 cm depth in a 30×30×30 cm{sup 3} solid water phantom and then irradiating on a linac with 6 MV photon beam at 10×10 cm{sup 2} field size to doses ranging from 1 to 5 Gy. Electron paramagnetic resonance (EPR) spectroscopy was used to determine the absorbed dose in alanine. An IMRT treatment plan was designed for a commercial heterogeneous CIRS thorax phantom and the dose values were calculated at three different points located in tissue, lung, and bone equivalent materials. A set of dose measurements was carried out to compare measured and calculated dose values by placing the alanine dosimeters at those selected locations inside the thorax phantom and delivering the IMRT to the phantom. Results: The alanine dose measurements and the IMRT plan dose calculations were found to be in agreement within ±2%. Specifically, the deviations were −0.5%, 1.3%, and −1.7% for tissue, lung, and bone; respectively. The slightly large deviations observed for lung and bone may be attributed to tissue inhomogeneity, steep dose gradients in these regions, and uncontrollable changes in spectrometer conditions. Conclusion: The results described herein confirmed that pure alanine dosimeter was suitable for in-phantom dosimetry of IMRT beams because of its high sensitivity and acceptable accuracy. This makes the dosimeter a promising option for quality control of the therapeutic beams, complementing the commonly used ionization chambers, TLDs, and films.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Computational Assessment of the Opening Mode Facesheet-to-Core Fracture Toughness of a Thin-Gage Facesheet Sandwich Structure

The single-cantilevered beam test (SCB) is currently being considered for standardization for measuring the opening mode fracture toughness of sandwich structures subjected to facesheet-to-core separation. At present, the Area method is considered the benchmark method for determining energy required to propagate the crack. This method is problematic since it requires periodic unloading and cannot be evaluated instantaneously. Therefore, a method is needed to accurately and instantaneously calculate opening mode fracture toughness. In previous work, SCB tests with thin-gauge facesheets were performed to compare the J-integral approach to data reduction methods including the Modified Beam Theory (MBT) and the Area method. Large deviations (over 10%) in the calculated fracture toughness were observed between the J-integral and MBT for large displacements. In this study, a simulated SCB experiment using a virtual crackclosure technique was performed to assess the previously observed discrepancies. The predicted load displacement response, the predicted angle at the load application point, and the predicted crack growth results are in good agreement (within 5%) with the corresponding experimental data. Using the predicted load and angles, the predicted J-integral is in good agreement (within 10%) with the experimentally obtained solution. Using the MBT approach, the previously observed deviations in the fracture toughness were determined to be associated with the empirical extrapolation procedure to estimate the crack-root rotation correction from nonlinear compliance curves. The nonlinear compliance is determined to be associated with large displacement effects observed within the predicted and experimental results.

Sandwich Composite↗

MO-F-CAMPUS-T-03: Data Driven Approaches for Determination of Treatment Table Tolerance Values for Record and Verification Systems

Purpose: To determine appropriate couch tolerance values for modern radiotherapy linac R&V systems with indexed patient setup. Methods: Treatment table tolerance values have been the most difficult to lower, due to many factors including variations in patient positioning and differences in table tops between machines. We recently installed nine linacs with similar tables and started indexing every patient in our clinic. In this study we queried our R&V database and analyzed the deviation of couch position values from the acquired values at verification simulation for all patients treated with indexed positioning. Mean and standard deviations of daily setup deviations were computed in the longitudinal, lateral and vertical direction for 343 patient plans. The mean, median and standard error of the standard deviations across the whole patient population and for some disease sites were computed to determine tolerance values. Results: The plot of our couch deviation values showed a gaussian distribution, with some small deviations, corresponding to setup uncertainties on non-imaging days, and SRS/SRT/SBRT patients, as well as some large deviations which were spot checked and found to be corresponding to indexing errors that were overriden. Setting our tolerance values based on the median + 1 standard error resulted in tolerance values of 1cm lateral and longitudinal, and 0.5 cm vertical for all non- SRS/SRT/SBRT cases. Re-analizing the data, we found that about 92% of the treated fractions would be within these tolerance values (ignoring the mis-indexed patients). We also analyzed data for disease site based subpopulations and found no difference in the tolerance values that needed to be used. Conclusion: With the use of automation, auto-setup and other workflow efficiency tools being introduced into radiotherapy workflow, it is very essential to set table tolerances that allow safe treatments, but flag setup errors that need to be reassessed before treatments.

60 APPLIED LIFE SCIENCES↗

Moderate deviations-based importance sampling for stochastic recursive equations

Subsolutions to the Hamilton–Jacobi–Bellman equation associated with a moderate deviations approximation are used to design importance sampling changes of measure for stochastic recursive equations. Analogous to what has been done for large deviations subsolution-based importance sampling, these schemes are shown to be asymptotically optimal under the moderate deviations scaling. Here, we present various implementations and numerical results to contrast their performance, and also discuss the circumstances under which a moderate deviation scaling might be appropriate.

97 MATHEMATICS AND COMPUTING↗

SU-F-T-588: Asymmetries in the CyberKnife Iris 2 Collimator

Purpose: The Iris 2 dodecahedral collimator in the CyberKnife unit can demonstrate visible asymmetry in the shaped radiation field, as seen on routine films. Specifically some edges of the collimator project longer than others. PDDs and output factors remain unaffected by this asymmetry. An in-house program was written to analyze the impact of this asymmetry on profiles. Methods: Two years of routine chromodynamic films were retrospectively analyzed using an in-house developed program. Films were obtained on a weekly basis in a plastic phantom on the vendor-provided mount. Profiles were obtained every degree. A dodecahedron was fitted to the 50% isodose line and the positions and sizes of the upper and lower hexagonal collimators were derived. The collimators’ profiles at every degree, symmetry and centering were calculated and trended over the months. Results: Asymmetries in the obtained profiles were apparent. The upper and lower collimators were off-centered (0.31 mm, p<10–8), they did not project the same field size at 80 cm (0.46 mm difference, p<10–31). Profiles as a function of angle demonstrate a periodicity, but short and long profiles are rarely 15° apart, but vary from 9° to 21°. In addition, a second, intermediate minimum is seen at every other short profile. A trend in the field size over time is noted and can be monitored to identify servicing needs. Conclusion: Despite its differences, the asymmetric Iris 2 collimator did not demonstrate large deviations in PDDs and output factors. Due to the asymmetry, the 15° apart profiles obtained at commissioning do not represent the true large and short blade profiles of the collimator. Averaging of the profiles in the commissioning process hides this difference.

60 APPLIED LIFE SCIENCES↗

Towards district scale flood simulations using conventional and anisotropic porosity shallow water models with high-resolution topographic information

Current topographic survey technology provides high-resolution (HR) datasets for urban environments. Incorporating this HR information in models aiming to provide flood risk assessment is desirable because the flood wave propagation is depending on the urban topographic features, i.e. buildings, bridges and street networks. Conceptual, numerical and practical challenges arise from the application of shallow water models to HR urban flood modeling. For example, numerical challenges are occurrence of wet-dry fronts, geometric discontinuities in the urban environment and discontinuous solutions, i.e. shock waves. These challenges can be overcome by using a Godunov-type scheme. However, the computational cost of this type of schemes is high, such that HR two-dimensional shallow water simulations with practical relevance have to be run on supercomputers. The porous shallow water model is an alternative approach that aims to reduce computational cost by using a coarse resolution and accounting for unresolved processes by means of the porosity terms. Usually, a speedup between two and three orders of magnitude in comparison to HR simulations can be obtained. This study presents preliminary results of a practical test case concerning pluvial flooding in a district of the city of Nice, France, caused by the intense rainfall event on October 3rd, 2015. HR topography data set on a 1 m resolution is available for the district, whereby street features of infra-metric dimensions have been included. A reference solution is calculated by a HR shallow water model on a 1 m by 1 m structured computational grid. The porous shallow water model is run on a 10 m by 10 m grid and the influence of the drag source term is studied. The model calculations show a large deviation, which is caused by the poor meshing strategy of the porous shallow water (AP) model. The study also summarizes practical challenges that arise during the application of the AP and HR models to a large urban catchment. The main difficulty is to obtain a good mesh. In smaller scale investigations, the mesh is currently constructed by hand such that the cell edges align with buildings. This approach is not feasible for large scale urban catchments with a large number of buildings. Future steps that have to be taken, such as a strategy for automatic mesh generation, are reported on.

42 ENGINEERING↗

TIDAL STIRRING OF SATELLITES WITH SHALLOW DENSITY PROFILES PREVENTS THEM FROM BEING TOO BIG TO FAIL

The “too big to fail” problem is revisited by studying the tidal evolution of populations of dwarf satellites with different density profiles. The high-resolution cosmological ΛCDM “ErisMod” set of simulations is used. These simulations can model both the stellar and dark matter components of the satellites, and their evolution under the action of the tides of a Milky Way (MW)-sized host halo at a force resolution better than 10 pc. The stronger tidal mass loss and re-shaping of the mass distribution induced in satellites with γ = 0.6 dark matter density distributions, as those resulting from the effect of feedback in hydrodynamical simulations of dwarf galaxy formation, are sufficient to bring the circular velocity profiles in agreement with the kinematics of MW’s dSphs. In contrast, in simulations in which the satellites retain cusps at z = 0 there are several “massive failures” with circular velocities in excess of the observational constraints. Various sources of deviations in the conventionally adopted relation between the circular velocity at the half-light radius and the one-dimensional line of sight velocity dispersions are found. Such deviations are caused by the response of circular velocity profiles to tidal effects, which also varies depending on the initially assumed inner density profile and by the complexity of the stellar kinematics, which include residual rotation and anisotropy. In addition, tidal effects naturally induce large deviations in the stellar mass–halo mass relation for halo masses below 10{sup 9} M {sub ⊙}, preventing any reliable application of the abundance matching technique to dwarf galaxy satellites.

79 ASTRONOMY AND ASTROPHYSICS↗

An Hp-Adaptive Minimum Action Method Based on a Posteriori Error Estimate

In this work, we develop an hp-adaptive minimum action method (MAM). MAM plays an important role to minimize the Freidlin-Wentzell action functional, which is the central object of the Freidlin-Wentzell theory of large deviations for transitions induced by small random perturbations. Because of the demanding computation cost, especially in spatially extended systems, numerical efficiency is a critical issue for MAM. Difficulties come from discretizations of both time and space. One hurdle for the application of MAM to large scale systems is the global reparametrization in time direction, which is required by most versions of MAM. We recently introduced a new version of MAM in [25], called tMAM, where we used some simple heuristic criteria to demonstrate that tMAM can be effectively coupled with h-adaptivity. The target of this paper is to integrate hp-adaptivity into tMAM using a posterior error estimation techniques. More specifically, we use the arc length constraint given by geometric minimum action method (gMAM) to define an indicator of the effect of linear time scaling in tMAM, and the derivative recovery technique to construct an error indicator and a regularity indicator for hp-refinement. Numerical results are presented.

Wan, Xiaoliang↗

Experimental and numerical study of a turbulent boundary layer with pressure gradients

Results are presented of an experimental and numerical study of a turbulent boundary layer with pressure gradients conducted using the recent 'fringe method' with its numerical advantages and good inflow quality. After an inflow transient good agreement is observed; the differences, of up to 13 percent, are discussed. Moderate deviations from the law of the wall are found in the velocity profiles of the simulation. They are fully correlated with the pressure gradient, are in fair quantitative agreement with the experimental results of Nagano et al. (1992), and are roughly the opposite of uncorrected mixing-length-model predictions. Large deviations from the wall scaling are observed for other quantities, notably for the turbulence dissipation rate. The a(1) structure parameter drops mildly in the upper layer with adverse pressure gradient.

Spalart, Philippe R.↗

Beam energy dependence of net- Λ fluctuations measured by the STAR experiment at the BNL Relativistic Heavy Ion Collider

The measurements of particle multiplicity distributions have generated considerable interest in understanding the fluctuations of conserved quantum numbers in the quantum chromodynamics (QCD) hadronization regime, in particular near a possible critical point and near the chemical freeze-out. Net-protons and net-kaons have been used as proxies for the net-baryon number and net-strangeness, respectively. We report the measurement of efficiency- and centrality-bin width-corrected cumulant ratios (C 2 /C 1 , C 3 /C 2 ) of net- Λ distributions, in the context of both strangeness and baryon number conservation, as a function of collision energy, centrality, and rapidity. The results are for Au + Au collisions at five beam energies ( √s NN = 19.6 , 27, 39, 62.4, and 200 GeV) recorded with the Solenoidal Tracker at RHIC (STAR). We compare our results to the Poisson and negative binomial (NBD) expectations, as well as to ultrarelativistic quantum molecular dynamics (UrQMD) and hadron resonance gas (HRG) model predictions. Both NBD and Poisson baselines agree with data within the statistical and systematic uncertainties. UrQMD describes the measured net- Λ C 1 and C 3 at 200 GeV reasonably well but deviates from C 2 , and the deviation increases as a function of collision energy. The ratios of the measured cumulants show no features of critical fluctuations. The chemical freeze-out temperatures extracted from a recent HRG calculation, which was successfully used to describe the net-proton, net-kaon, and net-charge data, indicate Λ freeze-out conditions similar to those of kaons. However, large deviations are found when comparing with temperatures obtained from net-proton fluctuations. The net- Λ cumulants show a weak but finite dependence on the rapidity coverage in the acceptance of the detector, which can be attributed to quantum number conservation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Long-term Ozone Changes and Associated Climate Impacts in CMIP5 Simulations

Ozone changes and associated climate impacts in the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations are analyzed over the historical (1960-2005) and future (2006-2100) period under four Representative Concentration Pathways (RCP). In contrast to CMIP3, where half of the models prescribed constant stratospheric ozone, CMIP5 models all consider past ozone depletion and future ozone recovery. Multimodel mean climatologies and long-term changes in total and tropospheric column ozone calculated from CMIP5 models with either interactive or prescribed ozone are in reasonable agreement with observations. However, some large deviations from observations exist for individual models with interactive chemistry, and these models are excluded in the projections. Stratospheric ozone projections forced with a single halogen, but four greenhouse gas (GHG) scenarios show largest differences in the northern midlatitudes and in the Arctic in spring (approximately 20 and 40 Dobson units (DU) by 2100, respectively). By 2050, these differences are much smaller and negligible over Antarctica in austral spring. Differences in future tropospheric column ozone are mainly caused by differences in methane concentrations and stratospheric input, leading to approximately 10DU increases compared to 2000 in RCP 8.5. Large variations in stratospheric ozone particularly in CMIP5 models with interactive chemistry drive correspondingly large variations in lower stratospheric temperature trends. The results also illustrate that future Southern Hemisphere summertime circulation changes are controlled by both the ozone recovery rate and the rate of GHG increases, emphasizing the importance of simulating and taking into account ozone forcings when examining future climate projections.

tropospheric ozone↗

Constructing auxiliary dynamics for nonequilibrium stationary states by variance minimization

We present a strategy to construct guiding distribution functions (GDFs) based on variance minimization. Auxiliary dynamics via GDFs mitigates the exponential growth of variance as a function of bias in Monte Carlo estimators of large deviation functions. The variance minimization technique exploits the exact properties of eigenstates of the tilted operator that defines the biased dynamics in the nonequilibrium system. We demonstrate our techniques in two classes of problems. In the continuum, we show that GDFs can be optimized to study the interacting driven diffusive systems where the efficiency is systematically improved by incorporating higher correlations into the GDF. On the lattice, we use a correlator product state ansatz to study the 1D weakly asymmetric simple exclusion process. We show that with modest resources, we can capture the features of the susceptibility in large systems that mark the phase transition from uniform transport to a traveling wave state. Our work extends the repertoire of tools available to study nonequilibrium properties in realistic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗