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 109 records · Page 6

Predicted electrothermal deicing of aircraft blades

A finite difference method is presented for the transient two-dimensional simulation of an electrothermal de-icer pad of an aircraft wing or blade. The irregular geometry of the composite ice laden blade is handled by use of a body fitted coordinate transformation. By this approach the various blade layers are mapped into a set of stacked rectangular strips in which the numerical solution takes place. Several heat conduction examples are presented in order to demonstrate the accuracy of the numerical procedure. Ice melting time predictions are made and compared to earlier predictions where possible. Finally, a new graphical presentation of thermal results is shown.

Keith, T. G., Jr.↗

Sonic-box method employing local Mach number for oscillating wings with thickness

A computer program was developed to account approximately for the effects of finite wing thickness in the transonic potential flow over an oscillating wing of finite span. The program is based on the original sonic-box program for planar wing which was previously extended to include the effects of the swept trailing edge and the thickness of the wing. Account for the nonuniform flow caused by finite thickness is made by application of the local linearization concept. The thickness effect, expressed in terms of the local Mach number, is included in the basic solution to replace the coordinate transformation method used in the earlier work. Calculations were made for a delta wing and a rectangular wing performing plunge and pitch oscillations, and the results were compared with those obtained from other methods. An input quide and a complete listing of the computer code are presented.

Ruo, S. Y.↗

Two-Dimensional Vernier Scale

Modified vernier scale gives accurate two-dimensional coordinates from maps, drawings, or cathode-ray-tube displays. Movable circular overlay rests on fixed rectangular-grid overlay. Pitch of circles nine-tenths that of grid and, for greatest accuracy, radii of circles large compared with pitch of grid. Scale enables user to interpolate between finest divisions of regularly spaced rule simply by observing which mark on auxiliary vernier rule aligns with mark on primary rule.

Juday, Richard D.↗

Boundary perturbation method for free boundary problem in convectively cooled continuous casting

Novel mathematical techniques are here used to treat general boundary conditions in convectively cooled continuous ingot casting problems, applying a constant convective heat transfer coefficient at the ingot sides (so that the temperature distribution along the sides is unknown). The problem is first inverted to treat the physical coordinates as dependent variables of temperature and heat flow, mapping the ingot into a region that is nearly rectangular. A boundary perturbation method is then used to obtain an analytical solution in this near-rectangular region. Solidification interface shapes depend on two dimensionless parameters: one is a function of the casting velocity, while the other is dependent on the convective heat transfer coefficient at the cooled wall.

Siegel, R.↗

Boundary-fitted curvilinear coordinate systems for solution of partial differential equations on fields containing any number of arbitrary two-dimensional bodies

A method is presented for automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multi-connected two-dimensional region containing any number of arbitrarily shaped bodies. No restrictions are placed on the shape of the boundaries, which may even be time-dependent, and the approach is not restricted in principle to two dimensions. With this procedure the numerical solution of a partial differential system may be done on a fixed rectangular field with a square mesh with no interpolation required regardless of the shape of the physical boundaries, regardless of the spacing of the curvilinear coordinate lines in the physical field, and regardless of the movement of the coordinate system in the physical plane. A number of examples of coordinate systems and application thereof to the solution of partial differential equations are given. The FORTRAN computer program and instructions for use are included.

Thompson, J. F.↗

Finite-volume solution of the Euler equations for steady three-dimensional transonic flow

A numerical finite-volume method is proposed for solving the Euler equations of a steady three-dimensional flow. The subsonic regions embedded in steady supersonic flow about aircraft are computed by introducing a time-dependent coordinate system that maps the fuselage and wing into one coordinate surface. The suitability of the method to a delta-winged aircraft is demonstrated, and its applicability to rectangular-winged aircraft and problems in internal transonic flow is pointed out.

Rizzi, A.↗

Magnetohydrodynamics of atmospheric transients. II - Two-dimensional numerical results for a model solar corona

