Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cylindrical coordinates”

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

Hybrid solid element with a traction-free cylindrical surface

An eight node solid element with two parallel faces and one traction-free cylindrical surface is derived using the assumed stress hybrid method. Cylindrical coordinates are used so that the assumed stresses satisfy the equilibrium equations as well as the traction-free condition over the cylindrical boundary. In the limiting case of plane stress conditions the assumed stresses also satisfy the compatibility conditions. Example solutions have demonstrated the advantage of using this special element for analyzing solids with circular holes.

Pian, T. H. H.↗

A shock capturing application of the finite element method

The paper is concerned with the development of finite element algorithms for the solution of viscous compressible flow problems with possible embedded shocks and recirculation regions. As a first step, the calculation in Cartesian coordinates of uniform flow on a rectangular region which encounters an embedded oblique shock with known turning angle is considered. A code is developed which is then used for computation of the boattail plume simulator problem in cylindrical coordinates. Oblique shock calculations are performed and the results are compared with a known finite difference solution.

Cooke, C. H.↗

HOTCFGM-2D: A Coupled Higher-Order Theory for Cylindrical Structural Components with Bi-Directionally Components with Bi-Directionally Graded Microstructures

The objective of this two-year project was to develop and deliver to the NASA-Glenn Research Center a two-dimensional higher-order theory, and related computer codes, for the analysis and design of cylindrical functionally graded materials/structural components for use in advanced aircraft engines (e.g., combustor linings, rotor disks, heat shields, brisk blades). To satisfy this objective, two-dimensional version of the higher-order theory, HOTCFGM-2D, and four computer codes based on this theory, for the analysis and design of structural components functionally graded in the radial and circumferential directions were developed in the cylindrical coordinate system r-Theta-z. This version of the higher-order theory is a significant generalization of the one-dimensional theory, HOTCFGM-1D, developed during the FY97 for the analysis and design of cylindrical structural components with radially graded microstructures. The generalized theory is applicable to thin multi-phased composite shells/cylinders subjected to steady-state thermomechanical, transient thermal and inertial loading applied uniformly along the axial direction such that the overall deformation is characterized by a constant average axial strain. The reinforcement phases are uniformly distributed in the axial direction, and arbitrarily distributed in the radial and circumferential direction, thereby allowing functional grading of the internal reinforcement in the r-Theta plane. The four computer codes fgmc3dq.cylindrical.f, fgmp3dq.cylindrical.f, fgmgvips3dq.cylindrical.f, and fgmc3dq.cylindrical.transient.f are research-oriented codes for investigating the effect of functionally graded architectures, as well as the properties of the multi-phase reinforcement, in thin shells subjected to thermomechanical and inertial loading, on the internal temperature, stress and (inelastic) strain fields. The reinforcement distribution in the radial and circumferential directions is specified by the user. The thermal and inelastic properties of the individual phases can vary with temperature. The inelastic phases are presently modeled by the power-law creep model generalized to multi-directional loading (within fgmc3dq.cylindrical.f and fgmc3dq.cylindrical.transient.f for steady-state and transient thermal loading, respectively), and incremental plasticity and GVIPS unified viscoplasticity theories (within the steady-state loading versions fgmp3dq.cylindrical.f and fgmgvips3dq.cylindrical.f).

Pindera, Marek-Jerzy↗

Polyatomic molecule vibrations

Polyatomic molecule vibrations are analyzed as harmonic vibrations along normal coordinates. The energy eigenvalues are found for linear and nonlinear symmetric triatomic molecules for valence bond models of the potential function with arbitrary coupling coefficients; such models can usually be fitted to observed energy levels with reasonably good accuracy. Approximate normal coordinates for the H2O molecule are discussed. Degenerate vibrational modes such as occur in CO2 are analyzed and expressions for Fermi resonance between close-lying states of the same symmetry are developed. The bending modes of linear triatomic molecules are expressed in terms of Laguerre polynomials in cylindrical coordinates as well as in terms of Hermite polynomials in Cartesian coordinates. The effects of large-amplitude bending such as occur in the C3 molecule are analyzed, along with anharmonic effects, which split the usually degenerate bending mode energy levels. Finally, the vibrational frequencies, degeneracies, and symmetry properties of XY3, X2Y2, and XY4 type molecules are discussed.

Source record↗

Navier-Stokes analysis of solid propellant rocket motor internal flows

