SEARCH · Search NASA
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.
In this paper we present a solution to the formation initialization (FI) problem for N distributed spacecraft located in deep space. Our solution to the FI problem is based on a three-stage sky search procedure that reduces the FI problem for N spacecraft to the simpler problem of initializing a set of sub-formations.
As NASA works to extend human exploration of the solar system, unprecedented capacity for electrical and thermal energy will be needed to support sustained human presence and the beginning of industrial activity on the Moon. Solar energy is abundant on the lunar surface but extended night hours (up to 350 consecutive hours) and extreme temperature change from daylight to nighttime operation creates complexity for solar power use. Analogous issues arise on Earth, where demand for additional renewable energy generation, including solar, is rising, but better power management, distribution, and energy storage solutions are needed to address intermittency and resiliency, among other issues. To accelerate the development of technologies for power distribution and energy storage that can operate in extreme lunar conditions, NASA’s Centennial Challenges Program launched the Watts on the Moon Challenge in September 2020. Submissions for Phase 1 of the challenge were due in March 2021 and winners were announced in May 2021. This paper provides an overview of the challenge, key technical requirements, outcomes from Phase 1, and an update on Phase 2 of the challenge, which is currently ongoing.
Numerical methods incorporated into a computer program to provide estimates of the subsonic aerodynamic performance of twisted and cambered wings of arbitrary planform with attainable thrust and vortex lift considerations are described. The computational system is based on a linearized theory lifting surface solution which provides a spanwise distribution of theoretical leading edge thrust in addition to the surface distribution of perturbation velocities. The approach used relies on a solution by iteration. The method also features a superposition of independent solutions for a cambered and twisted wing and a flat wing of the same planform to provide, at little additional expense, results for a large number of angles of attack or lift coefficients. A previously developed method is employed to assess the portion of the theoretical thrust actually attainable and the portion that is felt as a vortex normal force.
The study of the inversion algorithm for the single scatter lidar equation, for quantitative determination of cloud (or aerosol) optical properties, has received much attention over the last thirty years. Some of the difficulties associated with the solution of this equation are not yet solved. One problem is that a single scatter lidar equation has two unknowns. Because of this, the determination of the far-end boundary value, in the case of Klett's algorithm, is a problem if the atmosphere is optically inhomogeneous. Another difficulty concerns multiple scattering. There is a large error in the extinction distribution solution, in many cases, if only the single scattering component is considered, while neglecting the multiple scattering component. However, the use of multiple scattering in the remote sensing of aerosol or cloud optical properties is promising. In our early study, an inversion method for simultaneous determination of the cloud (or aerosol) Extinction Coefficient Distribution (ECD) and its Forward Scattering Phase Function (FSPF) was proposed according to multiply scattered lidar returns with two fields of view for the receiver. The method is based on a parameterized multiple scatter lidar equation. This paper is devoted to further numerical tests and an experimental study of lidar measurements of cloud ECD and FSPF using this method.
The extent, and some of the consequences, of nonequilibrium effects in the inviscid orbiter flow field are examined. Calculations have been made for the reacting flow along selected streamlines in the windward plane of symmetry of a blunted delta shape representative of a shuttle orbiter. An approximate model of the pressure field in the windward plane of symmetry of a blunted delta wing was constructed from wind tunnel measurements of the body surface pressure distribution and approximate calculations of the shock shape and streamline locations. Numerical solutions for the flow along streamlines were obtained for this pressure distribution, and solutions were obtained for both equilibrium and finite-reaction rate, nonequilibrium chemical reactions. Since the streamtube solutions are much easier to compute than fully three-dimensional, reacting flows, results were obtained for a wide range of conditions. The validity of this approach is based on the assumptions that the pressure field is insensitive to nonequilibrium effects and that the flow along streamlines is quasi-one dimensional. Numerical solutions have been used to verify this approach for other blunted slender shapes such as spherically blunted cones.
The new method presented for finding distribution functions, which depend only on the classical integrals of energy and angular momentum for stellar systems with known axisymmetric densities, is the analog for the axisymmetric case of Eddington's classical solution for the isotropic distribution function, depending only on energy, of a known spherical density. It is required that density be expressed as a function of the potential and of a radial coordinate. Our solution is also an integral which is derived directly from the density, and hence can be used with complicated densities. A numerical quadrature is generally required to evaluate this solution, but contour integrals can be computed accurately by numerical quadrature; this is preferable to an explicit evaluation if the latter is an infinite series, such as is obtained using Fricke's method. We give several examples, including some for which our distribution functions are new. Our method can be extended simply to the related problems of finding anisotropic distribution functions for spherical or disk systems.
Turbulent Couette flow between parallel plates was studied from a statistical mechanics approach utilizing a model equation, similar to the Boltzmann equation of kinetic theory, which was proposed by Lundgren from the velocity distribution of fluid elements. Solutions to this equation are obtained numerically, employing the discrete ordinate method and finite differences. Two types of boundary conditions on the distribution function are considered, and the results of the calculations are compared to available experimental data. The research establishes that Lundgren's equation provides a very good description of turbulence for the flow situation considered and that it offers an analytical tool for further study of more complex turbulent flows. The present work also indicates that modelling of the boundary conditions is an area where further study is required.
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 Boundary Force Method (BFM) was formulated for the three fundamental problems of elasticity: the stress boundary value problem, the displacement boundary value problem, and the mixed boundary value problem. Because the BFM is a form of an indirect boundary element method, only the boundaries of the region of interest are modeled. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate is used as the fundamental solution. Thus, unlike other boundary element methods, here the crack face 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 a center-cracked specimen subjected to mixed boundary conditions and a three-hole cracked configuration subjected to traction boundary conditions. The results obtained are in good agreement with accepted numerical solutions. The method is then used to generate stress-intensity solutions for two common cracked configurations: an edge crack emanating from a semi-elliptical notch, and an edge crack emanating from a V-notch. The BFM is a versatile technique that can be used to obtain very accurate stress intensity factors for complex crack configurations subjected to stress, displacement, or mixed boundary conditions. The method requires a minimal amount of modeling effort.
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.
Solutions of the radiative transfer equation describing photon interactions with vegetation canopies are important in remote sensing since they provide the canopy reflectance distribution required in the interpretation of satellite acquired information. The general one-dimensional two-angle transport problem for a finite copy of arbitrary leaf angle distribution is considered. Analytical solutions are obtained in terms of generalized Chandrasekhar's X- and Y-functions by invoking the principles of invariance. A critical step in the formulation involves the decomposition of the integral of the scattering phase function into a product of known functions of the incident and scattered photon directions. Several simplified cases previously considered in the literature are derived from the generalized solution. Various symmetries obeyed by the scattering operator and reciprocity relations are formally proved.
Theoretical analysis of pressure distribution, load capacity, load angle, and friction force for single-step concentric journal bearing
The dynamic normal-load distribution across a strip that is required to maintain a plane progressive wave along its length is studied for the case where the strip is of infinite length and lies on the surface of a homogeneous isotropic elastic half space. This configuration is proposed as a preliminary idealized model for analyzing the dynamic interaction between soils and flexible foundations. The surface load distribution across the strip and the motion of the strip are related by a pair of dual integral equations. An asymptotic solution is obtained for the limiting case of small wavelength. The nature of this solution depends importantly on the propagation velocity of the strip-traveling wave in comparison with the Rayleigh wave speed, the shear wave speed and the dilatational wave speed. When the strip-traveling wave propagates faster than the Rayleigh wave speed, a pattern of trailing Rayleigh waves is shed from the strip. The limiting amplitude of the trailing waves is provided by the asymptotic solution.
Heat transfer in cavity flow is numerically analyzed by a new numerical method called the finite-analytic method. The basic idea of the finite-analytic method is the incorporation of local analytic solutions in the numerical solutions of linear or nonlinear partial differential equations. In the present investigation, the local analytic solutions for temperature, stream function, and vorticity distributions are derived. When the local analytic solution is evaluated at a given nodal point, it gives an algebraic relationship between a nodal value in a subregion and its neighboring nodal points. A system of algebraic equations is solved to provide the numerical solution of the problem. The finite-analytic method is used to solve heat transfer in the cavity flow at high Reynolds number (1000) for Prandtl numbers of 0.1, 1, and 10.
Flat plate and stagnation flow heat transfer coefficients, similarity solutions of the laminar boundary layer for air in dissociation equilibrium and calculation of laminar heat-transfer distribution on blunt three-dimensional bodies in high speed flow
Three-dimensional numerical computations of natural convection in a cylindrical ampoule (L/R = 4) are presented for a Ga-doped Ge melt. The ampoule is maintained at isothermal end conditions, and the g vector is oriented at an angle gamma with the cylinder axis. Detailed velocity, temperature, and concentration field distributions are presented for different gravity levels and a range of inclination angles (gamma between 0 and 180 deg). For terrestrial conditions, complex, multicellular flow is found to occur for gamma between 0 and 180 deg. The strong convection results in significant isotherm distortions and enhanced heat transfer. As the gravity level is reduced to 0.001 g(0), the overall convection strength decays but the 3D multicellular flow persists, causing appreciable mass transfer. For a further reduction in the gravity level to 0.00001 g(0), the system thermo/solutal characteristics are close to a purely diffusive behavior. Flow and thermal and solutal distributions are presented in the principal, orthogonal, and cross sectional planes, and the average system heat and mass transfer are calculated.
A procedure based on the method of similar solutions is presented by which the skin friction, heat transfer, and boundary-layer thickness in a laminar hypersonic flow with pressure gradient may be rapidly evaluated if the pressure distribution is known. This solution, which at present is. restricted to power-law variations of pressure with surface distance, is presented for a wide range of exponents in the power law corresponding to both favorable and adverse pressure gradients. This theory has been compared to results from heat-transfer experiments on blunt-nose flat plates and a hemisphere cylinder at free-stream Mach numbers of 4 and 6.8. The flat-plate experiments included tests made at a Mach number of 6.8 over a range of angle of attack of +/- 10 deg. Reasonable agreement of the experimental and theoretical heat-transfer coefficients has been obtained as well as good correlation of the experimental results over the entire range of angle of attack studied. A similar comparison of theory with experiment was not feasible for boundary-layer-thickness data; however, the hypersonic similarity theory was found to account satisfactorily for the variation in boundary-layer thickness due to local pressure distribution for several sets of measurements.
The measured macroscopic growth rates of the (110) and (101) faces of tetragonal lysozyme show an unexpectedly complex dependence on the supersaturation. The growth rates decay asymptotically to zero when the supersaturation is lowered to zero and increase rapidly when the supersaturation is increased. When supersaturations are increased still further the growth rates attain a maximum before starting to decrease. However, growth of these crystals is known to proceed by the classical dislocation and 2D nucleation growth mechanisms. This anomaly can be explained if growth is assumed to occur not by monomer units but by lysozyme aggregates. Analysis of the molecular packing of these crystals revealed that they were constructed of strongly bonded 4(sub 3) helices, while weaker bonds were responsible for binding the helices to each other. It follows that during crystal growth the stronger bonds are formed before the weaker ones. Thus, the growth of these crystals could be viewed as a two step process: aggregate growth units corresponding to the 4(sub 3) helix are first formed in the bulk solution by stronger intermolecular bonds and then attached to the crystal face by weaker bonds on dislocation hillocks or 2D islands. This will lead to a distribution of aggregates in the solution with monomers and lower order aggregates being predominant at low supersaturations and higher order aggregates being predominant at high supersaturations. If the crystal grows mostly by higher order aggregates, such as tetramers and octamers, it would explain the anomalous dependence of the growth rates on the supersaturation. Besides the analysis of molecular packing, a comprehensive analysis of the measured (110) and (101) growth rates was also undertaken in this study. The distribution of aggregates in lysozyme nutrient solutions at various solution conditions were determined from reversible aggregation reactions at equilibrium. The supersaturation was defined for each aggregate species with respect to its concentration at saturation in order to apply growth rate models to this process. The measured growth rates were then compared with the predicted ones from several dislocation and 2D nucleation growth models, employing tetramer and octamer growth units in polydisperse solutions and monomer units in monodisperse solutions. For the (110) face, the calculations consistently showed that the measured growth rates followed the expected model relations with octamer growth units. For the (101) face, it is not possible to obtain a clear agreement between the predicted and measured growth rates for a single growth unit as done for the (110) face. However, the calculations do indicate that the average size of the growth unit is between a tetramer and an octamer. This suggests that tetramers, octamers and other intermediate size growth units all participate in the growth process for this face. These calculations show that it is possible to model the macroscopic protein crystal growth rates if the molecular level processes can be account for, particularly protein aggregation processes in the bulk solution. Our recent investigations of tetragonal lysozyme crystals employing high resolution atomic force microscopy scans have further confirmed the growth of these crystals by aggregate growth units corresponding to 4(sub 3) helices.