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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1987

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

NASA NTRS1988

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

NASA NTRS1988

A simple, analytic 3-dimensional downburst model based on boundary layer stagnation flow

A simple downburst model is developed for use in batch and real-time piloted simulation studies of guidance strategies for terminal area transport aircraft operations in wind shear conditions. The model represents an axisymmetric stagnation point flow, based on velocity profiles from the Terminal Area Simulation System (TASS) model developed by Proctor and satisfies the mass continuity equation in cylindrical coordinates. Altitude dependence, including boundary layer effects near the ground, closely matches real-world measurements, as do the increase, peak, and decay of outflow and downflow with increasing distance from the downburst center. Equations for horizontal and vertical winds were derived, and found to be infinitely differentiable, with no singular points existent in the flow field. In addition, a simple relationship exists among the ratio of maximum horizontal to vertical velocities, the downdraft radius, depth of outflow, and altitude of maximum outflow. In use, a microburst can be modeled by specifying four characteristic parameters, velocity components in the x, y and z directions, and the corresponding nine partial derivatives are obtained easily from the velocity equations.

Oseguera, Rosa M.↗

NASA NTRS1988

Numerical prediction of flow in slender vortices

The slender vortex approximation was investigated using the Navier-Stokes equations written in cylindrical coordinates. It is shown that, for free vortices without external pressure gradient, the breakdown length is proportional to the Reynolds number. For free vortices with adverse pressure gradients, the breakdown length is inversely proportional to the value of its gradient. For low Reynolds numbers, the predictions of the simplified system agreed well with the ones obtained from solutions of the full Navier-Stokes equations, whereas for high Reynolds numbers, the flow became quite sensitive to pressure fluctuations; it was found that the failure of the slender vortex equations corresponded to the critical condition as identified by Benjamin (1962) for inviscid flows. The predictions obtained from the approximating system were compared with available experimental results. For low swirl, a good agreement was obtained; for high swirl, on the other hand, upstream effects on the pressure gradient produced by the breakdown bubble caused poor agreement.

Reyna, Luis G.↗

NASA NTRS1989

Exponential Finite-Difference Technique

Report discusses use of explicit exponential finite-difference technique to solve various diffusion-type partial differential equations. Study extends technique to transient-heat-transfer problems in one dimensional cylindrical coordinates and two and three dimensional Cartesian coordinates and to some nonlinear problems in one or two Cartesian coordinates.

Handschuh, Robert F.↗

NASA NTRS1988

Investigation of single-event upset (SEU) in an advanced bipolar process

An extensive analytical and experimental study SEU in an advanced silicon bipolar process was made. The modeling used process and device parameters to model the SEU charge, collection, and circuit response derived from a special version of PISCES in cylindrical coordinates and SPICE, respectively. Data are reported for test cells of various sizes.

Zoutendyk, John A.↗

NASA NTRS1988

Lateral charge transport from heavy-ion tracks in integrated circuit chips

A 256K DRAM has been used to study the lateral transport of charge (electron-hole pairs) induced by direct ionization from heavy-ion tracks in an IC. The qualitative charge transport has been simulated using a two-dimensional numerical code in cylindrical coordinates. The experimental bit-map data clearly show the manifestation of lateral charge transport in the creation of adjacent multiple-bit errors from a single heavy-ion track. The heavy-ion data further demonstrate the occurrence of multiple-bit errors from single ion tracks with sufficient stopping power. The qualitative numerical simulation results suggest that electric-field-funnel-aided (drift) collection accounts for single error generated by an ion passing through a charge-collecting junction, while multiple errors from a single ion track are due to lateral diffusion of ion-generated charge.

Zoutendyk, J. A.↗

NASA NTRS1988

Thermal analysis of heat storage canisters for a solar dynamic, space power system