A multidimensional implicit Navier-Stokes analysis that uses numerical solution of the ensemble-averaged Navier-Stokes equations in a nonorthogonal, body-fitted, cylindrical coordinate system has been applied to the simulation of the steady mean flow in solid propellant rocket motor chambers. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. The code was validated by comparing computed results with the experimental data obtained in cylindrical-port cold-flow tests. The agreement between the computed and experimentally measured mean axial velocities is excellent. The axial location of transition to turbulent flow predicted by the two-equation (k-epsilon) turbulence model used in the computations also agrees well with the experimental data. Computations performed to simulate the axisymmetric flowfield in the vicinity of the aft field joint in the Space Shuttle solid rocket motor using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot.

Sabnis, J. S.↗

Developing Procedures to Implement Geometric Imperfections Beyond Right Circular Cylindrical Shells in Finite Element Method Models

Analysis of aerospace structures is frequently conducted using nominal dimensions and frequently assumes ideal conditions in loading, contact, constraints, et cetera. Off-nominal dimensions and nonideal conditions, however, are present in all structures. These are the result of widely ranging causes from coefficient of thermal expansion mismatches, manufacturing tooling anomalies, to assembly procedures that inadvertently alter the structure. Specifically, geometric imperfections can have potentially significant influence on the response of a structural test article observed in an experiment versus the response given by a numerical simulation. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS) was previously presented as a set of Python scripts to calculate and implement as-manufactured geometric midsurface and thickness imperfections into finite element method (FEM) shell models of nominally right circular cylinders. By taking advantage of the simple shape of a right circular cylinder, interpolations of the measured data points were able to be performed along directions that aligned to the cylindrical coordinate system axes of the entire structure. By taking advantage of the shell representation of the real structure as opposed to modeling using a continuum representation, the thickness variation was able to be implemented by shell section definitions instead of having to modify the position of multiple nodes in the thickness direction. Py_TIGIRS is a useful tool that established a procedural example on how to implement geometric imperfections in right circular cylindrical shell structures. Three new procedures, each expanded from concepts established in Py_TIGIRS, are proposed for various test-article designs and are intended to broaden the range of structures that can be modeled with measured geometric imperfections in the structural analysis community. Each test-article design introduces new challenges to successfully implement geometric imperfections into a FEM model. The first test-article design consists of a carbon fiber reinforced polymer square plate with a hat-shaped stiffener co-cured on one side. This test-article design was for a novel seven-point bend test that was also previously presented. Manufacturing and cure-cycle imperfections are observed using digital image correlation (DIC) techniques. As thermal expansion coefficient mismatches between the plate and stiffener materials were anticipated, a thermal analysis study with continuum shell and solid elements was conducted to capture the global shape observed prior to testing. The second test-article design is of a similar hat-stiffened plate configuration, but with a side length ratio near 3:1 with elongation in the stiffener direction. The test article was used to characterize the response to uniaxial compressive loading in the direction of the stiffener. Due to differing manufacturing steps, a thermal analysis like the one developed for the seven-point bend configuration was unable to mimic the observed geometric imperfections. Instead, a strategy based on applying deformations directly to the structure during analysis was developed for continuum shell and solid element representation of a stiffened panel.

Geometric imperfections↗

Linear magnetic spring and spring/motor combination

A magnetic spring, or a spring and motor combination, providing a linear spring force characteristic in each direction from a neutral position, in which the spring action may occur for any desired coordinate of a typical orthogonal coordinate system. A set of magnets are disposed, preferably symmetrically about a coordinate axis, poled orthogonally to the desired force direction. A second set of magnets, respectively poled opposite the first set, are arranged on the sprung article. The magnets of one of the sets are spaced a greater distance apart than those of the other, such that an end magnet from each set forms a pair having preferably planar faces parallel to the direction of spring force, the faces being offset so that in a neutral position the outer edge of the closer spaced magnet set is aligned with the inner edge of the greater spaced magnet set. For use as a motor, a coil can be arranged with conductors orthogonal to both the magnet pole directions and the direction of desired spring force, located across from the magnets of one set and fixed with respect to the magnets of the other set. In a cylindrical coordinate system having axial spring force, the magnets are radially poled and motor coils are concentric with the cylinder axis.

Patt, Paul J.↗

Steady temperature and density distributions in a gas containing heat sources

Computer program, STADDIG, is based on steady state, one dimensional heat transfer calculation using cylindrical coordinates. Program allows for conduction across gas and container walls. Heat is dissipated from walls by forced convection cooling with incompressible coolant. Heat sources are included in coolant, gas, and walls.

Davison, H. W.↗

Computer Program for Calculating the Temperature Field of Face Seals

FORTRAN 4 program for use on IBM 7094 computer for calculating temperature field of shaft seals is described. Shaft seals are composed of basically axisymmetric bodies. Cylindrical coordinate system is used as basis for analysis. Program is quite general and can be applied to variety of axisymmetric body problems.

