Search NASA⌕ Search

SEARCH · Search NASA

Results for “POLAR COORDINATE”

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 73 records · Page 4

Duct flows with swirl

The physics of the flow interaction between swirl and secondary flow was studied in duct bends relevant to the design of advanced aircraft nozzle systems. Both laminar and turbulent subsonic flows were investigated in generic duct bends for different amounts of swirl. The flow calculations are based on an economical three-dimensional spatial marching method employed in an existing computer code (PEPSIG). The computational method and code were extended to allow azimuthal periodicity and solutions in which the polar coordinate singularity occurs in the interior of the flow field. These extensions are needed to address swirling flow and twisted centerlines arising in out-of-plane bends. It was found that appropriate amounts of swirl can reduce total pressure loss relative to nonswirling cases. This conclusion was found to be insensitive to computational mesh.

Tsai, Tommy M.↗

Navier-Stokes analysis of two- and three-dimensional flow field in solid rocket motors with segment joints

A multidimensional implicit Navier-Stokes analysis which uses numerical solution of ensemble-averaged Navier-Stokes equations in a nonorthogonal bodyfitted cylindrical-polar coordinate system has been applied to simulation of the internal flow field in solid-propellant rocket motor chambers with segment joints. 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. Computations performed to simulate the axisymmetric flow field in the vicinity of the aft field joint in the Space Shuttle SRB using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot. Calculations were also performed for two cases involving asymmetric three-dimensional flow in the vicinity of the aft field joint in the SRB using 721,525 grid points to estimate circumferential velocities and pressure gradients at the joint.

Sabnis, J. S.↗

Optical Alignment Device For Laser Communication

Optical alignment device under development enables continuous tracking and coalignment of two beams of light. Intended primarily for laser-communication station, in which transmitted beam must be aligned with received beam to ensure transmitted beam falls on receiver at other station. Expected to consume less power and be smaller and less complicated than alignment shutter and drive previously used. Prism and filter separate two overlapping collimated light beams of different wavelength or polarization. Coordinates of two beams tracked on charge-coupled device to determine degree of directional misalignment between two beams.

Casey, William L.↗

Latitude dependence of eddy variances

The eddy variance of a meteorological field must tend to zero at high latitudes due solely to the nature of spherical polar coordinates. The zonal averaging operator defines a length scale: the circumference of the latitude circle. When the circumference of the latitude circle is greater than the correlation length of the field, the eddy variance from transient eddies is the result of differences between statistically independent regions. When the circumference is less than the correlation length, the eddy variance is computed from points that are well correlated with each other, and so is reduced. The expansion of a field into zonal Fourier components is also influenced by the use of spherical coordinates. As is well known, a phenomenon of fixed wavelength will have different zonal wavenumbers at different latitudes. Simple analytical examples of these effects are presented along with an observational example from satellite ozone data. It is found that geometrical effects can be important even in middle latitudes.

Bowman, Kenneth P.↗

Rarefied gas flow through two-dimensional nozzles

A kinetic theory analysis is made of the flow of a rarefied gas from one reservoir to another through two-dimensional nozzles with arbitrary curvature. The Boltzmann equation simplified by a model collision integral is solved by means of finite-difference approximations with the discrete ordinate method. The physical space is transformed by a general grid generation technique and the velocity space is transformed to a polar coordinate system. A numerical code is developed which can be applied to any two-dimensional passage of complicated geometry for the flow regimes from free-molecular to slip. Numerical values of flow quantities can be calculated for the entire physical space including both inside the nozzle and in the outside plume. Predictions are made for the case of parallel slots and compared with existing literature data. Also, results for the cases of convergent or divergent slots and two-dimensional nozzles with arbitrary curvature at arbitrary knudsen number are presented.

De Witt, Kenneth J.↗

The glideslope approach

The problems associated with optimum approaches to an orbiting spacecraft at the completion of rendezvous are exacerbated when the maneuvering spacecraft is to dock with such large and extended structures as the NASA Space Station, which also has a torque equilibrium attitude due to which the docking port does not point directly into the orbital track. Attention is presently given to a candidate 'operationally optimum' approach to a local vertical-local horizontal stabilized target; this glideslope approach is derived on the basis of Hill's relative-motion equations, which have been expanded to encompass constant relative accelerations in a closed form solution. The motion is described in polar coordinates.

Pearson, Don J.↗

On a class of unsteady three-dimensional Navier Stokes solutions relevant to rotating disc flows: Threshold amplitudes and finite time singularities

