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 271 records · Page 15

Log-normal model for microbial survival in heat sterilization

An analytical model is developed for the survival times of organisms in heat sterilization in which the probability of inactivation as a function of exposure time is log-normally distributed. Experimental data is examined relative to this model and it is concluded that the model is valid except during the initial period of heating when an additional interaction between the organisms and their surrounding medium appears to be present. At long heating times, the log-normal model appears more accurate for extrapolating to low survival probabilities than the usual logarithmic survivor curves and is generally more conservative.

SURVIVAL↗

Particle Morphology and Size Results from the Smoke Aerosol Measurement Experiment-2

Results are presented from the Reflight of the Smoke Aerosol Measurement Experiment (SAME-2) which was conducted during Expedition 24 (July-September 2010). The reflight experiment built upon the results of the original flight during Expedition 15 by adding diagnostic measurements and expanding the test matrix. Five different materials representative of those found in spacecraft (Teflon, Kapton, cotton, silicone rubber and Pyrell) were heated to temperatures below the ignition point with conditions controlled to provide repeatable sample surface temperatures and air flow. The air flow past the sample during the heating period ranged from quiescent to 8 cm/s. The smoke was initially collected in an aging chamber to simulate the transport time from the smoke source to the detector. This effective transport time was varied by holding the smoke in the aging chamber for times ranging from 11 to 1800 s. Smoke particle samples were collected on Transmission Electron Microscope (TEM) grids for post-flight analysis. The TEM grids were analyzed to observe the particle morphology and size parameters. The diagnostics included a prototype two-moment smoke detector and three different measures of moments of the particle size distribution. These moment diagnostics were used to determine the particle number concentration (zeroth moment), the diameter concentration (first moment), and the mass concentration (third moment). These statistics were combined to determine the diameter of average mass and the count mean diameter and, by assuming a log-normal distribution, the geometric mean diameter and the geometric standard deviations can also be calculated. Overall the majority of the average smoke particle sizes were found to be in the 200 nm to 400 nm range with the quiescent cases producing some cases with substantially larger particles.

Urban, David L.↗

An analysis of fracture trace patterns in areas of flat-lying sedimentary rocks for the detection of buried geologic structure

Two study areas in a cratonic platform underlain by flat-lying sedimentary rocks were analyzed to determine if a quantitative relationship exists between fracture trace patterns and their frequency distributions and subsurface structural closures which might contain petroleum. Fracture trace lengths and frequency (number of fracture traces per unit area) were analyzed by trend surface analysis and length frequency distributions also were compared to a standard Gaussian distribution. Composite rose diagrams of fracture traces were analyzed using a multivariate analysis method which grouped or clustered the rose diagrams and their respective areas on the basis of the behavior of the rays of the rose diagram. Analysis indicates that the lengths of fracture traces are log-normally distributed according to the mapping technique used. Fracture trace frequency appeared higher on the flanks of active structures and lower around passive reef structures. Fracture trace log-mean lengths were shorter over several types of structures, perhaps due to increased fracturing and subsequent erosion. Analysis of rose diagrams using a multivariate technique indicated lithology as the primary control for the lower grouping levels. Groupings at higher levels indicated that areas overlying active structures may be isolated from their neighbors by this technique while passive structures showed no differences which could be isolated.

Podwysocki, M. H.↗

A hill-sliding strategy for initialization of Gaussian clusters in the multidimensional space

A hill sliding technique was devised to extract Gaussian clusters from the multivariate probability density estimate of sample data for the first step of iterative unsupervised classification. Each cluster was assumed to posses a unimodal normal distribution. A clustering function proposed distinguished elements of a cluster under formation from the rest in the feature space. Initial clusters were extracted one by one according to the hill sliding tactics. A dimensionless cluster compactness parameter was proposed as a universal measure of cluster goodness and used satisfactorily in test runs with LANDSAT multispectral scanner data. The normalized divergence, defined by the cluster divergence divided by the entropy of the entire sample data, was utilized as a general separability measure between clusters. An overall clustering objective function was set forth in terms of cluster covariance matrices, from which the cluster compactness measure could be deduced. Minimal improvement of initial data partitioning was evaluated by this objective function in eliminating scattered sparse data points. The hill sliding clustering technique developed herein has the potential applicability to decomposition any multivariate mixture distribution into a number of unimodal distributions when an appropriate distribution function to the data set is employed.

Park, J. K.↗

