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 577 records · Page 32

Missing matter in the vicinity of the sun

The Poisson and Vlasov equations are solved numerically for realistic Galaxy models which include multiple disk components, a Population II spheroid, and an unseen massive halo. The total amount of matter in the vicinity of the sun is determined by comparing the observed distributions of tracer stars, samples of F dwarfs, and K giants with the predictions of the Galaxy models. Results are obtained for a number of different assumed distributions of the unseen disk mass. For all the observed samples, typical models imply that about half of the mass in the solar vicinity must be in the form of unobserved matter. The volume density of unobserved material near the sun is about 0.1 solar mass/cu pc; the corresponding column density is about 30 solar mass/sq pc. This so far unseen material must be in a disk with an exponential scale height of less than 0.7 kpc.

Bahcall, John N.↗

An Empirical State Error Covariance Matrix Orbit Determination Example

State estimation techniques serve effectively to provide mean state estimates. However, the state error covariance matrices provided as part of these techniques suffer from some degree of lack of confidence in their ability to adequately describe the uncertainty in the estimated states. A specific problem with the traditional form of state error covariance matrices is that they represent only a mapping of the assumed observation error characteristics into the state space. Any errors that arise from other sources (environment modeling, precision, etc.) are not directly represented in a traditional, theoretical state error covariance matrix. First, consider that an actual observation contains only measurement error and that an estimated observation contains all other errors, known and unknown. Then it follows that a measurement residual (the difference between expected and observed measurements) contains all errors for that measurement. Therefore, a direct and appropriate inclusion of the actual measurement residuals in the state error covariance matrix of the estimate will result in an empirical state error covariance matrix. This empirical state error covariance matrix will fully include all of the errors in the state estimate. The empirical error covariance matrix is determined from a literal reinterpretation of the equations involved in the weighted least squares estimation algorithm. It is a formally correct, empirical state error covariance matrix obtained through use of the average form of the weighted measurement residual variance performance index rather than the usual total weighted residual form. Based on its formulation, this matrix will contain the total uncertainty in the state estimate, regardless as to the source of the uncertainty and whether the source is anticipated or not. It is expected that the empirical error covariance matrix will give a better, statistical representation of the state error in poorly modeled systems or when sensor performance is suspect. In its most straight forward form, the technique only requires supplemental calculations to be added to existing batch estimation algorithms. In the current problem being studied a truth model making use of gravity with spherical, J2 and J4 terms plus a standard exponential type atmosphere with simple diurnal and random walk components is used. The ability of the empirical state error covariance matrix to account for errors is investigated under four scenarios during orbit estimation. These scenarios are: exact modeling under known measurement errors, exact modeling under corrupted measurement errors, inexact modeling under known measurement errors, and inexact modeling under corrupted measurement errors. For this problem a simple analog of a distributed space surveillance network is used. The sensors in this network make only range measurements and with simple normally distributed measurement errors. The sensors are assumed to have full horizon to horizon viewing at any azimuth. For definiteness, an orbit at the approximate altitude and inclination of the International Space Station is used for the study. The comparison analyses of the data involve only total vectors. No investigation of specific orbital elements is undertaken. The total vector analyses will look at the chisquare values of the error in the difference between the estimated state and the true modeled state using both the empirical and theoretical error covariance matrices for each of scenario.

Frisbee, Joseph H., Jr.↗

Harnessing distributed GPU computing for generalizable graph convolutional networks in power grid reliability assessments

