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 37 records · Page 2

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

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↗

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↗

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

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

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

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↗

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

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

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

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

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

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

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↗

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↗

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

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

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