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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Software to model AXAF image quality

This draft final report describes the work performed under this delivery order from May 1992 through June 1993. The purpose of this contract was to enhance and develop an integrated optical performance modeling software for complex x-ray optical systems such as AXAF. The GRAZTRACE program developed by the MSFC Optical Systems Branch for modeling VETA-I was used as the starting baseline program. The original program was a large single file program and, therefore, could not be modified very efficiently. The original source code has been reorganized, and a 'Make Utility' has been written to update the original program. The new version of the source code consists of 36 small source files to make it easier for the code developer to manage and modify the program. A user library has also been built and a 'Makelib' utility has been furnished to update the library. With the user library, the users can easily access the GRAZTRACE source files and build a custom library. A user manual for the new version of GRAZTRACE has been compiled. The plotting capability for the 3-D point spread functions and contour plots has been provided in the GRAZTRACE using the graphics package DISPLAY. The Graphics emulator over the network has been set up for programming the graphics routine. The point spread function and the contour plot routines have also been modified to display the plot centroid, and to allow the user to specify the plot range, and the viewing angle options. A Command Mode version of GRAZTRACE has also been developed. More than 60 commands have been implemented in a Code-V like format. The functions covered in this version include data manipulation, performance evaluation, and inquiry and setting of internal parameters. The user manual for these commands has been formatted as in Code-V, showing the command syntax, synopsis, and options. An interactive on-line help system for the command mode has also been accomplished to allow the user to find valid commands, command syntax, and command function. A translation program has been written to convert FEA output from structural analysis to GRAZTRACE surface deformation file (.dfm file). The program can accept standard output files and list files from COSMOS/M and NASTRAN finite analysis programs. Some interactive options are also provided, such as Cartesian or cylindrical coordinate transformation, coordinate shift and scale, and axial length change. A computerized database for technical documents relating to the AXAF project has been established. Over 5000 technical documents have been entered into the master database. A user can now rapidly retrieve the desired documents relating to the AXAF project. The summary of the work performed under this contract is shown.

Ahmad, Anees↗
NASA NTRS1992

Radiated noise of ducted fans

The differences in the radiated acoustic fields of ducted and unducted propellers of the same thrust operating under similar conditions are investigated. An FEM model is created for the generation, propagation, and radiation of steady, rotor alone noise and exit guide vane interaction noise of a ducted fan. For a specified number of blades, angular mode harmonic, and rotor angular velocity, the acoustic field is described in a cylindrical coordinate system reduced to only the axial and radial directions. It is found that, contrary to the usual understanding of the Tyler and Sofrin (1962) result, supersonic tip speed rotor noise can be cut off if the tip Mach number is only slightly in excess of unity and if the number of blades is relatively small. If there are many blades, the fundamental angular mode number is large, and the Tyler and Sofrin result for thin annuli becomes more relevant. Shrouding of subsonic tip speed propellers is a very effective means of controlling rotor alone noise.

Eversman, Walter↗
NASA NTRS1993

Linear theory of fast reconnection at an X-type neutral point

A complete analytic theory of dynamic linear reconnection at an X-type neutral point is developed. An eigenmode analysis, using cylindrical coordinates centered on the neutral point, extends the work of Craig & McClymont (1991) to include nonazimuthally symmetric perturbations on the flux function. It is demonstrated that all physically significant disturbances, both reconnective and nonreconnective, decay resistively on a 'fast' time scale. The significance of the linear theory is discussed within the general context of steady state and dynamic reconnection studies. Disturbances can be expected to focus explosively in the vicinity of the neutral point. This suggests the formation of a 'flux pile-up' current layer in which the bulk of the magnetic energy is released as heat rather than kinetic energy of mass motion.

Craig, I. J. D.↗
NASA NTRS1992

A comparative study of dynamically expanding force-free, constant-alpha magnetic configurations with applications to magnetic clouds

We contrast two different solutions of the constant alpha, force-free MHD equation, both of which have been suggested as models for magnetic clouds: a solution in cylindrical coordinates and one in spherical coordinates. In line with the observation that magnetic clouds expand, we generalize these static models and construct their expanding counterparts. We find that expansion introduces in both cases a large asymmetry in the field strength signature which is in the same sense as that seen the the data, i.e. towards the leading edge of the cloud. We then do a least squares fit of the respective models to one-spacecraft data on a magnetic cloud. We find that the fitting routine converges in both cases. However, while purely formally we cannot distinguish between the two models using data from one spacecraft, the field components in the 'spherical' model have features not compatible with data on magnetic clouds.