A thermal analysis was performed of a thermal energy storage canister of a type suggested for use in a solar receiver for an orbiting Brayton cycle power system. Energy storage for the eclipse portion of the cycle is provided by the latent heat of a eutectic mixture of LiF and CaF2 contained in the canister. The chief motivation for the study is the prediction of vapor void effects on temperature profiles and the identification of possible differences between ground test data and projected behavior in microgravity. The first phase of this study is based on a two-dimensional, cylindrical coordinates model using an interim procedure for describing void behavor in 1-g and microgravity. The thermal analysis includes the effects of solidification front behavior, conduction in liquid/solid salt and canister materials, void growth and shrinkage, radiant heat transfer across the void, and convection in the melt due to Marangoni-induced flow and, in 1-g, flow due to density gradients. A number of significant differences between 1-g and o-g behavior were found. This resulted from differences in void location relative to the maximum heat flux and a significantly smaller effective conductance in 0-g due to the absence of gravity-induced convection.

Wichner, R. P.↗

NASA NTRS1989

Application of spectral collocation techniques to the stability of swirling flows

The linearized stability equations in cylindrical coordinates of a Chebyshev spectral collocation method for the temporal and spatial stability of swirling flows are presently solved with the eigenvalues obtained through the use of the QZ routine. The algorithm thus created is robust and easily adaptable to a range of flow configurations encompassing internal and external flows with minor boundary condition application modifications. Accuracy and efficiency tests of the method are made for the cases of plane Poiseulle, rotating-pipe, and trailing line vortex flows.

Khorrami, Mehdi R.↗

NASA NTRS1989

Three-dimensional multigrid Navier-Stokes computations for turbomachinery applications

The fully three-dimensional, time-dependent compressible Navier-Stokes equations in cylindrical coordinates are presently used, in conjunction with the multistage Runge-Kutta numerical integration scheme for solution of the governing flow equations, to simulate complex flowfields within turbomechanical components whose pertinent effects encompass those of viscosity, compressibility, blade rotation, and tip clearance. Computed results are presented for selected cascades, emphasizing the code's capabilities in the accurate prediction of such features as airfoil loadings, exit flow angles, shocks, and secondary flows. Computations for several test cases have been performed on a Cray-YMP, using nearly 90,000 grid points.

Subramanian, S. V.↗

NASA NTRS1990

Evolution of midplate hotspot swells: Numerical solutions

The evolution of midplate hotspot swells on an oceanic plate moving over a hot, upwelling mantle plume is numerically simulated. The plume supplies a Gaussian-shaped thermal perturbation and thermally-induced dynamic support. The lithosphere is treated as a thermal boundary layer with a strongly temperature-dependent viscosity. The two fundamental mechanisms of transferring heat, conduction and convection, during the interaction of the lithosphere with the mantle plume are considered. The transient heat transfer equations, with boundary conditions varying in both time and space, are solved in cylindrical coordinates using the finite difference ADI (alternating direction implicit) method on a 100 x 100 grid. The topography, geoid anomaly, and heat flow anomaly of the Hawaiian swell and the Bermuda rise are used to constrain the models. Results confirm the conclusion of previous works that the Hawaiian swell can not be explained by conductive heating alone, even if extremely high thermal perturbation is allowed. On the other hand, the model of convective thinning predicts successfully the topography, geoid anomaly, and the heat flow anomaly around the Hawaiian islands, as well as the changes in the topography and anomalous heat flow along the Hawaiian volcanic chain.

Liu, Mian↗

NASA NTRS1990

An analytic model of axisymmetric mantle plume due to thermal and chemical diffusion

