SEARCH · Engineering Papers

Results for “coordinate system”

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 973 records · Page 54

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Numerical solution of the Navier-Stokes equations for arbitrary 2-dimensional multi-element airfoils

Numerical solutions of the Navier-Stokes equations, with an algebraic turbulence model, for time-dependent two dimensional flow about multi-element airfoils were developed. Fundamental to these solutions was the use of numerically-generated boundary-conforming curvilinear coordinate systems to allow bodies of arbitrary shape to be treated. A general two dimensional grid generation code for multiple-body configuration was written as a part of this project and made available through the COSMIC code library.

Thompson, J. F.↗

NASA NTRS1984

Surface grid generation for wing-fuselage bodies

In the application of finite-difference methods to obtain numerical solutions of viscous compressible fluid flow about wing-fuselage bodies, it is advantageous to transform the governing equations to an idealized boundary-fitted coordinate system. The advantages are reduced computational complexity and added accuracy in the application of boundary conditions. The solution process requires that a grid be superimposed on the physical solution domain which corresponds to a uniform grid on a rectangular computational domain (uniform rectangular parallel-epiped). Grid generation is the determination of a one to one relationship between grid points in the physical domain and grid points in the computational domain. A technique for computing wing-fuselage surface grids using the Harris geometry and software for smooth-surface representation is described. Grid spacing control concepts which govern the relationship between the wing-fuselage surface and the computational grid are also presented.

Smith, R. E.↗

NASA NTRS1984

Interactive modeling using the PATRAN-G program

PATRAN-G is an interactive finite element generating and graphics display program. The program is currently executable on the VAX or PRIME computers and is compatible with 19 commercially available graphics terminals. The PATRAN program can be used in two modes: generation of a finite element model, or postprocessing of output data from a finite element analysis program. The finite element generating mode is divided into two phases. Phase I is the geometry construction. During this phase, the geometric boundaries of the model are defined and the material properties are input. Mirror imaging, translation and scaling of lines, surfaces, and volumes, and rotation of regions permit the geometry to be described efficiently. Points may be defined in alternate rectangular, cylindrical, or spherical coordinate systems. Also, PATRAN-G can calculate the intersection of lines or surfaces. This capability allows a user to easily model, for example, the intersection of a cylinder which cuts through a sphere at an angle. Once the geometry of the model has been described in terms of points, lines, surfaces, and volumes. Phase II is used to generate the computer finite element model.

Stalnaker, W. A.↗

NASA NTRS1984

The Z3 model of Saturns magnetic field and the Pioneer 11 vector helium magnetometer observations

Magnetic field observations obtained by the Pioneer 11 vector helium magnetometer are compared with the Z(sub 3) model magnetic field. These Pioneer 11 observations, obtained at close-in radial distances, constitute an important and independent test of the Z(sub 3) zonal harmonic model, which was derived from Voyager 1 and Voyager 2 fluxgate magnetometer observations. Differences between the Pioneer 11 magnetometer and the Z(sub 3) model field are found to be small (approximately 1%) and quantitatively consistent with the expected instrumental accuracy. A detailed examination of these differences in spacecraft payload coordinates shows that they are uniquely associated with the instrument frame of reference and operation. A much improved fit to the Pioneer 11 observations is obtained by rotation of the instrument coordinate system about the spacecraft spin axis by 1.4 degree. With this adjustment, possibly associated with an instrumental phase lag or roll attitude error, the Pioneer 11 vector helium magnetometer observations are fully consistent with the Voyager Z(sub 3) model.

Connerney, J. E. P.↗

NASA NTRS1984

Three-dimensional transonic potential flow about complex 3-dimensional configurations

An analysis has been developed and a computer code written to predict three-dimensional subsonic or transonic potential flow fields about lifting or nonlifting configurations. Possible condfigurations include inlets, nacelles, nacelles with ground planes, S-ducts, turboprop nacelles, wings, and wing-pylon-nacelle combinations. The solution of the full partial differential equation for compressible potential flow written in terms of a velocity potential is obtained using finite differences, line relaxation, and multigrid. The analysis uses either a cylindrical or Cartesian coordinate system. The computational mesh is not body fitted. The analysis has been programmed in FORTRAN for both the CDC CYBER 203 and the CRAY-1 computers. Comparisons of computed results with experimental measurement are presented. Descriptions of the program input and output formats are included.

Reyhner, T. A.↗

NASA NTRS1984

Three-dimensional unsteady viscous flow analysis over airfoil sections

A three-dimensional solution procedure for the approximate form of the Navier-Stokes equation was exercised in the two- and three-dimensional modes to compute the unsteady turbulent boundary layer on a flat plate corresponding to the data of Karlsson. The procedure is based on the use of a consistently split Linearized Block Implicit technique in conjunction with a QR operator scheme. New time-dependent upstream boundary conditions were developed that yielded realistic solutions for the interior in the vicinity of the upstream boundary. Comparisons of the computation employing these boundary conditions with the data indicate that both qualitative and quantitative agreement was obtained for the mean velocity and the in phase and out of phase components of the first harmonic of the velocity. In addition, the calculation gave results for the skin friction phase angle that had expected physical behavior for large distances downstream of the inflow boundary. For the three-dimensional case, the two-dimensional data of Karlsson was considered, but in a coordinate system skewed at 45 deg to the free stream direction. The results of the calculations were in excellent agreement with the data and the two-dimensional computations.

