Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Poisson factorization”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

IPACS (Integrated Probabilistic Assessment of Composite Structures): Code development and applications

A methodology and attendant computer code have been developed and are described to computationally simulate the uncertain behavior of composite structures. The uncertain behavior includes buckling loads, stress concentration factors, displacements, stress/strain etc., which are the consequences of the inherent uncertainties (scatter) in the primitive (independent random) variables (constituent, ply, laminate and structural) that describe the composite structures. The computer code, IPACS (Integrated Probabilistic Assessment of Composite Structures), can handle both composite mechanics and composite structures. Application to probabilistic composite mechanics is illustrated by its uses to evaluate the uncertainties in the major Poisson's ratio and in laminate stiffness and strength. IPACS application to probabilistic structural analysis is illustrated by its use to evaluate the uncertainties in the buckling of a composite plate, in the stress concentration factor in a composite panel and in the vertical displacement and ply stress in a composite aircraft wing segment.

Chamis, C. C.↗
DOE OSTI2025

Parameterized anomalous transport model for current-carrying collisionless plasmas in pulsed power inertial confinement fusion

Current delivery in pulsed power inertial confinement fusion is influenced by collisionless current-carrying microturbulent plasmas, which are sourced from electrode surfaces. In this setting, the lower hybrid drift instability—triggered by plasma acceleration—is a leading candidate driver of difficult-to-predict momentum and energy transport. To characterize the nonlinear state of the microturbulent plasma, a parameterized anomalous transport model is developed for the instability, with analytic formulas for anomalous collision frequency, resistivity, and species heating rates. The formulas are expressed in terms of linear-theory variables and four dimensionless parameters that characterize the macroscopic plasma state. The model is built on linear theory analysis, power law analysis, and quasilinear theory analysis, and is validated using a series of nonlinear continuum kinetic Vlasov–Poisson simulations. The theoretical and computational investigation demonstrates that the anomalous collision frequency associated with the instability can be reliably approximated, within about a factor of five or better, by the unscaled linear theory growth rate of the fastest-growing wavenumber mode. This finding enables efficient calculation of anomalous resistivity and species heating rates over a wide range of plasma conditions, resulting in improved predictive capabilities.

Complex functions↗
NASA NTRS1990

Calculation of the binomial survivor function

A method is presented for calculating the binomial SF (cumulative binomial distribution), binfc(k;p,n), especially for a large n, beyond the range of existing tables, where conventional computer programs fail because of underflow and overflow, and Gaussian or Poisson approximations yield insufficient accuracy for the purpose at hand. This method is used to calculate and sum the individual binomial terms while using multiplication factors to avoid underflow; the factors are then divided out of the partial sum whenever it has the potential to overflow. A computer program uses this technique to calculate the binomial SF for arbitrary inputs of k, p, and n. Two other algorithms are presented to determine the value of p needed to yield a specified SF for given values of k and n and calculate the value where p = SF for a given k and n. Reliability applications of each algorithm/program are given, e.g., the value of p needed to achieve a stated k-out-of-n:G system reliability and the value of p for which k-out-of-n:G system reliability equals p.

Bowerman, Paul N.↗
NASA NTRS1945

The Elastic Constants for Wrought Aluminum Alloys

There are several constants which have been devised as numerical representations of the behavior of metals under the action of loadings which stress the metal within the range of elastic action. Some of these constants, such as Young's modulus of elasticity in tension and compression, shearing modulus of elasticity, and Poisson's ratio, are regularly used in engineering calculations. Precise tests and experience indicate that these elastic constants are practically unaffected by many of the factors which influence the other mechanical properties of materials and that a few careful determinations under properly controlled conditions are more useful and reliable than many determinations made under less favorable conditions. It is the purpose of this paper to outline the methods employed by the Aluminum Research Laboratories for the determination of some of these elastic constants, to list the values that have been determined for some of the wrought aluminum alloys, and to indicate the variations in the values that may be expected for some of the commercial products of these alloys.

ALUMINUM ALLOYS, WROUGHT-ELASTICITY↗
NASA NTRS2006

Application of p-Multigrid to Discontinuous Galerkin Formulations of the Poisson Equation

