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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Heat transfer of phase-change materials in two-dimensional cylindrical coordinates

Two-dimensional phase-change problem is numerically solved in cylindrical coordinates (r and z) by utilizing two Taylor series expansions for the temperature distributions in the neighborhood of the interface location. These two expansions form two polynomials in r and z directions. For the regions sufficiently away from the interface the temperature field equations are numerically solved in the usual way and the results are coupled with the polynomials. The main advantages of this efficient approach include ability to accept arbitrarily time dependent boundary conditions of all types and arbitrarily specified initial temperature distributions. A modified approach using a single Taylor series expansion in two variables is also suggested.

Labdon, M. B.↗
NASA NTRS2000

A New Method for Accurate Treatment of Flow Equations in Cylindrical Coordinates Using Series Expansions

The motivation of this work is the ongoing effort at the Center for Turbulence Research (CTR) to use large eddy simulation (LES) techniques to calculate the noise radiated by jet engines. The focus on engine exhaust noise reduction is motivated by the fact that a significant reduction has been achieved over the last decade on the other main sources of acoustic emissions of jet engines, such as the fan and turbomachinery noise, which gives increased priority to jet noise. To be able to propose methods to reduce the jet noise based on results of numerical simulations, one first has to be able to accurately predict the spatio-temporal distribution of the noise sources in the jet. Though a great deal of understanding of the fundamental turbulence mechanisms in high-speed jets was obtained from direct numerical simulations (DNS) at low Reynolds numbers, LES seems to be the only realistic available tool to obtain the necessary near-field information that is required to estimate the acoustic radiation of the turbulent compressible engine exhaust jets. The quality of jet-noise predictions is determined by the accuracy of the numerical method that has to capture the wide range of pressure fluctuations associated with the turbulence in the jet and with the resulting radiated noise, and by the boundary condition treatment and the quality of the mesh. Higher Reynolds numbers and coarser grids put in turn a higher burden on the robustness and accuracy of the numerical method used in this kind of jet LES simulations. As these calculations are often done in cylindrical coordinates, one of the most important requirements for the numerical method is to provide a flow solution that is not contaminated by numerical artifacts. The coordinate singularity is known to be a source of such artifacts. In the present work we use 6th order Pade schemes in the non-periodic directions to discretize the full compressible flow equations. It turns out that the quality of jet-noise predictions using these schemes is especially sensitive to the type of equation treatment at the singularity axis. The objective of this work is to develop a generally applicable numerical method for treating the singularities present at the polar axis, which is particularly suitable for highly accurate finite-differences schemes (e.g., Pade schemes) on non-staggered grids. The main idea is to reinterpret the regularity conditions developed in the context of pseudo-spectral methods. A set of exact equations at the singularity axis is derived using the appropriate series expansions for the variables in the original set of equations. The present treatment of the equations preserves the same level of accuracy as for the interior scheme. We also want to point out the wider utility of the method, proposed here in the context of compressible flow equations, as its extension for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical coordinates with finite-differences schemes of various level of accuracy is straightforward. The robustness and accuracy of the proposed technique is assessed by comparing results from simulations of laminar forced-jets and turbulent compressible jets using LES with similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r = 0.

Constantinescu, G.S.↗
NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗
NASA NTRS1979

A computational fluid dynamic technique valid at the centerline for non-axisymmetric problems in cylindrical coordinates

A technique is described for the numerical solution of non-axisymmetric flow problems posed in cylindrical coordinates when the z-axis is included in the flowfield. The highlight of the technique is the manner in which the singularities at the centerline are handled. Specifically, the governing flowfield equations at r = 0 are put in a special form by applying L'Hospital's Rule. The required radial derivatives are evaluated using a one-sided, second-order accurate, first-difference. This leads to a smooth, convergent calculation of the flowfield at the centerline. This appears to be the first generally applicable numerical method for avoiding coordinate system singularities in the context of a finite-difference scheme, and could have application to many nonaxisymmetric flows. The technique is illustrated by specific results for the time-dependent flowfield inside an internal combustion engine.

Griffin, M. D.↗
NASA NTRS1974

Liquefaction dynamics of n-octadecane in cylindrical coordinates

The study presents one method of mathematically modeling the effects of gravity-induced free convection upon the melting phenomena of a finite paraffin slab. The theoretical model is based on the numerical computer solution of the solid phase and liquid phase energy equations, coupled with an imposed ideal-viscous flow approximation. An implicit alternating-direction technique is used as the finite-difference approximation for the basic energy equations in the numerical solution. The theoretical predictions are compared to experimental data from an annular test cell, in which n-octadecane, C18H38, was used as the test material. The results show that the model is a good first approximation to the solution of solid-liquid phase change involving gravity-induced free convection.

