Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rectangular 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 127 records · Page 7

Calculation of supersonic three-dimensional free-mixing flows using the parabolic-elliptic Navier-Stokes equations

A numerical method is presented which is valid for integration of the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to the three-dimensional supersonic flow of a jet issuing into a supersonic free stream. Difficulties associated with the imposition of free-stream boundary conditions are noted, and a coordinate transformation, which maps the point at infinity onto a finite value, is introduced to alleviate these difficulties. Results are presented for calculations of a square jet and varying-aspect-ratio rectangular jets. The solution behavior varies from axisymmetry for the square jet to nearly two-dimensional for the high-aspect-ratio rectangle, although the computation always calculates the flow as though it were truly three-dimensional.

Hirsh, R. S.↗

Studies of the effects of curvature on dilution jet mixing

An analytical program was conducted using both three-dimensional numerical and empirical models to investigate the effects of transition liner curvature on the mixing of jets injected into a confined crossflow. The numerical code is of the TEACH type with hybrid numerics; it uses the power-law and SIMPLER algorithms, an orthogonal curvilinear coordinate system, and an algebraic Reynolds stress turbulence model. From the results of the numerical calculations, an existing empirical model for the temperature field downstream of single and multiple rows of jets injected into a straight rectangular duct was extended to model the effects of curvature. Temperature distributions, calculated with both the numerical and empirical models, are presented to show the effects of radius of curvature and inner and outer wall injection for single and opposed rows of cool dilution jets injected into a hot mainstream flow.

Holdeman, James D.↗

A variational approach to three-dimensional elasticity solutions of laminated composite plates

The displacements in a laminated composite are represented as products of two sets of unknown functions, one of which is only a function of the thickness coordinate and the other is a function of the in-plane coordinates (i.e., separation of variables approach), and the minimization of the total potential energy is reduced to a sequence of iterative linear problems. Analytical solutions are developed for cross-ply and angle-ply laminated composite rectangular plates. The solution for simply-supported cross-ply plates under sinusoidal transverse load reduces to that of Pagano. Numerical results for stresses and is placements for antisymmetric angle-ply laminates are presented. The three-dimensional elasticity solutions developed are important because they can be used to study the behavior of composite laminates, in addition to serving as reference for approximate solutions by numerical methods and two-dimensional theories.

Savoia, M.↗

Computation of transonic flow about lifting wing-cylinder combinations

Some results are described for a lifting rectangular wing centrally located on a circular-cylindrical body. This simple configuration has been utilized in order to assess the merits of a mapping technique for wing-body configurations. The procedure employed makes use of a coordinate transformation to simplify specification of the surface boundary condition in the computation of the flow about the wing. The method can be extended to incorporate wing sweep, finite length body of noncircular cross section, and arbitrary wing placement; however, these extensions involve a considerable increase in complexity of the problem.

Klunker, E. B.↗

Development and applications of algorithms for calculating the transonic flow about harmonically oscillating wings

A finite difference method to solve the unsteady transonic flow about harmonically oscillating wings was investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The differential equation for the unsteady velocity potential is linear with spatially varying coefficients and with the time variable eliminated by assuming harmonic motion. An alternating direction implicit procedure was investigated, and a pilot program was developed for both two and three dimensional wings. This program provides a relatively efficient relaxation solution without previously encountered solution instability problems. Pressure distributions for two rectangular wings are calculated. Conjugate gradient techniques were developed for the asymmetric, indefinite problem. The conjugate gradient procedure is evaluated for applications to the unsteady transonic problem. Different equations for the alternating direction procedure are derived using a coordinate transformation for swept and tapered wing planforms. Pressure distributions for swept, untaped wings of vanishing thickness are correlated with linear results for sweep angles up to 45 degrees.

Ehlers, F. E.↗

Generation of three-dimensional body-fitted coordinates using hyperbolic partial differential equations

An efficient numerical mesh generation scheme capable of creating orthogonal or nearly orthogonal grids about moderately complex three dimensional configurations is described. The mesh is obtained by marching outward from a user specified grid on the body surface. Using spherical grid topology, grids have been generated about full span rectangular wings and a simplified space shuttle orbiter.

Steger, J. L.↗

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

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

Navier-Stokes calculations of hovering rotor flowfields

Unsteady, thin-layer Navier-Stokes equations written in rotor coordinates are solved using a partially flux-split, implicit numerical algorithm to calculate the flowfields of a hovering rotor blade at subsonic and transonic conditions. Numerical results are in good agreement with experimental data for both nonlifting and lifting rectangular blades. For the lifting case, the rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blades. Alternate methods of calculating the hovering rotor flowfields in blade-fixed mode that have the same circulation distribution as a hovering blade are explored. All of the results presented in this paper were computed on a Cray-2 supercomputer.

Srinivasan, G. R.↗

Numerical solution of supersonic three-dimensional free-mixing flows using the parabolic-elliptic Navier-Stokes equations

