Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Truncated 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗
NASA NTRS1983

Evaluation of truncation error and adaptive grid generation for the transonic full potential flow calculations

The effects of truncation error on the numerical solution of transonic flows using the full potential equation are studied. The effects of adapting grid point distributions to various solution aspects including shock waves is also discussed. A conclusion is that a rapid change of grid spacing is damaging to the accuracy of the flow solution. Therefore, in a solution adaptive grid application an optimal grid is obtained as a tradeoff between the amount of grid refinement and the rate of grid stretching.

Nakamura, S.↗
DOE OSTI2022

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Explosion Discrimination Using Seismic Gradiometry and Spectral Filtering of Data

Here, we present a new method to discriminate between earthquakes and buried explosions using observed seismic data. The method is different from previous seismic discrimination algorithms in two main ways. First, we use seismic spatial gradients, as well as the wave attributes estimated from them (referred to as gradiometric attributes), rather than the conventional three-component seismograms recorded on a distributed array. The primary advantage of this is that a gradiometer is only a fraction of a wavelength in aperture compared with a conventional seismic array or network. Second, we use the gradiometric attributes as input data into a machine learning algorithm. The resulting discrimination algorithm uses the norms of truncated principal components obtained from the gradiometric data to distinguish the two classes of seismic events. Using high-fidelity synthetic data, we show that the data and gradiometric attributes recorded by a single seismic gradiometer performs as well as a conventional distributed array at the event type discrimination task.

58 GEOSCIENCES↗
NASA NTRS1974

Investigation of lineaments on Skylab and ERTS images of Peninsular Ranges, Southwestern California

The author has identified the following significant results. Northwest trending faults such as the Elsinore and San Jacinto are prominently displayed on Skylab and ERTS images of the Peninsular Ranges, southern California. Northeast, north-south, and west-north-west trending lineaments and faults are also apparent on satellite imagery. Several of the lineaments represent previously unmapped faults. Other lineaments are due to erosion along foliation directions and sharp bends in basement rock contacts rather than faulting. The northeast trending Thing Valley fault appears to be offset by the south branch of the Elsinore fault near Agua Caliente Hot Springs. Larger horizontal displacement along the Elsinore fault further northwest may be distributed along several faults which branch from the Elsinore fault in the Peninsular Ranges. The northeast and west-northwest trending faults are truncated by the major northwest trending faults and appear to be restricted to basement terrane. Limited data on displacement direction suggests that the northeast and west-northwest trending faults formed in response to an earlier period of east-northeast, west-southwest crustal shortening. Such a stress system is consistent with the plate tectonic model of a subduction zone parallel to the continental margin suggested in the late Mesozoic and early Tertiary.

Merifield, P. M.↗
DOE OSTI2023

Flares, Warps, Truncations, and Satellite: The Ultra-thin Galaxy UGC 11859

The structure of the outskirts of galaxies provides valuable information about their past and evolution. Due to their projected orientation, edge-on isolated galaxies effectively serve as test labs in which to study the three-dimensional structures of galaxies, including warps and flares, and to explore the possible sources of such distortions. We analyzed the structure of the apparently isolated edge-on ultra-thin galaxy UGC 11859 to look for the presence of distortions. The deep optical imaging observations (${\mu }_{\mathrm{lim}}=30.6$ and 30.0 $\mathrm{mag}\ {\mathrm{arcsec}}^{-2}$ in the g- and r-bands, respectively) we acquired with the 10.4 m Gran Telescopio Canarias are used to derive the radial and vertical surface brightness profiles and g - r color radial profile. We find that UGC 11859's disk displays a significant gravitational distortion. A warp is clearly detected on one side of the disk, and the galactic plane on both sides of the center shows increasing scale height with an increasing galactocentric radius, indicating the presence of a flare in the stellar distribution. The surface brightness profile of the disk shows a sharp break at 24 kpc galactocentric radius, and a steep decline to larger radii, an "edge-on truncation," which we associate with the presence of the flare. The present study is the first observational support for a connection between truncations and flares. Just beyond the warped side of the disk, a faint galaxy is observed within a small angular distance, identified as a potential interacting companion. Based on ultradeep g and r photometry we estimate that if the potential companion is at the same distance as UGC 11859, the stellar mass of the satellite galaxy is log 10 (M ⊙ ) = ${6.33}_{-0.02}^{+0.02}$ .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

