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.
We present a method of conversion of the lunar neutron counting rate measured by the Lunar Reconnaissance Orbiter (LRO) Lunar Exploration Neutron Detector (LEND) instrument collimated neutron detectors, to water equivalent hydrogen (WEH) in the top approximately 1 m layer of lunar regolith. Polar maps of the Moon’s inferred hydrogen abundance are presented and discussed.
The modeling of the behavior of fluids in low- or microgravity is a field with little validated data but of great interest to space agencies and companies in the face of humanity’s expanding presence beyond low Earth orbit. Efforts to refuel and repair satellites in geostationary orbit require an understanding of the behavior of sloshing fuel in microgravity. The stability of remote sensing equipment and satellites depends on robust fluid slosh dampening or very accurate modeling, and the Artemis program includes the refueling of large landing systems in its architecture. However, current methodologies are quite binary: either use a relatively low-accuracy equivalent mechanical model or use computational fluid dynamics which are not practical for onboard systems due to their relatively high computational cost. In this work, a novel method of tracking the bubble rather than tracking the fuel itself is presented. Computational fluid dynamics simulations are used to validate the model based on available experimental data, and the dynamics of a refueler spacecraft based on Gateway’s Logistics Module are simulated both for the sloshing case and without sloshing considered.
A promising path for developing a low-boom configuration is a multifidelity approach that (1) starts from a low-fidelity low-boom design, (2) refines the low-fidelity design with computational fluid dynamics (CFD) equivalent-area (Ae) analysis, and (3) improves the design with sonic-boom analysis by using CFD off-body pressure distributions. The focus of this paper is on the third step of this approach, in which the design is improved with sonic-boom analysis through the use of CFD calculations. A new inverse design process for off-body pressure tailoring is formulated and demonstrated with a low-boom supersonic configuration that was developed by using the mixed-fidelity design method with CFD Ae analysis. The new inverse design process uses the reverse propagation of the pressure distribution (dp/p) from a mid-field location to a near-field location, converts the near-field dp/p into an equivalent-area distribution, generates a low-boom target for the reversed equivalent area (Ae,r) of the configuration, and modifies the configuration to minimize the differences between the configuration s Ae,r and the low-boom target. The new inverse design process is used to modify a supersonic demonstrator concept for a cruise Mach number of 1.6 and a cruise weight of 30,000 lb. The modified configuration has a fully shaped ground signature that has a perceived loudness (PLdB) value of 78.5, while the original configuration has a partially shaped aft signature with a PLdB of 82.3.
Here, we present a scalable approach to solve a class of partial differential equation (PDE)‐constrained optimization problems with bound constraints. This approach utilizes a robust full‐space interior‐point (IP)‐Gauss–Newton optimization method. To cope with the poorly‐conditioned IP‐Gauss–Newton saddle‐point linear systems that need to be solved approximately, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE‐constrained optimization problems. A block Gauss–Seidel preconditioner is proposed for the GMRES‐based solution of the IP‐Gauss–Newton linear systems. It is shown, for a large‐class of PDE‐ and bound‐constrained optimization problems, that the spectrum of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix is asymptotically independent of discretization and is not impacted by the ill‐conditioning that notoriously plagues interior‐point methods. We exploit symmetry of the IP‐Gauss–Newton linear systems and propose a regularization and log‐barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)‐based solution of the equivalent IP‐Gauss–Newton–Schur complement linear systems. The eigenvalues of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix, that are not equal to one, are identical to the eigenvalues of the regularization and log‐barrier Hessian preconditioned Schur complement matrix. The scalability of the approach is demonstrated on two example problems. The numerical solution of these optimization problems is shown to require a discretization independent number of IP‐Gauss–Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill‐conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of the preconditioner, achieved via algebraic multigrid component solvers when applicable, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE‐ and bound‐constrained problems.
A rigorous mathematical framework is provided for a substructuring-based domain-decomposition approach for nonlocal problems that feature interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal domain-decomposition methods.
A rigorous mathematical framework is provided for a substructuring-based domain-decomposition approach for nonlocal problems that feature interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal domain-decomposition methods.
Materials for space applications must have low outgassing rates at normal operating temperatures, and the outgassing products should include a minimum of condensables at the temperatures of nearby surfaces. A screening method, developed several years ago and used at many space laboratories, consists of holding a material sample at 398 K (125 C) for 24 hours and measuring its percentage total mass loss (TML) and the percentage volatile condensable mass (VCM) accreted on a 298-K (25-C) collector. In general, the material is acceptable if the TML is less than 1 percent and the VCM is less than 0.1 percent. An analysis of the test and its results is presented.
When designing a Mesosphere-Stratosphere-Troposphere (MST) radar antenna one has to trade between the choices to optimize the effective aperture or to optimize the sidelobe suppression. An optimization of the aperture increases the sensitivity. Suppression of side-lobes by tapering attenuates undesirable signals which spoil the estimates of reflectivity and velocity. Generally, any sidelobe effects are equivalent to a broadening of the antenna beam. The return signal is due to a product of the antenna pattern with the varying atmospheric reflectivity structures. Thus, knowing the antenna pattern, it is in principle possible to find the signal spectra, which, however, may be a tedious computational and ambiguous procedure. For vertically pointing main beams the sidelobe effects are efficiently suppressed because of the aspect sensitivity. It follows that sidelobes are a minor problem for spaced antenna methods. However, they can be crucial for Doppler methods, which need off-vertical beams. If a sidelobe is pointing towards the zenith a larger power may be received from the vertical than off-vertical directions, but quantitative estimates of this effect are not yet known. To get an error estimate of sidelobe effects with an off-vertical main beam a 1-dimensional example is considered.
The pairwise kinematic Sunyaev–Zel’dovich (kSZ) effect measures both the pairwise motion between galaxy groups and clusters and the amount of gas within them, providing a tracer for cosmic growth. To interpret the cosmological information in the kSZ measurements, it is crucial to understand the optical-cluster-selection bias on the kSZ observables. Line-of-sight structures that contribute to both the optical observable (e.g., richness) and the cosmological signal can induce a correlation between these two quantities at a fixed cluster mass. The selection bias arising from this correlation is a key systematic effect for cosmological analyses. For cosmological observables such as cluster abundance and weak lensing, controlling this selection bias may help explain the tension between the DES-Y1 results and the Planck constraints. In order to test for a kSZ effect equivalent of such a bias, we adopted an alternative mock richness based on galaxy counts within cylindrical volumes along the line of sight. We applied the cylindrical count method to hydrodynamical simulations across a wide range of galaxy-selection criteria, assigning richness consistent with DES-Y1 to the mock clusters. When comparing optically selected clusters to mass-selected halos, we find no significant bias on pairwise kSZ signals, pairwise velocities, or optical depth within our uncertainty limits of approximately 16, 10, and 8%, respectively.
A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.
Radiotelescopic observations of an impulsive outburst of the sun on Mar. 27, 1969 are discussed in terms of the location and size of the source region. Two different methods have been used for this purpose. It is concluded that the outburst region was 30 seconds or less in diameter. Based on this diameter, a peak source equivalent temperature of 7 million K is indicated.
Specifications for shock testing of components that will be used on the Space Shuttle vehicles require very high acceleration levels. A special shock machine was built for testing of rocket components to determine if they can meet the specified accelerations. Calibrations of transducers and methods to monitor the shock tests raised several signature-analysis questions. In this report, calibration capabilities of shock accelerometers are found to be limited to 10,000g. Equivalency of the mechanical shock test and the rocket pyrotechnic shock are examined, and two simple relationships for equivalency are proposed. Five different pulse signature-analysis techniques are tested on analytical and experimental pulse data and recommendations are made for the signature technique which most clearly identifies the magnitude of the impulse applied to the test specimen.
Laser speckle velocimetry, the measurement of fluid velocity by measuring the translation of speckle pattern or individual particles that are moving with the fluid, is described. The measurement is accomplished by illuminating the fluid with consecutive pulses of Laser Light and recording the images of the particles or the speckles on a double exposed photographic plate. The plate contains flow information throughout the image plane so that a single double exposure may provide data at hundreds or thousands of points in the illuminated region of the fluid. Conventional interrogation of the specklegram involves illuminating the plate to form Young's fringes, whose spacing is inversely proportional to the speckle separation. Subsequently the fringes are digitized and analyzed in a computer to determine their frequency and orientation, yielding the velocity magnitude and orientation. The Young's fringe technique is equivalent to performing a 2-D spatial correlation of the double exposed specklegram intensity pattern, and this observation suggests that correlation should be considered as an alternative processing method. The principle of the correlation technique is examined.
A method based on the Clebsch velocity decomposition is presented to solve the steady, inviscid, supersonic flow field about arbitrary conical geometries. The system of equations developed, although formally equivalent to the Euler equations, retains the computational efficiency of type-dependent potential flow solutions. Accurate rotational-flow solutions are developed using shock-fitting procedures at the bow and imbedded waves along with special treatment of the vortical layer. Solutions are presented for several high Mach number conical flows and compared with existing Euler solutions and experimental data.
The derivation, logic, and performance of a simplified atmospheric guidance algorithm for aeroassist orbital-transfer vehicles (AOTVs) are presented. The algorithm was developed to meet the demands for an aerobraking trajectory guidance technique that was uncomplicated, easily integrated into existing trajectory programs, adaptable to a range of vehicle aerodynamic configurations, capable of performance equivalent to currently available guidance programs in compensating for dispersions in entry conditions, atmospheric conditions, and aerodynamic characteristics. The result was a hybrid lifting guidance algorithm combining the method of reference-profile generation with the method of predictor/corrector schemes. The resulting performance is good (less than 3 n.m. error from desired apogee despite uncertainties of + or - 50 percent atmospheric density, + or - 0.2 deg entry flight-path angle, or + or - 50 percent L/D. Combinations of these same dispersions with lesser magnitudes have also been successful, although performance with density 'pockets' within the atmosphere requires more analysis.
Rotating drive improves quality of holes made by electrical-arc machining. Mechanism (Uni-tek, rotary head, or equivalent) attached to electrical-arc system. Drive rotates electrode as though it were mechanical drill, while an arc disintegrates metal in workpiece, thereby creating hole. Rotating electrode method often used in electric-discharge machining. NASA innovation is application of technique to electrical-arc machining.
A computer code for fast automatic generation of quasi-three-dimensional grid systems for aerospace configurations is described. The code employs a homotopic method to algebraically generate two-dimensional grids in cross-sectional planes, which are stacked to produce a three-dimensional grid system. Implementation of the algebraic equivalents of the homotopic relations for generating body geometries and grids are explained. Procedures for controlling grid orthogonality and distortion are described. Test cases with description and specification of inputs are presented in detail. The FORTRAN computer program and notes on implementation and use are included.
The last decade has seen tremendous growth in cloud dynamical and microphysical models that are able to simulate storms and storm systems with very high spatial resolution, typically of the order of a few kilometers. The fairly realistic distributions of cloud and hydrometeor properties that these models generate has in turn led to a renewed interest in the three-dimensional microwave radiative transfer modeling needed to understand the effect of cloud and rainfall inhomogeneities upon microwave observations. Monte Carlo methods, and particularly backwards Monte Carlo methods have shown themselves to be very desirable due to the quick convergence of the solutions. Unfortunately, backwards Monte Carlo methods are not well suited to treat polarized radiation. This study reviews the existing Monte Carlo methods and presents a new polarized Monte Carlo radiative transfer code. The code is based on a forward scheme but uses aliasing techniques to keep the computational requirements equivalent to the backwards solution. Radiative transfer computations have been performed using a microphysical-dynamical cloud model and the results are presented together with the algorithm description.