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 concept of dissipative Hamiltonian system is converted to Lorentz-covariant form, with evolution generated jointly by two scalar functionals, the Lagrangian action and the global entropy. A bracket formulation yields the local covariant laws of energy-momentum conservation and of entropy production. The formalism is illustrated by a derivation of the covariant Landau kinetic equation.
Experimental data have been obtained in a three-dimensional, shear-driven, turbulent boundary layer subjected to sudden transverse strain. Measurements made with a newly developed, three-dimensional laser Doppler velocimeter include all components of mean flow, turbulent Reynolds stresses, and triple-product correlations. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) are extracted from the data. It was shown experimentally that eddy viscosity is strongly direction-dependent (not a scalar), indicating a need for full Reynolds stress modeling. Turbulence models for eddy viscosity, stress diffusion, pressure rate-of-strain, and anisotropy of Reynolds stresses are tested using the data. Calculations using a Reynolds stress equation model are compared with the data.
This paper describes a simple alternate approach to the difficult problem of modeling material behavior. Starting from a general representation for a rate-tpe constitutive equation, it is shown by example how sets of test data may be used to derive restrictions on the scalar functions appearing in the representation. It is not possible to determine these functions from experimental data, but the aforementioned restrictions serve as a guide in their eventual definition. The implications are examined for hypo-elastic, isotropically hardening plastic, and kinematically hardening plastic materials. A simple model for the evolution of the 'back-stress,' in a kinematic-hardening plasticity theory, that is entirely analogous to a hypoelastic stress-strain relation is postulated and examined in detail in modeling finitely plastic tension-torsion test. The implementation of rate-type material models in finite element algorithms is also discussed.
A time-accurate multiple-grid algorithm is described. The algorithm allows one to take much larger time steps with an explicit time-marching scheme than would otherwise be the case. Sample calculations of a scalar advection equation and the Euler equations for an oscillating airfoil are shown. For the oscillating airfoil, time steps an order of magnitude larger than the single-grid algorithm are possible.
The Harten (1983, 1984) total variation-diminishing (TVD) schemes, constituting a one-parameter explicit and implicit, second-order-accurate family, have the property of not generating spurious oscillations when applied to one-dimensional, nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. These methods are presently extended to the multidimensional hyperbolic conservation laws in curvilinear coordinates. Means by which to linearize the implicit operator and solution strategies, in order to improve the computation efficiency of the implicit algorithm, are discussed. Numerical experiments with steady state airfoil calculations indicate that the proposed linearized implicit TVD schemes are accurate and robust.
The impact of upstream-biased corrections for third-order spatial truncation error on the stability and phase error of the two-dimensional Crowley combined advective scheme with the cross-space term included is analyzed, putting primary emphasis on phase error reduction. The various versions of the Crowley scheme are formally defined, and their stability and phase error characteristics are intercompared using a linear Fourier component analysis patterned after Fromm (1968, 1969). The performances of the schemes under prototype simulation conditions are tested using time-dependent numerical experiments which advect an initially cone-shaped passive scalar distribution in each of three steady nondivergent flows. One such flow is solid rotation, while the other two are diagonal uniform flow and a strongly deformational vortex.
Data from the MAGSAT mission were used to derive a global scalar magnetic anomaly map at an average altitude of about 400 km. It was possible to work with 2 data sets corresponding to dawn and dusk. The anomalies which were repeatable at dawn and at dusk was identified and the error limits of these anomalies were estimated. The repeatable anomalies were downward continued to about 10 km altitude. The anomalies over Canada were correlated quantitatively with bandpass filtered magnetic anomalies derived from aeromagnetic surveys. The close correlation indicates that the repeatable anomalies detected from orbit are due to geological causes. This correlation supports the geological significance of the global anomaly map.
Crustal anomaly detection with MAGSAT data is frustrated by inherent resolving power of the data and by contamination from external and core fields. Quality of the data might be tested by modeling specific tectonic features which produce anomalies that fall within proposed resolution and crustal amplitude capabilities of MAGSAT fields. To test this hypothesis, north African hotspots associated with Ahaggar, Tibesti and Darfur were modeled as magnetic induction anomalies. MAGSAT data were reduced by subtracting external and core fields to isolate scalar and vertical component crustal signals. Of the three volcanic areas, only the Ahaggar region had an associated anomaly of magnitude above error limits of the data. Hotspot hypothesis was tested for Ahaggar by seeing if predicted magnetic signal matched MAGSAT anomaly. Predicted model magnetic signal arising from surface topography of the uplift and the Curie isothermal surface was calculated at MAGSAT altitudes by Fourier transform technique modified to allow for variable magnetization. Curie isotherm surface was calculated using a method for temperature distribution in a moving plate above a fixed hotspot. Magnetic signal was calculated for a fixed plate as well as a number of plate velocities and directions.
A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations.
Crustal anomaly detection with MAGSAT data is frustrated by the inherent resolving power of the data and by contamination from the external and core fields. The quality of the data might be tested by modeling specific tectonic features which produce anomalies that fall within the proposed resolution and crustal amplitude capabilities of the MAGSAT fields. To test this hypothesis, the north African hotspots associated with Ahaggar, Tibestia and Darfur have been modeled as magnetic induction anomalies due solely to shallower depth to the Curie isotherm surface beneath these features. The MAGSAT data were reduced by subtracting the external and core fields to isolate the scalar and vertical component crustal signals. The predicted model magnetic signal arising from the surface topography of the uplift and the Curie isotherm surface was calculated at MAGSAT altitudes by the Fourier transform technique modified to allow for variable magnetization. In summary it is suggested that the region beneath Ahaggar is associated with a strong thermal anomaly and the predicted anomaly best fits the associated MAGSAT anomaly if the African plate is moving in a northeasterly direction.
An algorithm for vectorized computation of stiffness matrices of an 8 noded isoparametric hexahedron element for geometric nonlinear analysis was developed. This was used in conjunction with the earlier 2-D program GAMNAS to develop the new program NAS3D for geometric nonlinear analysis. A conventional, modified Newton-Raphson process is used for the nonlinear analysis. New schemes for the computation of stiffness and strain energy release rates is presented. The organization the program is explained and some results on four sample problems are given. The study of CPU times showed that savings by a factor of 11 to 13 were achieved when vectorized computation was used for the stiffness instead of the conventional scalar one. Finally, the scheme of inputting data is explained.
The design of scalar and multi-variable feedback control systems for the GET700 turboshaft engine coupled to a helicopter rotor system is examined. A series of linearized models are presented and analyzed. Robustness and performance specifications are posed in the frequency domain. The linear-quadratic-Gaussian with loop-transfer-recovery (LQG/LTR) methodology is used to obtain a sequence of three feedback designs. Even in the single-input/single-output case, comparison of the current control system with that derived from the LQG/LTR approach shows significant performance improvement. The multi-variable designs, evaluated using linear and nonlinear simulations, show even more potential for performance improvement.
Using laser-induced fluorescence to probe nitrogen flows seeded with small amounts of nitric oxide, simultaneous measurements of all three thermodynamic scalar quantities temperature, density, and pressure, were demonstrated in a supersonic turbulent boundary layer. Instrumental uncertainty is 1% for temperature and 2% for density and pressure, making the techniques suitable for measurements of turbulent fluctuations. This technology is currently being transferred to an experimental program designed to use these optical techniques in conjunction with traditional methods to make measurements in turbulent flowfields that were not possible before. A detailed descritpion of the research progress and pertinent results are presented.
The Galerkin finite element solutions for the scalar homogeneous Helmholtz equation are presented for no reflection, hard wall, and potential relief exit terminations with a variety of triangular element orientations. For this group of problems, the correlation between the accuracy of the solution and the orientation of the linear triangle is examined. Nonsymmetric element patterns are found to give generally poor results in the model problems investigated, particularly for cases where standing waves exist. For a fixed number of vertical elements, the results showed that symmetric element patterns give much better agreement with corresponding exact analytical results. In laminated wave guide application, the symmetric pyramid pattern is convenient to use and is shown to give excellent results.
The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.
The response of beam and plate members to pulse and impact loading is investigated, and numerical methods are analyzed and compared. The equations of motion, the finite difference and finite element methods, the time integration/Runge-Kutta techniques, and material modeling are discussed in detail. It is found that both the finite difference and the finite element methods could be accurately employed to discretize the spatial variation in the displacements. Central differences or a fourth order Runge-Kutta algorithm could be used for the time integration. The total energy of the system would give the stability and accuracy of the solution. Results showed that the finite element method provided better efficiency in obtaining accurate solution than the finite difference method when a scalar processor is used. However, the finite difference method was more efficient on the vector processor. Therefore, the most efficient method of solution depends on the type of computer present for the analysis.
The Ahaggar volcanic province of North Central Africa is considered a region of excess heat flow (hot spot) and hence elevated Curie isotherm. Using a modified version of the Parker FFT potential field representation, magnetic signals were calculated at Magsat altitudes for models in which the African Plate is both fixed and moving. The moving-plate model extends the Curie isotherm anomaly in the direction of plate motion and provides a satisfactory match to vertical component anomaly data when the magnitude of plate velocity is 0.75 cm/yr. Although the signal levels are marginal for the scalar component anomalies of this region, the same model provides an adequate match to this data set and is clearly preferable to a fixed-plate model.
The errors in ship wind reports of light wings tend to significantly bias their mean. This occurs because wind speed is a scalar quantity that is constrained to zero or positive values. Therefore, observations tend to overestimate the light winds because of the one-sided distribution of errors, but the bias disappears under stronger winds. A method for removing this bias from ship data is presented. In particular, the method is applied to interpreting the ratios of wind speeds observed by ships to those obtained from tracking low level clouds. Corrected ratios allow low cloud speeds to serve as proxy data for surface based observations.