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 initial effort was concentrated on developing the quasi-analytical approach for two-dimensional transonic flow. To keep the problem computationally efficient and straightforward, only the two-dimensional flow was considered and the problem was modeled using the transonic small perturbation equation.
The authors treat the robust stability issue using the characteristic polynomial, for two different cases: first in coefficient space with respect to perturbations in the coefficient of the characteristic polynomial; and then for a control system containing perturbed parameters in the transfer function description of the plant. In coefficient space, a simple expression is first given for the l-(squared) stability margin for both the monic and nonmonic cases. Following this, a method is given to find the l(infinity) margin, and the method is extended to reveal much larger stability regions. In parameter space the authors consider all single-input (multi-output) or single-output (multi-input) systems with a fixed controller and a plant described by a set of transfer functions which are ratios of polynomials with variable coefficients. A procedure is presented to calculate the radius of the largest stability ball in the space of these variable parameters. The calculation serves as a stability margin for the control system. The formulas that result are quasi-closed-form expressions for the stability margin and are computationally efficient.
A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over the airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small perturbation about a steady inviscid flow over a NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 degrees. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time marching scheme are also presented. This procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of terms for artificial viscosity.
NASA's Mar Exploration Rover (MER) robotic prospectors were produced in an environment of unusually challenging schedule, volume, and mass restrictions. The technical challenges pushed the system s design towards extensive integration of function, which resulted in complex system engineering issues. One example of the system's integrated complexity can be found in the deployment process for the rover. Part of this process, rover "standup", is outlined in this paper. Particular attention is given to the Rover Lift Mechanism's (RLM) role and its design. Analysis methods are presented and compared to test results. It is shown that because prudent design principles were followed, a robust mechanism was created that minimized the duration of integration and test, and enabled recovery without perturbing related systems when reasonably foreseeable problems did occur. Examples of avoidable, unnecessary difficulty are also presented.
A Matlab-based toolbox has been developed for the wavefront control and optimization of segmented optical surfaces to correct for possible misalignments of James Webb Space Telescope (JWST) using influence functions. The toolbox employs both iterative and non-iterative methods to converge to an optimal solution by minimizing the cost function. The toolbox could be used in either of constrained and unconstrained optimizations. The control process involves 1 to 7 degrees-of-freedom perturbations per segment of primary mirror in addition to the 5 degrees of freedom of secondary mirror. The toolbox consists of a series of Matlab/Simulink functions and modules, developed based on a "wrapper" approach, that handles the interface and data flow between existing commercial optical modeling software packages such as Zemax and Code V. The limitations of the algorithm are dictated by the constraints of the moving parts in the mirrors.
A symmetrized version of the recently developed refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] is proposed. This model takes into account line coupling effects and hence allows the calculation of the off-diagonal elements of the relaxation matrix, without neglecting the rotational structure of the perturbing molecule. The formalism is applied to the isotropic Raman spectra of autoperturbed N2 for which a benchmark quantum relaxation matrix has recently been proposed. The consequences of the classical path approximation are carefully analyzed. Methods correcting for effects of inelasticity are considered. While in the right direction, these corrections appear to be too crude to provide off diagonal elements which would yield, via the sum rule, diagonal elements in good agreement with the quantum results. In order to overcome this difficulty, a re-normalization procedure is applied, which ensures that the off-diagonal elements do lead to the exact quantum diagonal elements. The agreement between the (re-normalized) semi-classical and quantum relaxation matrices is excellent, at least for the Raman spectra of N2, opening the way to the analysis of more complex molecular systems.
NASA’s Mar Exploration Rover (MER) robotic prospectors were produced in an environment of unusually challenging schedule, volume, and mass restrictions. The technical challenges pushed the system’s design towards extensive integration of function, which resulted in complex system engineering issues. One example of the system’s integrated complexity can be found in the deployment process for the rover. Part of this process, rover “standup”, is outlined in this paper. Particular attention is given to the Rover Lift Mechanism’s (RLM) role and its design. Analysis methods are presented and compared to test results. It is shown that because prudent design principles were followed, a robust mechanism was created that minimized the duration of integration and test, and enabled recovery without perturbing related systems when reasonably foreseeable problems did occur. Examples of avoidable, unnecessary difficulty are also presented.
A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.
Water with even a small amount of dissolved solids has an electrical conductivity orders of magnitude higher than dry rock and is therefore a near-ideal exploration target on Mars for low frequency, diffusive electromagnetic methods. Models of the temperature- and frequency-dependent electrical properties of rock-ice-water mixtures are used to predict the electromagnetic response of the Martian subsurface. Detection of ice is difficult unless it is massively segregated. In contrast, liquid water profoundly affects soundings, and even a small amount of adsorbed water in the cryosphere can be detected. Subcryospheric water is readily distinguishable at frequencies as low as 100 Hz for fresh water to 10 mHz for brines. These responses can be measured using either natural or artificial sources. Ultra low frequency signals from solar wind and diurnal-heating perturbations of the ionosphere are likely, and disturbances of regional crustal magnetic fields may also be observable. Spherics, or extremely to very low frequency signals from lightning discharge, would provide optimal soundings; however, lightning may be the least likely of the possible natural sources. Among the active techniques, only the time-domain electromagnetic (TDEM) method can accommodate a closely spaced transmitter and receiver and sound to depths of hundreds of meters or more. A ground- or aircraft-based TDEM system of several kilograms can detect water to a depth of several hundred meters, and a system of tens of kilograms featuring a large, fixed, rover- or ballistically deployed loop can detect water to several kilometers depth.
A state-space approach to parameter estimation in linear functional differential equations is developed using the theory of linear evolution equations. A locally convergent quasi-Newton type algorithm is applied to distributed systems with particular emphasis on parameters that induce unbounded perturbations of the state. The algorithm is computationally implemented on several functional differential equations, including coefficient and delay estimation in linear delay-differential equations.
A method is presented in satellite altimetry that attempts to simultaneously determine the geoid and sea surface topography with minimum wavelengths of about 500 km and to reduce the radial orbit error caused by geopotential errors. The modeling of the radial orbit error is made using the linearized Lagrangian perturbation theory. Secular and second order effects are also included. After a rather extensive validation of the linearized equations, alternative expressions of the radial orbit error are derived. Numerical estimates for the radial orbit error and geoid undulation error are computed using the differences of two geopotential models as potential coefficient errors, for a SEASAT orbit. To provide statistical estimates of the radial distances and the geoid, a covariance propagation is made based on the full geopotential covariance. Accuracy estimates for the SEASAT orbits are given which agree quite well with already published results. Observation equations are develped using sea surface heights and crossover discrepancies as observables. A minimum variance solution with prior information provides estimates of parameters representing the sea surface topography and corrections to the gravity field that is used for the orbit generation. The simulation results show that the method can be used to effectively reduce the radial orbit error and recover the sea surface topography.
The one-electron bond-orbital model of Harrison (1973) for tetrahedrally coordinated solids is modified to a two-electron model using a generalization of Falicov and Harris' (1969) method for treating the hydrogen molecule. Expressions for nuclear-exchange and pseudodipolar interactions are obtained in terms of bond-orbital model parameters. The eigenfunctions and eigenvalues of the Hamiltonian are employed along with standard perturbation theory to calculate the energy shifts resulting from the application of an external electric field and induced by the magnetic interaction between the nuclear spins and the bond-electron spins. The quantitative predictions of the derived formulas are examined for real semiconductors involving elements and compounds from periodic groups III, IV, and V. It is concluded that the present two-electron model makes it possible to assess the quantitative effects of electron correlation in the anion-cation bond for both nonmagnetic properties (electric susceptibility, dielectric constant) and magnetic properties (nuclear-exchange and pseudodipolar interactions) of various materials.
An instrument capable of observing the natural electron flux in the energy range from 0.1 to 12.0 kiloelectron volts is discussed for use in an experiment intended as a forerunner of a method that will utilize artificially accelerated electrons as tracer particles for electron fields parallel to the magnetic field. Effects that are of importance either as means of detecting the echo beam or as causes of beam perturbations (e.g., spacecraft charging effects and electron background) are to be studied. The use of electron accelerators as a tool to probe magnetospheric processes rather than to modify them is planned.
A method is presented that permits the thermal-acoustic efficiency spectrum in a long turbulent burner to be recovered from the corresponding far-field sound spectrum. An acoustic source/propagation model is used based on the perturbation solution of the equations describing the unsteady one-dimensional flow of an inviscid ideal gas with a distributed heat source. The technique is applied to a long cylindrical hydrogen-flame burner operating over power levels of 4.5-22.3 kW. The results show that the thermal-acoustic efficiency at a given frequency, defined as the fraction of the total burner power converted to acoustic energy at that frequency, is rather insensitive to burner power, having a maximum value on the order of 10 to the -4th at 150 Hz and rolling off steeply with increasing frequency. Evidence is presented that acoustic agitation of the flame at low frequencies enhances the mixing of the unburned fuel and air with the hot products of combustion. The paper establishes the potential of the technique as a useful tool for characterizing the acoustic source structure in any burner, such as a gas turbine combustor, for which a reasonable acoustic propagation model can be postulated.
Active and passive microwave sensors were applied effectively to the problem of estimating the surface soil moisture in a variety of environmental conditions. Research to date has shown that both types of sensors are also sensitive to the surface roughness and the vegetation cover. In estimating the soil moisture, the effect of the vegetation and roughness are often corrected either by acquiring multi-configuration (frequency and polarization) data or by adjusting the surface parameters in order to match the model predictions to the measured data. Due to the limitations on multi-configuration spaceborne data and the lack of a priori knowledge of the surface characteristics for parameter adjustments, it was suggested that the synergistic use of the sensors may improve the estimation of the soil moisture over the extreme range of naturally occurring soil and vegetation conditions. To investigate this problem, the backscattering and emission from a bare soil surface using the classical rough surface scattering theory were modeled. The model combines the small perturbation and the Kirchhoff approximations in conjunction with the Peak formulation to cover a wide range of surface roughness parameters with respect to frequency for both active and passive measurements. In this approach, the same analytical method was used to calculate the backscattering and emissivity. Therefore, the active and passive simulations can be combined at various polarizations and frequencies in order to estimate the soil moisture more actively. As a result, it is shown that (1) the emissivity is less dependent on the surface correlation length, (2) the ratio of the backscattering coefficient (HH) over the surface reflectivity (H) is almost independent of the soil moisture for a wide range of surface roughness, and (3) this ratio can be approximated as a linear function of the surface rms height. The results were compared with the data obtained by a multi-frequency radiometer-scatterometer system working at frequencies between 3.0 GHz to 11.0 GHz. The data were acquired over bare soil surfaces with moisture variations due to freezing and thawing and roughness changes due to rain and erosion. The model predictions are shown to be in reasonable agreement with the data. In addition, it was shown that the same ratio when calculated from the data shows almost no dependence on the soil moisture. Finally, a simple technique which combines the backscattering coefficient at HH polarization (active sensing) with the emissivity at H polarization (passive sensing) is suggested for retrieving the soil moisture from bare soil surfaces.
Two autopilots providing autonomous flight capabilities are presented herein. The first is the Pico-Pilot, demonstrated for the 12-inch size class of micro air vehicles. The second is the MicroPilot MP2028(sup g), where its integration into a 36-inch Zagi airframe (tailless, elevons only configuration) is investigated and is the main focus of the report. Analytical methods, which include the use of the Advanced Aircraft Analysis software from DARCorp, were used to determine the stability and control derivatives, which were then validated through wind tunnel experiments. From the aerodynamic data, the linear, perturbed equations of motion from steady-state flight conditions may be cast in terms of these derivatives. Using these linear equations, transfer functions for the control and navigation systems were developed and feedback control laws based on Proportional, Integral, and Derivative (PID) control design were developed to control the aircraft. The PID gains may then be programmed into the autopilot software and uploaded to the microprocessor of the autopilot. The Pico-Pilot system was flight tested and shown to be successful in navigating a 12-inch MAV through a course defined by a number of waypoints with a high degree of accuracy, and in 20 mph winds. The system, though, showed problems with control authority in the roll and pitch motion of the aircraft: causing oscillations in these directions, but the aircraft maintained its heading while following the prescribed course. Flight tests were performed in remote control mode to evaluate handling, adjust trim, and test data logging for the Zagi with integrated MP2028(sup g). Ground testing was performed to test GPS acquisition, data logging, and control response in autonomous mode. Technical difficulties and integration limitations with the autopilot prevented fully autonomous flight from taking place, but the integration methodologies developed for this autopilot are, in general, applicable for unmanned air vehicles within the 36-inch size class or larger that use a PID control based autopilot.
A microwave oven and microwave heating method generates microwaves within a cavity in a predetermined mode such that there is a known region of uniform microwave field. Samples placed in the region will then be heated in a relatively identical manner. Where perturbations induced by the samples are significant, samples are arranged in a symmetrical distribution so that the cumulative perturbation at each sample location is the same.
The Morbiter software numerically averages an osculating orbit s equations of motion (EOM) to arrive at the mean orbit s EOMs, which are then numerically propagated to obtain the long-term orbital ephemerides. The long-term evolution characteristics, and stability, of an orbit are best characterized using a mean element propagation of the perturbed, two-body variational equations of motion. The average process eliminates short period terms, leaving only secular and long period effects. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods.