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

Galactic gamma rays and the cosmic-ray halo problem

Detailed models of diffused halos of various sizes are considered. In such models, the scale perpendicular to the plane has a strong effect in determining the distribution of cosmic rays. Radial distributions were calculated for cylindrical coordinate models. The implied gamma-ray longitude distributions were then calculated and compared with the SAS-2 data. Assuming the sources to be supernova remnants or pulsars, only cosmic-ray nucleon halo models with an upper limit scale height of about 3 kpc provide a good fit to the gamma ray data. Consideration of the gamma-ray latitude data gives a half thickness of 2 + or - 2 kpc for the cosmic ray electron halo.

Stecker, F. W.↗

An extension of MacCormack's method for flows with higher-order equations and in different configurations

The numerical scheme for the computation of a shock discontinuity developed by MacCormack has been extended to solve a number of differential equations, including cases explicitly containing higher-order derivatives: (1) Korteweg-de Vries equation with a term of third-order derivative, (2) a system of nonlinear equations governing nonsteady one-dimensional plasma flow in cylindrical coordinate, (3) equations of solar wind. Comparisons with previous results are made, if available, to illustrate the advantages of the present method. The question of convergence of the numerical calculation is discussed.

Ying, S. J.↗

The influence of the Martian polar caps on the diurnal tide

The large horizontal heating gradients that exist near the edge of the Martian polar caps during spring are shown to be capable of exciting large oscillations in the diurnal tide. To a lesser extent, the daily mass cycling between cap and atmosphere can also contribute. The calculations which demonstrate this are based on classical tidal theory as applied to the cylindrical coordinate system. This is done to facilitate the representation of the heating function. Results are presented for the horizontal surface winds only. They indicate a circulation at the cap edge somewhat analogous to the smaller scale terrestrial sea breeze. The amplitude of the zonal component is largest and is increased from 1 to 10 m/sec by the modeled influence of the polar cap. When coupled with the basic flow these cap-edge tides can produce strong surface winds during spring. Such a mechanism may contribute to the ability of the south polar cap winds to generate the local dust storms observed near the cap edge at this season.

Haberle, R. M.↗

A computer program for thermal radiation from gaseous rocket exhuast plumes (GASRAD)

A computer code is presented for predicting incident thermal radiation from defined plume gas properties in either axisymmetric or cylindrical coordinate systems. The radiation model is a statistical band model for exponential line strength distribution with Lorentz/Doppler line shapes for 5 gaseous species (H2O, CO2, CO, HCl and HF) and an appoximate (non-scattering) treatment of carbon particles. The Curtis-Godson approximation is used for inhomogeneous gases, but a subroutine is available for using Young's intuitive derivative method for H2O with Lorentz line shape and exponentially-tailed-inverse line strength distribution. The geometry model provides integration over a hemisphere with up to 6 individually oriented identical axisymmetric plumes, a single 3-D plume, Shading surfaces may be used in any of 7 shapes, and a conical limit may be defined for the plume to set individual line-of-signt limits. Intermediate coordinate systems may specified to simplify input of plumes and shading surfaces.

Reardon, J. E.↗

Numerical simulation of reentry flow around the Shuttle Orbiter including real gas effects

The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.

Li, C. P.↗

The tailward magnetopause field beyond 10 RE

The geomagnetic field in the near-earth magnetotail is made up of three major contributions related to the dipole field, the magnetotail current system, and the magnetopause current system exclusive of the return currents in the magnetotail current system. The present investigation is concerned with the third contribution. The calculation of the magnetospheric field is discussed, taking into account the calculation of the field inside the magnetopause on the basis of the numerical integration of the Biot Savart integral. The representation of the field for distances greater than 10 earth radii on the dark side of the magnetosphere is considered. The field internal to the magnetopause is expressed as the negative gradient of a scalar potential which is expressed as a sum of solutions to Laplace's equation in cylindrical coordinates. The origin of coordinates is the earth.

Beard, D. B.↗

'Coriolis resonance' within a rotating duct

An investigation of the unsteady disturbances of a fixed frequency within a radial duct rotating at a set speed is presented. The flow is assumed to be compressible, inviscid, and of a fluid which is a perfect gas. Equations are developed for the steady and the unsteady parts of the flow in cylindrical coordinates. The unsteady disturbances are expressed by Fourier decomposition in angular position, distance into the duct, and in time. It is found that a resonance is possible when the frequency of flow disturbances is twice the shaft-rotation frequency, considering only the radial and tangential disturbances and not the radial and circumferential disturbances. The particular point at which the resonance occurs indicates the occurrence is due to the Coriolis force, which is only present in the radial and tangential directions. It is noted that the Coriolis force can only be present in open-ended ducts, such as those found in centrifugal compressors.

Kurosaka, M.↗

Investigation of compositional segregation during unidirectional solidification of solid solution semiconducting alloys

Compositional segregation of solid solution semiconducting alloys in the radial direction during unidirectional solidification was investigated by calculating the effect of a curved solid liquid interface on solute concentration at the interface on the solid. The formulation is similar to that given by Coriell, Boisvert, Rehm, and Sekerka except that a more realistic cylindrical coordinate system which is moving with the interface is used. Analytical results were obtained for very small and very large values of beta with beta = VR/D, where V is the velocity of solidification, R the radius of the specimen, and D the diffusivity of solute in the liquid. For both very small and very large beta, the solute concentration at the interface in the solid C(si) approaches C(o) (original solute concentration) i.e., the deviation is minimal. The maximum deviation of C(si) from C(o) occurs for some intermediate value of beta.

Wang, J. C.↗

Prediction of high speed propeller flow fields using a three-dimensional Euler analysis

To overcome the limitations of classical propeller theory, a computer program, NASPROP-E, was developed which solves for the flow field surrounding a multibladed propeller and axisymmetric nacelle combination using a finite difference method. The governing equations are the three dimensional unsteady Euler equations written in a cylindrical coordinate system. They are marched in time until a steady state solution is obtained. The Euler equations require no special treatment to model the blade work vorticity. The equations are solved using an implicit approximate factorization method. Numerical results are presented which have greatly increased the understanding of high speed propeller flow fields. Numerical results for swirl angle downstream of the propeller and propeller power coefficient are higher than experimental results. The radial variation of coefficient are higher than experimental results. The radial variation of swirl angle, however, is in reasonable agreement with the experimental results. The predicted variation of power coefficient with blade angle agrees very well with data.

Bober, L. J.↗

Stress-derivative control of keystoning deformation in finite element codes

The stress derivative technique for control of keystoning deformation in under-integrated finte elements is based on expansion of the stress in a Taylor series about the element center and retention of additional terms beyond the constant stress term. It has the advantage over other control techniques that keystoning resistance is provided by actual rather than artificial material properties. Application of this technique to the quadrilateral ring elements used for modelling solids of revolution subjected to axisymmetric loads is described. In a cylindrical coordinate system additional terms appear in the formulation which must be dealt with in arriving at a workable keystoning control scheme.

Schulz, J. C.↗

Effect of boundary layers on solid walls in three-dimensional subsonic wind tunnels

A solution for the tunnel wall boundary layer effects for three-dimensional subsonic tunnels is presented. The model potentials are represented with simple singularities placed on the centerline of the tunnel and Laplace's equation in cylindrical coordinates is solved for either the conventional homogeneous slotted-wall boundary condition, the solid-wall viscous boundary condition, or a combination of them. The most pronounced wall boundary layer effect is on solid blockage for completely closed wind tunnels. Boundary layers on the wall reduce the blockage from the solid-wall, no-boundary-layer case in a manner similar to opening slots in a solid wall. Additionally, for solid-wall tunnel configurations, the streamline curvature interference factor is reduced by a significant amount, whereas the lift interference factor at the model station does not depend on the boundary layer parameter. For combination wall configurations, the slot effect of the horizontal walls dominates the viscous effect of the solid sidewalls.

Adcock, J. B.↗

Relativistic electron in curved magnetic fields

Making use of the perturbation method based on the nonlinear differential equation theory, the author investigates the classical motion of a relativistic electron in a class of curved magnetic fields which may be written as B=B(O,B sub phi, O) in cylindrical coordinates (R. phi, Z). Under general astrophysical conditions the author derives the analytical expressions of the motion orbit, pitch angle, etc., of the electron in their dependence upon parameters characterizing the magnetic field and electron. The effects of non-zero curvature of magnetic field lines on the motion of electrons and applicabilities of these results to astrophysics are also discussed.

An, S.↗

Difference equation technique for determining the evolution of radially symmetric ultrasonic fields