Bain, R. L.↗
NASA NTRS1992

Radiative transfer solutions for cylindrical coordinates with emitting, absorbing, and anisotropic scattering medium

The radiative transfer equation is solved by the S-N discrete ordinates method in one-dimensional radial and two-dimensional r-z coordinates systems. The walls of the enclosures are diffuse, and the participating medium absorbs, emits and anisotropically scatters the radiative energy. Diffuse wall incidence, isothermal medium emission and collimated incidence problems are considered. Effects of the scattering phase functions on average incident radiation and net radiative heat fluxes are studied. In addition, the effects of scattering albedo, optical thickness and the wall emissivity are briefly discussed.

Jendoubi, Slah↗
NASA NTRS1992

A collection of edge-based elements

Edge-based elements have proved useful in solving electromagnetic problems since they are nondivergent. Previous authors have presented several two and three dimensional elements. Herein, we present four types of elements which are suitable for modeling several types of three dimensional geometries. Distorted brick and triangular prism elements are given in cartesian coordinates as well as the specialized cylindrical shell and pie-shaped prism elements which are suitable for problems best described in polar cylindrical coordinates.

Kempel, Leo C.↗
NASA NTRS2001

Simple Analytic Expressions for the Magnetic Field of a Circular Current Loop

Analytic expressions for the magnetic induction (magnetic flux density, B) of a simple planar circular current loop have been published in Cartesian and cylindrical coordinates [1,2], and are also known implicitly in spherical coordinates [3]. In this paper, we present explicit analytic expressions for B and its spatial derivatives in Cartesian, cylindrical, and spherical coordinates for a filamentary current loop. These results were obtained with extensive use of Mathematica "TM" and are exact throughout all space outside of the conductor. The field expressions reduce to the well-known limiting cases and satisfy V · B = 0 and V x B = 0 outside the conductor. These results are general and applicable to any model using filamentary circular current loops. Solenoids of arbitrary size may be easily modeled by approximating the total magnetic induction as the sum of those for the individual loops. The inclusion of the spatial derivatives expands their utility to magnetohydrodynamics where the derivatives are required. The equations can be coded into any high-level programming language. It is necessary to numerically evaluate complete elliptic integrals of the first and second kind, but this capability is now available with most programming packages.

Simpson, James C.↗
NASA NTRS1996

Marching Cubes in Cylindrical and Spherical Coordinates

Isosurface extraction is a common analysis and visualization technique for three-dimensional scalar data. Marching Cubes is the most commonly-used algorithm for finding polygonal representations of isosurfaces in such data. We extend Marching Cubes to produce geometry for data sets that lie in spherical and cylindrical coordinate systems as well as show the steps for derivation of transformations for other coordinate systems.

techniques marching cubes isosurfaces coordinate s↗
NASA NTRS1973

Relaxation solutions for inviscid axisymmetric transonic flow over blunt or pointed bodies.

A finite-difference relaxation method is presented for numerical solution of the full potential equation and exact boundary conditions for general axisymmetric bodies is inviscid, steady transonic flow. Body-normal coordinates are used in the nose region and sheared cylindrical coordinates are used on the afterbody to accommodate corners such as boattails and flares. An improved difference scheme is used which does not require that the flow be nearly alined with a coordinate direction in supersonic regions, and which treats either subsonic or supersonic free streams. Numerical results are illustrated for some simple classical shapes such as spheres and ellipsoids, and for more practical shapes like tangent-ogives with boattails. Special attention is given to bodies which have been studied for area-rule applications. Agreement with available experimental results is good in cases where viscous effects and wind-tunnel wall interference are not important.

South, J. C., Jr.↗
NASA NTRS1982

Computation of transonic potential flow about 3 dimensional inlets, ducts, and bodies

An analysis was developed and a computer code, P465 Version A, written for the prediction of transonic potential flow about three dimensional objects including inlet, duct, and body geometries. Finite differences and line relaxation are used to solve the complete potential flow equation. The coordinate system used for the calculations is independent of body geometry. Cylindrical coordinates are used for the computer code. The analysis is programmed in extended FORTRAN 4 for the CYBER 203 vector computer. The programming of the analysis is oriented toward taking advantage of the vector processing capabilities of this computer. Comparisons of computed results with experimental measurements are presented to verify the analysis. Descriptions of program input and output formats are also presented.

Reyhner, T. A.↗
NASA NTRS1982

GRID3C: Computer program for generation of C type multilevel, three dimensional and boundary conforming periodic grids