Although machine learning (ML) has emerged as a powerful tool for rapidly assessing grid contingencies, prior studies have largely considered a static grid topology in their analyses. This limits their application, since they need to be re-trained for every new topology. Here, this paper explores the development of generalizable graph convolutional network (GCN) models by pre-training them across a range of grid topologies and contingency types. We found that a GCN model with auto-regressive moving average (ARMA) layers with a line graph representation of the grid offered the best predictive performance in predicting voltage magnitudes (VM) and voltage angles (VA). We introduced the concept of phantom nodes to consider disparate grid topologies with a varying number of nodes and lines. For pre-training the GCN ARMA model across a variety of topologies, distributed graphics processing unit (GPU) computing afforded us significant training scalability. The predictive performance of this model on grid topologies that were part of the training data is substantially better than the direct current (DC) approximation. Although direct application of the pre-trained model to topologies that are not part of the grid is not particularly satisfactory, fine-tuning with small amounts of data from a specific topology of interest significantly improves predictive performance. In general, this paper highlights the feasibility of training large-scale GNN models to assess the reliability of power grids by considering a wide variety of grid topologies and contingency types. With the advent of foundational models in ML and the exponential increase in GPU computing clusters, generalizable ML models will significantly enhance how utilities manage power systems and make decisions in real-time or near-real-time.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. I. Stochastic Fragmentation with Optical-depth-dependent Rate and Universal Fragment Mass

Fragmentation in a gravitationally unstable accretion disk can be an important pathway for forming stellar/planetary companions. To characterize quantitatively the condition for and outcome of fragmentation under realistic thermodynamics, we perform global 3D simulations of gravitationally unstable disks at various cooling rates and cooling types, including the first global simulations of gravitational instability that employ full radiation transport. We find that fragmentation is a stochastic process, with the fragment generation rate per disk area p frag showing an exponential dependence on the parameter β ≡ Ω K t cool , where Ω K is the Keplerian rotation frequency, and t cool is the average cooling timescale. Compared to a prescribed constant β, radiative cooling in the optically thin/thick regime makes p frag decrease slower/faster in β; the critical β corresponding to ∼1 fragment per orbit is ≈3, 5, and 2 for constant β, optically thin, and optically thick cooling, respectively. The distribution function of the initial fragment mass is remarkably insensitive to disk thermodynamics. Regardless of cooling rate and optical depth, the typical initial fragment mass is m frag ≈ 40M tot h 3 , with M tot being the total (star+disk) mass and h = H/R being the disk aspect ratio. Applying this result to typical Class 0/I protostellar disks, we find m frag ∼ 20M J , suggesting that fragmentation more likely forms brown dwarfs. Given the finite width of the m frag distribution, forming massive planets is also possible.

gravitational instability↗

Superbubble blowout dynamics

Multiple supernovae and stellar winds from OB associations carve large holes filled with hot gas in the galactic disk. These superbubbles sweep up H I into cold, thin, dense shells and eventually grow large enough to blow completely out of the galactic H I disk. When superbubbles blow out of the disk, they vent hot gas and supernova energy into the galactic corona. In this paper ZEUS, a two-dimensional hydrodynamics code, is used to model the blowout of a superbubble from exponential and Gaussian models for the vertical density stratification. The results are compared to those from the Kompaneets (thin-shell) approximation. It is found that this approximation works very well, and that most of the mass of the shell remains in the plane, with 5 percent of it accelerating upward. The venting of the hot gas and the stability of the shell depends strongly on the model of the density distribution. It is suggested that the low galactic halo actually consists of a froth of merged superbubbles.

Mac Low, Mordecai-Mark↗

Real-Time System Verification by Kappa-Induction

We report the first formal verification of a reintegration protocol for a safety-critical, fault-tolerant, real-time distributed embedded system. A reintegration protocol increases system survivability by allowing a node that has suffered a fault to regain state consistent with the operational nodes. The protocol is verified in the Symbolic Analysis Laboratory (SAL), where bounded model checking and decision procedures are used to verify infinite-state systems by k-induction. The protocol and its environment are modeled as synchronizing timeout automata. Because k-induction is exponential with respect to k, we optimize the formal model to reduce the size of k. Also, the reintegrator's event-triggered behavior is conservatively modeled as time-triggered behavior to further reduce the size of k and to make it invariant to the number of nodes modeled. A corollary is that a clique avoidance property is satisfied.

Pike, Lee S.↗

The analytic two-loop soft function for leading-jet p T