A simplified difference equation technique has been developed to calculate the three-dimensional evolution of radially symmetric, complex ultrasonic fields. The parabolic approximation to the Helmholtz equation is obtained by assuming that localized changes in the magnitude of the field with respect to distance along the axis of propagation are small compared to those with respect to radial displacement. The Helmholtz equation in cylindrical coordinates may then be approximated, after separating variables, using truncated series expansions. The resulting difference equation determines complex field values at discrete locations in a plane using corresponding adjacent locations in the preceding plane. Stability criteria which determine acceptable values for various parameters are discussed. The technique has been applied to the propagation of ultrasonic fields generated by two-dimensional focusing and non-focusing, uniform and Gaussian profile, piezoelectric transducers in both profiled and non-profiled media.

Dockery, G. D.↗

A two-layer diagnostic model of the long-term physical evolution of warm-core ring 82B

The present shipboard data in the 1982 time series on a single Gulf Stream warm core ring are composited on a cylindrical coordinate system following ring motion; 10 C-isotherm depth measurements are used in a two-layer model of the ring's structure to compute gradient current, kinetic energy, available potential energy (APE), and potential vorticity. The ring's evolution can be divided into two periods: April-late June, while the ring is isolated from strong Gulf Stream interaction, during which the ring loses APE at a low rate, and from July onward, following ring interactions with the Gulf Stream and topography with much higher energy loss. Within measurement errors, the potential vorticity at the center of the ring is conserved from April through August, implying nearly inviscid dynamics.

Olson, D. B.↗

Computation of three-dimensional, rotational flow through turbomachinery blade rows for improved aerodynamic design studies

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.

Subramanian, S. V.↗

Gamma-ray constraints on Na-22 yields in nova explosions

The longitudinal distribution of diffuse 1.275 MeV emission produced by the decay of interstellar Na-22 generated by a limited class of energetic novae-accreting ONeMg white dwarfs is modeled. The probability distributions for the locations of disk novae are derived as functions of Galactocentric cylindrical coordinates. The occurrence rate for novae of all classes in the Galactic disk is about 30/yr. An occurrence rate of about 16/yr is found for spheroid novae of all classes. By comparing these simulation results with HEAO 3 upper limits on diffuse 1.275 MeV emission, upper limits of (3-5.6) x 10 to the -7th solar mass are determined for the mean Na-22 yield for energetic novae generated by accreting ONeMg white dwarfs. These Na-22 limits constrain future models of explosive nucleosynthesis in energetic nova explosions.

Higdon, J. C.↗

Acoustic absorption by sunspots

The paper presents the initial results of a series of observations designed to probe the nature of sunspots by detecting their influence on high-degree p-mode oscillations in the surrounding photosphere. The analysis decomposes the observed oscillations into radially propagating waves described by Hankel functions in a cylindrical coordinate system centered on the sunspot. From measurements of the differences in power between waves traveling outward and inward, it is demonstrated that sunspots appear to absorb as much as 50 percent of the incoming acoustic waves. It is found that for all three sunspots observed, the amount of absorption increases linearly with horizontal wavenumber. The effect is present in p-mode oscillations with wavelengths both significantly larger and smaller than the diameter of the sunspot umbrae. Actual absorption of acoustic energy of the magnitude observed may produce measurable decreases in the power and lifetimes of high-degree p-mode oscillations during periods of high solar activity.

Braun, D. C.↗

Fully-coupled analysis of jet mixing problems. Three-dimensional PNS model, SCIP3D

Numerical procedures formulated for the analysis of 3D jet mixing problems, as incorporated in the computer model, SCIP3D, are described. The overall methodology closely parallels that developed in the earlier 2D axisymmetric jet mixing model, SCIPVIS. SCIP3D integrates the 3D parabolized Navier-Stokes (PNS) jet mixing equations, cast in mapped cartesian or cylindrical coordinates, employing the explicit MacCormack Algorithm. A pressure split variant of this algorithm is employed in subsonic regions with a sublayer approximation utilized for treating the streamwise pressure component. SCIP3D contains both the ks and kW turbulence models, and employs a two component mixture approach to treat jet exhausts of arbitrary composition. Specialized grid procedures are used to adjust the grid growth in accordance with the growth of the jet, including a hybrid cartesian/cylindrical grid procedure for rectangular jets which moves the hybrid coordinate origin towards the flow origin as the jet transitions from a rectangular to circular shape. Numerous calculations are presented for rectangular mixing problems, as well as for a variety of basic unit problems exhibiting overall capabilities of SCIP3D.

Wolf, D. E.↗