Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

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

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

At least 1,495 records · Page 83

TU-E-213-03: Tools for Ensuring Quality and Safety

Purpose: To develop the first 4D robust optimization (RO) method accounting for respiratory motion and evaluate its potential to improve plan robustness and optimality compared to 3D RO and PTV-based optimization. Methods: A set of 4D CT images are used to track respiratory motion and deformation of tumors and organs. For each of 10 respiration phases, dose distributions for nine different uncertainty scenarios including the nominal one, those incorporating ±5mm setup uncertainties long x, y and z directions and ±3.5% range uncertainties are calculated. All 90 dose distributions are simultaneously optimized to achieve full dose coverage of 10 CTVs and sparing of normal structures. ITV-based 3D RO and PTV-based optimization based on the average CT are also carried for the same patient using same dose volume constrains. After optimization, 4D robustness evaluation was performed for all resulting plans. The CTV coverage and the sparing of normal tissue in 10 phases are evaluated and compared among the three methods. The widths of DVH bands represent the robustness of dose distributions in the structures. Results: For one patient studied so far, the worst case CTV coverage by the prescription dose among all 90 scenarios is: 99% for 4D RO; 88.9% for 3D RO, and 85.2% for PTV based optimization. 4D RO also results in best robustness with the narrowest DVH’ bandwidths for the CTV. 4D and 3D RO have similar organ sparing while PTV based optimization results in worst organ sparing. Conclusion: 4D robust optimization which accounts for anatomy motion and deformation in the optimization process, significantly improves plan robustness and achieves higher quality treatment plans for lung cancer patients. The method is being evaluated for multiple patients with different tumor and motion characteristics.

60 APPLIED LIFE SCIENCES↗

Potential benefits of dosimetric VMAT tracking verified with 3D film measurements

Purpose: To evaluate three different plan adaptation strategies using 3D film-stack dose measurements of both focal boost and hypofractionated prostate VMAT treatments. The adaptation strategies (a couch shift, geometric tracking, and dosimetric tracking) were applied for three realistic intrafraction prostate motions. Methods: A focal boost (35 × 2.2 and 35 × 2.7 Gy) and a hypofractionated (5 × 7.25 Gy) prostate VMAT plan were created for a heterogeneous phantom that allows for internal prostate motion. For these plans geometric tracking and dosimetric tracking were evaluated by ionization chamber (IC) point dose measurements (zero-D) and measurements using a stack of EBT3 films (3D). The geometric tracking applied translations, rotations, and scaling of the MLC aperture in response to realistic prostate motions. The dosimetric tracking additionally corrected the monitor units to resolve variations due to difference in depth, tissue heterogeneity, and MLC-aperture. The tracking was based on the positions of four fiducial points only. The film measurements were compared to the gold standard (i.e., IC measurements) and the planned dose distribution. Additionally, the 3D measurements were converted to dose volume histograms, tumor control probability, and normal tissue complication probability parameters (DVH/TCP/NTCP) as a direct estimate of clinical relevance of the proposed tracking. Results: Compared to the planned dose distribution, measurements without prostate motion and tracking showed already a reduced homogeneity of the dose distribution. Adding prostate motion further blurs the DVHs for all treatment approaches. The clinical practice (no tracking) delivered the dose distribution inside the PTV but off target (CTV), resulting in boost dose errors up to 10%. The geometric and dosimetric tracking corrected the dose distribution’s position. Moreover, the dosimetric tracking could achieve the planned boost DVH, but not the DVH of the more homogeneously irradiated prostate. A drawback of both the geometric and dosimetric tracking was a reduced MLC blocking caused by the rotational component of the MLC aperture corrections. Because of the used CTV to PTV margins and the high doses in the considered fractionation schemes, the TCP differed less than 0.02 from the planned value for all targets and all correction methods. The rectal NTCP constraints, however, could not be realized using any of these methods. Conclusions: The geometric and dosimetric tracking use only a limited input, but they deposit the dose distribution with higher geometric accuracy than the clinical practice. The latter case has boost dose errors up to 10%. The increased accuracy has a modest impact [Δ(NT)CP < 0.02] because of the applied margins and the high dose levels used. To allow further margin reduction tracking methods are vital. The proposed methodology could further be improved by implementing a rotational correction using collimator rotations.

60 APPLIED LIFE SCIENCES↗

Large-scale motions in the universe: Using clusters of galaxies as tracers