A class of exact steady and unsteady solutions of the Navier Stokes equations in cylindrical polar coordinates is given. The flows correspond to the motion induced by an infinite disc rotating with constant angular velocity about the z-axis in a fluid occupying a semi-infinite region which, at large distances from the disc, has velocity field proportional to (x,-y,O) with respect to a Cartesian coordinate system. It is shown that when the rate of rotation is large, Karman's exact solution for a disc rotating in an otherwise motionless fluid is recovered. In the limit of zero rotation rate a particular form of Howarth's exact solution for three-dimensional stagnation point flow is obtained. The unsteady form of the partial differential system describing this class of flow may be generalized to time-periodic equilibrium flows. In addition the unsteady equations are shown to describe a strongly nonlinear instability of Karman's rotating disc flow. It is shown that sufficiently large perturbations lead to a finite time breakdown of that flow whilst smaller disturbances decay to zero. If the stagnation point flow at infinity is sufficiently strong, the steady basic states become linearly unstable. In fact there is then a continuous spectrum of unstable eigenvalues of the stability equations but, if the initial value problem is considered, it is found that, at large values of time, the continuous spectrum leads to a velocity field growing exponentially in time with an amplitude decaying algebraically in time.

Hall, Philip↗

Theoretical hypervelocity ballistic limit for single or double plates using nonlinear modal analysis

The original research on the use on nonlinear vibration technique to solve for the hypervelocity ballistic limit for double plates is examined. Such structure is commonly found in typical Space Station design where the incoming space or man-made debris would be fragmented upon hitting the outer plate (shield) and the subsequent impact on the main wall would result in a much reduced damge of the space station or spacecraft. The existing few theoretical impact equations do not agree well with each other. The existing computer code bumper used at NASA-Johnson Space Center appears to predict an unconservative ballistic limit when compared with experimental data where the velocity ranges from 3 km/s to 8 km/s. Such unconservative prediction is unacceptable from a practical safe design point of view. The bumper code is based on Wilkinson's (1968) paper and his equations have not been improved nor modified even though they are viewed with suspicion due to lack of agreement with experiments. To make matters worse, there is not another theory which is better than Wilkinson's equation and the designers are forced to use purely empirical Nysmith (1969) or semiempirical equations developed by Cour-Palais in 1969. The Cour-Palais equations were later modified empirically in 1989. The purpose of the present investigation is to examine the many assumptions of the Wilkinson equation. An attempt is made to present design charts based on the modified-Wilkinson equation so that the designer can get a feel of the ranges of the parameters which are of interest and discard a huge range of parameters, thus, significantly reducing the number of test shots required. The analysis is based on a solution of the governing nonlinear differential equations for a plate, assuming axisymmetric behavior using polar coordinates.

Hui, David↗

The migration of a compound drop due to thermocapillarity

The quasistatic thermocapillary motion of a compound drop in an unbounded fluid possessing a uniform temperature gradient is analyzed. For completeness, gravitational effects are included in the treatment. The general model is formulated, and the equations for the concentric case are solved using spherical polar coordinates, while the eccentric case is handled using bispherical coordinates. Results are given for the velocity of the drop as well as that of the droplet with respect to the drop, along with useful approximations. Illustrative results are presented graphically for the thermocapillary migration of a compound drop in the special case when the droplet is a gas bubble. In addition to the velocities of the drop and the bubble, representative isotherms and streamlines also are presented which display interesting qualitative features.

Morton, David S.↗

Hypersonic merged layer blunt body flows with wakes

An attempt is made here to understand the basic physics of the flowfield with wake on a blunt body of revolution under hypersonic rarefied conditions. A merged layer model of flow is envisioned. Full steady-state Navier-Stokes equations in spherical polar coordinate system are computed from the surface with slip and temperature jump conditions to the free stream by the Accelerated Successive Replacement method of numerical integration. Analysis is developed for bodies of arbitrary shape, but actual computations have been carried out for a sphere and sphere-cone body. Particular attention is paid to set the limit of the onset of separation, wake closure, shear-layer impingement, formation and dissipation of the shocks in the flowfield. Validity of the results is established by comparing the present results for sphere with the corresponding results of the SOFIA code in the common region of their validity and with the experimental data.

Jain, Amolak C.↗

Noise Radiation by Instability Waves in Coaxial Jets