A systematic study of dynamic response of the inner solar corona is made within the context of two-dimensional, time-dependent plane hydromagnetics. The governing equations are written in r-phi coordinates (i.e., in the solar equatorial plane), and numerical solutions are obtained by introducing an impulsive temperature enhancement within a rectangular region (i.e., 'box') in an initially isothermal corona in magnetohydrostatic equilibrium. Effects of the magnetic field configuration are illustrated for initially open (radial) and closed (azimuthal) magnetic fields by comparison with the nonmagnetic response. The channeling or blocking effects by the magnetic fields on the mass motion of solar plasma as the consequence of the evolution of fast and slow mode MHD shock waves are demonstrated. Some physically significant applications of the results, useful for the interpretation of observations, are discussed. It is noted, for example, that coronal transients observed in white light probably occur within essentially radial field topologies.

Wu, S. T.↗

Analysis and computation of three-dimensional flow in strongly curved ducts

Numerical solutions are presented for three-dimensional laminar and turbulent flow in curved ducts of rectangular cross section and significant curvature. The analysis is based on a primary-secondary velocity decomposition in a given coordinate system, and leads to approximate governing equations which correct an a priori inviscid solution for viscous effects, secondary flows, total pressure distortion, heat transfer, and internal flow blockage and losses. Solution of the correction equations is accomplished as an initial-value problem in space using an implicit forward-marching technique. The overall solution procedure requires significantly less computational effort than Navier-Stokes algorithms. The present solution procedure is effective even with the extreme local mesh resolution which is necessary to resolve near-wall sublayer regions in turbulent flow calculations. Computed solutions for both laminar and turbulent flow compare very favorably with available analytical and experimental results.

Kreskovsky, J. P.↗

Three-dimensional solutions for the thermal buckling and sensitivity derivatives of temperature-sensitive multilayered angle-ply plates

Analytic three-dimensional thermoelasticity solutions are presented for the thermal buckling of multilayered angle-ply composite plates with temperature-dependent thermoelastic properties. Both the critical temperatures and the sensitivity derivatives are computed. The sensitivity derivatives measure the sensitivity of the buckling response to variations in the different lamination and material parameters of the plate. The plates are assumed to have rectangular geometry and an antisymmetric lamination with respect to the middle plane. The temperature is assumed to be independent of the surface coordinates, but has an arbitrary symmetric variation through the thickness of the plate. The prebuckling deformations are accounted for. Numerical results are presented, for plates subjected to uniform temperature increase, showing the effects of temperature-dependent material properties on the prebuckling stresses, critical temperatures, and their sensitivity derivatives.

Noor, A. K.↗

Structural Dynamic/CFD Interaction for Computation of Vertical Tail Buffet

The multi-disciplinary problem of tail buffeting is solved using three sets of equations. The first set is the unsteady, compressible, full Navier Stokes equations which are used for obtaining the flow-field vector and the aerodynamic loads. The second set is the coupled aeroelastic equations which are used for obtaining the bending and torsional deflections of the tail. The third set is the grid-displacement equations which are used for updating the grid coordinates due to the tail deflections. For the computational applications a sharp-edged delta wing of aspect ratio one and a rectangular vertical tail of aspect ratio one placed in the plane of geometric symmetry behind the wing are considered. The configuration is pitched at a critical angle of attack (alpha = 38 deg.) which produces asymmetric, vortex-breakdown flow of the delta wing primary vortices. The results show the effects of coupled and uncoupled bending-torsional responses and the effects of the Reynolds number.

Kandil, Osama A.↗

Digital array scanned interferometers for astronomy

Investigations are reported of digital array scanned interferometers (DASI) with silicon CCD array detectors to define the operational capabilities of a mapping (polarimetric) spectrometer for astronomical applications based on these instruments. For spectral mapping, spatially resolved spectra using a cylindrical lens to image in the interferometers's redundant coordinate are given. The signal-to-noise characteristics of the Fourier transformed data are demonstrated with regard to the effects of a rectangular sampling function, spectral multiplexing and the pixel-to-pixel variation of the CCD array. These data indicate that DASIs can offer simple, versatile (polarization) mapping spectrometers suitable for spectral mapping observations from the ultraviolet to the infrared of extended sources at variable spatial resolution, particularly where long term stable operation is essential, as for spacecraft instruments.

Smith, Wm. Hayden↗

Illustration of airfoil shape effect on forward-swept wing divergence