Can clusters of galaxies be used to trace the large-scale peculiar velocity field of the universe? We answer this question by using large-scale cosmological simulations to compare the motions of rich clusters of galaxies with the motion of the underlying matter distribution. Three models are investigated: Omega = 1 and Omega = 0.3 cold dark matter (CDM), and Omega = 0.3 primeval baryonic isocurvature (PBI) models, all normalized to the Cosmic Background Explorer (COBE) background fluctuations. We compare the cluster and mass distribution of peculiar velocities, bulk motions, velocity dispersions, and Mach numbers as a function of scale for R greater than or = 50/h Mpc. We also present the large-scale velocity and potential maps of clusters and of the matter. We find that clusters of galaxies trace well the large-scale velocity field and can serve as an efficient tool to constrain cosmological models. The recently reported bulk motion of clusters 689 +/- 178 km/s on approximately 150/h Mpc scale (Lauer & Postman 1994) is larger than expected in any of the models studied (less than or = 190 +/- 78 km/s).

Gramann, Mirt↗

Sub- and super-Maxwellian evaporation of simple gases from liquid water

Non-Maxwellian evaporation of light atoms and molecules (particles) such as He and H{sub 2} from liquids has been observed experimentally. In this work, we use simulations to study systematically the evaporation of Lennard-Jones particles from liquid water. We find instances of sub- and super-Maxwellian evaporation, depending on the mass of the particle and the particle-water interaction strength. The observed trends are in qualitative agreement with experiment. We interpret these trends in terms of the potential of mean force and the effectiveness and frequency of collisions during the evaporation process. The angular distribution of evaporating particles is also analyzed, and it is shown that trends in the energy from velocity components tangential and normal to the liquid surface must be understood separately in order to interpret properly the angular distributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distributed Load Restoration in Unbalanced Active Distribution Systems

We present that the resilience of modern power systems has been constantly threatened by various natural disasters and manmade attacks. Under those threats, fast and automatic restoration actions are critical for restoring impacted systems to normal operating conditions and preventing from additional disastrous consequences. Smart grid technologies, such as distributed energy resources (DERs) and microgrids, provide both opportunities and challenges for distribution system restoration. In this paper, distributed load restoration (DLR), in unbalanced active distribution networks is developed using the alternating direction method of multipliers (ADMM) algorithm. The step-by-step restoration actions are provided to schedule the output of DERs and distribution system control devices, as well as the power consumption from transmission system restoration. The nonlinearities from three-phase unbalanced power flow and distribution components modeling are relaxed into a convex quadratic programming model. Then, the problem is decomposed into subproblems for each node by applying ADMM-DLR; and solves through each agent by exchanging limited information with neighboring nodes in an iterative procedure. The developed models and algorithms are validated and demonstrated through testing of the IEEE 13-node and 123-node test feeders. Lastly, simulation results demonstrate the effectiveness of restoring unbalanced distribution system with the hybrid bottom-up and top-down restoration strategies.

14 SOLAR ENERGY↗

i- flow: High-dimensional integration and sampling with normalizing flows

In many fields of science, high-dimensional integration is required. Numerical methods have been developed to evaluate these complex integrals. We introduce the code i-flow, a python package that performs high-dimensional numerical integration utilizing normalizing flows. Normalizing flows are machine-learned, bijective mappings between two distributions. i-flow can also be used to sample random points according to complicated distributions in high dimensions. We compare i-flow to other algorithms for high-dimensional numerical integration and show that i-flow outperforms them for high dimensional correlated integrals. The i-flow code is publicly available on gitlab at https://gitlab.com/i-flow/i-flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Z-checker: A framework for assessing lossy compression of scientific data