A Gaussian measure of quantum phase noise

We study the width of the semiclassical phase distribution of a quantum state in its dependence on the average number of photons (m) in this state. As a measure of phase noise, we choose the width, delta phi, of the best Gaussian approximation to the dominant peak of this probability curve. For a coherent state, this width decreases with the square root of (m), whereas for a truncated phase state it decreases linearly with increasing (m). For an optimal phase state, delta phi decreases exponentially but so does the area caught underneath the peak: all the probability is stored in the broad wings of the distribution.

Schleich, Wolfgang P.↗
DOE OSTI2022

Generative Ensemble Regression: Learning Particle Dynamics from Observations of Ensembles with Physics-informed Deep Generative Models

Here, we propose a new method for inferring the governing stochastic ordinary differential equations (SODEs) by observing particle ensembles at discrete and sparse time instants, i.e., multiple “snapshots.” Particle coordinates at a single time instant, possibly noisy or truncated, are recorded in each snapshot but are unpaired across the snapshots. By training a physics-informed generative model that generates “fake” sample paths, we aim to fit the observed particle ensemble distributions with a curve in the probability measure space, which is induced from the inferred particle dynamics. We employ different metrics to quantify the differences between distributions, e.g., the sliced Wasserstein distances and the adversarial losses in generative adversarial networks. We refer to this method as generative “ensemble-regression” (GER), in analogy to the classic “point-regression,” where we infer the dynamics by performing regression in the Euclidean space. We illustrate the GER by learning the drift and diffusion terms of particle ensembles governed by SODEs with Brownian motions and Lévy processes up to 100 dimensions. We also discuss how to treat cases with noisy or truncated observations. Apart from systems consisting of independent particles, we also tackle nonlocal interacting particle systems with unknown interaction potential parameters by constructing a physics-informed loss function. Finally, we investigate scenarios of paired observations and discuss how to reduce the dimensionality in such cases by proving a convergence theorem that provides theoretical support.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

Complex Dynamics of Air Traffic Flow

Air traffic in the United States has continued to grow at a steady pace since 1980, except for a dip immediately after the tragic events of September 11, 2001. There are different growth scenarios associated both with the magnitude and the composition of the future air traffic. The Terminal Area Forecast (TAF), prepared every year by the FAA, projects the growth of traffic in the United States. Both Boeing and Airbus publish market outlooks for air travel annually. Although predicting the future growth of traffic is difficult, there are two significant trends: heavily congested major airports continue to see an increase in traffic, and the emergence of regional jets and other smaller aircraft with fewer passengers operating directly between non-major airports. The interaction between air traffic demand and the ability of the system to provide the necessary airport and airspace resources can be modeled as a network. The size of the resulting network varies depending on the choice of its nodes. It would be useful to understand the properties of this network to guide future design and development. Many questions, such as the growth of delay with increasing traffic demand and impact of the en route weather on future air traffic, require a systematic understanding of the properties of the air traffic network. There has been a major advance in the understanding of the behavior of networks with a large number of components. Several theories have been advanced about the evolution of large biological and engineering networks by authors in diversified disciplines like physics, mathematics, biology and computer science. Several networks exhibit a scale-free property in the sense that the probabilistic distribution of their nodes as a function of connections decreases slower than an exponential. These networks are characterized by the fact that a small number of components have a disproportionate influence on the performance of the network. Scale-free networks are tolerant to random failure of components, but are vulnerable to selective attack on components. This paper examines two network representations for the baseline air traffic system. A network defined with the 40 major airports as nodes and with standard flight routes as links has a characteristic scale: all nodes have 60 or more links and no node has more than 460 links. Another network is defined with baseline aircraft routing structure exhibits an exponentially truncated scale-free behavior. Its degree ranges from 2 connections to 2900 connections, and 225 nodes have more than 250 connections. Furthermore, those high-degree nodes are homogeneously distributed in the airspace. A consequence of this scale-free behavior is that the random loss of a single node has little impact, but the loss of multiple high-degree nodes (such as occurs during major storms in busy airspace) can adversely impact the system. Two future scenarios of air traffic growth are used to predict the growth of air traffic in the United States. It is shown that a three-times growth in the overall traffic may result in a ten-times impact on the density of traffic in certain parts of the United States.

