Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Exponential 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Evolution of the age-included nearest pair distribution in disperse multiphase flows

The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Sizes of pure and doped helium droplets from single shot x-ray imaging

Advancements in x-ray free-electron lasers on producing ultrashort, ultrabright, and coherent x-ray pulses enable single-shot imaging of fragile nanostructures, such as superfluid helium droplets. Furthermore, this imaging technique gives unique access to the sizes and shapes of individual droplets. In the past, such droplet characteristics have only been indirectly inferred by ensemble averaging techniques. Here, we report on the size distributions of both pure and doped droplets collected from single-shot x-ray imaging and produced from the free-jet expansion of helium through a 5 μm diameter nozzle at 20 bars and nozzle temperatures ranging from 4.2 to 9 K. This work extends the measurement of large helium nanodroplets containing 10 9 –10 11 atoms, which are shown to follow an exponential size distribution. Additionally, we demonstrate that the size distributions of the doped droplets follow those of the pure droplets at the same stagnation condition but with smaller average sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Biophotonic rogue waves in red blood cell suspensions

Rogue waves are ubiquitous in nature, appearing in a variety of physical systems ranging from acoustics, microwave cavities, optical fibers, and resonators to plasmas, superfluids, and Bose–Einstein condensates. Unlike nonlinear solitary waves, rogue waves are extreme events that can occur even without nonlinearity by, for example, spontaneous synchronization of waves with different spatial frequencies in a linear system. Here, we report the observation of rogue-wave-like events in human red blood cell (RBC) suspensions under weak light illumination, characterized by an abnormal L-shaped probability distribution. Such biophotonic extreme events arise mostly due to the constructive interference of Mie-scattered waves from the suspended RBCs, whose biconcave shape and mutable orientation give rise to a time-dependent random phase modulation to an incident laser beam. We numerically simulate the beam propagation through the colloidal suspensions with added disorder in both spatial and temporal domains to mimic random scattering due to Brownian motion. In addition, at high power levels, nonlinear beam self-focusing is also observed, leading to a dual-exponential probability distribution associated with the formation of multiple soliton-like spots. Such rogue wave events should also exist in environments with cells of other species such as swimming bacteria, and understanding of their underlying physics may lead to unexpected biophotonic applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Perturbative readout-error mitigation for near-term quantum computers

Readout errors on near-term quantum computers can introduce significant error to the empirical probability distribution sampled from the output of a quantum circuit. These errors can be mitigated by classical postprocessing given the access of an experimental response matrix that describes the error associated with the measurement of each computational basis state. However, the resources required to characterize a complete response matrix and to compute the corrected probability distribution scale exponentially with the number of qubits, n . In this work, we modify standard matrix inversion techniques using perturbative approximations with significantly reduced complexity and bounded error when the likelihood of high-order bit-flip events is strongly suppressed. Given a characteristic error rate q , we discuss a method to recover the probability of the all-zeros bit string p 0 by sampling only a small subspace of the response matrix before inverting readout error, resulting in a relative speedup of poly [ 2 n / ( n w ) ] , which we motivate using a simplified error model for which the approximation incurs only O ( q w ) error for some integer w . We then provide a generalized technique to efficiently recover full output distributions with O ( q w ) error in the perturbative limit. These approximate techniques for readout-error correction may greatly accelerate near-term quantum computing applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Counting generations in birth and death processes with competing Erlang and exponential waiting times

Lymphocyte populations, stimulated in vitro or in vivo, grow as cells divide. Stochastic models are appropriate because some cells undergo multiple rounds of division, some die, and others of the same type in the same conditions do not divide at all. If individual cells behave independently, then each cell can be imagined as sampling from a probability density of times to division and death. The exponential density is the most mathematically and computationally convenient choice. It has the advantage of satisfying the memoryless property, consistent with a Markov process, but it overestimates the probability of short division times. With the aim of preserving the advantages of a Markovian framework while improving the representation of experimentally-observed division times, we consider a multi-stage model of cellular division and death. We use Erlang-distributed (or, more generally, phase-type distributed) times to division, and exponentially distributed times to death. We classify cells into generations, using the rule that the daughters of cells in generation n are in generation n + 1. In some circumstances, our representation is equivalent to established models of lymphocyte dynamics. We find the growth rate of the cell population by calculating the proportions of cells by stage and generation. The exponent describing the late-time cell population growth, and the criterion for extinction of the population, differs from what would be expected if N steps with rate λ were equivalent to a single step of rate λ/N. We link with a published experimental dataset, where cell counts were reported after T cells were transferred to lymphopenic mice, using Approximate Bayesian Computation. In the comparison, the death rate is assumed to be proportional to the generation and the Erlang time to division for generation 0 is allowed to differ from that of subsequent generations. The multi-stage representation is preferred to a simple exponential in posterior distributions, and the mean time to first division is estimated to be longer than the mean time to subsequent divisions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