A numerical method is presented for solving the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to three-dimensional supersonic laminar jet flow issuing parallel with a supersonic free stream. A coordinate transformation is introduced which maps the boundaries at infinity into a finite computational domain in order to eliminate difficulties associated with the imposition of free-stream boundary conditions. Results are presented for an approximate circular jet, a square jet, varying aspect ratio rectangular jets, and interacting square jets. The solution behavior varies from axisymmetric to nearly two-dimensional in character. For cases where comparisons of the present results with those obtained from shear layer calculations could be made, agreement was good.

Hirsh, R. S.↗

Computation of Thin-Walled Prismatic Shells

We consider a prismatic shell consisting of a finite number of narrow rectangular plates and having in the cross-section a finite number of closed contours (fig. 1(a)). We shall assume that the rectangular plates composing the shell are rigidly joined so that there is no motion of any kind of one plate relative to the others meeting at a given connecting line. The position of a point on the middle prismatic surface is considered to be defined by the coordinate z, the distance to a certain initial cross-section z = O, end the coordinate s determining its position on the contour of the cross-section.

Vlasov, V. Z.↗

Numerical generation of two-dimensional grids by the use of Poisson equations with grid control at boundaries

A method for generating boundary-fitted, curvilinear, two dimensional grids by the use of the Poisson equations is presented. Grids of C-type and O-type were made about airfoils and other shapes, with circular, rectangular, cascade-type, and other outer boundary shapes. Both viscous and inviscid spacings were used. In all cases, two important types of grid control can be exercised at both inner and outer boundaries. First is arbitrary control of the distances between the boundaries and the adjacent lines of the same coordinate family, i.e., stand-off distances. Second is arbitrary control of the angles with which lines of the opposite coordinate family intersect the boundaries. Thus, both grid cell size (or aspect ratio) and grid cell skewness are controlled at boundaries. Reasonable cell size and shape are ensured even in cases wherein extreme boundary shapes would tend to cause skewness or poorly controlled grid spacing. An inherent feature of the Poisson equations is that lines in the interior of the grid smoothly connect the boundary points (the grid mapping functions are second order differentiable).

Sorenson, R. L.↗

Flow in serpentine coolant passages with trip strips

Under the subject contract, an effort is being conducted at Scientific Research Associates, Inc. (SRA) to obtain flow field measurements in the coolant passage of a rotating turbine blade with ribbed walls, both in the stationary and rotating frames. The data obtained will be used for validation of computational tools and assessment of turbine blade cooling strategies. The configuration of the turbine blade passage model is given, and the measuring plane locations are given. The model has a four-pass passage with three 180 turns. This geometry was chosen to allow analyses of the velocity measurements corresponding to the heat transfer results obtained by Wagner. Two passes of the passage have a rectangular cross-section of 1.0 in x 0.5 in. Another two passes have a square cross-section of 0.5 in x 0.5 in. Trips with a streamwise pitch to trip height (P/e) = 5 and trip height to coolant passage width (e/Z) = 0.1, were machined along the leading and trailing walls. These dimensions are typical of those used in turbine blade coolant passages. The trips on these walls are staggered by the half-pitch. The trips are skewed at +/- 45 deg, and this allows the effect of trip orientation to be examined. Experiments will be conducted with flow entering the model through the 1.0 in x 0.5 in rectangular passage (Configuration C) and the 0.5 in x 0. 5 in square passage (Configuration D) to examine the effect of passage aspect ratio. Velocity measurements were obtained with a Reynolds number (Re) of 25,000, based on the hydraulic diameter of and bulk mean velocity in the half inch square passage. The coordinate system used in presenting the results for configurations C and D, respectively, is shown. The first, second and third passes of the passage will be referred to as the first, second and third passages, respectively, in later discussion. Streamwise distance (x) from the entrance is normalized by the hydraulic diameter (D). Vertical (y) and tangential (z) distances are normalized by the half passage height (H) and width (Z), respectively. The x coordinate and U component are positive in the streamwise direction. The y coordinate and V component are positive against gravity. The z coordinate and W component are positive in the direction of rotation. The velocities are normalized by the bulk mean velocity (Ub) of 3.44 m/s based on the half-inch square passage. The contours of the 1.0 in x 0.5 in and 0.5 in x 0.5 in passages were evaluated from 11 x 30 and 9 x 30 measurement grids, respectively.

Tse, D. G.-N.↗

An approximate buckling analysis for rectangular orthotropic plates with centrally located cutouts

An approximate analysis for predicting buckling of rectangular orthotropic composite plates with centrally located cutouts is presented. In this analysis, prebuckling and buckling problems are converted from a two-dimensional to a one-dimensional system of linear differential equations with variable coefficients. The conversion is accomplished by expressing the displacements as series with each element containing a trigonometric function of one coordinate and a coefficient that is an arbitrary function of the other coordinate. Ordinary differential equations are then obtained from a variational principle. Analytical results obtained from the approximate analysis are compared with finite element analyses for isotropic plates and for specially orthotropic plates with central circular cutouts of various sizes. Experimental results for the specially orthotropic plates are also presented. In nearly all cases, the approximate analysis predicts the buckling mode shapes correctly and predicts the buckling loads to within a few percent of the finite element and experimental results.

