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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Development of a fractional-step method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A fractional step method is developed for solving the time-dependent three-dimensional incompressible Navier-Stokes equations in generalized coordinate systems. The primitive variable formulation uses the pressure, defined at the center of the computational cell, and the volume fluxes across the faces of the cells as the dependent variables, instead of the Cartesian components of the velocity. This choice is equivalent to using the contravariant velocity components in a staggered grid multiplied by the volume of the computational cell. The governing equations are discretized by finite volumes using a staggered mesh system. The solution of the continuity equation is decoupled from the momentum equations by a fractional step method which enforces mass conservation by solving a Poisson equation. This procedure, combined with the consistent approximations of the geometric quantities, is done to satisfy the discretized mass conservation equation to machine accuracy, as well as to gain the favorable convergence properties of the Poisson solver. The momentum equations are solved by an approximate factorization method, and a novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two- and three-dimensional laminar test cases are computed and compared with other numerical and experimental results to validate the solution method. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗

DOE OSTI2024

Evaluating county-level lung cancer incidence from environmental radiation exposure, PM 2.5 , and other exposures with regression and machine learning models

Characterizing the interplay between exposures shaping the human exposome is vital for uncovering the etiology of complex diseases. For example, cancer risk is modified by a range of multifactorial external environmental exposures. Environmental, socioeconomic, and lifestyle factors all shape lung cancer risk. However, epidemiological studies of radon aimed at identifying populations at high risk for lung cancer often fail to consider multiple exposures simultaneously. For example, moderating factors, such as PM 2.5 , may affect the transport of radon progeny to lung tissue. This ecological analysis leveraged a population-level dataset from the National Cancer Institute’s Surveillance, Epidemiology, and End-Results data (2013–17) to simultaneously investigate the effect of multiple sources of low-dose radiation (gross γ activity and indoor radon) and PM 2.5 on lung cancer incidence rates in the USA. County-level factors (environmental, sociodemographic, lifestyle) were controlled for, and Poisson regression and random forest models were used to assess the association between radon exposure and lung and bronchus cancer incidence rates. Tree-based machine learning (ML) method perform better than traditional regression: Poisson regression: 6.29/7.13 (mean absolute percentage error, MAPE), 12.70/12.77 (root mean square error, RMSE); Poisson random forest regression: 1.22/1.16 (MAPE), 8.01/8.15 (RMSE). The effect of PM 2.5 increased with the concentration of environmental radon, thereby confirming findings from previous studies that investigated the possible synergistic effect of radon and PM 2.5 on health outcomes. In summary, the results demonstrated (1) a need to consider multiple environmental exposures when assessing radon exposure’s association with lung cancer risk, thereby highlighting (1) the importance of an exposomics framework and (2) that employing ML models may capture the complex interplay between environmental exposures and health, as in the case of indoor radon exposure and lung cancer incidence.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

DOE OSTI2024

Enhancing scalability and accuracy of quantum poisson solver

The Poisson equation has many applications across the broad areas of science and engineering. Most quantum algorithms for the Poisson solver presented so far either suffer from lack of accuracy and/or are limited to very small sizes of the problem and thus have no practical usage. In this regard, our previous work showed a proof-of-concept demonstration in advancing quantum Poisson solver algorithm and validated preliminary results for a simple case of 3 x 3 problem. In this work, we delve into comprehensive research details, presenting the results on up to 15 x 15 problems that include step-by-step improvements in Poisson equation solutions, scaling performance, and experimental exploration. In particular, we demonstrate the implementation of eigenvalue amplification by a factor of up to 2 8 , achieving a significant improvement in the accuracy of our quantum Poisson solver and comparing that to the exact solution. Additionally, we present success probability results, highlighting the reliability of our quantum Poisson solver. Moreover, we explore the scaling performance of our algorithm against the circuit depth and width, demonstrating how our approach scales with larger problem sizes and thus further solidifies the practicality of easy adaptation of this algorithm in real-world applications. We also discuss a multilevel strategy for how this algorithm might be further improved to explore much larger problems with greater performance. Finally, through our experiments on the IBM quantum hardware, we conclude that though overall results on the existing NISQ hardware are dominated by the error in the CNOT gates, this work opens a path to realizing a multidimensional Poisson solver on near-term quantum hardware.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1992

Elevated temperature axial and torsional fatigue behavior of Haynes 188