A static aeroelastic analysis is presented of the divergence of untapered wings with conventional and supercritical airfoil sections at sweep angles of zero and -15 deg. One bending and one torsion mode were employed for a uniform rectangular cantilevered beam with the elastic axis at midchord, and calculations were based on a two-dimensional differential equations formulation in the structural coordinate system and in simple strip theory. A minimum divergence speed in the transonic range is obtained which is associated with the rearward shift of the aerodynamic center, and a 17% difference in minimum divergence dynamic pressure is found between a supercritical and a conventional wing. It is noted that although the strip method employed allows the assessment of the sensitivity of airfoil shapes to divergence, three-dimensional transonic aerodynamic methods should be used to predict wing divergence characteristics.

Bland, S. R.↗

Accelerated iterative calculation of transonic nacelle flowfields

A method is presented for the calculation of inviscid, supercritical flowfields about axisymmetric inlet cowls. A finite-difference calculation is performed in a simple, rectangular domain obtained from the nacelle geometry by a nearly-conformal mapping procedure. Type-dependent finite-differences are constructed using a coordinate-independent, 'rotated' differencing scheme. Methods of accelerating convergence of the iterative solution are demonstrated including a hybrid fast-Poisson-solver/relaxation scheme and an extrapolated relaxation procedure. Calculated pressure distributions are compared with experimental data for a variety of Mach numbers and mass-flow ratios, and show generally good agreement.

Caughey, D. A.↗

Vibrations of blades with variable thickness and curvature by shell theory

A procedure for analyzing the vibrations of rotating turbomachinery blades has been previously developed. This procedure is based upon shallow shell theory, and utilizes the Ritz method to determine frequencies and mode shapes. However, it has been limited heretofore to blades of uniform thickness, uniform curvature, and/or twist and rectangular planform. The present work shows how the procedure may be generalized to eliminate the aforementioned restrictions. Nonrectangular planforms are dealt with by a suitable coordinate transformation. This, as well as variable thickness, curvature and twist, require using numerical integration. The procedure is demonstrated on four examples of cantilevered blades for which theoretical and experimental data have been previously published: (1) flat plate with spanwise taper, (2) flat plate with chordwise taper, (3) twisted plate with chordwise taper, and (4) cylindrical shell with chordwise taper.

Lee, J. K.↗

Streamwise computation of three-dimensional parabolic flows

A new approach to calculate three-dimensional parabolic flows is presented. The flow field is computed by calculating velocity along a set of streamlines. The dependent variables commonly used in the computation of three-dimensional flows are the three velocity components. In contrast, the dependent variables in the present approach are the streamwise velocity and the two coordinates, in the cross-stream plane, of the chosen streamlines. The streamwise velocity is calculated from the finite difference equations obtained by applying Euler's momentum theorem to streamtubes constructed around the chosen streamlines; the streamline coordinates are calculated from the conservation of mass. Results of the calculations, based on the present approach, are compared with the experimental data for flow through rectangular ducts; the agreement is satisfactory.

Greywall, M. S.↗

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗

Geometric and structural properties of a rectangular supercritical wing oscillated in pitch for measurement of unsteady transonic pressure distributions

Wind-tunnel tests to measure unsteady aerodynamic data in the transonic region have been completed on an aspect ratio 2.0 rectangular wing with a supercritical airfoil. The geometric and structural properties of the wing are presented. (Other references contain the measured aerodynamic data.) Both measured and design airfoil coordinates are presented and compared. In addition, measured wing bending and torsional stiffness distributions and some trailing-edge flexibility influence coefficients are presented.

Ricketts, R. H.↗

Reference manual for the Thermal Analyst's Help Desk Expert System

This document provides technical information and programming guidance for the maintenance and future development of the Thermal Analyst's Help Desk. Help Desk is an expert system that operates within the EXSYSTM expert system shell, and is used to determine first approximations of thermal capacity for spacecraft and instruments. The five analyses supported in Help Desk are: (1) surface area required for a radiating surface, (2) equilibrium temperature of a surface, (3) enclosure temperature and heat loads for a defined position in orbit, (4) enclosure temperature and heat loads over a complete orbit and, (5) selection of appropriate surface properties. The two geometries supported by Help Desk are a single flat plate and a rectangular box enclosure. The technical information includes the mathematical approach and analytical derivations used in the analyses such as: radiation heat balance, view factor calculation, and orbit determination with coordinate transformation. The programming guide for developers describes techniques for enhancement of Help Desk. Examples are provided showing the addition of new features, user interface development and enhancement, and external program interfaces.

Ormsby, Rachel A.↗