Farrugia, C. J.↗
NASA NTRS1992

Three-dimensional calculations of rotor-airframe interaction in forward flight

A method for analyzing the mutual aerodynamic interaction between a rotor and an airframe model has been developed. This technique models the rotor implicitly through the source terms of the momentum equations. A three-dimensional, incompressible, laminar, Navier-Stokes solver in cylindrical coordinates was developed for analyzing the rotor/airframe problem. The calculations are performed on a simplified model at an advance ratio of 0.1. The airframe surface pressure predictions are found to be in good agreement with wind tunnel test data. Results are presented for velocity and pressure field distributions in the wake of the rotor.

Zori, Laith A. J.↗
NASA NTRS1993

Burnett solutions along the stagnation line of a cooled cylinder in low-density hypersonic flows

The Burnett equations in the cylindrical coordinate system have been applied to the hypersonic flow over a cylinder in low density environments. Through a local similarity transformation, the Burnett equations are degenerated into nonlinear third-order ordinary differential equations which are only valid along the stagnation line. These equations are solved by the multi-staged Runge-Kutta integration technique. Steady state solutions are asymptotically reached by the time-marching procedure. A newly-developed computer code is used to solve both the Navier-Stokes and the Burnett equations within the shock layer. The computed heat transfer coefficients are consistent with wind-tunnel and STS flight data. Comparisons with the DSMC solutions and the Boltzmann solutions show that the Burnett equations are more appropriate as the governing equations than the Navier-Stokes equations in the transitional flow regime.

Liaw, G. S.↗
NASA NTRS1992

Comparison of radiated noise from shrouded and unshrouded propellers

The ducted propeller in a free field is modeled using the finite element method. The generation, propagation, and radiation of sound from a ducted fan is described by the convened wave equation with volumetric body forces. Body forces are used to introduce the blade loading for rotating blades and stationary exit guide vanes. For an axisymmetric nacelle or shroud, the problem is formulated in cylindrical coordinates. For a specified angular harmonic, the angular coordinate is eliminated, resulting in a two-dimensional representation. A finite element discretization based on nine-node quadratic isoparametric elements is used.

Eversman, Walter↗
NASA NTRS1990

Skyline based terrain matching

Skyline-based terrain matching, a new method for locating the vantage point of stereo camera or laser range-finding measurements on a global map previously prepared by satellite or aerial mapping is described. The orientation of the vantage is assumed known, but its translational parameters are determined by the algorithm. Skylines, or occluding contours, can be extracted from the sensory measurements taken by an autonomous vehicle. They can also be modeled from the global map, given a vantage estimate from which to start. The two sets of skylines, represented in cylindrical coordinates about either the true or the estimated vantage, are employed as 'features' or reference objects common to both sources of information. The terrain matching problem is formulated in terms of finding a translation between the respective representations of the skylines, by approximating the two sets of skylines as identical features (curves) on the actual terrain. The search for this translation is based on selecting the longest of the minimum-distance vectors between corresponding curves from the two sets of skylines. In successive iterations of the algorithm, the approximation that the two sets of curves are identical becomes more accurate, and the vantage estimate continues to improve. The algorithm was implemented and evaluated on a simulated terrain. Illustrations and examples are included.

Page, Lance A.↗
NASA NTRS1995

Solutions of the benchmark problems by the dispersion-relation-preserving scheme

The 7-point stencil Dispersion-Relation-Preserving scheme of Tam and Webb is used to solve all the six categories of the CAA benchmark problems. The purpose is to show that the scheme is capable of solving linear, as well as nonlinear aeroacoustics problems accurately. Nonlinearities, inevitably, lead to the generation of spurious short wave length numerical waves. Often, these spurious waves would overwhelm the entire numerical solution. In this work, the spurious waves are removed by the addition of artificial selective damping terms to the discretized equations. Category 3 problems are for testing radiation and outflow boundary conditions. In solving these problems, the radiation and outflow boundary conditions of Tam and Webb are used. These conditions are derived from the asymptotic solutions of the linearized Euler equations. Category 4 problems involved solid walls. Here, the wall boundary conditions for high-order schemes of Tam and Dong are employed. These conditions require the use of one ghost value per boundary point per physical boundary condition. In the second problem of this category, the governing equations, when written in cylindrical coordinates, are singular along the axis of the radial coordinate. The proper boundary conditions at the axis are derived by applying the limiting process of r approaches 0 to the governing equations. The Category 5 problem deals with the numerical noise issue. In the present approach, the time-independent mean flow solution is computed first. Once the residual drops to the machine noise level, the incident sound wave is turned on gradually. The solution is marched in time until a time-periodic state is reached. No exact solution is known for the Category 6 problem. Because of this, the problem is formulated in two totally different ways, first as a scattering problem then as a direct simulation problem. There is good agreement between the two numerical solutions. This offers confidence in the computed results. Both formulations are solved as initial value problems. As such, no Kutta condition is required at the trailing edge of the airfoil.

