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 91 records · Page 5

Influence of Viscosity on Dynamic Magnetization of Thermally Blocked Iron Oxide Nanoparticles Characterized by a Sensitive AC Magnetometer

In this work, we show that the viscosity of carrier liquid affects the dynamic magnetization of thermally blocked multi-core iron oxide nanoparticles. The core size of the nanoparticles was determined from the magnetization curve measured by a specially developed high-T{sub c} SQUID magnetometer and calculated to be 11.7 nm. Using an AC magnetometer developed based on induction coils, the dynamic magnetization of the multi-core iron oxide nanoparticle solution was measured from 3 Hz to 10.48 kHz. Later, we reconstructed of the hydrodynamic size distribution of the particles by assuming a log-normal distribution of particle size in an AC susceptibility model by Shliomis and Stepanov, which accounts for anisotropic directions of the easy axes of magnetic nanoparticles with respect to the excitation field direction. The reconstructed hydrodynamic sizes showed an average diameter of 130 nm and agreed with the size determined by dynamic light scattering method. In the case of increasing viscosity of the carrier liquids from 0.89 to 8.11 mPa s, the dynamic magnetization peaks of the imaginary component have shifted to a lower frequency region. We showed that the harmonics ratio and phase delay upon the magnetic field excitation at 30 Hz could also be used to determine the viscosity of carrier liquid independently.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

The thermal pressure distribution of a simulated cold neutral medium

We numerically study the thermal pressure distribution in a gas with thermal properties similar to those of the cold neutral interstellar gas by analyzing three-dimensional hydrodynamic models in boxes with sides of 100 pc with turbulent compressible forcing at 50 pc and different Mach numbers. We find that at high pressures and for large Mach numbers, both the volume-weighted and the density-weighted distributions can be appropriately described by a log-normal distribution, whereas for small Mach numbers they are better described by a power law. Thermal pressure distributions resulting from similar simulations but with self-gravity differ only for low Mach numbers; in this case, they develop a high pressure tail.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic particle imaging of blood coagulation

We investigated the feasibility of visualizing blood coagulation using a system for magnetic particle imaging (MPI). A magnetic field-free line is generated using two opposing neodymium magnets and transverse images are reconstructed from the third-harmonic signals received by a gradiometer coil, using the maximum likelihood-expectation maximization algorithm. Our MPI system was used to image the blood coagulation induced by adding CaCl{sub 2} to whole sheep blood mixed with magnetic nanoparticles (MNPs). The “MPI value” was defined as the pixel value of the transverse image reconstructed from the third-harmonic signals. MPI values were significantly smaller for coagulated blood samples than those without coagulation. We confirmed the rationale of these results by calculating the third-harmonic signals for the measured viscosities of samples, with an assumption that the magnetization and particle size distribution of MNPs obey the Langevin equation and log-normal distribution, respectively. We concluded that MPI can be useful for visualizing blood coagulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

FIRST DIRECT MEASUREMENTS OF TRANSVERSE WAVES IN SOLAR POLAR PLUMES USING SDO/AIA

There is intense interest in determining the precise contribution of Alfvénic waves propagating along solar structures to the problems of coronal heating and solar wind acceleration. Since the launch of SDO/AIA, it has been possible to resolve transverse oscillations in off-limb solar polar plumes and recently McIntosh et al. concluded that such waves are energetic enough to play a role in heating the corona and accelerating the fast solar wind. However, this result is based on comparisons to Monte Carlo simulations and confirmation via direct measurements is still outstanding. Thus, this Letter reports on the first direct measurements of transverse wave motions in solar polar plumes. Over a four hour period, we measure the transverse displacements, periods, and velocity amplitudes of 596 distinct oscillations observed in the 171 Å channel of SDO/AIA. We find a broad range of non-uniformly distributed parameter values which are well described by log-normal distributions with peaks at 234 km, 121 s, and 8 km s{sup –1}, and mean and standard deviations of 407 ± 297 km, 173 ± 118 s, and 14 ± 10 km s{sup –1}. Within standard deviations, our direct measurements are broadly consistent with previous results. However, accounting for the whole of our observed non-uniform parameter distribution we calculate an energy flux of 9-24 W m{sup –2}, which is 4-10 times below the energy requirement for solar wind acceleration. Hence, our results indicate that transverse magnetohydrodynamic waves as resolved by SDO/AIA cannot be the dominant energy source for fast solar wind acceleration in the open-field corona.