Scale-free Networks↗
NASA NTRS1995

Back pain during 6 degrees head-down tilt approximates that during actual microgravity

Astronauts often experience back pain during spaceflight. Retrospectively, Wing et al. (14) found that during spaceflight, 14 of 19 Shuttle crewmembers experienced back pain, which they described as dull (62%), localized to the lower back (50%), and with an intensity of 2 on a 5-point scale. Further, the spine lengthens 4-7 cm in microgravity. Our objective was to compare back pain and spinal lengthening (body height increase) during simulated microgravity (6 degrees head-down tilt, HDT) with the same parameters during actual microgravity. Eight male subjects completed a modified McGill pain questionnaire with intensity graded from zero (no pain) to five (intense and incapacitating pain) each day at 7:00 pm during 2 d pre-HDT control, 16 d HDT, and 1 d post-HDT recovery periods. Also, the subjects' heights were measured each day while supine (control and recovery) and during HDT. Back pain increased from zero (pre-tilt control period) to 2.3 +/- 0.4 at days 1 to 3 of HDT, and was categorized as dull and/or burning pain in subjects' lower backs. Only 2 subjects reported any pain after day 9 of HDT and during recovery. Heights increased 2.1 +/- 0.5 cm by day 3 of HDT and remained at that level until the end of the HDT period. Although spinal lengthening in space is greater than that during HDT, the HDT model approximates the level, type, distribution, and time course of back pain associated with actual microgravity. In the HDT model, pain subsides in intensity when spinal lengthening stops.(ABSTRACT TRUNCATED AT 250 WORDS).