The results of high-temperature axial and torsional low-cycle fatigue experiments performed on Haynes 188, a wrought cobalt-base superalloy, are reported. Fatigue tests were performed at 760 C in air on thin-walled tubular specimens at various ranges under strain control. Data are also presented for coefficient of thermal expansion, elastic modulus, and shear modulus at various temperatures from room to 1000 C, and monotonic and cyclic stress-strain curves in tension and in shear at 760 C. The data set is used to evaluate several multiaxial fatigue life models (most were originally developed for room temperature multiaxial life prediction) including von Mises equivalent strain range (ASME boiler and pressure vessel code), Manson-Halford, Modified Multiaxiality Factor (proposed here), Modified Smith-Watson-Topper, and Fatemi-Socie-Kurath. At von Mises equivalent strain ranges (the torsional strain range divided by the square root of 3, taking the Poisson's ratio to be 0.5), torsionally strained specimens lasted, on average, factors of 2 to 3 times longer than axially strained specimens. The Modified Multiaxiality Factor approach shows promise as a useful method of estimating torsional fatigue life from axial fatigue data at high temperatures. Several difficulties arose with the specimen geometry and extensometry used in these experiments. Cracking at extensometer probe indentations was a problem at smaller strain ranges. Also, as the largest axial and torsional strain range fatigue tests neared completion, a small amount of specimen buckling was observed.

Bonacuse, Peter J.↗

NASA NTRS1995

Elevated temperature axial and torsional fatigue behavior of Haynes 188

The results are reported for high-temperature axial and torsional low-cycle fatigue experiments performed at 760 C in air on thin-walled tubular specimens of Haynes 188, a wrought cobalt-based superalloy. Data are also presented for mean coefficient of thermal expansion, elastic modulus, and shear modulus at various temperatures from room to 1000 C, and monotonic and cyclic stress-strain curves in tension and in shear at 760 C. This data set is used to evaluate several multiaxial fatigue life models (most were originally developed for room temperature multiaxial life prediction) including von Mises equivalent strain range (ASME Boiler and Pressure Code), Manson-Halford, modified multiaxiality factor (proposed in this paper), modified Smith-Watson-Topper, and Fatemi-Socie-Kurath. At von Mises equivalent strain ranges (the torsional strain range divided by the square root of 3, taking the Poisson's ratio to be 0.5), torsionally strained specimens lasted, on average, factors of 2 to 3 times longer than axially strained specimens. The modified multiaxiality factor approach shows promise as a useful method of estimating torsional fatigue life from axial fatigue data at high temperatures. Several difficulties arose with the specimen geometry and extensometry used in these experiments. Cracking at extensometer probe indentations was a problem at smaller strain ranges. Also, as the largest axial and torsional strain range fatigue tests neared completion, a small amount of specimen buckling was observed.

Bonacuse, Peter J.↗

DOE OSTI2021

Power spectrum and form factor in random diagonal matrices and integrable billiards

Triggered by a controversy surrounding a universal behavior of the power spectrum in quantum systems exhibiting regular classical dynamics, we focus on a model of random diagonal matrices (RDM), often associated with the Poisson spectral universality class, and examine how the power spectrum and the form factor get affected by two-sided truncations of RDM spectra. Having developed a nonperturbative description of both statistics, we perform their detailed asymptotic analysis to demonstrate explicitly how a traditional assumption (lying at the heart of the controversy) – that the power spectrum is merely determined by the spectral form factor – breaks down for truncated spectra. This observation has important consequences as we further argue that bounded quantum systems with integrable classical dynamics are described by heavily truncated rather than complete RDM spectra. High-precision numerical simulations of semicircular and irrational rectangular billiards lend independent support to these conclusions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1987

Finite-element analysis of corner cracks in rectangular bars

Stress intensity factors for a wide range of quarter-elliptical corner cracks in rectangular bars are presented. Cracked configurations were subjected to remote tension, in-plane bending, or out-of-plane bending. Results were compared with other numerical results for a quarter-circular corner crack configuration. The present results along the interior of the crack generally agreed within 3 percent with those from the literature. Some larger differences (3 to 13 percent) were observed near the intersection of the crack front and the free surfaces in the boundary layer region. Analyses were performed to study the effect of mesh refinement in the boundary layer region and the influence of Poisson's ratio on the distribution of stress-intensity factors.

Raju, L. S.↗

NASA NTRS1992

Algorithm Calculates Cumulative Poisson Distribution

Algorithm calculates accurate values of cumulative Poisson distribution under conditions where other algorithms fail because numbers are so small (underflow) or so large (overflow) that computer cannot process them. Factors inserted temporarily to prevent underflow and overflow. Implemented in CUMPOIS computer program described in "Cumulative Poisson Distribution Program" (NPO-17714).

Bowerman, Paul N.↗

NASA NTRS1999

Direct Detection Doppler Lidar for Spaceborne Wind Measurement

Aerosol and molecular based versions of the double-edge technique can be used for direct detection Doppler lidar spaceborne wind measurement. The edge technique utilizes the edge of a high spectral resolution filter for high accuracy wind measurement using direct detection lidar. The signal is split between an edge filter channel and a broadband energy monitor channel. The energy monitor channel is used for signal normalization. The edge measurement is made as a differential frequency measurement between the outgoing laser signal and the atmospheric backscattered return for each pulse. As a result the measurement is insensitive to laser and edge filter frequency jitter and drift at a level less than a few parts in 10(exp 10). We have developed double edge versions of the edge technique for aerosol and molecular-based lidar measurement of the wind. Aerosol-based wind measurements have been made at Goddard Space Flight Center and molecular-based wind measurements at the University of Geneva. We have demonstrated atmospheric measurements using these techniques for altitudes from 1 to more than 10 km. Measurement accuracies of better than 1.25 m/s have been obtained with integration times from 5 to 30 seconds. The measurements can be scaled to space and agree, within a factor of two, with satellite-based simulations of performance based on Poisson statistics. The theory of the double edge aerosol technique is described by a generalized formulation which substantially extends the capabilities of the edge technique. It uses two edges with opposite slopes located about the laser frequency at approximately the half-width of each edge filter. This doubles the signal change for a given Doppler shift and yields a factor of 1.6 improvement in the measurement accuracy compared to the single edge technique. The use of two high resolution edge filters substantially reduces the effects of Rayleigh scattering on the measurement, as much as order of magnitude, and allows the signal to noise ratio to be substantially improved in areas of low aerosol backscatter. We describe a method that allows the Rayleigh and aerosol components of the signal to be independently determined using the two edge channels and an energy monitor channel. The effects of Rayleigh scattering may then subtracted from the measurement and we show that the correction process does not significantly increase the measurement noise for Rayleigh to aerosol ratios up to 10. We show that for small Doppler shifts a measurement accuracy of 0.4 m/s can be obtained for 5000 detected photon, 1.2 m/s for 1000 detected photons, and 3.7 m/s for 50 detected photons for a Rayleigh to aerosol ratio of 5. Methods for increasing the dynamic range of the aerosol-based system to more than +/- 100 m/s are given.

Korb, C. Laurence↗

NASA NTRS2012

Role of Temperature, Humidity and Rainfall on Influenza Transmission in Guatemala, El Salvador and Panama

Worldwide, seasonal influenza causes about 500,000 deaths and 5 million severe illnesses per year. The environmental drivers of influenza transmission are poorly understood especially in the tropics. We aimed to identify meteorological factors for influenza transmission in tropical Central America. We gathered laboratory-confirmed influenza case-counts by week from Guatemala City, San Salvador Department (El Salvador) and Panama Province from 2006 to 2010. The average total cases per year were: 390 (Guatemala), 99 (San Salvador) and 129 (Panama). Meteorological factors including daily air temperature, rainfall, relative and absolute humidity (RH, AH) were obtained from ground stations, NASA satellites and land models. For these factors, we computed weekly averages and their deviation from the 5-yr means. We assessed the relationship between the number of influenza case-counts and the meteorological factors, including effects lagged by 1 to 4 weeks, using Poisson regression for each site. Our results showed influenza in San Salvador would increase by 1 case within a week of every 1 day with RH>75% (Relative Risk (RR)= 1.32, p=.001) and every 1C increase in minimum temperature (RR=1.29, p=.007) but it would decrease by 1 case for every 1mm-above mean weekly rainfall (RR=0.93,p<.001) (model pseudo-R2=0.55). Within 2 weeks, influenza in Panama was increased by 1 case for every 1% increase in RH (RR=1.04, p=.003), and it was increased by 2 cases for every 1C increase of minimum temperature (RR=2.01, p<.001) (model pseudo-R2=0.4). Influenza counts in Guatemala had 1 case increase for every 1C increase in minimum temperature in the previous week (RR=1.21, p<.001), and for every 1mm/day-above normal increase of rainfall rate (RR=1.03, p=.03) (model pseudo-R2=0.54). Our findings that cases increase with temperature and humidity differ from some temperate-zone studies. But they indicate that climate parameters such as humidity and temperature could be predictive of influenza activity and should be incorporated into country-specific influenza transmission models

Soebiyanto, Radina P.↗

NASA NTRS1989

Fast, adaptive summation of point forces in the two-dimensional Poisson equation

A comparatively simple procedure is presented for the direct summation of the velocity field introduced by point vortices which significantly reduces the required number of operations by replacing selected partial sums by asymptotic series. Tables are presented which demonstrate the speed of this algorithm in terms of the mere doubling of computational time in dealing with a doubling of the number of vortices; current methods involve a computational time extension by a factor of 4. This procedure need not be restricted to the solution of the Poisson equation, and may be applied to other problems involving groups of points in which the interaction between elements of different groups can be simplified when the distance between groups is sufficiently great.

Van Dommelen, Leon↗

DOE OSTI2024

Convolutional Non-Homogeneous Poisson Process and its Application to Wildfire Ignition Risk Quantification for Power Delivery Networks

To quantify wildfire ignition risks on power delivery networks, the current practice predominantly relies on the empirically calculated fire danger indices, which may not well capture the effects of dynamically changing environmental factors. This article proposes a spatio-temporal point process model, known as the Convolutional Non-homogeneous Poisson Process (cNHPP), and applies the model to quantify wildfire ignition risks for power delivery networks. The proposed model captures both the current (i.e., instantaneous) and cumulative (i.e., historical) effects of key environmental processes (i.e., covariates) on wildfire risks, as well as the spatio-temporal dependency among different segments of the power delivery network. The computation and interpretation of the intensity function are thoroughly investigated. We apply the proposed approach to estimate wildfire ignition risks on major transmission lines in California, using historical fire data, meteorological and vegetation data obtained from the National Oceanic and Atmospheric Administration and National Aeronautics and Space Administration. Here, a comprehensive comparison study is performed to show the applicability and predictive capability of the proposed approach.

Non-homogeneous Poisson Process↗

NASA NTRS1986

Fourier analysis of the SOR iteration

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. It is shown that for the standard model problem of Poisson's equation on a rectangle, the optimal omega and corresponding convergence rate can be rigorously obtained by Fourier analysis. The trick is to tilt the space-time grid so that the SOR stencil becomes symmetrical. The tilted grid also gives insight into the relation between convergence rates of several variants.

Leveque, R. J.↗

NASA NTRS1998

The Double Edge Aerosol and Molecular Techniques for Doppler Lidar Wind Measurement

We have developed the theory for aerosol- and molecular-based lidar measurements of the wind using double edge versions of the edge technique. Aerosol-based wind measurements have been made at Goddard Space Flight Center and molecular-based wind measurements at the University of Geneva. We have demonstrated atmospheric measurements using these techniques for altitudes from 1 to more than 10 km. Measurement accuracies of better than 1.25 m/s have been obtained with integration times from 5 to 30 seconds. The measurements can be scaled to space and agree, within a factor of two, with satellite-based simulations of performance based on Poisson statistics.

Korb, C. Laurence↗

NASA NTRS1988

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe↗

DOE OSTI2020

Spectral fluctuations in the Sachdev-Ye-Kitaev model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N ) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q= 4 or higher even q-body interactions, but also for the q= 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Direct Observations of Ion Densities at Functionalized Interfaces to Test Hypotheses Regarding the Origin of Specific Ion Effects

Interactions of anions with protonatable groups were investigated using X-ray fluorescence near total reflection (XFNTR) on floating monolayers at the surface of water. The number of ions attracted to the interfacial region, which XFNTR measures directly, is ion-specific as well as monolayer-specific. Our observation of the distinctly different behaviors of ClO 4 – and ReO 4 – , two ions with the same tetrahedral structure and almost the same sizes and hydration enthalpies, challenges current theories of ion specificity. Our observations are inconsistent with not only the Gouy–Chapman model (as expected) but also size-modified Poisson–Boltzmann theory and the “law of matching water affinity”. Furthermore, we suggest that factors other than ion size and ion–water interactions, including possibly ion–ion interactions and lateral ordering at the interface, must be considered to account for specific ion effects.

36 MATERIALS SCIENCE↗

NASA NTRS1982

The crack problem for a nonhomogeneous plane

The plane elasticity problem for a nonhomogeneous medium containing a crack is considered. It is assumed that the Poisson's ratio of the medium is constant and the Young's modulus E varies exponentially with the coordinate parallel to the crack. First the half plane problem is formulated and the solution is given for arbitrary tractions along the boundary. Then the integral equation for the crack problem is derived. It is shown that the integral equation having the derivative of the crack surface displacement as the density function has a simple Cauchy type kernel. Hence, its solution and the stresses around the crack tips have the conventional square root singularity. The solution is given for various loading conditions. The results show that the effect of the Poisson's ratio and consequently that of the thickness constraint on the stress intensity factors are rather negligible.

Delale, F.↗