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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1983

An implicit, transonic, full-potential code for cascade flow on H-grid topology

A transonic, full-potential code is developed for computing the flow through two-dimensional cascades using an H-type grid topology that employs an implicit approximate-factorization scheme. The body-conforming H-grid is generated numerically by solving Poisson's equation. The flow-solution algorithm at the coordinate mapping singularity associated with this grid is investigated using two different types of finite-difference schemes. The grid-geometry effect on these schemes is also studied by noting free-stream capturing properties. It is found that by implementing a consistent spatial differencing scheme, the mapping singularities can be resolved numerically, and the grid-geometry-induced error minimized. The code is verified by computing model cascade flow problems.

Kwak, D.↗

NASA NTRS1983

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.↗

NASA NTRS1991

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

The time-dependent, three-dimensional incompressible Navier-Stokes equations are presently solved in generalized coordinate systems by means of a fractional-step method whose primitive variable formulation uses as dependent variables, in place of the Cartesian components of the velocity: (1) pressure (defined at the center of the computational cell), and (2) volume fluxes across the faces of the cells. The momentum equations are solved by means of an approximate factorization method. A novel 'ZEBRA' scheme incorporating four-color ordering efficiently solves the Poisson equation. Illustrative two- and three-dimensional laminar flow test cases are computed and evaluated relative to extant numerical and experimental results, and good agreement is obtained.

Rosenfeld, Moshe↗

NASA NTRS1999

Structural Evaluation of the RSRM Nozzle Replacement Adhesive

