SEARCH · Engineering Papers
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.
The BAT (Batten-Augmented Triangular) BEAM is characterized by battens which are buckled in the deployed state, thus preloading the truss. The preload distribution is determined, and the effects of various external loading conditions are investigated. The conceptual design of a deployer is described and loads are predicted. The influence of joint imperfections on effective member stiffness is investigated. The beam is assessed structurally.
A niobium alloy with 1 percent zirconium and 0.063 percent carbon by weight was investigated in the as-rolled and annealed conditions, and after high-temperature (1350 and 1400 K) exposure with and without an applied stress. In the as-rolled and annealed conditions, large metastable carbides were observed in addition to a regular distribution of small particles. During the high-temperature exposure, the majority of the large carbides were dissolved and a more stable carbide phase formed. This finely dispersed phase had a composition determined to be approximately 70 percent ZrC and approximately 30 percent NbC and showed some evidence of an orientation relationship with the matrix. The precipitates appeared to coarsen slightly after approximately 5000 hr exposure in the presence of an applied stress resulted in a decrease in the size and in the interparticle spacing of the stable precipitates. However, the composition of the precipitate phase and its ability to pin dislocations were not affected by the temperature or stress conditions.
The Variable Sweep Transition Flight Experiment (VSTFE) was initiated to establish a boundary-layer transition data base for laminar flow wing design. For this experiment, full-span upper-surface gloves will be fitted to a variable sweep F-14 aircraft. The results of two initial tasks are documented: a parametric pressure distribution/boundary-layer stability study and the design of an upper-surface glove for Mach 0.8. The first task was conducted to provide a data base from which wing-glove pressure distributions could be selected for glove designs. Boundary-layer stability analyses were conducted on a set of pressure distributions for various wing sweep angles, Mach numbers, and Reynolds number in the range of those anticipated for the flight-test program. The design procedure for the Mach 0.8 glove is described, and boundary-layer stability calculations and pressure distributions are presented both at design and off-design conditions. Also included is the analysis of the clean-up glove (smoothed basic wing) that will be flight-tested initially and the analysis of a Mach 0.7 glove designed at the NASA Langley Research Center.
A method is described for the analysis of the unsteady, incompressible potential flow associated with a helicopter rotor and it's wake in forward flight. This method is particularly useful in low advance ratio flight due to the major contribution, in the near field, of the deformed wake. The rotor geometry is prescribed and the unsteady wake geometry is computed from the local flow perturbation velocities. The wake is modeled as a full vortex lattice. The rotor geometry is arbitrary and several rotor blades can be represented. The unsteady airloads on the rotor blades are computed in the presence of the deformed rotor wake by a time-stepping technique. Solution for the load distribution on the blade surfaces is found by prescribing boundary conditions in a reference system which rotates with the blade tips. Transformation tensors are used to describe the contribution of the wake in the inertial system to the rotor in the rotating reference system. The effects of blade cyclic pitch variation are computed using a rotation tensor. The deformation of the wake is computed in the inertial frame. The wake is started impulsively from rest, allowing a natural convection of the wake with time.
The Boundary Force Method (BFM) was formulated for the two-dimensional stress analysis of complex crack configurations. In this method, only the boundaries of the region of interest are modeled. The boundaries are divided into a finite number of straight-line segments, and at the center of each segment, concentrated forces and a moment are applied. This set of unknown forces and moments is calculated to satisfy the prescribed boundary conditions of the problem. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate are used as the fundamental solution. Thus, the crack need not be modeled as part of the boundary. The formulation of the BFM is described and the accuracy of the method is established by analyzing several crack configurations for which accepted stress-intensity factor solutions are known. The crack configurations investigated include mode I and mixed mode (mode I and II) problems. The results obtained are, in general, within + or - 0.5 percent of accurate numerical solutions. The versatility of the method is demonstrated through the analysis of complex crack configurations for which limited or no solutions are known.
The reflection from and transmission through a semiinfinite chiral medium are analyzed by obtaining the Fresnel equations in terms of parallel- and perpendicular-polarized modes, and a comparison is made with results reported previously. The chiral medium is described electromagnetically by the constitutive relations D = (epsilon)E+i(gamma)B and H = i(gamma)E+(1/mu)B. The constants epsilon, mu and gamma are real and have values that are fixed by the size, the shape, and the spatial distribution of the elements that collectively compose the medium. The conditions are obtained for the total internal reflection of the incident wave from the interface and for the existence of the Brewster angle. The effects of the chirality on the polarization and the intensity of the reflected wave from the chiral half-space are discussed and illustrated by using the Stokes parameters. The propagation of electromagnetic wave through an infinite slab of chiral medium is formulated for oblique incidence and solved analytically for the case of normal incidence.
Seasat A scatterometer measurements of surface wind components made under normal weather conditions during August and September 1978 are examined. The longitudinal distributions of the monthly mean zonal component for each season are described. The average monthly mean slope of the wave number spectra throughout the 550- to 2200-km wavelength band was -1.7. The spectra levels of the zonal wind, but not the meridional component, are substantially different in each equatorial ocean.
The current sheet and plasma sheet boundary layer (PSBL) in a magnetotail model are discussed. A test particle code is used to study the response of ensembles of particles to a two-dimensional, time-dependent model of the geomagnetic tail, and test the proposition (Coroniti, 1985a, b; Buchner and Zelenyi, 1986; Chen and Palmadesso, 1986; Martin, 1986) that the stochasticity of the particle orbits in these fields is an important part of the physical mechanism for magnetospheric substorms. The realistic results obtained for the fluid moments of the particle distribution with this simple model, and their insensitivity to initial conditions, is consistent with this hypothesis.
An analytical study of linear-amplifying instabilities of a laminar boundary layer as found in the experimental data of the LaRC/8-foot laminar-flow control (LFC) experiment was completed and the results are presented. The LFC airfoil used for this experiment was a swept, supercritical design which removed suction air through spanwise slots. The amplification of small disturbances by linear processes on a swept surface such as this can be due to either Tollmien-Schlichting (TS) and/or crossflow (CF) mechanisms. This study consists of the examination of these two instabilities by both the commonly used incompressible (SALLY and MARIA) analysis and the more involved compressible (COSAL) analysis. A wide range of experimental test conditions with variations in Mach number, Reynolds number, and suction distributions were available for this study. Experimentally determined transition locations were found from thin-film techniques and were used to correlate the n-factors at transition for the range of test cases.
A wind-tunnel interference assessment method applicable to test sections with discrete finite-length wall slots is described. The method is based on high order panel method technology and uses mixed boundary conditions to satisfy both the tunnel geometry and wall pressure distributions measured in the slotted-wall region. Both the test model and its sting support system are represented by distributed singularities. The method yields interference corrections to the model test data as well as surveys through the interference field at arbitrary locations. These results include the equivalent of tunnel Mach calibration, longitudinal pressure gradient, tunnel flow angularity, wall interference, and an inviscid form of sting interference. Alternative results which omit the direct contribution of the sting are also produced. The method was applied to the National Transonic Facility at NASA Langley Research Center for both tunnel calibration tests and tests of two models of subsonic transport configurations.
A numerical solution has been developed for the anisotropic multiple scattering of light in a spherical shell comet atmosphere. The code has been run for a spherically symmetric coma distribution, benchmarked against past studies, and then run for the conditions of an axisymmetric dust jet at the subsolar point of the comet. The radiant flux impinging on the nucleus surface and the mean intensity of light throughout the coma were investigated.
We have completed computer simulations that model shock effects in the venusian atmosphere caused during the passage of two cometlike bodies 100 m and 1000 m in diameter and an asteroidlike body 10 km in diameter. Our objective is to examine hypervelocity-generated shock effects in the venusian atmosphere for bodies of different types and sizes in order to understand the following: (1) their deceleration and depth of penetration through the atmosphere; and (2) the onset of possible ground-surface shock effects such as splotches, craters, and ejecta formations. The three bodies were chosen to include both a range of general conditions applicable to Venus as well as three specific cases of current interest. These calculations use a new multiphase computer code (DICE-MAZ) designed by California Research & Technology for shock-dynamics simulations in complex environments. The code was tested and calibrated in large-scale explosion, cratering, and ejecta research. It treats a wide range of different multiphase conditions, including material types (vapor, melt, solid), particle-size distributions, and shock-induced dynamic changes in velocities, pressures, temperatures (internal energies), densities, and other related parameters, all of which were recorded in our calculations.
A probabilistic global/local method is proposed to reduce the computational requirements of probabilistic structural analysis. A coarser global model is used for most of the computations with a local more refined model used only at key probabilistic conditions. The global model is used to establish the cumulative distribution function (cdf) and the Most Probable Point (MPP). The local model then uses the predicted MPP to adjust the cdf value. The global/local method is used within the advanced mean value probabilistic algorithm. The local model can be more refined with respect to the g1obal model in terms of finer mesh, smaller time step, tighter tolerances, etc. and can be used with linear or nonlinear models. The basis for this approach is described in terms of the correlation between the global and local models which can be estimated from the global and local MPPs. A numerical example is presented using the NESSUS probabilistic structural analysis program with the finite element method used for the structural modeling. The results clearly indicate a significant computer savings with minimal loss in accuracy.
Tabular and plotted results are presented (as a function of lamination angle and temperature) for the first two vibration frequencies and mode shapes of a heated, single-layer Boron/epoxy laminate. A Ritz formulation is used to determine the thermally induced in-plane forces and the out-of-plane vibration frequencies. This approach is formulated in terms of the out-of-plane displacement and the in-plane stress function; both the displacement and the stress function are expressed as simple algebraic polynomials modified to satisfy the essential boundary conditions. Heating is applied as a steady-state temperature distribution; material properties are considered to be temperature-dependent. To establish the correctness and effectiveness of the method, the in-plane forces are compared with Fourier series analytical results.
A two-dimensional finite element model was used to determine the freeze/thaw characteristics of an external artery heat pipe. During startup, the working fluid, which was located in the liquid channel and the circumferential wall grooves, experienced a phase transformation from a solid to a liquid state. The transient heat conduction equations with moving interfacial conditions were solved using the appropriate initial boundary conditions. The modelling results include the cross-sectional temperature distribution and the interfacial or melt front position as a function of time. A fixed grid approach was adopted in the model for the phase-change process during thawing of frozen working fluid. The interfacial position between the liquid and solid regions was found by balancing the latent heat caused by interfacial movement with the heat addition or extraction at the related grid points.
The performance of the turbine component of an NACA research jet engine was investigated with cold air. The interaction and the matching of the turbine with the NACA eight-stage compressor were computed with the combination considered as a jet engine. The over-all performance of the engine was then determined. The internal aerodynamics were studied to the extent of investigating the performance of the first stator ring and its influence on the turbine performance. For this ring, the stream-filament method for computing velocity distribution permitted efficient sections to be designed, but the design condition of free-vortex flow with uniform axial velocities was not obtained.
The external wave drag of bodies of revolution moving at supersonic speeds can be expressed either in terms of the geometry of the body, or in terms of the body-simulating axial source distribution. For purposes of deriving optimum bodies under various given conditions, it is found that the second of the methods mentioned is the more tractable. By use of a quasi-cylindrical theory, that is, the boundary conditions are applied on the surface of a cylinder rather than on the body itself, the variational problems of the optimum bodies having prescribed volume or caliber are solved. The streamline variations of cross-sectional area and drags of the bodies are exhibited, and some numerical results are given.
An advanced higher harmonic control (HHC) analysis has been developed and applied to investigate its effect on vibration reduction levels, blade and control system fatigue loads, rotor performance, and power requirements of servo-actuators. The analysis is based on a finite element method in space and time. A nonlinear time domain unsteady aerodynamic model, based on the indicial response formulation, is used to calculate the airloads. The rotor induced inflow is computed using a free wake model. The vehicle trim controls and blade steady responses are solved as one coupled solution using a modified Newton method. A linear frequency-domain quasi-steady transfer matrix is used to relate the harmonics of the vibratory hub loads to the harmonics of the HHC inputs. Optimal HHC is calculated from the minimization of the vibratory hub loads expressed in term of a quadratic performance index. Predicted vibratory hub shears are correlated with wind tunnel data. The fixed-gain HHC controller suppresses completely the vibratory hub shears for most of steady or quasi-steady flight conditions. HHC actuator amplitudes and power increase significantly at high forward speeds (above 100 knots). Due to the applied HHC, the blade torsional stresses and control loads are increased substantially. For flight conditions where the blades are stalled considerably, the HHC input-output model is quite nonlinear. For such cases, the adaptive-gain controller is effective in suppressing vibratory hub loads, even though HHC may actually increase stall areas on the rotor disk. The fixed-gain controller performs poorly for such flight conditions. Comparison study of different rotor systems indicates that a soft-inplane hingeless rotor requires less actuator power at high speeds (above 130 knots) than an articulated rotor, and a stiff-inplane hingeless rotor generally requires more actuator power than an articulated or a soft-inplane hingeless rotor. Parametric studies for a hingeless rotor operating in a transition flight regime and for an articulated rotor operating at the level-flight boundary (high speed and high thrust conditions) indicate that blade parameters including flap, lag, torsion stiffness distributions, linear pretwist, chordwise offset of center-of-mass from elastic axis and chordwise offset of elastic axis from aerodynamic center can be selected to minimize the actuator power requirements for HHC.