79 ASTRONOMY AND ASTROPHYSICS↗

Model-based Dynamic Control of Speculative Forays in Parallel Computation

In simulations running in parallel, the processors would have to synchronize with other processors to maintain correct global order of computations. This can be done either by blocking computation until correct order is guaranteed, or by \textit{speculatively} proceeding with the best guess (based on local information) and later correcting errors if/as necessary. Since the gainful lengths of speculative forays depend on the dynamics of the application software and hardware at runtime, an online control system is necessary to dynamically choose and/or switch between the blocking and speculative strategies. In this paper, we formulate the reversible speculative computing in large-scale parallel computing as a dynamic linear feedback control (optimization) system model and evaluate its performance in terms of time and cost savings as compared to the traditional (forward) computing. We illustrate with an exact analogy in the form of vehicular travel under dynamic, delayed route information. The objective is to assist in making the optimal decision on what computational approach is to be chosen, by predicting the amount of time and cost savings (or losing) under different environments represented by different parameters and probability distribution functions. We consider the cases of Gaussian, exponential and log-normal distribution functions. The control system is intended for incorporating into speculative parallel applications such as optimistic parallel discrete event simulations to decide at runtime when and to what extent speculative execution can be performed gainfully.

Perumalla, Kalyan↗

Grid Frequency Extreme Event Analysis and Modeling: Preprint

Sudden losses of generation or load can lead to instantaneous changes in electric grid frequency and voltage. Extreme frequency events pose a major threat to grid stability. As renewable energy sources supply power to grids in increasing proportions, it becomes increasingly important to examine when and why extreme events occur to prevent destabilization of the grid. To better understand frequency events, including extrema, historic data were analyzed to fit probability distribution functions to various frequency metrics. Results showed that a standard Cauchy distribution fit the difference between the frequency nadir and prefault frequency (f_(C-A)) metric well, a standard Cauchy distribution fit the settling frequency (f_B) metric well, and a standard normal distribution fit the difference between the settling frequency and frequency nadir (f_(B-C)) metric very well. Results were inconclusive for the frequency nadir (f_C) metric, meaning it likely has a more complex distribution than those tested. This probabilistic modeling should facilitate more realistic modeling of grid faults.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Prediction of interest rate using CKLS model with stochastic parameters

The Chan, Karolyi, Longstaff and Sanders (CKLS) model is a popular one-factor model for describing the spot interest rates. In this paper, the four parameters in the CKLS model are regarded as stochastic. The parameter vector φ{sup (j)} of four parameters at the (J+n)-th time point is estimated by the j-th window which is defined as the set consisting of the observed interest rates at the j′-th time point where j≤j′≤j+n. To model the variation of φ{sup (j)}, we assume that φ{sup (j)} depends on φ{sup (j−m)}, φ{sup (j−m+1)},…, φ{sup (j−1)} and the interest rate r{sub j+n} at the (j+n)-th time point via a four-dimensional conditional distribution which is derived from a [4(m+1)+1]-dimensional power-normal distribution. Treating the (j+n)-th time point as the present time point, we find a prediction interval for the future value r{sub j+n+1} of the interest rate at the next time point when the value r{sub j+n} of the interest rate is given. From the above four-dimensional conditional distribution, we also find a prediction interval for the future interest rate r{sub j+n+d} at the next d-th (d≥2) time point. The prediction intervals based on the CKLS model with stochastic parameters are found to have better ability of covering the observed future interest rates when compared with those based on the model with fixed parameters.

97 MATHEMATICS AND COMPUTING↗

Collection strategy, inner morphology, and size distribution of dust particles in ASDEX Upgrade