Thermal pressure in the cold neutral medium of nearby galaxies

Dynamic and thermal processes regulate the structure of the multiphase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure (P{sub th}) of the diffuse, neutral gas. By combining [C ii] 158 μm, H i, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of ∼1 kiloparsec. We find a distribution of thermal pressures in the P{sub th}/k∼10{sup 3}--10{sup 5} K cm{sup −3} range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of P{sub th}/k ≈ 3600 K cm{sup −3}. These results are consistent with thermal pressure measurements using other observational methods. We find that P{sub th} increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral media could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about ∼30% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical analyses of reservoir and fracturing parameters for a multifractured shale oil reservoir in Mississippi

The use of multifractured horizontal wells has improved the efficiency of hydrocarbon extraction from shale gas and oil plays. It is highly desirable to estimate the characteristics of the reservoir and well fracturing through production data analysis. Production rate transient data from 16 wells in the Mississippi sections were analyzed to estimate local reservoir permeability and hydraulic fracture parameters. Key factors affecting well productivity have been identified for optimizing future well completion. Based on the theory of distance of investigation during transient linear flow, the reservoir matrix permeability of TMS was estimated to be between 53 nd and 210 nd, averaging at 116 nd. The matrix permeability follows the lognormal distribution and is considered homogeneous according to the coefficient of variation. The fracture half-length was estimated using the matrix permeability data from the rate transient analysis. The fracture half-length was found to have a mean value of 234 ft with a standard deviation of 66 ft in a normal distribution. The fracture conductivity was back-calculated by matching pseudosteady production rate data to Guo et al's (SPE Reserve Eval Eng.12, 2009, 879) productivity model for boundary-dominated flow. The fracture conductivity was estimated to range from 0.4 md-ft to 3.2 md-ft, averaging at 1.4 md-ft. Based on Guo et al's (SPE Reserve Eval Eng.12, 2009, 879) well productivity model applied to TMS condition, well productivity can be improved significantly by increasing fracture length and conductivity. Future TMS wells should be completed with conductivity values of greater than 10 md-ft.

04 OIL SHALES AND TAR SANDS↗

Evaluating quantitative 3-D image analysis as a design tool for low enriched uranium fuel compacts for the transient reactor test facility: A preliminary study

In this study, 3-D image analysis when combined with a non-destructive examination technique such as X-ray computed tomography (CT) provides a highly quantitative tool for the investigation of a material’s structure. In this investigation 3-D image analysis and X-ray CT were combined to analyze the microstructure of a preliminary subsized fuel compact for the Transient Reactor Test Facility’s low enriched uranium conversion program to assess the feasibility of the combined techniques for use in the optimization of the fuel compact fabrication process. The quantitative image analysis focused on determining the size and spatial distribution of the surrogate fuel particles and the size, shape, and orientation of voids within the compact. Additionally, the maximum effect of microstructural features on heat transfer through the carbonaceous matrix of the preliminary compact was estimated. The surrogate fuel particles occupied 0.8% of the compact by volume with a log-normal distribution of particle sizes with a mean diameter of 39 μm and a standard deviation of 16 μm. Roughly 39% of the particles had a diameter greater than the specified maximum particle size of 44 μm suggesting that the particles agglomerate during fabrication. The local volume fraction of particles also varies significantly within the compact although uniformities appear to be evenly dispersed throughout the analysed volume. The voids produced during fabrication were on average plate-like in nature with their major axis oriented perpendicular to the compaction direction of the compact. Finally, the microstructure, mainly the large preferentially oriented voids, may cause a small degree of anisotropy in the thermal diffusivity within the compact. α∥/α⊥, the ratio of thermal diffusivities parallel to and perpendicular to the compaction direction are expected to be no less than 0.95 with an upper bound of 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probabilistic SSME blades structural response under random pulse loading

The purpose is to develop models of random impacts on a Space Shuttle Main Engine (SSME) turbopump blade and to predict the probabilistic structural response of the blade to these impacts. The random loading is caused by the impact of debris. The probabilistic structural response is characterized by distribution functions for stress and displacements as functions of the loading parameters which determine the random pulse model. These parameters include pulse arrival, amplitude, and location. The analysis can be extended to predict level crossing rates. This requires knowledge of the joint distribution of the response and its derivative. The model of random impacts chosen allows the pulse arrivals, pulse amplitudes, and pulse locations to be random. Specifically, the pulse arrivals are assumed to be governed by a Poisson process, which is characterized by a mean arrival rate. The pulse intensity is modelled as a normally distributed random variable with a zero mean chosen independently at each arrival. The standard deviation of the distribution is a measure of pulse intensity. Several different models were used for the pulse locations. For example, three points near the blade tip were chosen at which pulses were allowed to arrive with equal probability. Again, the locations were chosen independently at each arrival. The structural response was analyzed both by direct Monte Carlo simulation and by a semi-analytical method.