High Resolution Data Analysis: Plans and Prospects [Book Chapter]

Herein, a report on the progress on the high resolution data analysis of the ADMX experimental results is presented. In this paper, tools are developed and tested on a blind injection mimicking a Maxwellian like signal in the frequency domain. This blind injection will be used as a test bed which can be later implemented on all the high resolution data. The high resolution data is stored in the Fermilab server. In this analysis a PostgreSQL query was made to ensure the blind injection is in the middle of the frequency spectrum and 19 such files were found. The time series data is read using a c++ program. An apodization function is applied on the time series data and zero filled to reduce the frequency spacing in order to achieve a better interpolation. A FFTW header is used to compute the Fourier transform of the time series data. A Savitzky–Golay filter is applied on the unnormalized power which then can be used to remove the spectral shape. Each frequency spectrum has a bandwidth of 50 kHz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Reducing dielectric loss and enhancing electrical insulation for multilayer polymer films by nanoconfined ion transport under high poling electric fields

High temperature polar polymers have demonstrated potential for good thermal stability and high dielectric constant at the same time. However, polarization of contaminated impurity ions in polar polymers, even at the ppm level, can significantly increase the dielectric loss at high temperature and low frequencies. One effective strategy to mitigate this problem is to multilayer them with a high temperature nonpolar dielectric polymer to confine impurity ion transport at the nanometer scale. In this study, confined ion transport in high temperature polycarbonate (HTPC)/poly(vinylidene fluoride) (PVDF) multilayer films under high AC electric fields was studied using a direct analytical simulation method. Different from the ion transport under low fields, the ion diffusion model failed to describe the ion transport under high electric fields. An exponential ion distribution profile, which was observed for the DC poling situation, was employed to implement the direct analytical simulation. Confined impurity ion transport under high AC fields was quantitatively understood. As the AC field increased, the mobile ion concentration decreased whereas the diffusion coefficient increased. The decrease of mobile ion concentration was explained by the blockage of impurity ions by the HTPC layers. This knowledge helped in the determination of optimal conditions to polarize impurity ions from the PVDF layers into the HTPC layers. After cooling below the glass transition temperature of HTPC, polarized impurity ions were locked inside the HTPC layers. As a result, increased discharge efficiency and enhanced electrical insulation (i.e., increased dielectric breakdown strength) were achieved for the polarized multilayer films.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A quantum algorithm to efficiently sample from interfering binary trees

Quantum computers provide an opportunity to efficiently sample from probability distributions that include non-trivial interference effects between amplitudes. Using a simple process wherein all possible state histories can be specified by a binary tree, we construct an explicit quantum algorithm that runs in polynomial time to sample from the process once. The corresponding naive Markov Chain algorithm does not produce the correct probability distribution and an explicit classical calculation of the full distribution requires exponentially many operations. Finally, the problem can be reduced to a system of two qubits with repeated measurements, shedding light on a quantum-inspired efficient classical algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The Cosmic Evolution of C IV Absorbers at 1.4 < z < 4.5: Insights from 100,000 Systems in DESI Quasars