Because of the vast volume of data being produced by today’s scientific simulations and experiments, lossy data compressor allowing user-controlled loss of accuracy during the compression is a relevant solution for significantly reducing the data size. However, lossy compressor developers and users are missing a tool to explore the features of scientific data sets and understand the data alteration after compression in a systematic and reliable way. To address this gap, we have designed and implemented a generic framework called Z-checker. On the one hand, Z-checker combines a battery of data analysis components for data compression. On the other hand, Z-checker is implemented as an open-source community tool to which users and developers can contribute and add new analysis components based on their additional analysis demands. In this article, we present a survey of existing lossy compressors. Then, we describe the design framework of Z-checker, in which we integrated evaluation metrics proposed in prior work as well as other analysis tools. Specifically, for lossy compressor developers, Z-checker can be used to characterize critical properties (such as entropy, distribution, power spectrum, principal component analysis, and autocorrelation) of any data set to improve compression strategies. For lossy compression users, Z-checker can detect the compression quality (compression ratio and bit rate) and provide various global distortion analysis comparing the original data with the decompressed data (peak signal-to-noise ratio, normalized mean squared error, rate–distortion, rate-compression error, spectral, distribution, and derivatives) and statistical analysis of the compression error (maximum, minimum, and average error; autocorrelation; and distribution of errors). Z-checker can perform the analysis with either coarse granularity (throughout the whole data set) or fine granularity (by user-defined blocks), such that the users and developers can select the best fit, adaptive compressors for different parts of the data set. Z-checker features a visualization interface displaying all analysis results in addition to some basic views of the data sets such as time series. To the best of our knowledge, Z-checker is the first tool designed to assess lossy compression comprehensively for scientific data sets.

97 MATHEMATICS AND COMPUTING↗

Nonplanar wing load-line and slender wing theory

Nonplanar load line, slender wing, elliptic wing, and infinite aspect ratio limit loading theories are developed. These are quasi two dimensional theories but satisfy wing boundary conditions at all points along the nonplanar spanwise extent of the wing. These methods are applicable for generalized configurations such as the laterally nonplanar wing, multiple nonplanar wings, or wing with multiple winglets of arbitrary shape. Two dimensional theory infers simplicity which is practical when analyzing complicated configurations. The lateral spanwise distribution of angle of attack can be that due to winglet or control surface deflection, wing twist, or induced angles due to multiwings, multiwinglets, ground, walls, jet or fuselage. In quasi two dimensional theory the induced angles due to these extra conditions are likewise determined for two dimensional flow. Equations are developed for the normal to surface induced velocity due to a nonplanar trailing vorticity distribution. Application examples are made using these methods.

Deyoung, J.↗

Complete quantification of parametric uncertainties in ( d , p ) transfer reactions

Deuteron-induced transfer reactions are a popular probe in nuclear structure and nuclear astro physics studies. The interpretation of these transfer measurements relies on reaction theory that takes as input effective interactions between the nucleons and the target nucleus. Previous work quantified the uncertainty associated with the optical potentials between the nucleons and the target. In this study, we extend that work by also including the parameters of the mean field associated with the overlap function of the final bound state, thus obtaining the full parametric uncertainty on transfer observables. We use Bayesian Markov Chain Monte Carlo simulations to obtain parameter posterior distributions. We use elastic-scattering cross sections to constrain the optical potential parameters and use the asymptotic normalization coefficient of the final state to constrain the bound state interaction. We then propagate these posteriors to the transfer angular distributions and obtain confidence intervals for this observable. We study (d,p) reactions on 14 C, 16 O, and 48 Ca at energies in the range $E$ $d$ = 7 - 24 MeV. Our results show a strong reduction in uncertainty by using the asymptotic normalization coefficient as a constraint, particularly for those reactions most sensitive to ambiguities in the mean field. For those reactions, the importance of constraining the bound state interaction is equal to that of constraining the optical potentials. Furthermore, the case of 14 C is an outlier because the cross section is less sensitive to the nuclear interior. When minimal constraints are used on the parameters of the nucleon-target interaction, the 1σ uncertainties on the differential cross sections are large (~ 500-100%). However, if elastic-scattering data and the asymptotic normalization coefficient are used in the analysis, with an error of 10% (5%), this uncertainty reduces to ~ 30% (~ 15%).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗

Measurement of top quark polarization in t t ¯ lepton + jets final states

We present a measurement of top quark polarization in t ¯ t pair production in p ¯ p collisions at √ s = 1.96 TeV using data corresponding to 9.7 fb -1 of integrated luminosity recorded with the D0 detector at the Fermilab Tevatron Collider. We consider final states containing a lepton and at least three jets. The polarization is measured through the distribution of lepton angles along three axes: the beam axis, the helicity axis, and the transverse axis normal to the t ¯ t production plane. This is the first measurement of top quark polarization at the Tevatron using lepton + jet final states and the first measurement of the transverse polarization in t ¯ t production. The observed distributions are consistent with standard model predictions of nearly no polarization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Minibeam radiation therapy for the management of osteosarcomas: A Monte Carlo study