Nemeth, M. P.↗

Hippocampal place-cell firing during movement in three-dimensional space

"Place" cells of the rat hippocampus are coupled to "head direction" cells of the thalamus and limbic cortex. Head direction cells are sensitive to head direction in the horizontal plane only, which leads to the question of whether place cells similarly encode locations in the horizontal plane only, ignoring the z axis, or whether they encode locations in three dimensions. This question was addressed by recording from ensembles of CA1 pyramidal cells while rats traversed a rectangular track that could be tilted and rotated to different three-dimensional orientations. Cells were analyzed to determine whether their firing was bound to the external, three-dimensional cues of the environment, to the two-dimensional rectangular surface, or to some combination of these cues. Tilting the track 45 degrees generally provoked a partial remapping of the rectangular surface in that some cells maintained their place fields, whereas other cells either gained new place fields, lost existing fields, or changed their firing locations arbitrarily. When the tilted track was rotated relative to the distal landmarks, most place fields remapped, but a number of cells maintained the same place field relative to the x-y coordinate frame of the laboratory, ignoring the z axis. No more cells were bound to the local reference frame of the recording apparatus than would be predicted by chance. The partial remapping demonstrated that the place cell system was sensitive to the three-dimensional manipulations of the recording apparatus. Nonetheless the results were not consistent with an explicit three-dimensional tuning of individual hippocampal neurons nor were they consistent with a model in which different sets of cells are tightly coupled to different sets of environmental cues. The results are most consistent with the statement that hippocampal neurons can change their "tuning functions" in arbitrary ways when features of the sensory input or behavioral context are altered. Understanding the rules that govern the remapping phenomenon holds promise for deciphering the neural circuitry underlying hippocampal function.

Non-NASA Center↗

Three-dimensional solutions for coupled thermoelectroelastic response of multilayered plates

Analytic three-dimensional solutions are presented for the coupled thermoelectroelastic response of multilayered hybrid composite plates. The plates consist of a combination of fiber-reinforced cross-ply and piezothermoelastic layers. Both the thermoelectroelastic static response and its sensitivity coefficients are computed. The sensitivity coefficients measure the sensitivity of the response to variations in different mechanical, thermal and piezoelectric material properties of the plate. A linear constitutive model is used, and the material properties are assumed to be independent of the temperature and the electric field. The plates are assumed to have rectangular geometry and special material symmetries. A mixed formulation is used with the fundamental unknowns consisting of the three transverse stress components; three displacement components; transverse component of the electric displacement field; electric potential; transverse heat flux component, and temperature change. Each of the fundamental unknowns is expressed in terms of a double Fourier series in the Cartesian surface coordinates. A state space approach is used to generate the static response and to evaluate the sensitivity coefficients. Extensive numerical results are presented showing the effects of variation in the geometric parameters of the plate on the different response quantities and their sensitivity coefficients.

Xu, Kangming↗

Effects of Geometric Azimuthal Asymmetries of the PPM Stack on Electron Beam Characteristics

The effects of geometric azimuthally asymmetric properties of a periodic permanent magnet (PPM) focusing stack on electron beam characteristics obtained using a fully three dimensional (3D) particle-in-cell (PIC) code will be presented. The simulation model, using MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm), incorporates 3D behavior of the beam immersed in static fields calculated directly from the exact geometry and material properties of the 3D magnetic focusing structure. The Hughes 8916H, 18-40 GHz helical TWT for the millimeter-wave power module (MMPM) was used as a prototype. Firstly, the effects of C-magnets used at the input and output of the TWT to allow for coupling of the RF signal into and out of the tube are considered. The 8916H input and output C-magnets differ because coaxial couplers are used at the input and waveguide couplers are used at the output The repositioning of the beam from its central axis due to the inclusion of the output C-magnet was found to be most significant. The modeled output C-magnet and its orientation in the Cartesian coordinate system is shown, and a two-dimensional beam profile including the output C-magnet is also shown. A table presents the shift of the beam center off the central axis relative to the average radius of the beam at the longitudinal points A, B and C designated on an enclosed figure. Secondly, the addition of shunts, or rectangular iron pieces applied manually by a skilled technician in order to improve beam transmission, is considered. The shunts are applied to the top of the tube; thus, azimuthal symmetry of the focusing stack is interrupted. Although shunts are typically added during RF focusing, they are also typically added at the input section of the tube where RF forces are minimal, making an electron optics analysis meaningful. Because several shunts are usually applied to one pole piece, the simulations have been simplified by modeling a half washer with the same radius and longitudinal length as a shunt over the entire x, positive-y half of the transverse plane. A modeled pole piece and shunt as described are shown. Lastly, in order to study the effects of magnet misalignments, a magnet in the PPM stack was arbitrarily chosen and adjusted so that its central axis was shifted both 0.7 percent and 1.0 percent of the magnet outer diameter in the positive-y direction. In practice, positioning the magnets so that their central axis is accurately aligned with the central axis of the tube is challenging. Thus, it is a strong possibility that one or more magnets will be misaligned relative to the tube central axis.

Kory, Carol L.↗