A fast computer program, GRID3C, was developed for accurately generating periodic, boundary conforming, three dimensional, consecutively refined computational grids applicable to realistic axial turbomachinery geometries. The method is based on using two functions to generate two dimensional grids on a number of coaxial axisymmetric surfaces positioned between the centerbody and the outer radial boundary. These boundary fitted grids are of the C type and are characterized by quasi-orthogonality and geometric periodicity. The built in nonorthogonal coordinate stretchings and shearings cause the grid clustering in the regions of interest. The stretching parameters are part of the input to GRID3C. In its present version GRID3C can generate and store a maximum of four consecutively refined three dimensional grids. The output grid coordinates can be calculated either in the Cartesian or in the cylindrical coordinate system.

Dulikravich, D. S.↗
NASA NTRS1994

A 3D-CFD code for accurate prediction of fluid flows and fluid forces in seals

Current and future turbomachinery requires advanced seal configurations to control leakage, inhibit mixing of incompatible fluids and to control the rotodynamic response. In recognition of a deficiency in the existing predictive methodology for seals, a seven year effort was established in 1990 by NASA's Office of Aeronautics Exploration and Technology, under the Earth-to-Orbit Propulsion program, to develop validated Computational Fluid Dynamics (CFD) concepts, codes and analyses for seals. The effort will provide NASA and the U.S. Aerospace Industry with advanced CFD scientific codes and industrial codes for analyzing and designing turbomachinery seals. An advanced 3D CFD cylindrical seal code has been developed, incorporating state-of-the-art computational methodology for flow analysis in straight, tapered and stepped seals. Relevant computational features of the code include: stationary/rotating coordinates, cylindrical and general Body Fitted Coordinates (BFC) systems, high order differencing schemes, colocated variable arrangement, advanced turbulence models, incompressible/compressible flows, and moving grids. This paper presents the current status of code development, code demonstration for predicting rotordynamic coefficients, numerical parametric study of entrance loss coefficients for generic annular seals, and plans for code extensions to labyrinth, damping, and other seal configurations.

Athavale, M. M.↗
NASA NTRS1984

The computation of mesh configurations for three-dimensional flow analysis

A geometry modeling method was developed to define three-dimensional mesh configurations, which were expressed by the set of surface intersections with a coordinate mesh. These configurations define bodies for transonic potential-flow computation in simple physical coordinates. Various configurations were prepared for 3-D potential flow analysis, including many bare inlet models in cylindrical coordinates. Trimmed surface models were used to prepare nacelle installation configurations. A variety of ducts and mixers were represented. Analytically defined mesh/surface intersections were compared with similar data obtained from the corresponding surface model. Parameters for an analytic surface definition from mesh-surface intersections obtained from a model of the surface were produced. In both studies the modeling accuracy was shown to be compatible with the flow analysis requirements. This technique fulfills the requirements of the associated three-dimensional potential-flow analysis. It preserves the flow analysis method's superior flexibility for adapting to unusual geometries. The geometry modeling technique is supported by a sufficient set of tools for effective use by general users, and sufficient accuracy was demonstrated.

Gibson, S. G.↗
NASA NTRS1988

Applications of an exponential finite difference technique

An exponential finite difference scheme first presented by Bhattacharya for one dimensional unsteady heat conduction problems in Cartesian coordinates was extended. The finite difference algorithm developed was used to solve the unsteady diffusion equation in one dimensional cylindrical coordinates and was applied to two and three dimensional conduction problems in Cartesian coordinates. Heat conduction involving variable thermal conductivity was also investigated. The method was used to solve nonlinear partial differential equations in one and two dimensional Cartesian coordinates. Predicted results are compared to exact solutions where available or to results obtained by other numerical methods.

Handschuh, Robert F.↗
NASA NTRS1987

An Exponential Finite Difference Technique for Solving Partial Differential Equations

An exponential finite difference algorithm, as first presented by Bhattacharya for one-dimensianal steady-state, heat conduction in Cartesian coordinates, has been extended. The finite difference algorithm developed was used to solve the diffusion equation in one-dimensional cylindrical coordinates and applied to two- and three-dimensional problems in Cartesian coordinates. The method was also used to solve nonlinear partial differential equations in one (Burger's equation) and two (Boundary Layer equations) dimensional Cartesian coordinates. Predicted results were compared to exact solutions where available, or to results obtained by other numerical methods. It was found that the exponential finite difference method produced results that were more accurate than those obtained by other numerical methods, especially during the initial transient portion of the solution. Other applications made using the exponential finite difference technique included unsteady one-dimensional heat transfer with temperature varying thermal conductivity and the development of the temperature field in a laminar Couette flow.

Handschuh, Robert F.↗