Russell, T. E.↗

Analytical solutions to the problem of transient heat transfer in living tissue.

An analytical model of transient heat transfer in living biological tissue is considered. The model includes storage, generation, conduction, and convective transport of heat in the tissue. Solutions for rectangular and cylindrical coordinates are presented and discussed. Transient times for reaching the ?locally fully developed' temperature profile were found to be of the order of 5 to 25 min. These transients are dominated by a geometrical parameters and, to a lesser extent, by a parameter representing the ratio of heat supplied by blood flow to heat conducted in the tissue.

Shitzer, A.↗

Numerical solutions of reactive fluid flows during postignition transients in hybrid rocket systems.

A computational method has been developed for the study of the post-ignition transients in hybrid rocket systems. The particular system chosen consisted of a gaseous oxidizer flowing within a tube of solid fuel, resulting in heterogeneous combustion. With the appropriate assumptions, two-dimensional, time-dependent conservation equations were derived for the reacting gas phase, and for the solid phase, in a cylindrical coordinate system. These were then programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and the Crank-Nicolson scheme for the solid phase. Appropriate initial and boundary conditions were represented, including heat and mass conservation at the interface between gas and solid. Initially, no attempt was made to relate the recession rate at the surface to the surface temperature, or to include heat transfer by radiation. A simple case was selected for preliminary calculations, with aluminum and oxygen as fuel and oxidizer, and aluminum oxide as the product.

Hung, W. S. Y.↗

Further analytical study of hybrid rocket combustion

Analytical studies of the transient and steady-state combustion processes in a hybrid rocket system are discussed. The particular system chosen consists of a gaseous oxidizer flowing within a tube of solid fuel, resulting in a heterogeneous combustion. Finite rate chemical kinetics with appropriate reaction mechanisms were incorporated in the model. A temperature dependent Arrhenius type fuel surface regression rate equation was chosen for the current study. The governing mathematical equations employed for the reacting gas phase and for the solid phase are the general, two-dimensional, time-dependent conservation equations in a cylindrical coordinate system. Keeping the simplifying assumptions to a minimum, these basic equations were programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and, the Crank-Nicolson scheme for the solid phase.

Hung, W. S. Y.↗

Hydrodynamic water impact

The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.

Kettleborough, C. F.↗

Evolution of rotating interstellar clouds. I - Numerical techniques

A method is described for the numerical calculation of the hydrodynamic evolution of a self-gravitating configuration in two space dimensions with assumed axial symmetry. The calculation is formulated in cylindrical coordinates with respect to a moving Eulerian grid and is solved using explicit hydrodynamics combined with implicit radiative transfer. The physics included is appropriate for calculation of the collapse of a rotating protostellar cloud. The gravitational field is obtained by means of an alternating-direction iterative technique. Numerical tests to demonstrate the correctness of the method are presented for special cases.

Black, D. C.↗

The diffusion approximation. An application to radiative transfer in clouds

It is shown how the radiative transfer equation reduces to the diffusion equation. To keep the mathematics as simple as possible, the approximation is applied to a cylindrical cloud of radius R and height h. The diffusion equation separates in cylindrical coordinates and, in a sample calculation, the solution is evaluated for a range of cloud radii with cloud heights of 0.5 km and 1.0 km. The simplicity of the method and the speed with which solutions are obtained give it potential as a tool with which to study the effects of finite-sized clouds on the albedo of the earth-atmosphere system.

Arduini, R. F.↗

The galactic halo question: New size constraints from galactic gamma-ray data

The SAS-2 gamma-ray data is analyzed making use of recent CO line emission and other data for determining the large-scale distribution of galactic gas. A nonuniform distribution of cosmic rays in the galaxy is implied. This fact rules out large trapping halo models and extragalactic origin models. Detailed models of diffusion halos of various sizes perpendicular to the galactic plane are considered. In such models, the scale perpendicular to the plane has a strong effect in determining the radial distribution of cosmic rays. Such radial distributions are calculated for cylindrical coordinate models. The implied gamma-ray longitude distributions are then calculated and compared with the SAS-2 data for goodness-of-fit. Assuming the sources to be supernova remnants or pulsars, cosmic ray nucleon halo models with scale heights greater than 3 kpc are found to provide a poor fit to the gamma-ray longitude data (probability of 6% or less). Thin halo, or source dominated diffusion models are found to provide a good fit to the gamma-ray data, with an upper limit scale height of approximately 3 kpc.

Stecker, F. W.↗