Shiao, Michael↗

Distributions of the interplanetary magnetic field revisited

The adequacy of the power spectrum to characterize the variations of a parameter depends on whether or not the parameter has a Gaussian distribution. We here perform very simple tests of Gaussianity on the distribution. We here perform very simple tests of Gaussianity on the distributions of the magnitudes of the interplanetary magnetic field, and on the distributions of the components; that is, we find the first four cumulants of the distributions (mean, variance, skewness, and kurtosis) and their solar cycle variations. We find, consistent with other recent analyses, that the traditional distributions of the 1-hour averaged magnitude are not distributed normally or lognomally as has often been assumed and the 1-hour averaged z component is found to have a nonzero kurtosis. Thus the power spectrum is insufficient to completely characterize these variations and polyspectra are needed. We have isolated variations in the 1/f frequency region of the spectrum and show that the distributions of the magnitudes have nonzero skewness and kurtosis, the magnitudes are not distributed lognormally, and the distributions of the components have nonzero kurtosis. Thus higher-order spectra are again needed for a full characterization.

Feynman, Joan↗

TU-H-CAMPUS-JeP3-05: Adaptive Determination of Needle Sequence HDR Prostate Brachytherapy with Divergent Needle-By-Needle Delivery

Purpose: To develop a new method which adaptively determines the optimal needle insertion sequence for HDR prostate brachytherapy involving divergent needle-by-needle dose delivery by e.g. a robotic device. A needle insertion sequence is calculated at the beginning of the intervention and updated after each needle insertion with feedback on needle positioning errors. Methods: Needle positioning errors and anatomy changes may occur during HDR brachytherapy which can lead to errors in the delivered dose. A novel strategy was developed to calculate and update the needle sequence and the dose plan after each needle insertion with feedback on needle positioning errors. The dose plan optimization was performed by numerical simulations. The proposed needle sequence determination optimizes the final dose distribution based on the dose coverage impact of each needle. This impact is predicted stochastically by needle insertion simulations. HDR procedures were simulated with varying number of needle insertions (4 to 12) using 11 patient MR data-sets with PTV, prostate, urethra, bladder and rectum delineated. Needle positioning errors were modeled by random normally distributed angulation errors (standard deviation of 3 mm at the needle’s tip). The final dose parameters were compared in the situations where the needle with the largest vs. the smallest dose coverage impact was selected at each insertion. Results: Over all scenarios, the percentage of clinically acceptable final dose distribution improved when the needle selected had the largest dose coverage impact (91%) compared to the smallest (88%). The differences were larger for few (4 to 6) needle insertions (maximum difference scenario: 79% vs. 60%). The computation time of the needle sequence optimization was below 60s. Conclusion: A new adaptive needle sequence determination for HDR prostate brachytherapy was developed. Coupled to adaptive planning, the selection of the needle with the largest dose coverage impact increases chances of reaching the clinical constraints. M. Borot de Battisti is funded by Philips Medical Systems Nederland B.V.; M. Moerland is principal investigator on a contract funded by Philips Medical Systems Nederland B.V.; G. Hautvast and D. Binnekamp are fulltime employees of Philips Medical Systems Nederland B.V.

60 APPLIED LIFE SCIENCES↗

Probability Distributions for Effective Permeability of Potentially Leaking Wells at CO 2 Sequestration Sites

Wells with defects originating from either poor completions or subsequent damage arising from chemical, thermal and/or mechanical stresses provide potential pathways for leakage of carbon dioxide (CO 2 ) and brine from a CO 2 sequestration site. This report provides a methodology for analyzing data on the frequency of occurrence of defective wells and the effective permeability of the damaged wells. The results include bimodal, log-normal distributions of the effective permeability of wells in addition to other distributions designed to provide a measure of the likelihood of encountering a leaking well in a CO 2 storage project.

58 GEOSCIENCES↗