Weinberg, B. C.↗

NASA NTRS1984

A conservative treatment of zonal boundaries for Euler equation calculations

Finite-difference calculations require the generation of a grid for the region of interest. A zonal approach, wherein the given region is subdivided into zones and the grid for each zone is generated independently, makes the grid-generation process for complicated topologies and for regions requiring selective grid refinement a fairly simple task. This approach results in new boundaries within the given region, that is, zonal boundaries at the interfaces of the various zones. The zonal-boundary scheme (the integration scheme used to update the points on the zonal boundary) for the Euler equations must be conservative, accurate, stable, and applicable to general curvilinear coordinate systems. A zonal-boundary scheme with these desirable properties is developed in this study. The scheme is designed for explicit, first-order-accurate integration schemes but can be modified to accommodate second-order-accurate explicit and implicit integration schemes. Results for inviscid flow, including supersonic flow over a cylinder, blast-wave diffraction by a ramp, and one-dimensional shock-tube flow are obtained on zonal grids. The conservative nature of the zonal-boundary scheme permits the smooth transition of the discontinuities associated with these flows from one zone to another. The calculations also demonstrate the continuity of contour lines across zonal boundaries that can be achieved with the present zonal scheme.

Rai, M. M.↗

NASA NTRS1984

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic, and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching step sizes with the implicit shock fitting procedure. A correlation is found between the marching step size, Reynolds number, and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗

NASA NTRS1984

Nonlinear computation of wing-body-vertical tail-wake flows at low supersonic speeds

A numerical method based on the conservation form of the full potential equation has been applied to the problem of three-dimensional supersonic flows with embedded subsonic regions. The governing equation is cast in a nonorthogonal coordinate system, and the theory of characteristics is used to accurately monitor the type-dependent flow field. A conservative switching scheme is employed to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. The newly developed computer program can handle arbitrary geometries with fuselage, wing, vertical tail and wake components at combined angle of attack and sideslip. Results are presented for a low supersonic Mach numbers flow over the Shuttle orbiter (including the OMS pods and vertical tail), and for flows over a realistic fighter type configuration. Comparisons with experimental data are shown to be in good agreement for various cases.

Szema, K.-Y.↗

NASA NTRS1984

A Newton multigrid method for the Euler equations

A multigrid method is used to apply Newton's method to the Euler equations in a two dimensional curvilinear coordinate system. The objective is to obtain rapid convergence for steady state problems. Solutions computed with the method evolve in a non-time-like manner. Stable pressure distributions typically develop in eight to ten Newton-multigrid steps, which is equivalent to the computational work of about 70 iterations with a factored implicit algorithm.

Childs, R. E.↗

NASA NTRS1983

A compressible model for the simulation of moist mountain waves

A two-dimensional, nonlinear, nonhydrostatic model is described which allows the calculation of moist airflow in mountainous terrain. The model is compressible, uses a terrain-following coordinate system, and employs lateral and upper boundary conditions which minimize wave reflections. The model's accuracy and sensitivity are examined. These tests suggest that in numerical simulations of vertically propagating, highly nonlinear mountain waves, a wave absorbing layer does not accurately mimic the effects of wave breakdown and dissipation at high levels in the atmosphere. In order to obtain a correct simulation, the region in which the waves are physically absorbed must generally be included in the computational domain (a nonreflective upper boundary condition should be used as well). The utility of the model is demonstrated in two examples (linear waves in a uniform atmosphere and the Jan. 11, 1972 Boulder windstorm) which illustrate how the presence of moisture can influence propagating waves. In both cases, the addition of moisture to the upstream flow greatly reduces the wave response.

Durran, D. R.↗

NASA NTRS1984

Barriers and dispersal surfaces in minimum-time interception

A method is proposed for mapping the barrier, dispersal, and control-level surfaces for a class of minimum-time interception and pursuit-evasion problems. Minimum-time interception of a target moving in a horizontal plane is formulated in a coordinate system whose origin is at the interceptor's terminal position and whose x-axis is along the terminal line of sight. This approach makes it possible to discuss the nature of the interceptor's extremals, using its extremal trajectory maps (ETMs), independently of target motion. The game surfaces are constructed by drawing sections of the isochrones, or constant minimum-time loci, from the interceptor and target ETMs. In this way, feedback solutions for the optimal controls are obtained. An example involving the interception of a target moving in a straight line at constant speed is presented.

Rajan, N.↗

NASA NTRS1984

The effect of convection upon charged particle transport in random magnetic fields