Non-NASA Center↗
DOE OSTI2020

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory (χ EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χ EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. Furthermore, we then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium χ EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

Hybrid-coordinate spacecraft dynamics using large-deformation modal coordinates

For a vehicle amenable to idealization as a set of n + 1 rigid bodies interconnected by n line hinges (implying tree topology), a set of minimum-dimension matrix equations of motion is recorded in terms of variables including the relative rotations of contiguous bodies. When a substructure or portion of the vehicle is subjected to large gross deformations but is restricted to 'small' strains or local deformations, the n-body model is (for sufficiently large n) assumed to experience 'small' relative rotations of contiguous bodies within the substructure, and the equations of motion are linearized in these variables. The application of coordinate transformations and truncations to the linearized subset of the system equations of motion then produces a hybrid-coordinate formulation of the problem involving both discrete and distributed coordinates.

Likins, P. W.↗
NASA NTRS1982

Estimation of distributed parameter systems

The estimation and control of large flexible space structures pose significant new control technology problems. One major problem area is the closed-loop stability of a structure that has been modeled with truncated modal dynamics. Stability problems arise due to the control design process beginning with a deficient model. This work provides the necessary conditions for the optimal estimation of infinite dimensional (partial-differential equation) systems. This approach can be particularly useful for initial control studies and for gaining considerable insight into what the optimal estimator for a truly infinite dimensional system should be. A detailed example of the estimation of the continuous shape of a string in tension is presented.

Schaechter, D. B.↗
NASA NTRS1987

Leading-edge vortex solutions with large total pressure losses

Computations are presented for a Lambda = 75 deg delta wing in a supersonic freestream under two conditions which lead to leading-edge vortices. For one condition, analysis of the computed vortical flow reveals a closed streamline in the core. From varying computational parameters, it appears that this is due to truncation error of the convective derivatives. For the other condition, comparisons are made with wind-tunnel data, and good agreement is noted for pitot pressure distributions, flow angles on the symmetry plane, and the position of an embedded shock. Many of the aerodynamic parameters are shown to be insensitive to grid spacing.

Murman, Earll M.↗
NASA NTRS1993

The distribution of flare parameters and implications for coronal heating

Assuming that solar hard X-ray flares are caused by bremsstrahlung from a distribution of accelerated electrons impinging on a thick target, we determine the relationship between the total flare energy released in the corona and observable flare parameters. It is shown that for the more numerous smaller flares to be energetically more important than larger ones, there must be an anticorrelation between observed flux and one or more of the photon spectral index, the total flare duration, the ratio of the total energy that goes into heating the corona to the energy that goes into acceleration of electrons, and the reciprocal of the low-energy cutoff of the accelerated electron spectrum. Using a new statistical test which properly takes account of data truncation due to observational selection effects, we find that neither spectral index nor duration show an anticorrelation with flux in ISEE 3/ICE hard X-ray data. Using a nonparametric method, we obtain single distributions of peak flux, spectral index, duration, and fluence, corrected for selection biases. The functional forms of these distributions do not support coronal heating by numerous small flares without very specific behavior by yet unobserved parameters.

Lee, Theodore T.↗
DOE OSTI2021

Quantifying and attributing time step sensitivities in present-day climate simulations conducted with EAMv1

Abstract. This study assesses the relative importance of time integration error in present-day climate simulations conducted with the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) at 1∘ horizontal resolution. We show that a factor-of-6 reduction of time step size in all major parts of the model leads to significant changes in the long-term mean climate. Examples of changes in 10-year mean zonal averages include the following: up to 0.5 K of warming in the lower troposphere and cooling in the tropical and subtropical upper troposphere, 1 %–10 % decreases in relative humidity throughout the troposphere, and 10 %–20 % decreases in cloud fraction in the upper troposphere and decreases exceeding 20 % in the subtropical lower troposphere. In terms of the 10-year mean geographical distribution, systematic decreases of 20 %–50 % are seen in total cloud cover and cloud radiative effects in the subtropics. These changes imply that the reduction of temporal truncation errors leads to a notable although unsurprising degradation of agreement between the simulated and observed present-day climate; to regain optimal climate fidelity in the absence of those truncation errors, the model would require retuning. A coarse-grained attribution of the time step sensitivities is carried out by shortening time steps used in various components of EAM or by revising the numerical coupling between some processes. Our analysis leads to the finding that the marked decreases in the subtropical low-cloud fraction and total cloud radiative effect are caused not by the step size used for the collectively subcycled turbulence, shallow convection, and stratiform cloud macrophysics and microphysics parameterizations but rather by the step sizes used outside those subcycles. Further analysis suggests that the coupling frequency between the subcycles and the rest of EAM significantly affects the subtropical marine stratocumulus decks, while deep convection has significant impacts on trade cumulus. The step size of the cloud macrophysics and microphysics subcycle itself appears to have a primary impact on cloud fraction in the upper troposphere and also in the midlatitude near-surface layers. Impacts of step sizes used by the dynamical core and the radiation parameterization appear to be relatively small. These results provide useful clues for future studies aiming at understanding and addressing the root causes of sensitivities to time step sizes and process coupling frequencies in EAM. While this study focuses on EAMv1 and the conclusions are likely model-specific, the presented experimentation strategy has general value for weather and climate model development, as the methodology can help researchers identify and understand sources of time integration error in sophisticated multi-component models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A fast approximation to supershell partition functions

A formula for supershell partition functions, which play a major role in the Super Transition Array approach to radiative-opacity calculations, is derived as a functional of the distribution of energies within the supershell. It consists in an alternative expansion for an arbitrary number of electrons or holes which also allows for quick approximate evaluations with truncated number of terms in the expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