In this paper predictions are made for the noise radiation from supersonic coaxial jets. The noise in the downstream arc of a supersonic jet is dominated by highly directional radiation from the supersonically convecting large scale structures in the jet mixing layer. Since the mean flow is not described easily in terms of simple analytic functions, a numerical prediction is made for its development. The compressible Reynolds-averaged boundary layer equations in cylindrical polar coordinates are solved. A mixing length turbulence model is used. Empirical correlations are developed the effects of velocity and temperature ratio and Mach number. Both normal and inverted velocity profiles are considered. Comparisons with measurements for both single and coaxial jets show good agreement. The large scale structures are modeled as instability waves. The noise radiation generated by the instability waves is determined by a matching between the inner instability wave solution and the outer acoustic solution. Predictions are made for the differences between the noise radiated by coaxial jets with different operating conditions and a single equivalent jet with the same exit area, thrust, and mass-flow.

Dahl, Milo D.↗

Three-Dimensional Orbits of Earth Satellites, Including Effects of Earth Oblateness and Atmospheric Rotation

The principal purpose of the present paper is to present sets of equations which may be used for calculating complete trajectories of earth satellites from outer space to the ground under the influence of air drag and gravity, including oblateness effects, and to apply these to several examples of entry trajectories starting from a circular orbit. Equations of motion, based on an "instantaneous ellipse" technique, with polar angle as independent variable, were found suitable for automatic computation of orbits in which the trajectory consists of a number of revolutions. This method is suitable as long as the trajectory does not become nearly vertical. In the terminal phase of the trajectories, which are nearly vertical, equations of motion in spherical polar coordinates with time as the independent variable were found to be more suitable. In the first illustrative example the effects of the oblateness component of the earth's gravitational field and of atmospheric rotation were studied for equatorial orbits. The satellites were launched into circular orbits at a height of 120 miles, an altitude sufficiently high that a number of revolutions could be studied. The importance of the oblateness component of the earth's gravitational field is shown by the fact that a satellite launched at circular orbital speed, neglecting oblateness, has a perigee some 67,000 feet lower when oblateness forces are included in the equations of motion than when they are not included. Also, the loss in altitude per revolution is double that of a satellite following an orbit not subject to oblateness. The effect of atmospheric rotation on the loss of altitude per revolution was small. As might be surmised, the regression of the line of nodes as predicted by celestial mechanics is unchanged when drag is included. It is clear that the inclination of the orbital plane to the equator will be relatively unaffected by drag for no atmospheric rotation since the drag lies in the orbital plane in this case. With the inclusion of atmospheric rotation it was found that the inclination of the plane changed about one-millionth of a radian per revolution. Thus the prediction of the position of the orbital plane of an earth satellite is not complicated by the introduction of drag. The line of apsides, which without drag but with oblateness moves slowly in space, tends to move with the satellite when drag is included in the calculations. As a results, the usual linearized solutions based on oblateness alone must be basically altered when drag is included to take into account the rapid movement of the line of apsides. In the second illustrative example the final revolution was calculated to impact for a number of trajectories in an orbital plane inclined at 650 to the equator. Of particular interest is the large effect the oblateness gravitational field and atmospheric rotation can have on the impact point. For a value of CDA/m of unity, and for an initial downward angle at 80 miles altitude of 0.01 radian, such as might be utilized for manned re-entry, oblateness had an influence of about 300 miles in the impact point, and atmospheric rotation had about a 150-mile influence.

Nielsen, Jack N.↗

Morphology and the Strength of Intermolecular Contact in Protein Crystals

The strengths of intermolecular contacts (macrobonds) in four lysozyme crystals were estimated based on the strengths of individual intermolecular interatomic interaction pairs. The periodic bond chain of these macrobonds accounts for the morphology of protein crystals as shown previously. Further in this paper, the surface area of contact, polar coordinate representation of contact site, Coulombic contribution on the macrobond strength, and the surface energy of the crystal have been evaluated. Comparing location of intermolecular contacts in different polymorphic crystal modifications, we show that these contacts can form a wide variety of patches on the molecular surface. The patches are located practically everywhere on this surface except for the concave active site. The contacts frequently include water molecules, with specific intermolecular hydrogen-bonds on the background of non-specific attractive interactions. The strengths of macrobonds are also compared to those of other protein complex systems. Making use of the contact strengths and taking into account bond hydration we also estimated crystal-water interfacial energies for different crystal faces.

Matsuura, Yoshiki↗

Reflector Surface Error Compensation in Dual-Reflector Antennas