We present the largest catalog to date of triply ionized carbon (C IV ) absorbers detected in quasar spectra from the Dark Energy Spectroscopic Instrument. Using an automated matched-kernel convolution method with adaptive signal-to-noise thresholds, we identify 101,487 C IV systems in the redshift range 1.4 < z < 4.5 from 300,637 quasar spectra. Completeness is estimated via Monte Carlo simulations, and the catalog is 50% complete at EW C IV ≥ 0.4 Å. The differential equivalent width frequency distribution declines exponentially and shows weak redshift evolution. The absorber incidence per unit comoving path increases by a factor of 2–5 from z ≈ 4.5 to z ≈ 1.4, with stronger redshift evolution for strong systems. Using column densities derived from the apparent optical depth method, we constrain the cosmic mass density of C IV , Ω C IV , which increases by a factor of ∼3.8 from (0.82 ± 0.05) × 10 −8 at z ≈ 4.5 to (3.16 ± 0.2) × 10 −8 at z ≈ 1.4. From Ω C IV , we estimate a lower limit on intergalactic medium metallicity ${\mathrm{log}}({Z}_{{\rm{IGM}}}/{Z}_{\odot })\gtrsim -3.25$ at z ∼ 2.3, with a smooth decline at higher redshifts. These trends trace the cosmic star formation history and He II photoheating rate, suggesting a link between C IV enrichment, star formation, and UV background over ∼3 Gyr. The catalog also provides a critical resource for future studies connecting circumgalactic metals to galaxy evolution, especially near cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Thixotropic spectra and Ashby-style charts for thixotropy

There is no universal model for thixotropy, and comparing thixotropic effects between different fluids is a subtle yet challenging problem. We introduce a generalized (model-insensitive) framework for comparing thixotropic properties based on thixotropic spectra. A superposition of exponential stress modes distributed over thixotropic time scales is used to quantify buildup and breakdown times and mode strengths in response to step-change input. This mathematical framework is tested with several experimental step-shear rate data on colloidal suspensions. Low-dimensional metrics based on moments of the distribution reveal characteristic average thixotropic properties, which are visualized on Ashby-style diagrams. This method outlines a framework for describing thixotropy across a diverse range of microstructures, supporting scientific studies as well as material selection for engineering design applications.

Mechanics↗
DOE OSTI2020

Integrating high resolution drone imagery and forest inventory to distinguish canopy and understory trees and quantify their contributions to forest structure and dynamics

Tree growth and survival differ strongly between canopy trees (those directly exposed to overhead light), and understory trees. However, the structural complexity of many tropical forests makes it difficult to determine canopy positions. The integration of remote sensing and ground-based data enables this determination and measurements of how canopy and understory trees differ in structure and dynamics. Here we analyzed 2 cm resolution RGB imagery collected by a Remotely Piloted Aircraft System (RPAS), also known as drone, together with two decades of bi-annual tree censuses for 2 ha of old growth forest in the Central Amazon. We delineated all crowns visible in the imagery and linked each crown to a tagged stem through field work. Canopy trees constituted 40% of the 1244 inventoried trees with diameter at breast height (DBH) > 10 cm, and accounted for ~70% of aboveground carbon stocks and wood productivity. The probability of being in the canopy increased logistically with tree diameter, passing through 50% at 23.5 cm DBH. Diameter growth was on average twice as large in canopy trees as in understory trees. Growth rates were unrelated to diameter in canopy trees and positively related to diameter in understory trees, consistent with the idea that light availability increases with diameter in the understory but not the canopy. The whole stand size distribution was best fit by a Weibull distribution, whereas the separate size distributions of understory trees or canopy trees > 25 cm DBH were equally well fit by exponential and Weibull distributions, consistent with mechanistic forest models. The identification and field mapping of crowns seen in a high resolution orthomosaic revealed new patterns in the structure and dynamics of trees of canopy vs. understory at this site, demonstrating the value of traditional tree censuses with drone remote sensing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Dimension Reduction and Redundancy Removal through Successive Schmidt Decompositions