In a coordinate system moving with the plasma and random magnetic fields of a wind that blows with constant velocity in the direction of the guiding field, transport of energetic particles is described by a Boltzmann equation which is similar to the one that describes unconvected transport. Although this formulation is mathematically identical to that developed by Luhmann, which refers to the system where the guiding field is static, there are both practical and fundamental reasons to adopt the new approach. It leads to first-order approximate transport equations which are similar to those that apply in the absence of convection. However, these equations are more general than Parker's description of diffusion and convection, for they describe the coherent modes of transport that appear when the mean free path is large compared to the scale length for spatial variations of the guiding field, and they are valid for arbitrary wind velocity. The latter characteristic opens up new possibilities for analyzing particle transport in relativistic flows seen in some astronomical objects.

Earl, J. A.↗

NASA NTRS1984

Flap-lag damping of an elastic rotor blade with torsion and dynamic inflow in hover from symbolically generated equations

A program written in FORTRAN IV and coded into subroutines for subsequent numerical study is used to derive symbolically the governing equations of motion of the elastic rotor blade and the dynamic inflow equations. By rearranging the coefficients of the equations, the blade and dynamic inflow equations are converted into equations in a multiblade coordinated system. The final multiblade equations can accommodate any number of elastic blade modes. It is noted that with a given data set, the entire process, from the derivation of the equations to the numerical calculations, is computerized and requires only limited user interface. To demonstrate the usefulness of the approach, damping data in hover with and without dynamic inflow effects are presented for various rotor blade models, including flap-lag, flap-lag-torsion, quasi-static torsion, and rigid blade lag-flap.

Reddy, T. S. R.↗

NASA NTRS1984

Numerical simulation of the viscous flow fields over three-dimensional complicated geometries

A 'thin-layer' Navier-Stokes code capable of predicting steady state viscous flows is applied to complicated three-dimensional flow fields. The code is written in a generalized coordinate system, and a recently developed grid generation procedure is used for the flow-field discretization. Application is made to the vortical flow over a delta wing at high angle of attack, and the computed results are compared with experimental results. The results indicate that the present method can capture the physical phenomenon well. It is observed that a leading-edge separation vortex is formed over the wing, as is a secondary separation vortex near the leading edge. The flow field behind the trailing edge is also well described. Application is also made to the transonic flow over the Shuttle configuration. The result appears to be reasonable even though no experimental data are available for comparison. These results indicate that the present approach is capable of computing complicated three-dimensional flow fields.

Fujii, K.↗

NASA NTRS1984

Observations of the dynamics of the polar thermosphere

The dynamics of the polar thermosphere are examined by using observations made from the Dynamics Explorer 2 satellite. The results used in this study were obtained primarily from the Fabry-Perot interferometer (FPI) and the wind and temperature spectrometer (WATS) during the time period from September 1981 through January 1982. Two primary geophysical conditions were examined: these were the southern summer and the northern winter polar regions. The results support the conclusion that above 60 degrees of latitude the neutral winds are strongly controlled by ion/neutral frictional momentum transfer resulting from magnetospheric convection. This implies that the natural coordinate system within which to display the neutral winds in the high polar thermosphere is magnetic. The collected observations of this study were used to assess the validity of two of the large thermospheric general circulation models. The result of this assessment was that the models reasonably represent the vector winds at high altitudes but do not, at present, accurately simulate the thermodynamics of that regime.

Hays, P. B.↗

NASA NTRS1984

An analytic solution for the orbital perturbations of the Venus Radar Mapper due to gravitational harmonics

Hill's variational equations are solved analytically for the orbital perturbations of a spacecraft nominally in an elliptic orbit around a non-spherical body. The rotation of the central planet about its spin-axis is not considered in the analysis. The perturbations are restricted to the planetary gravitational harmonics only. An extremely simple algorithm is derived to transform the spherical harmonic potentials to the orbital coordinate system, and the resulting accelerations are shown to be simply trigonometric functions of the true anomaly. With the principal matrix solution for the differential equations of the adjoint system given in closed form, the orthogonality of the trigonometric functions makes it possible to obtain an analytic solution for the non-homogeneous problem, at intervals of 2 pi in true anomaly. The solution for orbital perturbations can be extended over several revolutions by applying well-known results from Floquet's theory. The technique is demonstrated with results presented on the spacecraft periapsis altitude for the forthcoming Venus Radar Mapper Mission.

Vijayaraghavan, A.↗

NASA NTRS1984

Penalty function finite element analysis of steady viscous incompressible flow in rotating coordinates

Finite element methods for incompressible viscous flow in turbomachines have not been presented in the literature previously. This paper develops a penalty function primitive variable method including Coriolis and centrifugal force terms for steady flow in a rotating coordinate system. Simplex elements are used with the result of solution times comparable to equivalent finite different solutions. Example cases considered are Couette flow, Poiseuile flow, flow over a step and flow in a rotating channel. Both laminar and turbulent flows are discussed. The accuracy of computed solutions compares well with theoretical solutions and experimenal measurements.

Rosen, M. C.↗