We investigate p-multigrid as a solution method for several different discontinuous Galerkin (DG) formulations of the Poisson equation. Different combinations of relaxation schemes and basis sets have been combined with the DG formulations to find the best performing combination. The damping factors of the schemes have been determined using Fourier analysis for both one and two-dimensional problems. One important finding is that when using DG formulations, the standard approach of forming the coarse p matrices separately for each level of multigrid is often unstable. To ensure stability the coarse p matrices must be constructed from the fine grid matrices using algebraic multigrid techniques. Of the relaxation schemes, we find that the combination of Jacobi relaxation with the spectral element basis is fairly effective. The results using this combination are p sensitive in both one and two dimensions, but reasonable convergence rates can still be achieved for moderate values of p and isotropic meshes. A competitive alternative is a block Gauss-Seidel relaxation. This actually out performs a more expensive line relaxation when the mesh is isotropic. When the mesh becomes highly anisotropic, the implicit line method and the Gauss-Seidel implicit line method are the only effective schemes. Adding the Gauss-Seidel terms to the implicit line method gives a significant improvement over the line relaxation method.

Helenbrook, B. T.↗
NASA NTRS2007

Managing Uncertainty Due to a Fundamental Error Source Arising from Scatterer Distribution Complexity in Radar Remote Sensing of Precipitation

The assumption that cloud and rain drops are spatially distributed according to a Poisson distribution within a scattering volume probed by a radar being used to estimate precipitation has represented bedrock theory in establishing 'rules of the game' for pulse averaging--the process needed to beat down noise to an acceptable level in the measurement of radar reflectivity factor. Based on relatively recent observations of 'realistic' spatial distributions of hydrometeor scatterers in a cloudy atmosphere motivates a renewed examination of the consequences of using a too simplified assumption underlying volume scattering--particularly in regards to the standard pulse averaging rule. Our investigation addresses two extremes, simple to complex, insofar as allowed for complexities in an underlying scatterer distribution. It is demonstrated that as the spatial distribution ranges from Poisson (a narrow distribution) to multi-fractal (much broader distribution), uncertainty in a measurement increases if the rule for pulse averaging goes unchanged from its Poisson distribution reference county. [A bounded cascade is used for the multi-fractal distribution, a regularly observed distribution vis-a-vis cloud liquid water content.] The resultant measurement uncertainty leads to a fundamental source of error in the estimation of rain rate from radar measurements, one that has been disregarded since the early 1950s when radar sets first began to be used for rainfall measuring. It is shown how this source of error can be 'managed'--under the assumption that number of data analysis experiments would be carried out, experiments involving pulse-by-pulse measurements obtained from a radar set modified to output individual pulses of reflectivity factor. For practical applications, a new parameter called normalized k-sample intensity invariance is developed to enable defining the required pulse average count according to a preferred degree of uncertainty.

Smith, Eric A.↗
DOE OSTI2022

Estimating the local dark matter density in a non-axisymmetric wobbling disc

The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Quantifying radiation quality for space relevant radiation types: Fitting excess risk models to outbred mice data

Accurately quantifying the differences in radiation quality between space and terrestrial environments is important for predicting health risks for astronauts.Recently, Edmundson et al. 2020[1] provided valuable new results from out-bred mice linking tumor induction and genetic background aer exposure to low and high-LET radiation. The goal of the current study is to more rigorously estimate a relative biological effectiveness (RBE) factor by leveraging the solid tumor data from Edmundson et al. 2020. Excess relative risk (ERR) models and excess absolute risk (EAR) models were fit using Poisson regression similar to the models that the Radiation Effects Research Foundation uses to fit atomic bomb survivor data. Linear ERR and EAR slopes were simulated usingBayesian analyses, and RBE values were calculated from the ratio of the heavy ion linear slope to the gamma linear slope using the full posterior distribution.

Lori J. Chappell↗
DOE OSTI2024

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Characterization of Ice for Return-to-Flight of the Space Shuttle: Soft Ice - Part 2