We present the calculation of the two-loop soft function for the transverse momentum distribution of the leading jet produced in association with any colour-singlet system (e.g. a Higgs or a Z boson). This constitutes a central ingredient for the resummation of the above distribution as well as the jet-vetoed cross section at the next-to-next-to-next-to-leading logarithmic order, both of which play an important role in the precision physics programme at the Large Hadron Collider. The calculation is performed in soft-collinear effective theory with an appropriate regularisation of the rapidity divergences that occur in the phase-space integrals. We obtain analytic results by employing an exponential regulator and by taking a Laurent expansion in the jet radius R . All expressions are presented as supplementary material attached to this article.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dark matter near the sun

The amount of dark matter in the disk of the Galaxy at the solar position is determined by comparing the observed distributions of tracer stars with the predictions obtained from different assumptions of how the unseen matter is distributed. The major uncertainties, observational and theoretical, are estimated. For all the observed samples, typical models imply that about half of the mass in the solar vicinity must be in the form of unobserved matter. The volume density of unobserved material near the sun is about 0.1 solar mass/cu pc; the corresponding column density is about 30 solar masses/cu pc. This, so far unseen, material must be in a disk with an exponential scale height of less than 0.7 kpc. All the existing observations are consistent with the unseen disk material being in the form of stars not massive enough to burn hydrogen. It is suggested that the unseen material that is required to hold up the rotation curves of galaxies and to satisfy the virial theorem for clusters of galaxies might also be in the form of low-mass stars.

Bahcall, J. N.↗

A New Approach to Distributed Hypothesis Testing and Non-Bayesian Learning: Improved Learning Rate and Byzantine-Resilience

Here, we study a setting where a group of agents, each receiving partially informative private signals, seek to collaboratively learn the true underlying state of the world (from a finite set of hypotheses) that generates their joint observation profiles. To solve this problem, we propose a distributed learning rule that differs fundamentally from existing approaches, in that it does not employ any form of “belief-averaging”. Instead, agents update their beliefs based on a min-rule. Under standard assumptions on the observation model and the network structure, we establish that each agent learns the truth asymptotically almost surely. As our main contribution, we prove that with probability 1, each false hypothesis is ruled out by every agent exponentially fast, at a network-independent rate that is strictly larger than existing rates. We then develop a computationally-efficient variant of our learning rule that is provably resilient to agents who do not behave as expected (as represented by a Byzantine adversary model) and deliberately try to spread misinformation.

42 ENGINEERING↗

An analytic model of the neutral cloud evolution in the earth's atmosphere

An analytic model for the expansion of a chemically reacting gas in the upper atmosphere is discussed. The spatial-temporal distribution of released particles in the initial regime is described on the basis of self-similar solutions of Euler's equations. In the case of transition to the diffusion regime an approximate solution of the kinetic equation with the collision integral in the form of BGK is performed. Gravitation and the atmospheric inhomogeneity are taken into account, The diffusion regime is described by an analytic solution of the diffusion equation in an exponential atmosphere taking account of possible losses of the gas due to chemical reactions. Some peculiarities of the expansion of various gases as well as the possibility of applying the model for describing ionospheric 'holes' are discussed.

Bruskin, L. G.↗

Laser Disdrometer Quantities (LDQUANTS) and Video Disdrometer Quantities (VDISQUANTS) Value-Added Products Report

Disdrometers are useful instruments that measure the distribution of raindrop sizes (drop size distribution: DSD) and the associated rainfall rates/accumulation as those raindrops fall to the ground. These DSDs are a quantity of interest to members of the modeler and observational communities. However, use of disdrometer data for model evaluation, radar monitoring, or other activities requires careful quality control and processing for key DSD properties of interest (e.g., the number concentration of drops) to ensure appropriate physical (scattering, fall speed) assumptions. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Laser Disdrometer Quantities and Video Disdrometer Quantities Value-Added Product (LDQUANTS/VDISQUANTS VAP) uses standard methods from Tokay et al. (e.g., 2013, 2014) to filter drops with unrealistic fall speeds. Further, it estimates several microphysical/geophysical quantities of parameterized DSDs (gamma or exponential assumption type fitting methods) as in previous disdrometer studies and ARM long-term efforts (e.g., Testud et al. 2001, Giangrande et al. 2014, Thompson et al. 2015). Disdrometers are also beneficial for cross-checks with other instrumentation, including rain gauges and radars. To support research interests and related radar monitoring activities, this product calculates radar equivalent quantities, including dual polarization radar quantities (e.g., Z, Differential Reflectivity ZDR, etc.), using T-Matrix scattering and additional wavelength, temperature, and drop shape assumptions (e.g., Thurai et al. 2007). All of these efforts allow disdrometer data sets to become more useful and easily handled by modeling and observational studies, or routine instrumentation checks.