Quantum computers are believed to have the ability to process huge data sizes, which can be seen in machine learning applications. In these applications, the data, in general, are classical. Therefore, to process them on a quantum computer, there is a need for efficient methods that can be used to map classical data on quantum states in a concise manner. On the other hand, to verify the results of quantum computers and study quantum algorithms, we need to be able to approximate quantum operations into forms that are easier to simulate on classical computers with some errors. Motivated by these needs, in this paper, we study the approximation of matrices and vectors by using their tensor products obtained through successive Schmidt decompositions. We show that data with distributions such as uniform, Poisson, exponential, or similar to these distributions can be approximated by using only a few terms, which can be easily mapped onto quantum circuits. The examples include random data with different distributions, the Gram matrices of iris flower, handwritten digits, 20newsgroup, and labeled faces in the wild. Similarly, some quantum operations, such as quantum Fourier transform and variational quantum circuits with a small depth, may also be approximated with a few terms that are easier to simulate on classical computers. Furthermore, we show how the method can be used to simplify quantum Hamiltonians: In particular, we show the application to randomly generated transverse field Ising model Hamiltonians. The reduced Hamiltonians can be mapped into quantum circuits easily and, therefore, can be simulated more efficiently.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantifying Quantum Chaos through Microcanonical Distributions of Entanglement

A characteristic feature of “quantum chaotic” systems is that their eigenspectra and eigenstates display universal statistical properties described by random matrix theory (RMT). However, eigenstates of local systems also encode structure beyond RMT. To capture this feature, we introduce a framework that allows us to compare the properties of eigenstates in local systems with those of pure random states. In particular, our framework defines a notion of distance between quantum state ensembles that utilizes the Kullback-Leibler divergence to compare the microcanonical distribution of entanglement entropy (EE) of eigenstates with a reference RMT distribution generated by pure random states (with appropriate constraints). This notion gives rise to a quantitative metric for quantum chaos that not only accounts for averages of the distributions but also higher moments. The differences in moments are compared on a highly resolved scale set by the standard deviation of the RMT distribution, which is exponentially small in system size. As a result, the metric can distinguish between chaotic and integrable behaviors and, in addition, quantify and compare the of chaos (in terms of proximity to RMT behavior) between two systems that are assumed to be chaotic. We implement our framework in local, minimally structured, Floquet random circuits, as well as a canonical family of many-body Hamiltonians, the mixed-field Ising model (MFIM). Importantly, for Hamiltonian systems, we find that the reference random distribution must be appropriately constrained to incorporate the effect of energy conservation in order to describe the ensemble properties of midspectrum eigenstates. The metric captures deviations from RMT across all models and parameters, including those that have been previously identified as strongly chaotic, and for which other diagnostics of chaos such as level spacing statistics look strongly thermal. In Floquet circuits, the dominant source of deviations is the second moment of the distribution, and this persists for all system sizes. For the MFIM, we find significant variation of the KL divergence in parameter space. Notably, we find a small region where deviations from RMT are minimized, suggesting that “maximally chaotic” Hamiltonians may exist in fine-tuned pockets of parameter space. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

A riming‐dependent parameterization of scattering by snowflakes using the self‐similar Rayleigh–Gans approximation

Abstract Riming is a key process of precipitation formation in ice‐containing clouds, but quantifying riming from observations is challenging, limiting our ability to evaluate the riming process in numerical weather models. One challenge for radar observations is that riming changes both the physical properties (mass, area cross‐section) and scattering properties of ice particles. These changes need to be implemented consistently as a function of riming in radar forward operators, which are required for retrievals and model evaluation in observation space. In this study, mass–size, cross‐section area–size, and backscattering cross‐section relations are developed as a function of the normalized rime mass for aggregates composed of various monomer types (columns, dendrites, needles, plates, and rosettes). The proposed framework allows us to simulate scattering properties of aggregated ice particles consistently as a function of riming in retrievals and radar forward operators. The parameterizations are developed from a large data set of simulated rimed aggregates of different sizes and monomer crystal types. The backscattering cross‐section parameterization (the “riming‐dependent parameterization”) is evaluated for radar frequencies of 35.6 and 94.0 GHz and is based on the Self‐Similar Rayleigh–Gans approximation (SSRGA), which is increasingly used to calculate microwave scattering of ice crystals and snowflakes. Compared with parameterizations from the literature that do not consider riming, the riming‐dependent parameterization leads to significantly smaller biases in terms of backscattering cross‐section. When using the particle masses and scattering properties of the individual particles simulated by the aggregation and riming model as a reference, the bias of our parameterization is below 1 dB when integrating over an exponential particle size distribution with sizes from 0.1–10 mm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling of transverse stimulated Raman scattering in KDP/DKDP in large-aperture plates suitable for polarization control