By probing the field on a small subreflector at a minimal number of points, the main reflector surface errors can be obtained and subsequently used to design a phase-correction subreflector that can compensate for main reflector errors. The compensating phase-error profile across the subreflector can be achieved either by a surface deformation or by the use of an array of elements such as patch antennas that can cause a phase shift between the incoming and outgoing fields. The second option is of primary interest here, but the methodology can be applied to either case. The patch array is most easily implemented on a planar surface. Therefore, the example of a flat subreflector and a parabolic main reflector (a Newtonian dual reflector system) is considered in this work. The subreflector is assumed to be a reflector array covered with patch elements. The phase variation on a subreflector can be detected by a small number of receiving patch elements (probes). By probing the phase change at these few selected positions on the subreflector, the phase error over the entire surface can be recovered and used to change the phase of all the patch elements covering the subreflector plane to compensate for main reflector errors. This is accomplished by using a version of sampling theorem on the circular aperture. The sampling is performed on the phase-error function on the circular aperture of the main reflector by a method developed using Zernike polynomials. This method is based upon and extended from a theory previously proposed and applied to reflector aperture integration. This sampling method provides for an exact retrieval of the coefficients of up to certain orders in the expansion of the phase function, from values on a specifically calculated set of points in radial and azimuthal directions in the polar coordinate system, on the circular reflector aperture. The corresponding points on the subreflector are then obtained and, by probing the fields at these points, a set of phase values is determined that is then transferred back to the main reflector aperture for recovering the phase function. Once this function is recovered, the corresponding phase function on the subreflector is calculated and used to compensate for main reflector surface errors. In going back and forth between sub and main reflectors, geometrical (ray) optics is employed, which even though it ignores edge diffraction and other effects, is shown to be accurate for phase recovery.

Jamnejad, Vahraz↗

The Eccentric Behavior of Nearly Frozen Orbits

Frozen orbits are orbits which have only short-period changes in their mean eccentricity and argument of periapse, so that they basically keep a fixed orientation within their plane of motion. Nearly frozen orbits are those whose eccentricity and argument of periapse have values close to those of a frozen orbit. We call them "nearly" frozen because their eccentricity vector (a vector whose polar coordinates are eccentricity and argument of periapse) will stay within a bounded distance from the frozen orbit eccentricity vector, circulating around it over time. For highly inclined orbits around the Earth, this distance is effectively constant over time. Furthermore, frozen orbit eccentricity values are low enough that these orbits are essentially eccentric (i.e., off center) circles, so that nearly frozen orbits around Earth are bounded above and below by frozen orbits.

spacecraft↗

Numerically Optimized Coronagraph Designs for the Habitable Exoplanet Imaging Mission (HabEx) Concept

The primary science goal of the Habitable Exoplanet Imaging Mission (HabEx), one of four candidate flagship missions under investigation, is to image and spectrally characterize Earth-like exoplanets. It is well known that pupil obscurations degrade coronagraphic performance and complicate coronagraph design, so HabEx is planned to have an off-axis, unobscured primary mirror. We utilize the circular symmetry of the aperture to investigate 1D-radial coronagraph optimization methods that are prohibitively time-consuming or intractable in 2D, such as diffractive pupil remapping and concurrent, multi-plane optimization. We also directly constrain sensitivities to dynamic, low-order Zernike aberrations, which are separable in polar coordinates and can thus be propagated as 1D-radial integrals. The mask technologies in our designs claim heritage from the extensive modeling and testbed experiments performed by the Wide-Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) project. In this paper, we detail our optimization methods and outline future work to complete our design survey.

Balasubramaian, Kunjithapatham↗

Computation of three-dimensional flow about aerobrake configurations

Ellipsoid, cone and cylinder aerobrake configurations are analyzed to provide comparison data between experimental and model predictions. An analytical model was devised to account for the shock layer ahead of the body and in the near-wake region in terms of the Navier-Stokes equations expressed in conformal polar and azimuthal-angle coordinates. Using polar coordinates simplified the equations by mapping the body onto a sphere, a procedure which also reduced the magnitude of the discretization errors. The equations are then solved using an alternating direction implicit (ADI) factorization technique. Computations were carried out for Mach 3-10 at various grid resolutions and compared with available wind tunnel data. The model generated pressure distributions, heat transfer coefficients and velocity profile data that agreed relatively well with experimental data at a reduced computational cost. Further work is necessary to identify the location of shocks and to model flows about asymmetric configurations.

Li, C. P.↗