47 OTHER INSTRUMENTATION↗

X-ray study of solar plage regions and a small flare.

The O VII and Ne IX resonance line and 5 keV continuum emission of a small solar flare are interpreted on the basis of a model in which the emission measure of the region is an exponentially decreasing function of temperature. Good agreement with all three sets of data is obtained - a result not possible on the basis of an isothermal approximation. The change in the emission of the flare over a two minute interval appears to be consistent with a simple cooling of the region. Computation of the contribution function versus temperature for each emission indicates that each predominantly comes from a physically separate volume within the feature. The emission measure distribution for the flare is distinctly flatter than that of five plage regions which tend to have similar slopes and which are in turn flatter than the general corona.

Acton, L. W.↗

Measurements and Modeling of Soot Formation and Radiation in Microgravity Jet Diffusion Flames

This is a computational and experimental study for soot formation and radiative heat transfer in jet diffusion flames under normal gravity (1-g) and microgravity (0-g) conditions. Instantaneous soot volume fraction maps are measured using a full-field imaging absorption technique developed by the authors. A compact, self-contained drop rig is used for microgravity experiments in the 2.2-second drop tower facility at NASA Lewis Research Center. On modeling, we have coupled flame structure and soot formation models with detailed radiation transfer calculations. Favre-averaged boundary layer equations with a k-e-g turbulence model are used to predict the flow field, and a conserved scalar approach with an assumed Beta-pdf are used to predict gaseous species mole fraction. Scalar transport equations are used to describe soot volume fraction and number density distributions, with formation and oxidation terms modeled by one-step rate equations and thermophoretic effects included. An energy equation is included to couple flame structure and radiation analyses through iterations, neglecting turbulence-radiation interactions. The YIX solution for a finite cylindrical enclosure is used for radiative heat transfer calculations. The spectral absorption coefficient for soot aggregates is calculated from the Rayleigh solution using complex refractive index data from a Drude- Lorentz model. The exponential-wide-band model is used to calculate the spectral absorption coefficient for H20 and C02. It is shown that when compared to results from true spectral integration, the Rosseland mean absorption coefficient can provide reasonably accurate predictions for the type of flames studied. The soot formation model proposed by Moss, Syed, and Stewart seems to produce better fits to experimental data and more physically sound than the simpler model by Khan et al. Predicted soot volume fraction and temperature results agree well with published data for a normal gravity co-flow laminar flames and turbulent jet flames. Predicted soot volume fraction results also agree with our data for 1-g and 0-g laminar jet names as well as 1-g turbulent jet flames.

Ku, Jerry C.↗

Entanglement structures in quantum field theories. II. Distortions of vacuum correlations through the lens of local observers

When observing a quantum field via detectors with access to only the mixed states of spatially separated, local regions—a ubiquitous experimental design—the capacity to access the full extent of distributed entanglement can be limited, shrouded by classical correlations. By performing projective measurements of the field external to two detection patches and classically communicating the results, underlying pure states may be identified for which entanglement quantification is clear. In the Gaussian continuous-variable states of the free scalar field vacuum, this protocol uncovers a disparity between the spacelike entanglement established within the field and that which is locally detectable. Further, this discrepancy is found to grow exponentially with the separation between observation regions. The protocol developed herein offers insight and practical guidance for clarifying the unavoidable distortion of quantum field correlations when viewed from the vantage of a pair of local observers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Widespread Nanoflare Variability Detected with Hinode/XRT in a Solar Active Region