In support of characterizing ice debris for return-to-flight (RTF) of NASA's space shuttle, we have determined the microstructure, density and compressive strength (at -10 C at approximately 0.3 per second) of porous or soft ice that was produced from both atmospheric water and consolidated snow. The study showed that the atmospheric material was generally composed of a mixture of very fine (0.1 to 0.3 millimeters) and coarser (5 to 10 millimeter) grains, plus air bubbles distributed preferentially within the more finely-grained part of the microstructure. The snow ice was composed of even finer grains (approximately 0.05 millimeters) and contained more pores. Correspondingly, the snow ice was of lower density than the atmospheric ice and both materials were significantly less dense than hard ice. The atmospheric ice was stronger (approximately 3.8 MPa) than the snow ice (approximately 1.9 MPa), but weaker by a factor of 2 to 5 than pore-free hard ice deformed under the same conditions. Zero Values are given for Young's modulus, compressive strength and Poisson's ratio that can be used for modeling soft ice from the external tank (ET).

Schulson, Erland M.↗
DOE OSTI2025

Fuzzy spheres in stringy matrix models: quantifying chaos in a mixed phase space

We consider a truncation of the BMN matrix model to a configuration of two fuzzy spheres, described by two coupled non-linear oscillators dependent on the mass parameter μ. The classical phase diagram of the system generically (μ ≠ 0) contains three equilibrium points: two centers and a center-saddle; as μ → 0 the system exhibits a pitchfork bifurcation. We demonstrate that the system is exactly integrable in quadratures for μ = 0, while for very large values of μ, it approaches another integrable point characterized by two harmonic oscillators. The classical phase space is mixed, containing both integrable islands and chaotic regions, as evidenced by the classical Lyapunov spectrum. At the quantum level, we explore indicators of early and late time chaos. The eigenvalue spacing is best described by a Brody distribution, which interpolates between Poisson and Wigner distributions; it dovetails, at the quantum level, the classical results and reemphasizes the notion that the quantum system is mixed. We also study the spectral form factor and the quantum Lyapunov exponent, as defined by out-of-time-ordered correlators. These two indicators of quantum chaos exhibit weak correlations with the Brody distribution. We speculate that the behavior of the system as μ → 0 dominates the spectral form factor and the quantum Lyapunov exponent, making these indicators of quantum chaos less effective in the context of a mixed phase space.

AdS-CFT correspondence↗
NASA NTRS1998

Correlations Between the Cosmic X-Ray and Microwave Backgrounds: Constraints on a Cosmological Constant

In universes with significant curvature or cosmological constant, cosmic microwave background (CMB) anisotropies are created very recently via the Rees-Sciama or integrated Sachs-Wolfe effects. This causes the CMB anisotropies to become partially correlated with the local matter density (z less than 4). We examine the prospects of using the hard (2- 10 keV) X-ray background as a probe of the local density and the measured correlation between the HEAO1 A2 X-ray survey and the 4-year COBE-DMR map to obtain a constraint on the cosmological constant. The 95% confidence level upper limit on the cosmological constant is OMega(sub Lambda) less than or equal to 0.5, assuming that the observed fluctuations in the X-ray map result entirely from large scale structure. (This would also imply that the X-rays trace matter with a bias factor of b(sub x) approx. = 5.6 Omega(sub m, sup 0.53)). This bound is weakened considerably if a large portion of the X-ray fluctuations arise from Poisson noise from unresolved sources. For example, if one assumes that the X-ray bias is b(sub x) = 2, then the 95% confidence level upper limit is weaker, Omega(sub Lambda) less than or equal to 0.7. More stringent limits should be attainable with data from the next generation of CMB and X-ray background maps.

Boughn, S. P.↗
DOE OSTI2025

Near-Efficient and Non-Asymptotic Multiway Inference