The dust collection and analysis strategy in ASDEX Upgrade (AUG) is described. During five consecutive operation campaigns (2007–2011), Si collectors were installed, which were supported by filtered vacuum sampling and collection with adhesive tapes in 2009. The outer and inner morphology (e.g. shape) and elemental composition of the collected particles were analysed by scanning electron microscopy. The majority of the ~50?000 analysed particles on the Si collectors of campaign 2009 contain tungsten—the plasma-facing material in AUG—and show basically two different types of outer appearance: spheroids and irregularly shaped particles. By far most of the W-dominated spheroids consist of a solid W core, i.e. solidified W droplets. A part of these particles is coated with a low-Z material; a process that seems to happen presumably in the far scrape-off layer plasma. In addition, some conglomerates of B, C and W appear as spherical particles after their contact with plasma. By far most of the particles classified as B-, C- and W-dominated irregularly shaped particles consist of the same conglomerate with varying fraction of embedded W in the B–C matrix and some porosity, which can exceed 50%. The fragile structures of many conglomerates confirm the absence of intensive plasma contact. Both the ablation and mobilization of conglomerate material and the production of W droplets are proposed to be triggered by arcing. The size distribution of each dust particle class is best described by a log-normal distribution allowing an extrapolation of the dust volume and surface area. The maximum in this distribution is observed above the resolution limit of 0.28 µm only for the W-dominated spheroids, at around 1 µm. The amount of W-containing dust is extrapolated to be less than 300 mg on the horizontal areas of AUG.

99 GENERAL AND MISCELLANEOUS↗

MENT-Flow: maximum-entropy phase space tomography using normalizing flows

Generative models can be trained to reproduce low-dimensional projections of high-dimensional phase space distributions. Normalizing flows are generative models that parameterize invertible transformations, allowing exact probability density evaluation and sampling. Consequently, flows are unbiased entropy estimators and could be used to solve the high-dimensional maximum-entropy tomography (MENT) problem. In this work, we evaluate a flow-based MENT solver (MENT-Flow) against exact maximum-entropy solutions and Minerbo's iterative MENT algorithm in two dimensions.

Hoover, Austin↗

Comparison of beam emission spectroscopy and gas puff imaging edge fluctuation measurements in National Spherical Torus Experiment

In this study, the close physical proximity of the Gas Puff Imaging (GPI) and Beam Emission Spectroscopy (BES) diagnostics on the National Spherical torus Experiment (NSTX) is leveraged to directly compare fluctuation measurements, and to study the local effects of the GPI neutral deuterium puff during H-mode plasmas without large Edge Localized Modes. The GPI and BES views on NSTX provide partially overlapping coverage of the edge and scrape-off layer (SOL) regions above the outboard midplane. The separation in the toroidal direction is 16°, and field lines passing through diagnostic views are separated by ∼20 cm in the direction perpendicular to the magnetic field. Strong cross-correlation is observed, and strong cross-coherence is seen for frequencies between 5 and 15 kHz. Also, probability distribution functions of fluctuations measured ∼3 cm inside the separatrix exhibit only minor deviations from a normal distribution for both diagnostics, and good agreement between correlation length estimates, decorrelation times, and structure velocities is found at the ±40% level. While the two instruments agree closely in many respects, some discrepancies are observed. Most notably, GPI normalized fluctuation levels exceed BES fluctuations by a factor of ∼9. BES mean intensity is found to be sensitive to the GPI neutral gas puff, and BES normalized fluctuation levels for frequencies between 1 and 10 kHz are observed to increase during the GPI puff.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The evolution of galaxy shapes in CANDELS: from prolate to discy

In this study, we model the projected b/a – log a distributions of CANDELS star-forming main-sequence galaxies, where a (b) is the half-light semimajor (semiminor) axis of the galaxy images measured by GALFIT . We find that smaller a galaxies are rounder at all stellar masses M* and redshifts, so we include a when analysing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate normal distribution of galaxy size and two shape parameters, we construct their intrinsic shape and size distributions to obtain the fractions of elongated (prolate), discy (oblate), and spheroidal galaxies in each redshift and mass bin. We find that galaxies tend to be prolate at low M* and high redshifts, and discy at high M* and low redshifts, qualitatively consistent with van der Wel et al., implying that galaxies tend to evolve from prolate to discy. These results are consistent with the predictions from simulations that the transition from prolate to oblate is caused by a compaction event at a characteristic mass range, making the galaxy centre baryon dominated. We give probabilities of a galaxy’s being elongated, discy, or spheroidal as a function of its M* , redshift, and projected b/a and a , which can facilitate target selections of galaxies with specific shapes at high redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗

An Exact and Near-Exact Distribution Approach to the Behrens–Fisher Problem