An analytic model of axisymmetric mantle plumes driven by either thermal diffusion or combined diffusion of both heat and chemical species from a point source is presented. The governing equations are solved numerically in cylindrical coordinates for a Newtonian fluid with constant viscosity. Instead of starting from an assumed plume source, constraints on the source parameters, such as the depth of the source regions and the total heat input from the plume sources, are deduced using the geophysical characteristics of mantle plumes inferred from modelling of hotspot swells. The Hawaiian hotspot and the Bermuda hotspot are used as examples. Narrow mantle plumes are expected for likely mantle viscosities. The temperature anomaly and the size of thermal plumes underneath the lithosphere can be sensitive indicators of plume depth. The Hawaiian plume is likely to originate at a much greater depth than the Bermuda plume. One suggestive result puts the Hawaiian plume source at a depth near the core-mantle boundary and the source of the Bermuda plume in the upper mantle, close to the 700 km discontinuity. The total thermal energy input by the source region to the Hawaiian plume is about 5 x 10(10) watts. The corresponding diameter of the source region is about 100 to 150 km. Chemical diffusion from the same source does not affect the thermal structure of the plume.

Liu, Mian↗

NASA NTRS1990

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗

NASA NTRS1989

Coolant side heat transfer with rotation: User manual for 3D-TEACH with rotation

This program solves the governing transport equations in Reynolds average form for the flow of a 3-D, steady state, viscous, heat conducting, multiple species, single phase, Newtonian fluid with combustion. The governing partial differential equations are solved in physical variables in either a Cartesian or cylindrical coordinate system. The effects of rotation on the momentum and enthalpy calculations modeled in Cartesian coordinates are examined. The flow of the fluid should be confined and subsonic with a maximum Mach number no larger than 0.5. This manual describes the operating procedures and input details for executing a 3D-TEACH computation.

Syed, S. A.↗

NASA NTRS1991

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗

NASA NTRS1991

Thermal modeling with solid/liquid phase change of the thermal energy storage experiment

A thermal model which simulates combined conduction and phase change characteristics of thermal energy storage (TES) materials is presented. Both the model and results are presented for the purpose of benchmarking the conduction and phase change capabilities of recently developed and unvalidated microgravity TES computer programs. Specifically, operation of TES-1 is simulated. A two-dimensional SINDA85 model of the TES experiment in cylindrical coordinates was constructed. The phase change model accounts for latent heat stored in, or released from, a node undergoing melting and freezing.

Skarda, J. Raymond Lee↗

NASA NTRS1991

A Local Condensation Analysis Representing Two-phase Annular Flow in Condenser/radiator Capillary Tubes

NASA's effort for the thermal environmental control of the Space Station Freedom is directed towards the design, analysis, and development of an Active Thermal Control System (ATCS). A two phase, flow through condenser/radiator concept was baselined, as a part of the ATCS, for the radiation of space station thermal load into space. The proposed condenser rejects heat through direct condensation of ATCS working fluid (ammonia) in the small diameter radiator tubes. Analysis of the condensation process and design of condenser tubes are based on the available two phase flow models for the prediction of flow regimes, heat transfer, and pressure drops. The prediction formulas use the existing empirical relationships of friction factor at gas-liquid interface. An attempt is made to study the stability of interfacial waves in two phase annular flow. The formulation is presented of a stability problem in cylindrical coordinates. The contribution of fluid viscosity, surface tension, and transverse radius of curvature to the interfacial surface is included. A solution is obtained for Kelvin-Helmholtz instability problem which can be used to determine the critical and most dangerous wavelengths for interfacial waves.

Karimi, Amir↗

NASA NTRS1991

On the position of the near-earth neutral sheet - A comparison of magnetic model predictions with empirical formulas

A detailed comparison of the near-earth neutral sheet position according to several approximations has been made. In particular, two empirical formulas for the neutral sheet location, recently derived from AMPTE/CCE data, have been compared with the corresponding positions derived from Tsyganenko's magnetic field models. Cylindrical coordinates referenced to the geodipole have been used, and the position of the neutral sheet has been identified by the reversal of the radial component of the magnetic field. The analysis reveals that, in the region X(GSM) = -9 to -5 earth radii, abs. value of Y(GSM) not greater than 5 earth radii, close agreement exists between the empirical formulas and the Tsyganenko models. Furthermore, the Tsyganenko models provide a representation for the neutral sheet position in the transition region between near-earth and distant tail models.

Peredo, Mauricio↗