Tam, Christopher K. W.↗
NASA NTRS1994

Gasdynamic models of the solar wind/interstellar medium interaction

The interaction between the solar wind and the interstellar medium is modeled self-consistently using numerical solutions of the time-dependent gasdynamic equations in spherical and cylindrical coordinates. For the results presented here it is assumed that the solar system moves through the surrounding medium with a supersonic velocity. After an initial (nonequilibrium) state has been specified, the numerical solution follows the evolution in time until the interaction relaxes to a dynamic equilibrium. As would be expected, the solutions show the formation of a bow shock upstream of the traveling solar system to deflect the interstellar plasma around the cavity created by the solar wind. A terminiation shock also forms to slow and compress the solar wind plasma. For the simulation in spherical coordinates, the downstream portion of the termination shock reaches equilibrium more than three times further from the Sun than the equilibrium distance to the termination shock on the upstream side.

Steinolfson, R. S.↗
NASA NTRS1995

Navier-Stokes calculations of rotor-airframe interaction in forward flight

A method for analyzing the mutual aerodynamic interaction between a rotor and an airframe model has been developed. This technique models the rotor implicitly through the source terms of the momentum equations. A three-dimensional, incompressible, laminar, Navier-Stokes solver in cylindrical coordinates was developed for analyzing the rotor-airframe problem. The calculations are performed on a simplified rotor-airframe model at an advance ratio of 0.1. The airframe surface pressure predictions are found to be in good agreement with wind tunnel test data. Results are also presented for velocity and pressure field distributions in the wake of the rotor.

Zori, Laith A. J.↗
NASA NTRS1995

Aspects of passive magnetic levitation based on high-T(sub c) superconducting YBCO thin films

Passive magnetic levitation systems reported in the past were mostly confined to bulk superconducting materials. Here we present fundamental studies on magnetic levitation employing cylindrical permanent magnets floating above high-T(sub c) superconducting YBCO thin films (thickness about 0.3 mu m). Experiments included free floating rotating magnets as well as well-established flexible beam methods. By means of the latter, we investigated levitation and drag force hysteresis as well as magnetic stiffness properties of the superconductor-magnet arrangement. In the case of vertical motion of the magnet, characteristic high symmetry of repulsive (approaching) and attractive (withdrawing) branches of the pronounced force-displacement hysteresis could be detected. Achievable force levels were low as expected but sufficient for levitation of permanent magnets. With regard to magnetic stiffness, thin films proved to show stiffness-force ratios about one order of magnitude higher than bulk materials. Phenomenological models support the measurements. Regarding the magnetic hysteresis of the superconductor, the Irie-Yamafuji model was used for solving the equation of force balance in cylindrical coordinates allowing for a macroscopic description of the superconductor magnetization. This procedure provided good agreement with experimental levitation force and stiffness data during vertical motion. For the case of (lateral) drag force basic qualitative characteristics could be recovered, too. It is shown that models, based on simple asymmetric magnetization of the superconductor, describe well asymptotic transition of drag forces after the change of the magnet motion direction. Virgin curves (starting from equilibrium, i.e. symmetric magnetization) are approximated by a linear approach already reported in literature only. This paper shows that basic properties of superconducting thin films allow for their application to magnetic levitation or - without need of levitation forces, e.g. microgravity - magnetic damping devices.

Schoenhuber, P.↗
NASA NTRS1995

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗
NASA NTRS1995

Rapid Numerical Simulation of Viscous Axisymmetric Flow Fields