Purpose: Minibeam radiation therapy (MBRT) exploits the well-established tissue-sparing effect provided by the combination of submillimetric field sizes and a spatial fractionation of the dose. The aim of this work is to evaluate the feasibility and potential therapeutic gain of MBRT, in comparison with conventional radiotherapy, for osteosarcoma treatments. Methods: Monte Carlo simulations (PENELOPE/PENEASY code) were used as a method to study the dose distributions resulting from MBRT irradiations of a rat femur and a realistic human femur phantoms. As a figure of merit, peak and valley doses and peak-to-valley dose ratios (PVDR) were assessed. Conversion of absorbed dose to normalized total dose (NTD) was performed in the human case. Several field sizes and irradiation geometries were evaluated. Results: It is feasible to deliver a uniform dose distribution in the target while the healthy tissue benefits from a spatial fractionation of the dose. Very high PVDR values (⩾20) were achieved in the entrance beam path in the rat case. PVDR values ranged from 2 to 9 in the human phantom. NTD{sub 2.0} of 87 Gy might be reached in the tumor in the human femur while the healthy tissues might receive valley NTD{sub 2.0} lower than 20 Gy. The doses in the tumor and healthy tissues might be significantly higher and lower than the ones commonly delivered used in conventional radiotherapy. Conclusions: The obtained dose distributions indicate that a gain in normal tissue sparing might be expected. This would allow the use of higher (and potentially curative) doses in the tumor. Biological experiments are warranted.

61 RADIATION PROTECTION AND DOSIMETRY↗

Normal-Force and Hinge-Moment Characteristics at Transonic Speeds of Flap-Type Ailerons at Three Spanwise Locations on a 4-Percent-Thick Sweptback-Wing-Body Model and Pressure-Distribution Measurements on an Inboard Aileron

An investigation has been conducted at the Langley 16-foot transonic tunnel to determine the loading characteristics of flap-type ailerons located at inboard, midspan, and outboard positions on a 45 deg. sweptback-wing-body combination. Aileron normal-force and hinge-moment data have been obtained at Mach numbers from 0.80 t o 1.03, at angles of attack up to about 27 deg., and at aileron deflections between approximately -15 deg. and 15 deg. Results of the investigation indicate that the loading over the ailerons was established by the wing-flow characteristics, and the loading shapes were irregular in the transonic speed range. The spanwise location of the aileron had little effect on the values of the slope of the curves of hinge-moment coefficient against aileron deflection, but the inboard aileron had the greatest value of the slope of the curves of hinge-moment coefficient against angle of attack and the outboard aileron had the least. Hinge-moment and aileron normal-force data taken with strain-gage instrumentation are compared with data obtained with pressure measurements.

Runckel, Jack F.↗

Highly Inclined Jets in Cross Flow

Results from an experimental investigation of flow field generated by pitched and yawed jets discharging from a flat plate into a cross-flow are presented. The circular jet was pitched at alpha = 20deg and 45deg and yawed between Beta = 0deg and 90deg in increments of 15deg. The measurements were performed with two ×-wires providing all three components of velocity and turbulent stresses. These data were obtained at downstream locations of x = 3, 5, 10 and 20, where the distance x, normalized by the jet diameter, is measured from the center of the orifice. Data for all configurations were acquired at a momentum-flux ratio J = 8. Additionally, for selected angles and locations, surveys were conducted for J = 1.5, 4, and 20. As expected, the jet penetration is found to be higher at larger alpha. With increasing beta the jet spreads more. The rate of reduction of peak streamwise vorticity, ω max, with the downstream distance is significantly less at higher Beta but is found to be practically independent of alpha. Thus, at the farthest measurement station x = 20, ωxmax is about five times larger for Beta = 75deg compared to the levels at Beta = 0deg. Streamwise velocity within the jet-vortex structure is found to depend on the parameter J. At J = 1.5 and 4, 'wake-like' velocity profiles are observed. In comparison, a 'jet-like' overshoot is present at higher J. Distributions of turbulent stresses for various cases are documented. Peak normal stresses are found to occur within the core of the streamwise vortices. With yaw, at lower values of J, high turbulence is also observed in the boundary layer underneath the jet-vortex structure

Milanovic, I. M.↗

Detection of scattering in the 2175 A interstellar band