Transverse stimulated Raman scattering (TSRS) in potassium dihydrogen phosphate (KDP) and deuterated potassium dihydrogen phosphate (DKDP) plates for large-aperture, inertial confinement fusion (ICF)-class laser systems is a well-recognized limitation giving rise to parasitic energy conversion and laser-induced damage. The onset of TSRS is manifested in plates exposed to the ultraviolet section of the beam. TSRS amplification is a coherent process that grows exponentially and is distributed nonuniformly in the crystal and at the crystal surfaces. To understand the growth and spatial distribution of TSRS energy in various configurations, a modeling approach has been developed to simulate the operational conditions relevant to ICF-class laser systems. Specific aspects explored in this work include (i) the behavior of TSRS in large-aperture crystal plates suitable for third-harmonic generation and use as wave plates for polarization control in current-generation ICF-class laser system configurations; (ii) methods, and their limitations, of TSRS suppression and (iii) optimal geometries to guide future designs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Parametric description of intermittent probability distribution functions in solar wind and magnetohydrodynamic turbulence

ABSTRACT In this work, we find empirical evidence that the scale-dependent statistical properties of solar wind and magnetohydrodynamic (MHD) turbulence can be described in terms of a family of parametric probability distribution functions (PDFs) known as Normal Inverse Gaussian (NIG). Understanding these PDFs is one of the most important goals in turbulence theory, as they are inherently connected to the intermittent properties of solar wind turbulence. We investigate the properties of PDFs of Elsasser increments based on a large statistical sample from solar wind observations and high-resolution numerical simulations of MHD turbulence. In order to measure the PDFs and their corresponding properties, three experiments are presented: fast and slow solar wind for experimental data and a simulation of reduced MHD (RMHD) turbulence. Conditional statistics on a 23-yr-long sample of WIND data near 1 au and high-resolution pseudo-spectral simulation of steadily driven RMHD turbulence on a $2048^3$ mesh are used to construct scale-dependent PDFs. The empirical PDFs are fitted to NIG distributions, which depend on four free parameters. Our analysis shows that NIG distributions accurately capture the evolution of the PDFs, with scale-dependent parameters, from large scales characterized by a Gaussian distribution, turning to exponential tails within the inertial range and stretched exponentials at dissipative scales. We also show that empirically-measured NIG parameters exhibit well-defined scaling properties that are similar across the three empirical data sets, which may be indicative of universal behaviour.

Astronomy & Astrophysics↗
DOE OSTI2022

Observation of Current-Driven Features of 2.5 Mev Ion Bunch With Complete and Efficient 5D Measurements at the SNS Beam Test Facility

The SNS Beam Test Facility research program is focused detailed studies of beam distributions for medium-energy ion beams, with the goal of reconstructing realistic 6D bunch distributions to enable halo prediction. For complete characterization of the initial distribution, scan time scales exponentially with scan dimension. Currently, a full 6D measurement with ~10 points across most dimensions requires 30 hours. However, measurement of the 5D distribution f(x, x’,y,y’,w) can be done very rapidly using a hybrid slit/screen method. This approach requires ~4 hours to obtain at least 32 points/dimension, with very high resolution (0.5 keV) in the energy distribution. This presentation reports on the approach and results for 5D characterization of the initial RFQ-formed bunch. This includes higher-resolution views of previously reported transverse-longitudinal dependence and additional interplane dependencies that were not previously reported.

Ruisard, Kiersten↗
DOE OSTI2024

Stress due to electric charge density distribution in a dielectric slab

The spatial distribution of electric field due to an imposed electric charge density profile in an infinite slab of dielectric material is derived analytically by integrating Gauss’s law. Various charge density distributions are considered, including exponential and power-law forms. Here, the Maxwell stress tensor is used to compute a notional static stress in the material due to the charge density and its electric field. Characteristics of the electric field and stress distributions are computed for example cases in polyethylene, showing that field magnitudes exceeding the dielectric strength would be required in order to achieve a stress exceeding the ultimate tensile strength.

42 ENGINEERING↗