This paper describes the structural performance evaluation of a replacement adhesive for the Reusable Solid Rocket Motor (RSRM) nozzle utilizing finite element analysis. Due to material obsolescence and industrial safety issues, the two current structural adhesives, EA 913 and EA 946 are to be replaced with a new adhesive. TIGA 321. The structural evaluation in support of the adhesive replacement effort includes residual stress, transportation, and flight analyses. Factors of safety are calculated using the stress response from each analysis. The factors of safety are used as the limiting criteria to compare the replacement adhesive against the current adhesives. Included in this paper are the analytical approach, assumptions and modeling techniques as well as the results of the evaluation. An important factor to the evaluation is the similarity in constitutive material properties (elastic modulus and Poisson's ratio) between TIGA 321 and EA 913. This similarity leads to equivalent material response from the two adhesives. However, TIGA 321 surpasses EA 913's performance due to higher material capabilities. Conversely, the change in stress response from EA 946 to TIGA 321 is more apparent: this is primarily attributed to the difference in the modulii of the two adhesives, which differ by two orders of magnitude. The results of the bondline evaluation indicate that the replacement adhesive provides superior performance than the current adhesives with only minor exceptions. Furthermore, TIGA 321 causes only a minor chance in the response of the phenolic and metal components.

Batista-Rodriguez, A.↗

NASA NTRS1989

An improved correction algorithm for number density measurements made with the forward scattering spectrometer probe

A correction factor to the number density measured by the Forward Scattering Spectrometer Probe (FSSP) which compensates for dead time and coincidence errors was determined by calculating the probabilities of, and the average number of particles in, the six possible types of dead time and coincidence events. These probabilities and averages were calculated by means of a probabilistic model based on Poisson statistics. A Monte Carlo computer simulation of the FSSP operation was also carried out and the number density correction factor was compared with the Monte Carlo data. For an actual number density of 2000/cu cm, it was found that the measured number density was of the order of 300/cu cm.

Lock, James A.↗

NASA NTRS2000

Nonlinear Local Bending Response and Bulging Factors for Longitudinal and Circumferential Cracks in Pressurized Cylindrical Shells

Results of a geometrically nonlinear finite element parametric study to determine curvature correction factors or bulging factors that account for increased stresses due to curvature for longitudinal and circumferential cracks in unstiffened pressurized cylindrical shells are presented. Geometric parameters varied in the study include the shell radius, the shell wall thickness, and the crack length. The major results are presented in the form of contour plots of the bulging factor as a function of two nondimensional parameters: the shell curvature parameter, lambda, which is a function of the shell geometry, Poisson's ratio, and the crack length; and a loading parameter, eta, which is a function of the shell geometry, material properties, and the applied internal pressure. These plots identify the ranges of the shell curvature and loading parameters for which the effects of geometric nonlinearity are significant. Simple empirical expressions for the bulging factor are then derived from the numerical results and shown to predict accurately the nonlinear response of shells with longitudinal and circumferential cracks. The numerical results are also compared with analytical solutions based on linear shallow shell theory for thin shells, and with some other semi-empirical solutions from the literature, and limitations on the use of these other expressions are suggested.

Young, Richard D.↗

DOE OSTI2026

Higher-order space-charge stability in anisotropic beams: Vlasov-Poisson derivation, refined dispersion relations, and stability charts

The Hofmann stability chart is used to screen working points in space-charge-dominated linacs. We identify two errors in its published higher-order dispersion relations: missing $(1\mp2\hatη^2/α)$ factors in the third-order $S^4$ coupling residues, and a sign error in the stated isotropic reduction of the fourth-order relation. Both corrections follow from Hofmann's Vlasov-Poisson equations without fitted parameters. They reproduce coherent tune-shift coefficients in the author's later monograph that the printed forms miss by 24% and 127%. Mode-resolved figures from a published application agree with the corrected relations and reject the printed forms, indicating an inconsistency between the 1998 equations and the calculations underlying those tested figures. We quantify the effect on the non-oscillatory stability chart. Inside the adopted $S^2\le10$ comparison domain, printed and corrected forms disagree on 0.73-2.11% of cells, with no preferred direction. Among excluded cells, disagreement reaches 22%, and the printed relation over-predicts instability at every sampled anisotropy. This concentration may help explain why the errors persisted, although it does not establish their historical cause. For PIP-II, the corrected chart flags four of thirty-two evaluable periods, including one on a third-order odd branch missed by a second-order screen. This count covers non-oscillatory modes only and remains conditional on an unresolved factor-five disagreement between two codes on transverse emittance growth.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

DOE OSTI2020

Experimental investigation on the effects of fabric architectures on mechanical and damage behaviors of carbon/epoxy woven composites

The mechanical behaviors and damage evolutions of carbon/epoxy woven fabric composites with three different geometries, i.e., one plain weave and two twill weave patterns with different areal densities, are studied under tensile loading. The effects of weave patterns on mechanical properties are investigated by monotonic and cyclic tension tests. Remarkable variations in stress–strain curve, Poisson’s ratio, residual strain and strain map exist in the three composites. Crimp ratio is found to be a critical factor to govern the mechanical properties. With smaller crimp ratio, a quasi-linear stress–strain curve with higher elastic modulus and strength is observed. The stress–strain curves of composites with higher crimp ratio contain transition stages with significant tangent modulus degradation. Elastic modulus, strength and damage initiation are all correlated with the crimp ratio linearly regardless of the fabric pattern. Dramatic nonlinear evolution in Poisson’s ratio occurs in the composite with higher crimp ratio. Cyclic tension results indicate that the residual strain is a more appropriate damage indicator than the unloading elastic modulus. Microstructure examination shows that damage developments are essentially related to the fabric geometry, and result in various mechanical behaviors. This work provides important insights into the geometry-deformation mechanism-mechanical property relationship of the woven composites.

36 MATERIALS SCIENCE↗

NASA NTRS1991

Charge collected by diffusion from an ion track under mixed boundary conditions

Charge-carrier diffusion from an ion track in a silicon substrate at least a few hundred microns thick is analyzed. The substrate upper surface is treated as reflective except for a small section, intended to represent a reverse-biased junction, which is treated as a sink. Total charge collected by the sink is calculated by assuming transport to be governed by an ambipolar diffusion equation with temporally constant and spatially uniform carrier lifetime and diffusion coefficient. Present results apply to a normally incident track but could easily be generalized to arbitrary track direction. The collected charge is found to depend on track length and on the electrostatic capacitance, rather than the area, of the sink. Theoretical predictions are compared to the results of a numerical simulation called the Poisson and Continuity Equation Solver (PISCES) for three cases and are found to agree within a factor of two in the worst case.

Edmonds, Larry D.↗

NASA NTRS2004

Rayleigh-Sommerfield Diffraction vs Fresnel-Kirchhoff, Fourier Propagation and Poisson's Spot

The boundary conditions imposed on the diffraction problem in order to obtain the Fresnel-Kirchhoff (FK) solution are well-known to be mathematically inconsistent and to be violated by the solution when the observation point is close to the diffracting screen 1-3. These problems are absent in the Rayleigh-Sommerfeld (RS) solution. The difference between RS and FK is in the inclination factor and is usually immaterial because the inclination factor is approximated by unity. But when this approximation is not valid, FK can lead to unacceptable answers. Calculating the on-axis intensity of Poisson s spot provides a critical test, a test passed by RS and failed by FK. FK fails because (a) convergence of the integral depends on how it is evaluated and (b) when the convergence problem is xed, the predicted amplitude at points near the obscuring disk is not consistent with the assumed boundary conditions.

RS (RAYLEIGH SOMMERFELD)↗

DOE OSTI2024

Bounds on galaxy stochasticity from halo occupation distribution modeling

The joint probability distribution of matter overdensity and galaxy counts in cells is a powerful probe of cosmology, and the extent to which variance in galaxy counts at fixed matter density deviates from Poisson shot noise is not fully understood. The lack of informed bounds on this stochasticity is currently the limiting factor in constraining cosmology with the galaxy–matter probability distribution function (PDF). We investigate stochasticity in the conditional distribution of galaxy counts along lines of sight with fixed matter density, and we present a halo occupation distribution (HOD)-based approach for obtaining plausible ranges for stochasticity parameters. To probe the high-dimensional space of possible galaxy–matter connections, we derive a set of HODs that conserve the galaxies’ linear bias and number density to produce RED M A G I C-like galaxy catalogs within the A BACUS S UMMIT suite of N -body simulations. We study the impact of individual HOD parameters and cosmology on stochasticity and perform a Monte Carlo search in HOD parameter space subject to the constraints on bias and density. In mock catalogs generated by the selected HODs, shot noise in galaxy counts spans both sub-Poisson and super-Poisson values, ranging from 80% to 133% of Poisson variance for cells with mean matter density. Nearly all of the derived HODs show a positive relationship between local matter density and stochasticity. For galaxy catalogs with higher stochasticity, modeling galaxy bias to second order is required for an accurate description of the conditional PDF of galaxy counts at fixed matter density. The presence of galaxy assembly bias also substantially extends the range of stochasticity in the super-Poisson direction. This HOD-based approach leverages degrees of freedom in the galaxy–halo connection to obtain informed bounds on nuisance model parameters and can be adapted to study other parametrizations of shot noise in galaxy counts, in particular to motivate prior ranges on stochasticity for cosmological analyses.

Britt, Dylan (ORCID:000000019905601X)↗

DOE OSTI2024

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems—like fixation of cells coupled with fluorescence microscopy—provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. In this study, these snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

On the Feynman-alpha method for reflected fissile assemblies

The Feynman-alpha method is a neutron noise technique that is used to estimate the prompt neutron period of fissile assemblies. The method and quantity are of widespread interest including in applications such as nuclear criticality safety, safeguards and nonproliferation, and stockpile stewardship; the prompt neutron period may also be used to infer the k eff multiplication factor. The Feynman-alpha method is predicated on time-correlated neutron detections that deviate from a Poisson random variable due to multiplication. Traditionally, such measurements are diagnosed with one-region point kinetics, but two-region models are required when the fissile assembly is reflected. This paper presents a derivation of the two-region point kinetics Feynman-alpha equations based on a double integration of the Rossi-alpha equations, develops novel propagation of measurement uncertainty, and validates the theory. Validation is achieved with organic scintillator measurements of weapons-grade plutonium reflected by various amounts of copper to achieve k eff values of 0.83–0.94 and prompt periods of 5–75 ns. The results demonstrate that Feynman-alpha measurements should use the two-region model instead of the one-region model. The simplified one-region model deviates from the validated two-region models by as much as 10% in the estimate of the prompt neutron period, and the two-region model reduces to the one-region model for small amounts of reflector. The Feynman-alpha estimates of the prompt neutron period are compared to those of the Rossi-alpha approach. The comparative results demonstrate that the Feynman-alpha method is more precise than the Rossi-alpha method and more accurate for k eff < 0.92, whereas the Rossi-alpha method is generally more accurate for higher multiplications. Here, the uncertainty propagation developed in this work should be used for all Feynman-alpha measurements and will therein improve fitting accuracy and appropriate precision estimates.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

NASA NTRS1982

Approximate factorization with an elliptic pressure solver for incompressible flow

Two-dimensional curvilinear coordinates are used to solve the incompressible Navier-Stokes equations, in conjunction with approximate factorization for the solution of the momentum equation and the successive overrelaxation by lines method for the solution of a Poisson equation for the pressure. The combined algorithm, although not fully explicit, is marginally stable at Reynolds numbers lower than 10,000 and time increments of 0.01. Pressure distributions calculated for attack angles of zero and 6 deg are of the same shape as the experimental curves, but are shifted to one side.

Bernard, R. S.↗

NASA NTRS1947

Elastic Properties in Tension and Shear of High Strength Nonferrous Metals and Stainless Steel - Effect of Previous Deformation and Heat Treatment

A resume is given of an investigation of the influence of plastic deformation and of annealing temperature on the tensile and shear elastic properties of high strength nonferrous metals and stainless steels in the form of rods and tubes. The data were obtained from earlier technical reports and notes, and from unpublished work in this investigation. There are also included data obtained from published and unpublished work performed on an independent investigation. The rod materials, namely, nickel, monel, inconel, copper, 13:2 Cr-Ni steel, and 18:8 Cr-Ni steel, were tested in tension; 18:8 Cr-Ni steel tubes were tested in shear, and nickel, monel, aluminum-monel, and Inconel tubes were tested in both tension and shear. There are first described experiments on the relationship between hysteresis and creep, as obtained with repeated cyclic stressing of annealed stainless steel specimens over a constant load range. These tests, which preceded the measurements of elastic properties, assisted in devising the loading time schedule used in such measurements. From corrected stress-set curves are derived the five proof stresses used as indices of elastic or yield strength. From corrected stress-strain curves are derived the secant modulus and its variation with stress. The relationship between the forms of the stress-set and stress-strain curves and the values of the properties derived is discussed. Curves of variation of proof stress and modulus with prior extension, as obtained with single rod specimens, consist in wavelike basic curves with superposed oscillations due to differences of rest interval and extension spacing; the effects of these differences are studied. Oscillations of proof stress and modulus are generally opposite in manner. The use of a series of tubular specimens corresponding to different amounts of prior extension of cold reduction gave curves almost devoid of oscillation since the effects of variation of rest interval and extension spacing were removed. Comparison is also obtained between the variation of the several properties, as measured in tension and in shear. The rise of proof stress with extension is studied, and the work-hardening rates of the various metals evaluated. The ratio between the tensile and shear proof stresses for the various annealed and cold-worked tubular metals is likewise calculated. The influence of annealing or tempering temperature on the proof stresses and moduli for the cold-worked metals and for air-hardened 13:2 Cr-Ni steel is investigated. An improvement of elastic strength generally is obtained, without important loss of yield strength, by annealing at suitable temperature. The variation of the proof stress and modulus of elasticity with plastic deformation or annealing temperature is explained in terms of the relative dominance of three important factors: namely, (a) internal stress, (b) lattice-expansion or work-hardening, and (c) crystal reorientation. Effective values of Poisson's ratio were computed from tensile and shear moduli obtained on tubular specimens. The variation of Poisson's ratio with plastic deformation and annealing temperature is explained in terms of the degree of anisotropy produced by changes of (a) internal stress and (b) crystal orientation.

METALS - ELASTICITY↗

NASA NTRS1996

Multigrid, Fractional-Step Computation Of Flow

Speed of computer code solving three-dimensional Navier-Stokes equations of flow of incompressible fluid by fractional-step method increased significantly by use of multigrid procedures. In method, equations solved on general nonorthogonal curvilinear coordinate grid, using volume fluxes. At each time step, computations performed in two fractional steps. In first step, equations of conservation of momentum solved by use of gradient of pressure from previous time step via explicit approximate-factorization method, yielding approximate flow field that does not satisfy equation of conservation of mass. In second step, discrete Poisson-like equation with Neumann-type boundary conditions, formed by combining equations of conservation of momentum and mass, solved iteratively.

Kwak, Dochan↗

DOE OSTI2024

Pion gravitational form factors in the QCD instanton vacuum. I

The pion form factors of the QCD energy-momentum tensor (EMT) are studied in the instanton liquid model (ILM) of the QCD vacuum. In this approach, the breaking of conformal symmetry is encoded in the form of stronger-than-Poisson fluctuations in the number of instantons. For the trace of the EMT, it is shown that the gluonic trace anomaly term contributes half the pion mass, with the other half coming from the quark-mass-dependent σ term. The Q 2 dependence of the form factors is governed by glueball and scalar meson exchanges. For the EMT, the spin-0 (trace) and spin-2 (traceless rank-2 tensor) form factors are computed at next-to-leading order in the instanton density using effective quark operators. Relations between the gluon and quark contributions to the EMT form factors are derived. The form factors are also expressed in terms of the pion light-front wave functions in the ILM. The results at the low resolution scale of the inverse instanton size are evolved to higher scales using the renormalization group equation. The ILM results compare well with those of recent lattice QCD calculations. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1993

Probabilistic assessment of composite structures

A methodology and attendant computer code were developed and are used 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 is IPACS (Integrated Probabilistic Assessment of Composite Structures). IPACS can simulate both composite mechanics and composite structural behavior. Application to probabilistic composite mechanics is illustrated by its use 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 used to evaluate the uncertainties in the buckling of a composite plate, the stress concentration factor in a composite panel, and the vertical displacement and ply stress in a composite aircraft wing segment. IPACS' application to probabilistic design is illustrated by its use to assess the thin composite shell (pipe).

Chamis, C. C.↗