Spectrophotometric observations of two reflection nebulae and their illuminating stars obtained with the International Ultraviolet Explorer Satellite have provided the first evidence for the presence of a scattering contribution to the 2175 A interstellar extinction band. Lower than normal far-UV extinction for the stars embedded in the nebulae indicates that the nebulae have a dust particle size distribution that is dominated by larger particles. The strength of the 2175 A band is larger than normal in both cases. The scattering is found to dominate the long-wavelength wing of the band, without shifting the central wavelength of the band by more than 20 A toward longer wavelengths. These observations are taken to indicate that the solid particles responsible for the 2175 A band can be considerably larger than the Rayleigh limit in some interstellar locations. The absence of a notable shift in the central wavelength of the band in such large particles presents a new severe constraint for models of interstellar grains.

Witt, A. N.↗

Volumetric Modulated Arc Therapy Planning for Primary Prostate Cancer With Selective Intraprostatic Boost Determined by {sup 18}F-Choline PET/CT

Purpose: This study evaluated expected tumor control and normal tissue toxicity for prostate volumetric modulated arc therapy (VMAT) with and without radiation boosts to an intraprostatically dominant lesion (IDL), defined by {sup 18}F-choline positron emission tomography/computed tomography (PET/CT). Methods and Materials: Thirty patients with localized prostate cancer underwent {sup 18}F-choline PET/CT before treatment. Two VMAT plans, plan{sub 79} {sub Gy} and plan{sub 100-105} {sub Gy}, were compared for each patient. The whole-prostate planning target volume (PTV{sub prostate}) prescription was 79 Gy in both plans, but plan{sub 100-105} {sub Gy} added simultaneous boost doses of 100 Gy and 105 Gy to the IDL, defined by 60% and 70% of maximum prostatic uptake on {sup 18}F-choline PET (IDL{sub suv60%} and IDL{sub suv70%}, respectively, with IDL{sub suv70%} nested inside IDL{sub suv60%} to potentially enhance tumor specificity of the maximum point dose). Plan evaluations included histopathological correspondence, isodose distributions, dose-volume histograms, tumor control probability (TCP), and normal tissue complication probability (NTCP). Results: Planning objectives and dose constraints proved feasible in 30 of 30 cases. Prostate sextant histopathology was available for 28 cases, confirming that IDL{sub suv60%} adequately covered all tumor-bearing prostate sextants in 27 cases and provided partial coverage in 1 case. Plan{sub 100-105} {sub Gy} had significantly higher TCP than plan{sub 79} {sub Gy} across all prostate regions for α/β ratios ranging from 1.5 Gy to 10 Gy (P<.001 for each case). There were no significant differences in bladder and femoral head NTCP between plans and slightly lower rectal NTCP (endpoint: grade ≥ 2 late toxicity or rectal bleeding) was found for plan{sub 100-105} {sub Gy}. Conclusions: VMAT can potentially increase the likelihood of tumor control in primary prostate cancer while observing normal tissue tolerances through simultaneous delivery of a steep radiation boost to a {sup 18}F-choline PET-defined IDL.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Scale Mixtures of Frechet Distributions as Asymptotic Approximations of Extreme Precipitation

This paper is a further development of the results of [20] where, based on the negative binomial model for the duration of wet periods measured in days [16], an asymptotic approximation was proposed for the distribution of the maximum daily precipitation volume within a wet period. This approximation has the form of a scale mixture of the Fr´echet distribution with the gamma mixing distribution and coincides with the distribution of a positive power of a random variable having the Snedecor–Fisher distribution. Here we extend this result to the mth extremes, m ∈ ℕ, and sample quantiles. The proof of this result is based on the representation of the negative binomial distribution as a mixed geometric (and hence, mixed Poisson) distribution [17] and limit theorems for extreme order statistics in samples with random sizes having mixed Poisson distributions [10]. Some analytic properties of the obtained limit distribution are described. In particular, it is demonstrated that under certain conditions the limit distribution of the maximum precipitation is mixed exponential and hence, is infinitely divisible. It is shown that under the same conditions this limit distribution can be represented as a scale mixture of stable or Weibull or Pareto or folded normal laws. The corresponding product representations for the limit random variable can be used for its computer simulation. The results of fitting this distribution to real data are presented illustrating high adequacy of the proposed model. It is also shown that the limit distribution of sample quantiles is the well-known Student distribution. Several methods are proposed for the estimation of the parameters of the asymptotic distributions. The obtained mixture representations for the limit laws and the corresponding asymptotic approximations provide better insight into the nature of mixed probability (“Bayesian”) models.

97 MATHEMATICS AND COMPUTING↗