The Behrens–Fisher problem occurs when testing the equality of means of two normal distributions without the assumption that the two variances are equal. This paper presents approaches based on the exact and near-exact distributions for the test statistic of the Behrens–Fisher problem, depending on different combinations of even or odd sample sizes. We present the exact distribution when both sample sizes are odd and the near-exact distribution when one or both sample sizes are even. The near-exact distributions are based on a finite mixture of generalized integer gamma (GIG) distributions, used as an approximation to the exact distribution, which consists of an infinite series. The proposed tests, based on the exact and the near-exact distributions, are compared with Welch’s t-test through Monte Carlo simulations, in particular for small and unbalanced sample sizes. The results show that the proposed approaches are competent solutions to the Behrens–Fisher problem, exhibiting precise sizes and better powers than Welch’s approach for those cases. Numerical studies show that the Welch’s t-test tends to be a bit more conservative than the test statistics based on the exact or near-exact distribution, in particular when sample sizes are small and unbalanced, situations in which the proposed exact or near-exact distributions obtain higher powers than Welch’s t-test.

Hong, Serim↗

Long- and short-term temporal variability in cloud condensation nuclei spectra over a wide supersaturation range in the Southern Great Plains site

Abstract. When aerosol particles seed the formation of liquid water droplets in the atmosphere, they are called cloud condensation nuclei (CCN). Different aerosols will act as CCN under different degrees of water supersaturation (relative humidity above 100 %), depending on their size and composition. In this work, we build and analyze a best-estimate CCN spectrum product, tabulated at ∼ 45 min resolution, generated using high quality data from seven independent instruments at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site. The data product spans a large supersaturation range, from 0.0001 % to ∼ 30 %, and time period of 5 years, from 2009–2013, and is available on the ARM data archive. We leverage this added statistical power to examine relationships that are unclear in smaller datasets. Our analysis is performed in three main areas. First, probability distributions of many aerosol and CCN metrics are found to exhibit skewed log-normal distribution shapes. Second, clustering analyses of CCN spectra reveal that the primary drivers of CCN differences are aerosol number size distributions, rather than hygroscopicity or composition, especially at supersaturations above 0.2 %, while also allowing for a simplified understanding of seasonal and diurnal variations in CCN behavior. The predictive ability of using limited hygroscopicity data with accurate number size distributions to estimate CCN spectra is investigated, and the uncertainties of this approach are estimated. Third, the dynamics of CCN spectral clusters and concentrations are examined with cross-correlation and autocorrelation analyses. We find that CCN concentrations change rapidly on the timescale of 1–3 h, with some conservation beyond that which is greatest for the lower supersaturation region of the spectrum.

54 ENVIRONMENTAL SCIENCES↗

Second-order statistical bootstrap for the uncertainty quantification of time-temperature-superposition analysis

Time-temperature superposition (TTS), which for decades has been a powerful method for long-term prediction from accelerated aging data, involves rigid-shifting isotherms in logarithmic time to produce a single master prediction curve. For simple thermo-rheological properties that accurately follow the TTS principle, the shifts can be easily determined, even manually by the eye. However, for many properties of interest, where the principle is obeyed only approximately, or the data is noisy, it is imperative to develop objective shifting techniques along with reliable uncertainty bounds. This work analyzes in detail the method of arclength-minimization as an unsupervised algorithm to determining optimum shifts and demonstrates that the method is nearly unbiased for all practical datasets with a variety of noise distributions. Moreover, if averaged over with-replacement (bootstrap) resamples, the predicted shifts follow a normal distribution, a fact that can be used to construct confidence interval for the master curve through a second-order bootstrap procedure.

36 MATERIALS SCIENCE↗

MO-FG-CAMPUS-JeP3-01: A Statistical Model for Analyzing the Rotational Error of Single Iso-Center Technique

