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.
Explore the source record for details and available documents.
A feasibility study is presented concerning an application of a superconducting linear synchronous motor (LSM) to a large-scale rocket launcher, whose acceleration guide tube of LSM armature windings is constructed 1,500 meters under the ground. The rocket is released from the linear launcher just after it gets to a peak speed of about 900 kilometers per hour, and it flies out of the guide tube to obtain the speed of 700 kilometers per hour at the height of 100 meters above ground. The linear launcher is brought to a stop at the ground surface for a very short time of 5 seconds by a quick control of deceleration. Very large current variations in the single-layer windings of the LSM armature, which are produced at the higher speed region of 600 to 900 kilometers per hour, are controlled successfully by adopting the double-layer windings. The proposed control method makes the rocket launcher ascend stably in the superconducting LSM system, controlling the Coriolis force.
A method for generating finite-dimensional approximations to the solutions of optimal control problems is introduced. By employing a description of the dynamical system in terms of its attainable sets in favor of using differential equations, the controls are completely eliminated from the system model. Besides reducing the dimensionality of the discretized problem compared to state-of-the-art collocation methods, this approach also alleviates the search for initial guesses from where standard gradient search methods are able to converge. The mechanics of the new method are illustrated on a simple double integrator problem. The performance of the new algorithm is demonstrated on a 1-D rocket ascent problem (`Goddard Problem') in presence of a dynamic pressure constraint.
This report contains a theoretical discussion of the load distribution in bolted or riveted joints in light-alloy structures which is applicable not only for loads below the limit of proportionality but also for loads above this limit. The theory is developed for double and single shear joints. The methods given are illustrated by numerical examples and the values assumed for the bolt (or rivet) stiffnesses are based partly on theory and partly on known experimental values. It is shown that the load distribution does not vary greatly with the bolt (or rivet) stiffnesses and that for design purposes it is usually sufficient to know their order of magnitude. The theory may also be directly used for spot-welded structures and, with small modifications, for seam-welded structures, The computational work involved in the methods described is simple and may be completed in a reasonable time for most practical problems. A summary of earlier theoretical and experimental investigations on the subject is included in the report.
Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.
A three part paper to describe the results of several gear drive types with a high gear ratio is presented. A single stage planetary gear train with double helical gears is described with methods to reduce transmission errors and improve load distribution by regulating backlash during assembly. A new arrangement for face gear is also described. This new mechanism can perform rotations between axes that are collinear and intersected. Finally the design and simulation of an isostatic planetary gear train is presented. Conditions that can lead to noise and vibration of the planetary gear drive are described.
Recent analysis suggests that the effectiveness of stratospheric aerosol climate engineering through emission of non-condensable vapors such as SO2 is limited because the slow conversion to H2SO4 tends to produce aerosol particles that are too large; SO2 injection may be so inefficient that it is difficult to counteract the radiative forcing due to a CO2 doubling. Here we describe an alternate method in which aerosol is formed rapidly in the plume following injection of H2SO4, a condensable vapor, from an aircraft. This method gives better control of particle size and can produce larger radiative forcing with lower sulfur loadings than SO2 injection. Relative to SO2 injection, it may reduce some of the adverse effects of geoengineering such as radiative heating of the lower stratosphere. This method does not, however, alter the fact that such a geoengineered radiative forcing can, at best, only partially compensate for the climate changes produced by CO2.
Nephele is a small atmospheric probe mission concept to determine whether the Venus clouds contain organic matter. Nephele complements larger missions that seek to conduct Venus atmospheric analyses, such as DAVINCI, Venera-D, and the Morning Star missions, by directly sampling and analyzing cloud and haze particles. AERACEPT is an aerosol sampling technology designed to minimize the mass, volume, and complexity required for in situ planetary atmosphere characterization. A single aeroshell functions as a combined entry vehicle, descent probe, and instrumented sonde, combining recent advances in thermal protection materials with heritage technology from planetary and airborne science. The velocity of descent drives aerosol through inlets embedded in the aeroshell’s thermal protection material; inertial separation isolates the particles from the gas stream; and the particle capture surface doubles as the analysis substrate. This method can work without heat shield separation, descent control, or active sample handling; these advantages trade against a smaller sample volume, higher risk of sample heating, and shorter time for analysis and data return. AERACEPT further benefits from modeling of the internal and external flow to understand the particle capture efficiency and bias for size and concentration distributions. AERACEPT is well suited for a Venus mission such as Nephele, where a passive descent trajectory is both subsonic over the altitudes of interest and provides sufficient velocity and pressure differential for the particle sizes of interest. A thermal and flow toolchain developed in-house (Traj→ DPLR→ PATO→ openFOAM) was used to model an 80-cm aeroshell with science operations 63–39 km. These models indicate that AERACEPT will provide ~24 µL accumulated sample of particles ≥ 0.2 µm, substantially above the 2 µL threshold requirement. The worst-case particle heating will be ~14 K above ambient, for the smallest particles at the lowest altitude; most particles will experience ≤ 1 K rise, at which almost all analytes of interest will be unaffected. Plasma jet testing of a small-scale aeroshell-material inlet was recently conducted to validate the thermal model parameters; wind tunnel testing to validate the flow model parameters is underway.
HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.
HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.
A simple approach for the strain energy release rate computations based on the finite element method and a singular fitting model is presented. The model uses the stress and displacement distributions at the delamination front. The method is applied to a mixed-mode double cracked-lap-shear composite configuration. The strain energy release rate components predicted by the model are compared with the finite element crack-closure method. The effect of the mesh size on the stress and displacement distribution is isolated. The strain energy release rates predicted by relatively coarse mesh sizes are in good agreement with the finite element crack closure method.
Using a primitive variable formulation of the finite-element method, the differences in the surface observables, such as topography and geoid, produced by single- and double-layered thermal convection, were compared. Both constant and depth-dependent viscosities have been considered. For the same Rayleigh number, larger surface perturbations are produced by single-cell convection. For the same Nusselt number, the magnitudes of the surface observables are greater for double-layered convection. For the same surface heat-flux, surface topographies have similar magnitudes, but the relative amplitudes of geoid anomalies depend greatly on the style of viscosity stratification. This difference in the geoid between the two systems increases with greater surface heat-flow, regardless of viscosity structure.
Common methods for frequency stabilizing diode lasers systems employ gratings, etalons, optical electric double feedback, atomic resonance, and a Faraday cell with low magnetic field. Our method, the Faraday Anomalous Dispersion Optical Transmitter (FADOT) laser locking, is much simpler than other schemes. The FADOT uses commercial laser diodes with no antireflection coatings, an atomic Faraday cell with a single polarizer, and an output coupler to form a compound cavity. This method is vibration insensitive, thermal expansion effects are minimal, and the system has a frequency pull in range of 443.2 GHz (9A). Our technique is based on the Faraday anomalous dispersion optical filter. This method has potential applications in optical communication, remote sensing, and pumping laser excited optical filters. We present the first theoretical model for the FADOT and compare the calculations to our experimental results.
Because passive microwave instruments are confined to polar-orbiting satellites, rainfall estimates must interpolate across long time periods, during which no measurements are available. In this paper the authors discuss a technique that allows one to partially overcome the sampling limitations by using frequent infrared observations from geosynchronous platforms. To accomplish this, the technique compares all coincident microwave and infrared observations. From each coincident pair, the infrared temperature threshold is selected that corresponds to an area equal to the raining area observed in the microwave image. The mean conditional rainfall rate as determined from the microwave image is then assigned to pixels in the infrared image that are colder than the selected threshold. The calibration is also applied to a fixed threshold of 235 K for comparison with established infrared techniques. Once a calibration is determined, it is applied to all infrared images. Monthly accumulations for both methods are then obtained by summing rainfall from all available infrared images. Two examples are used to evaluate the performance of the technique. The first consists of a one-month period (February 1988) over Darwin, Australia, where good validation data are available from radar and rain gauges. For this case it was found that the technique approximately doubled the rain inferred by the microwave method alone and produced exceptional agreement with the validation data. The second example involved comparisons with atoll rain gauges in the western Pacific for June 1989. Results here are overshadowed by the fact that the hourly infrared estimates from established techniques, by themselves, produced very good correlations with the rain gauges. The calibration technique was not able to improve upon these results.
A method of measuring the Mode I fracture toughness of core/facesheet bonds in sandwich Structures is desired, particularly with the widespread use of models that need this data as input. This study examined if a critical strain energy release rate, G(sub IC), can be obtained from the climbing drum peel (CDP) test. The CDP test is relatively simple to perform and does not rely on measuring small crack lengths such as required by the double cantilever beam (DCB) test. Simple energy methods were used to calculate G(sub IC) from CDP test data on composite facesheets bonded to a honeycomb core. Facesheet thicknesses from 2 to 5 plies were tested to examine the upper and lower bounds on facesheet thickness requirements. Results from the study suggest that the CDP test, with certain provisions, can be used to find the GIG value of a core/facesheet bond.
Without any reference to the theory of differential equations, the initial value problem of the nonlinear, nonconservative double pendulum system is solved by the application of the method of Ritz to the equation of Hamilton. Also shown is an example of the reduction of the traditional eigenvalue problem of linear, homogeneous, differential equations of motion to the solution of a set of nonhomogeneous algebraic equations. No theory of differential equations is used. Solution of the time-space path of the linear oscillator is demonstrated and compared to the exact solution.
The 'coupled cluster single, double, and triple' (CCSDT) excitation model has been used to ascertain electron correlation energies for 14 different molecules representing a variety of chemical bonds, in conjunction with several methods of this type 'CCSDT-x', which include only an approximate treatment of connected triple excitations; these methods encompass CCSDT-1a, -1b, -2, -3, and -4, as well as the novel CCSD(T). While all methods treat the effects of connected triple excitations iteratively, CCSD(T) approaches then perturbationally. For the 14 molecules considered, the CCSD(T) method's average error relative to CCSDT is substantially lower than any of the CCSDT-x methods.
Mechanically induced viscoelastic dissipation is difficult to compute. When the constitutive model is defined by history integrals, the formula for dissipation is a double convolution integral. Since double convolution integrals are difficult to approximate, coupled thermo-mechanical analyses of highly viscous rubber-like materials cannot be made with most commercial finite element software. In this study, we present a method to approximate the dissipation for history integral constitutive models that represent Maxwell-like materials without approximating the double convolution integral. The method requires that the total stress can be separated into elastic and viscous components, and that the relaxation form of the constitutive law is defined with a Prony series. Numerical data is provided to demonstrate the limitations of this approximate method for determining dissipation. Rubber cylinders with imbedded steel disks and with an imbedded steel ball are dynamically loaded, and the nonuniform heating within the cylinders is computed.