Non-Gaussian Distributions of TOA SW Flux as Observed by MISR and CERES

The Top of Atmosphere (TOA) shortwave (SW) flux, converted from Terra Multi angle Imaging‐SpectroRadiometer (MISR) narrow band albedos, is compared with that measured from Clouds and the Earth’s Radiant Energy System (CERES). We describe the probability density function (PDF) of the monthly TOA SW flux and how the statistical third moment, skewness, can impact the quantification of the flux. The PDF of the SW flux is not normally distributed but positively skewed. In both sets of observations, the near-global (80S-80N) median value of the SW flux is ≈3 W/m 2 less than the mean value, due to the positive skewness of the distribution. The near-global mean TOA SW flux converted from MISR is about 7 W/m 2 (≈7%) less than CERES measured flux during the last two decades. Surprisingly, hemispheric asymmetry exists with TOA SW observations from Terra platform. SH reflects 3.92 W/m 2 and 1.15 W/m 2 more mean SW flux than NH, from MISR and CERES Single Scanner Footprint products, respectively. We can infer that the offsetting by morning clouds in the SH is greater than the effect of hemispheric imbalance of SW flux caused by different land masses in two hemispheres. While the characteristics of the two SW fluxes broadly agree with each other, differences in the regional PDF from two different SW fluxes are substantially different over high cloud regions and high altitude regions. Our analysis shows that some parts of the different skewness from two measurements may be attributed to the different calibration of the radiance anisotropy over high cloud scenes.

satellite measurements↗

High-frequency wave normals in the solar wind

High-frequency (0.01-0.04 Hz) magnetic fluctuations in 506 ten-minute intervals of contemporaneous Explorer 35 and Apollo 12 measurements made in the solar wind near the morning side of the earth's bow shock show the presence of a large population of disturbances resembling Alfven waves. Each wavefront normal is systematically aligned (median deviation = 35 deg) with the associated ten-minute average of the magnetic field. Because of variability in the direction of this ten-minute average, from one interval to another, the coupled wavefront-normal distribution is nearly isotropic in solar ecliptic coordinates, in contrast with the results of other studies of waves at much lower frequency indicating outward propagation from the sun. Presumably the high frequency waves discussed here are stirred into isotropy (in solar ecliptic coordinates) by following the low frequency fluctuations. As these waves maintain their alignment of wavefront-normal with the field ten-minute average despite the great variation of the field ten-minute average, a strong physical alignment constraint is inferred.

Herbert, F.↗

Concrete mixture characterization. Cementitious barriers partnership

This report summarizes the characterization study performed on two concrete mixtures used for radioactive waste storage. Both mixtures were prepared with approximately 425 kg of binder. The testing protocol mostly focused on determining the transport properties of the mixtures; volume of permeable voids (porosity), diffusion coefficients, and water permeability were evaluated. Tests were performed after different curing durations. In order to obtain data on the statistical distribution of transport properties, the measurements after 2 years of curing were performed on 10+ samples. Overall, both mixtures exhibited very low tortuosities and permeabilities, a direct consequence of their low water-to-binder ratio and the use of supplementary cementitious materials. The data generated on 2-year old samples showed that porosity, tortuosity and permeability follow a normal distribution. Chloride ponding tests were also performed on test samples. They showed limited chloride ingress, in line with measured transport properties. These test results also showed that both materials react differently with chloride, a consequence of the differences in the binder chemical compositions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

THE MASS DISTRIBUTION FUNCTION OF PLANETS

The distribution of orbital period ratios of adjacent planets in extrasolar planetary systems discovered by the Kepler space telescope exhibits a peak near ∼1.5–2, a long tail of larger period ratios, and a steep drop-off in the number of systems with period ratios below ∼1.5. We find from these data that the dimensionless orbital separations have an approximately log-normal distribution. Using Hill’s criterion for the dynamical stability of two planets, we find an upper bound on planet masses such that the most common planet mass does not exceed 10{sup −3.2}m{sub ∗}, or about two-thirds of Jupiter’s mass for solar-mass stars. Assuming that the mass ratio and the dynamical separation (orbital spacings in units of mutual Hill radius) of adjacent planets are independent random variates, and adopting empirical distributions for these, we use Hill’s criterion in a statistical way to estimate the planet mass distribution function from the observed distribution of orbital separations. We find that the planet mass function is peaked in logarithm of mass, with a peak value and standard deviation of logm/M{sub ⊕} of ∼(0.6−1.0) and ∼(1.1−1.2), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Reconciling Chord Length Distributions and Area Distributions for Fields of Fractal Cumulus Clouds