Purpose: To develop a generalized statistical model that incorporates the treatment uncertainty from the rotational error of single iso-center technique, and calculate the additional PTV (planning target volume) margin required to compensate for this error. Methods: The random vectors for setup and additional rotation errors in the three-dimensional (3D) patient coordinate system were assumed to follow the 3D independent normal distribution with zero mean, and standard deviations σx, σy, σz, for setup error and a uniform σR for rotational error. Both random vectors were summed, normalized and transformed to the spherical coordinates to derive the chi distribution with 3 degrees of freedom for the radical distance ρ. PTV margin was determined using the critical value of this distribution for 0.05 significant level so that 95% of the time the treatment target would be covered by ρ. The additional PTV margin required to compensate for the rotational error was calculated as a function of σx, σy, σz and σR. Results: The effect of the rotational error is more pronounced for treatments that requires high accuracy/precision like stereotactic radiosurgery (SRS) or stereotactic body radiotherapy (SBRT). With a uniform 2mm PTV margin (or σx =σy=σz=0.7mm), a σR=0.32mm will decrease the PTV coverage from 95% to 90% of the time, or an additional 0.2mm PTV margin is needed to prevent this loss of coverage. If we choose 0.2 mm as the threshold, any σR>0.3mm will lead to an additional PTV margin that cannot be ignored, and the maximal σR that can be ignored is 0.0064 rad (or 0.37°) for iso-to-target distance=5cm, or 0.0032 rad (or 0.18°) for iso-to-target distance=10cm. Conclusions: The rotational error cannot be ignored for high-accuracy/-precision treatments like SRS/SBRT, particularly when the distance between the iso-center and target is large.

60 APPLIED LIFE SCIENCES↗

Slip band characteristics in the presence of grain boundaries in nickel-based superalloy

Shear strain profiles along slip bands in a modified Rolls-Royce nickel superalloy (RR1000) were analyzed for a tensile sample deformed by 2%. The strain increased with distance away from a grain boundary (GB), with maximum shear strain towards the center of the grain, indicating that dislocation nucleation generally occurred in the grain interior. The strain gradients in the neighborhood of the GBs were quantified and generally correlated with rotation about the active slip system line direction. This leads to an ability to determine the active slip system in these regions. The dislocation spacing and pileup stresses were inferred. The dislocation spacing closely follows an Eshelby analytical solution for a single ended pileup of dislocations under an applied stress. The distribution of pileup stress values for GBs of a given misorientation angle follows a log-normal distribution, with no correlation between the pileup stress and the GB misorientation angle. Furthermore, there is no observed correlation between various transmissivity factors and slip band pileup stress. Hence it appears that the obstacle strength of any of the observed GBs is adequate to facilitate the dislocation pileups present in the slip bands. However, slip band transmission does correlate with transmissivity factors, with the current study focusing on the Luster and Morris m’-factor. Observation of strain profiles of transmitted bands indicate dislocation nucleation locations.

36 MATERIALS SCIENCE↗

Conductance distribution near the Anderson transition

Using a modification of the Shapiro approach, we introduce the two-parameter family of conductance distributions W(g), defined by simple differential equations, which are in the one-to-one correspondence with conductance distributions for quasi-one-dimensional systems of size L{sup d–1} × L{sub z}, characterizing by parameters L/ξ and L{sub z}/L (ξ is the correlation length, d is the dimension of space). This family contains the Gaussian and log-normal distributions, typical for the metallic and localized phases. For a certain choice of parameters, we reproduce the results for the cumulants of conductance in the space dimension d = 2 + ϵ obtained in the framework of the σ-model approach. The universal property of distributions is existence of two asymptotic regimes, log-normal for small g and exponential for large g. In the metallic phase they refer to remote tails, in the critical region they determine practically all distribution, in the localized phase the former asymptotics forces out the latter. A singularity at g = 1, discovered in numerical experiments, is admissible in the framework of their calculational scheme, but related with a deficient definition of conductance. Apart of this singularity, the critical distribution for d = 3 is well described by the present theory. One-parameter scaling for the whole distribution takes place under condition, that two independent parameters characterizing this distribution are functions of the ratio L/ξ.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fitting FFT-derived spectra: Theory, tool, and application to solar radio spike decomposition

Spectra derived from fast Fourier transform (FFT) analysis of time-domain data intrinsically contain statistical fluctuations whose distribution depends on the number of accumulated spectra contributing to a measurement. The tail of this distribution, which is essential for separating the true signal from the statistical fluctuations, deviates noticeably from the normal distribution for a finite number of accumulations. In this paper, we develop a theory to properly account for the statistical fluctuations when fitting a model to a given accumulated spectrum. The method is implemented in software for the purpose of automatically fitting a large body of such FFT-derived spectra. We apply this tool to analyze a portion of a dense cluster of spikes recorded by our FASR Subsystem Testbed instrument during a record-breaking event that occurred on 2006 December 6. The outcome of this analysis is briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