It is generally agreed that small impulsive energy bursts called nanoflares are responsible for at least some of the Sun s hot corona, but whether they are the explanation for most of the multi-million degree plasma has been a matter of ongoing debate. We here present evidence that nanoflares are widespread in an active region observed by the X-Ray Telescope on-board the Hinode mission. The distributions of intensity fluctuations have small but important asymmetries, whether taken from individual pixels, multi-pixel subregions, or the entire active region. Negative fluctuations (corresponding to reduced intensity) are greater in number but weaker in amplitude, so that the median fluctuation is negative compared to a mean of zero. Using Monte Carlo simulations, we show that only part of this asymmetry can be explained by Poisson photon statistics. The remainder is explainable with a tendency for exponentially decreasing intensity, such as would be expected from a cooling plasma produced, e.g., from a nanoflare. We suggest that nanoflares are a universal heating process within active regions.

Reale, Fabio↗

Molecular mechanisms of thickness-dependent water desalination in polyamide reverse-osmosis membranes

Here, using non-equilibrium molecular simulations, we systematically elaborate the relationships among synthesis (membrane's thickness and morphology), atomic-scale transport mechanism, and separation performance (permeability and selectivity) for 3D-printed polyamide (PA) membranes. Results indicated that water diffusion, swelling ratio, water flux, and water permeance proportionally decrease with increasing membrane thickness (4.0–32.5 nm). PA membranes with different thicknesses can achieve almost 100% salt rejection over the simulation time. Importantly, water permeability exponentially decreases with increasing thickness, and 15 nm is identified as the critical membrane thickness for efficient water transport. The discontinuous water-available space spreads all over PA membranes with thicknesses greater than 5 nm, allowing water molecules to jump by way of the temporary open-and-close pores. However, the connected water-useable space exists in PA membranes with thicknesses below 5 nm, offering the continuous channels to dominate water transport. More significantly, pore distribution is more homogeneous as the thickness increases. The applied high pressures can lead to membrane compaction during reverse osmosis and the thicker membranes show a lower compression ratio. In short, these investigations provide molecular insights for effectively designing and manufacturing PA membranes for water desalination and treatment at the molecular level.

42 ENGINEERING↗

Analytic solutions of the DGLAP evolution and theoretical uncertainties

The energy dependence for the singlet sector of Parton Distributions Functions (PDFs) is described by an entangled pair of ordinary linear differential equations. Although there are no exact analytic solutions, it is possible to provide approximated results depending on the assumptions and the methodology adopted. These results differ in their sub-leading, neglected terms and ultimately they are associated with different treatments of the theoretical uncertainties. In this work, a novel analytic approach in Mellin space is presented and a new methodology for obtaining closed and exponentiated analytic solutions is devised. Different results for the DGLAP evolution at Next-Leading-Order are compared, discussing advantages and disadvantages for each solution. The generalizations to higher orders are addressed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Neutral hydrogen in the outer Galaxy

The H I 21 cm emission-line survey data of Stark et al. (1991) have been used to study the distribution of H I in the outer Galaxy and to improve understanding of the various systematic problems which influence all studies of Galactic H I. For n(H) greater than 0.01/cu cm, the distributions are similar to those found by Burton and te Lintel Hekkert (1986) but are 'better behaved' at substantially lower values of n(H). The z distribution of gas in the H I layer is best described as complicated. The behavior of the gaseous matter is consistent with the observed distribution of H I found in external spiral galaxies. In the directions of the maximum warp the average distance of the gas away from the Galactic plane reaches about 4 kpc at R about 24 kpc, while the flaring of the layer as measured by z(rms) increases to about 3 kpc. In directions where the warp layers particularly well measured, the Galaxy is followed to R about 30 kpc. The surface density of the atomic gas decreases exponentially.

Diplas, Athanassios↗