We establish non-asymptotic efficiency guarantees for tensor decomposition–based inference in count data models. Under a Poisson framework, we consider two related goals: (i) parametric inference , the estimation of the full distributional parameter tensor, and (ii) multiway analysis , the recovery of its canonical polyadic (CP) decomposition factors. Our main result shows that in the rank-one setting, a rank-constrained maximum-likelihood estimator achieves multiway analysis with variance matching the Cramér–Rao Lower Bound (CRLB) up to absolute constants and logarithmic factors. This provides a general framework for studying “near-efficient” multiway estimators in finite-sample settings. For higher ranks, we illustrate that our multiway estimator may not attain the CRLB; nevertheless, CP-based parametric inference remains nearly minimax optimal, with error bounds that improve on prior work by offering more favorable dependence on the CP rank. Numerical experiments corroborate near-efficiency in the rank-one case and highlight the efficiency gap in higher-rank scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Chandra measurements of gas homogeneity and turbulence at intermediate radii in the Perseus Cluster

ABSTRACT We present a Chandra study of surface brightness fluctuations in the diffuse intracluster medium of the Perseus Cluster. Our study utilizes deep, archival imaging of the cluster core as well as a new mosaic of 29 short $\, 5\,$-ks observations extending in eight different directions out to radii of $r_{500} \sim 2.2\, r_{2500}$. Under the assumption that the distribution of densities at a given radius is log-normally distributed, two important quantities can be derived from the width of the log-normal density distribution on a given spatial scale: the density bias, that is equal to the square root of the clumping factor C; and the one-component turbulent velocity, vk,1D. We forward-model all contributions to the measured surface brightness, including astrophysical and particle background components, and account for the Poisson nature of the measured signal. Measuring the distribution of surface brightness fluctuations in 1 arcmin2 regions, spanning the radial range $0.3\rm {-}2.2\, r_{2500}$ ($7.8\rm {-}57.3\,$ arcmin), we find a small to moderate average density bias of around $3{{\ \rm per\ cent}}$ at radii below $1.6\, r_{2500}$. We also infer an average turbulent velocity at these radii of V1D < 400 km s−1. Direct confirmation of our results on turbulent velocities inferred from surface brightness fluctuations should be possible using the X-ray calorimeter spectrometers to be flown aboard the XRISM and Athena observatories.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗
NASA NTRS1989

Operator induced multigrid algorithms using semirefinement

A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two and three dimensional model problems are presented, together with a two level analysis explaining these results.

Decker, Naomi↗
NASA NTRS1989

Operator induced multigrid algorithms using semirefinement

A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two- and three-dimensional model problems are presented, together with a two level analysis explaining these results.

Decker, Naomi Henderson↗
NASA NTRSDate not supplied

Evaluating the Feasibility of Data Simulation for Radiation Research: Estimating RBE Values for HZE Radiation

As the drive for deep space exploration intensifies, a comprehensive understanding of the health effects of radiation exposure becomes paramount to the future of human space flight. However, epidemiological data for radiation exposure, particularly to high-energy (HZE) ions, is limited, partially due to the financial and logistical costs of radiation studies. As an alternative, this study aims to assess the viability of data simulation strategies to accurately model potential study parameters prior to utilizing laboratory conditions. This study estimates a relative biological effectiveness (RBE) factor based on the solid tumor data for outbred mice provided by Edmundson et al. 2020. Excess relative risk (ERR) models for HZE-Fe ions and gamma radiation were estimated using Poisson regression with Weibull models to represent the background solid tumor hazard without radiation. RBE values were calculated from the ratio of the heavy ion linear slope to the gamma linear slope. The parameters from these models were then applied to simulate iterations of 300 datasets across HZE-Fe doses of 0.05, 0.2, 0.4, and 0.75 Gy; gamma radiation doses of 0.75, 2, and 3 Gy, as well as an unirradiated control group. Sample size per dose varied from 100-500 mice across simulations (800-4000 mice total per dataset). 1500 datasets were generated total (300 each for sample sizes 800, 1600, 3200, and 4000). For each dataset, ERR per radiation type and RBE for HZE-Fe were calculated. The RBE from the Edmondson data was calculated to be 5.55. The RBEs from the simulated data converged around this value across the increasing sample sizes. For n = 4000, the mean RBE across the 300 datasets was 5.71 (95% CI: 5.47-5.95). These results suggest that RBEs calculated from simulated data are credible. Based on this exploration, data simulation is a viable method of testing radiation studies. It provides a method of testing study assumptions and refining research questions prior to embarking on costly laboratory experiments.

Sneha Mehta↗