A two-dimensional Navier-Stokes code has been developed for rapid numerical simulation of axisymmetric flow fields, including flow fields with an azimuthal velocity component. The azimuthal-invariant Navier-Stokes equations in a cylindrical coordinate system are mapped to a general body-fitted coordinate system, with the streamwise viscous terms then neglected by applying the thin-layer approximation. Turbulence effects are modeled using an algebraic model, typically the Baldwin-Lomax turbulence model, although a modified Cebeci-Smith model can also be used. The equations are discretized using central finite differences and solved using a multistage Runge-Kutta algorithm with a spatially varying time step and implicit residual smoothing. Results are presented for calculations of supersonic flow over a waisted body-of-revolution, transonic flow through a normal shock wave in a straight circular duct of constant cross sectional area, swirling supersonic (inviscid) flow through a strong shock in a straight radial duct, and swirling subsonic flow in an annular-to-circular diffuser duct. Comparisons between computed and experimental results are in fair to good agreement, demonstrating that the viscous code can be a useful tool for practical engineering design and analysis work.

Tweedt, Daniel L.↗
NASA NTRS1996

Nonlinear Analysis and Modeling of Tires

The objective of the study was to develop efficient modeling techniques and computational strategies for: (1) predicting the nonlinear response of tires subjected to inflation pressure, mechanical and thermal loads; (2) determining the footprint region, and analyzing the tire pavement contact problem, including the effect of friction; and (3) determining the sensitivity of the tire response (displacements, stresses, strain energy, contact pressures and contact area) to variations in the different material and geometric parameters. Two computational strategies were developed. In the first strategy the tire was modeled by using either a two-dimensional shear flexible mixed shell finite elements or a quasi-three-dimensional solid model. The contact conditions were incorporated into the formulation by using a perturbed Lagrangian approach. A number of model reduction techniques were applied to substantially reduce the number of degrees of freedom used in describing the response outside the contact region. The second strategy exploited the axial symmetry of the undeformed tire, and uses cylindrical coordinates in the development of three-dimensional elements for modeling each of the different parts of the tire cross section. Model reduction techniques are also used with this strategy.

Noor, Ahmed K.↗
NASA NTRS1996

An MHD Study of the Interaction Between the Solar Wind and the Interstellar Medium

The overall objective of this research program is to obtain a better understanding of the interaction between the solar wind and the interstellar medium through the use of numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations. The simulated results have been compared with observations where possible and with the results from previous analytic and numerical studies. The primary accomplishment of this project has been the development of codes for 2-D models in both spherical and cylindrical coordinates and the application of the codes to the solar wind/interstellar medium interaction. Computations have been carried out for both a relatively simple gas-dynamic interaction and a flow-aligned interstellar magnetic field. The results have been shown to compare favorably with models that use more approximations and to modify and extend the previous results as would be expected. The simulations have also been used along with a data analysis study to provide a quantitative estimate of the distance to the termination and bow shocks. Some of the specific topics that have been studied are: (1) gas dynamic models of the solar wind/interstellar medium interaction, (2) termination shock response to large-scale solar wind fluctuations, and (3) distances to the termination shock and heliopause. The main results from each of these studies are summarized. The results were published in three papers which are included as attachments.

Steinolfson, R. S.↗
NASA NTRS1996

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗
NASA NTRS1997

Turbulence Measurements in the Near Field of a Wingtip Vortex

The roll-up of a wingtip vortex, at Reynolds number based on chord of 4.6 million was studied with an emphasis on suction side and near wake measurements. The research was conducted in a 32 in. x 48 in. low-speed wind tunnel. The half-wing model had a semi-span of 36 in. a chord of 48 in. and a rounded tip. Seven-hole pressure probe measurements of the velocity field surrounding the wingtip showed that a large axial velocity of up to 1.77 U(sub infinity) developed in the vortex core. This level of axial velocity has not been previously measured. Triple-wire probes have been used to measure all components of the Reynolds stress tensor. It was determined from correlation measurements that meandering of the vortex was small and did not appreciably contribute to the turbulence measurements. The flow was found to be turbulent in the near-field (as high as 24 percent RMS w - velocity on the edge of the core) and the turbulence decayed quickly with streamwise distance because of the nearly solid body rotation of the vortex core mean flow. A streamwise variation of the location of peak levels of turbulence, relative to the core centerline, was also found. Close to the trailing edge of the wing, the peak shear stress levels were found at the edge of the vortex core, whereas in the most downstream wake planes they occurred at a radius roughly equal to one-third of the vortex core radius. The Reynolds shear stresses were not aligned with the mean strain rate, indicating that an isotropic-eddy-viscosity based prediction method cannot accurately model the turbulence in the cortex. In cylindrical coordinates, with the origin at the vortex centerline, the radial normal stress was found to be larger than the circumferential.

Chow, Jim↗