While the total cover of broken cloud fields can in principle be obtained from one-dimensional measurements, the cloud size distribution normally differs between two-dimensional (area) and one-dimensional retrieval (chord length) methods. In this study, we use output from high-resolution Large Eddy Simulations to generate a transfer function between the two. We retrieve chord lengths and areas for many clouds, and plot the one as a function of the other, and vice versa. We find that the cloud area distribution conditional on the chord length behaves like a gamma distribution with well-behaved parameters, with a mean μ=1.1L and a shape parameter β=L−0.645. Using this information, we are able to generate a transfer function that can adjust the chord length distribution so that it comes much closer to the cloud area distribution. Our transfer function improves the error in predicting the mean cloud size, and is performs without strong biases for smaller sample sizes. However, we find that the method is still has difficulties in accurately predicting the frequency of occurrence of the largest cloud sizes.

54 ENVIRONMENTAL SCIENCES↗

MO-FG-202-06: Improving the Performance of Gamma Analysis QA with Radiomics- Based Image Analysis

Purpose: The use of gamma analysis for IMRT quality assurance has well-known limitations. Traditionally, a simple thresholding technique is used to evaluated passing criteria. However, like any image the gamma distribution is rich in information which thresholding mostly discards. We therefore propose a novel method of analyzing gamma images that uses quantitative image features borrowed from radiomics, with the goal of improving error detection. Methods: 368 gamma images were generated from 184 clinical IMRT beams. For each beam the dose to a phantom was measured with EPID dosimetry and compared to the TPS dose calculated with and without normally distributed (2mm sigma) errors in MLC positions. The magnitude of 17 intensity histogram and size-zone radiomic features were derived from each image. The features that differed most significantly between image sets were determined with ROC analysis. A linear machine-learning model was trained on these features to classify images as with or without errors on 180 gamma images.The model was then applied to an independent validation set of 188 additional gamma distributions, half with and half without errors. Results: The most significant features for detecting errors were histogram kurtosis (p=0.007) and three size-zone metrics (p<1e-6 for each). The sizezone metrics detected clusters of high gamma-value pixels under mispositioned MLCs. The model applied to the validation set had an AUC of 0.8, compared to 0.56 for traditional gamma analysis with the decision threshold restricted to 98% or less. Conclusion: A radiomics-based image analysis method was developed that is more effective in detecting error than traditional gamma analysis. Though the pilot study here considers only MLC position errors, radiomics-based methods for other error types are being developed, which may provide better error detection and useful information on the source of detected errors. This work was partially supported by a grant from the Agency for Healthcare Research and Quality, grant number R18 HS022244-01.

60 APPLIED LIFE SCIENCES↗

Evaluating Safeguards Statistical Assumptions via Stochastic Simulation

Herein, the authors built and tested a stochastic simulation to estimate achieved detection probabilities (DPs) on a stratum basis, over a tailorable range of diverted amounts from 0 to 2 SQ, using typical IAEA inspection data: i.e., SQ in stratum, number of items, number of gross/partial/bias defect measurements conducted, and realistic relative standard deviation (RSD) values for typical IAEA verification measurements. For bulk strata, the model calculates achieved DP at 0.01 SQ diversion increments; for item strata, the model calculates DP using the smallest realistic diversion increment (e.g., a plate, pin, or coupon). After successfully benchmarking against IAEA deterministic models, the simulation was used to test the sensitivity of DP to certain standard assumptions and selected input parameters. First, the equal defect assumption was tested; the results suggest significant complexity in the effectiveness of partial defect measurements. Next, the authors explored the sensitivity of DP to the assumed RSD of attribute tests. Then, the authors compared non-normal models for instrument performance (e.g., logistic, step, or arbitrary functions) to the typical results from a normal distribution (characterized by RSD). This last comparison was supplemented with experimentally derived performance data for an HM-5. The HM-5 was used to make enrichment measurements on both LEU and HEU MTR fuel elements as plates were removed, and the results fit with logistic and step curves and applied in the simulation. These stochastic DP results were compared to DP estimates from a deterministic model assuming a normal curve and typical RSD, yielding insights that could improve effectiveness in the field. These early results illustrate the potential of stochastic models to better understand achieved DP and to improve safeguards effectiveness.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