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.
A general purpose solver package for constructing and solving a range of sparse linear systems arising from discretization of PDEs on unstructured meshes is developed. Once the sparse symmetric complex matrix is constructed, it can be solved by either a preconditioned bi-conjugate gradient solver, a two-stage Cholesky LDLT factorization solver, or a hybrid solver combining the above two methods.
A general purpose solver package for constructing and solving a range of sparse linear systems arising from discretization of PDEs on unstructured meshes is developed. Once the sparse symmetric complex matrix is constructed, it can be solved by either a preconditioned bi-conjugate gradient solver, a two-stage Cholesky LDLT factorization solver, or a hybrid solver combining the above two methods. (More detailed than 95-0127).
The scattering of rotor noise is an area that has received little attention over the years, yet the limited work that has been done has shown that both the directivity and intensity of the acoustic field may be significantly modified by the presence of scattering bodies. One of the inputs needed to compute the scattered acoustic field is the acoustic pressure gradient on a scattering surface. Two new analytical formulations of the acoustic pressure gradient have been developed and implemented in the PSU-WOPWOP rotor noise prediction code. These formulations are presented in this paper. The first formulation is derived by taking the gradient of Farassat's retarded-time Formulation 1A. Although this formulation is relatively simple, it requires numerical time differentiation of the acoustic integrals. In the second formulation, the time differentiation is taken inside the integrals analytically. The acoustic pressure gradient predicted by these new formulations is validated through comparison with the acoustic pressure gradient determined by a purely numerical approach for two model rotors. The agreement between analytic formulations and numerical method is excellent for both stationary and moving observers case.
In this paper, we propose an adaptive shape servoing method to deform a soft object into a desired 3-D shape. The high dimensional representation and the unknown deformation properties of the soft object pose a challenge to actively manipulate its shape. To address this issue, we develop a method to compute the deformation Jacobian matrix in real-time. The Jacobian is estimated using a set of basis functions and its corresponding parameters to capture the dynamics of the system and relate the applied input motion to changes in the soft object's shape. An integral concurrent learning (ICL) based adaptive update law is derived using Lyapunov analysis to estimate the deformation parameters and prove its convergence. A physics-based simulation is used to validate the proposed method and controller by performing manipulation tasks with different desired configurations. The performance is compared with a standard gradient update law to demonstrate the accuracy and robustness of our approach.
A heritage wine-rack thermal/mechanical design for the nickel-hydrogen batteries was the baseline at the Landsat-7 Preliminary Design Review. An integrated thermal and power analysis of the batteries performed by the author in 1994 revealed that the maximum cell-to-cell gradient was 6.6 C. The author proposed modifying the heritage wine-rack design by enhancing heat conduction from cells to cells, and from cells to battery frame. At the 1995 Intersociety Energy Conversion Engineering Conference (IECEC), the author presented a paper on methods of modifying the wine-rack design. It showed that the modified wine-rack option, which uses a metallic filler, could reduce the maximum cell-to-cell temperature gradient to 1.30 C, and could also reduce the maximum cell temperature by as much as 80 C. That design concept was adopted by the Landsat7 Project Office, and a design change was made at the Critical Design Review. Results of the spacecraft thermal vacuum and thermal balance tests, and temperature data in flight show that the temperatures of the battery cells are very uniform. The maximum cell-to-cell gradient is 1.50 C. They validate the modified wine-rack thermal design.
A time-dependent extension of the first-order hyperbolic system method for advection-diffusion problems is introduced. Diffusive/viscous terms are written and discretized as a hyperbolic system, which recovers the original equation in the steady state. The resulting scheme offers advantages over traditional schemes: a dramatic simplification in the discretization, high-order accuracy in the solution gradients, and orders-of-magnitude convergence acceleration. The hyperbolic advection-diffusion system is discretized by the second-order upwind residual-distribution scheme in a unified manner, and the system of implicit-residual-equations is solved by Newton's method over every physical time step. The numerical results are presented for linear and nonlinear advection-diffusion problems, demonstrating solutions and gradients produced to the same order of accuracy, with rapid convergence over each physical time step, typically less than five Newton iterations.
A method of general perturbations, based on the use of Lie series to generate approximate canonical transformations, is applied to study the long-term effects of gravity-gradient torque and orbital evolution on the rotational motion of a triaxial, rigid satellite. The center of mass of the satellite is constrained to move in an elliptic orbit about an attracting point mass. The orbit, which has a constant inclination, is constrained to precess and spin with constant rates. The method of general perturbations is used to obtain the Hamiltonian for the nonresonant secular and long-period rotational motion of the satellite to second order in n/omega sub 0, where n is the orbital mean motion of the center of mass and omega sub 0 is a reference value of the magnitude of the satellite's rotational angular velocity. The differential equations derivable from the transformed Hamiltonian are integrable, and the solution for the long-term motion may be expressed in terms of Jacobian elliptic functions and elliptic integrals. Geometrical aspects of the long-term rotational motion are discussed, and a comparison of theoretical results with observations is made.
The simulation of the diffusion of a continuum field by the random walk (RW) displacement of a set of particles is considered. Elements of the gradients of the diffusive concentration are transported by computational particles. It is demonstrated that, by the use of concentration gradients in the RW process, statistical errors are reduced and each realization of the numerical solution is a representation of the exact solution. The algorithm is grid-free, and the computational elements move to follow the gradients; hence, the algorithm is self-adaptive, and uniform resolution is achieved for all times.
A significant reduction in the drag of transonic aircraft can be achieved by using an active method of boundary-layer control to maintain laminar flow on the aerodynamic surfaces. The stabilizing influence of suction for both favorable and adverse pressure gradients is demonstrated by means of the self-similar incompressible three-dimensional boundary layers on yawed wedges. Profile instability is measured by the maximum amplification rate of fixed-frequency disturbances computed according to linearized, locally parallel, spatial stability theory. Suction is found to be more effective in controlling Tollmien-Schlichting instability than stationary cross-flow disturbances. The effectiveness of surface cooling as a method of stabilization is compared with suction for the boundary layers on a transonic 35 deg swept wing of infinite span. It is found from compressible stability theory that when the surface is cooled to a uniform temperature such that the maximum cross-flow velocity is reduced by the same amount as with a uniform suction distribution, the stabilizing effect of cooling on cross-flow disturbances is less than with suction.
The effects of fin leading-edge radius and sweep angle on peak heating rates due to shock-shock interactions were investigated in the NASA Langley Research Center 20-inch Mach 6 Air Tunnel. The fin model leading edges, which represent cylindrical leading edges or struts on hypersonic vehicles, were varied from 0.25 inches to 0.75 inches in radius. A 9deg wedge generated a planar oblique shock at 16.7deg to the flow that intersected the fin bow shock, producing a shock-shock interaction that impinged on the fin leading edge. The fin angle of attack was varied from 0deg (normal to the free-stream) to 15deg and 25deg swept forward. Global temperature data was obtained from the surface of the fused silica fins using phosphor thermography. Metal oil flow models with the same geometries as the fused silica models were used to visualize the streamline patterns for each angle of attack. High-speed zoom-schlieren videos were recorded to show the features and temporal unsteadiness of the shock-shock interactions. The temperature data were analyzed using one-dimensional semi-infinite as well as one- and two-dimensional finite-volume methods to determine the proper heat transfer analysis approach to minimize errors from lateral heat conduction due to the presence of strong surface temperature gradients induced by the shock interactions. The general trends in the leading-edge heat transfer behavior were similar for the three shock-shock interactions, respectively, between the test articles with varying leading-edge radius. The dimensional peak heat transfer coefficient augmentation increased with decreasing leading-edge radius. The dimensional peak heat transfer output from the two-dimensional code was about 20% higher than the value from a standard, semi-infinite onedimensional method.
The Knudsen cell technique has been used for over a century and is a valuable tool for measurement of vapor pressures and thermodynamic properties. It is based on a small enclosure (~1 cm long x 1 cm diameter) in which a condensed phase/vapor equilibria can be established. A small (<1 mm) orifice on the cell allows sampling of the vapor via a variety of techniques including weight loss, torsion effusion, target collection, and mass spectrometry. Many excellent measurements based on these methods have been reported. However in order to obtain reliable measurements, a variety of factors must be considered. They include proper cell material selection, accurate and uniform temperature control and measurement, and proper sampling of the vapor. Each of these factors are discussed in detail in this chapter. Typically these studies are conducted at high temperatures and it is a challenge to select an inert container material. Recommended materials are discussed and in some cases the container may be used as part of the system under study. Temperature control and measurement is perhaps the most important issue. In most systems, the furnace must be compact yet there can be no temperature gradient in the cell. Temperatures are measured with either a thermocouple or pyrometer and the relative advantages of each are discussed. Sampling method considerations depend on the particular technique. It is essential that all of the vapor or a representative portion of the vapor be sampled. The distribution of the effusate from a Knudsen cell is discussed and sampling positions discussed. Mass spectrometry is often used to study the effusing vapor and the relations between ion current and vapor pressure are discussed.
The effects of fin-leading-edge radius and sweep angle on peak heating rates due to shock-shock interactions were investigated in the NASA Langley Research Center 20-inch Mach 6 Air Tunnel. The fin model leading edges, which represent cylindrical leading edges or struts on hypersonic vehicles, were varied from 0.25 inches to 0.75 inches in radius. A 9deg wedge generated a planar oblique shock at 16.7deg to the flow that intersected the fin bow shock, producing a shock-shock interaction that impinged on the fin leading edge. The fin angle of attack was varied from 0deg (normal to the free-stream) to 15deg and 25deg swept forward. Global temperature data was obtained from the surface of the fused silica fins through phosphor thermography. Metal oil flow models with the same geometries as the fused silica models were used to visualize the streamline patterns for each angle of attack. High-speed zoom-schlieren videos were recorded to show the features and temporal unsteadiness of the shock-shock interactions. The temperature data were analyzed using one-dimensional semi-infinite as well as one- and two-dimensional finite-volume methods to determine the proper heat transfer analysis approach to minimize errors from lateral heat conduction due to the presence of strong surface temperature gradients induced by the shock interactions. The general trends in the leading-edge heat transfer behavior were similar for the three shock-shock interactions, respectively, between the test articles with varying leading-edge radius. The dimensional peak heat transfer coefficient augmentation increased with decreasing leading-edge radius. The dimensional peak heat transfer output from the two-dimensional code was about 20% higher than the value from a standard, semi-infinite one-dimensional method.
The effects of fin-leading-edge radius and sweep angle on peak heating rates due to shock-shock interactions were investigated in the NASA Langley Research Center 20-Inch Mach 6 Air Tunnel. The fin model leading edges, which represent cylindrical leading edges or struts on hypersonic vehicles, were varied from 0.25 inches to 0.75 inches in radius. A 9° wedge generated a planar oblique shock at 16.7° to the flow that intersected the fin bow shock, producing a shock-shock interaction that impinged on the fin leading edge. The fin angle of attack was varied from 0° (normal to the free-stream) to 15° and 25° swept forward. Global temperature data were obtained from the surface of the fused silica fins through phosphor thermography. Metal oil flow models with the same geometries as the fused silica models were used to visualize the streamline patterns for each angle of attack. High-speed zoom-schlieren videos were recorded to show the features and temporal unsteadiness of the shock-shock interactions. The temperature data were analyzed using one-dimensional semi-infinite as well as one- and two-dimensional finite-volume methods to determine the proper heat transfer analysis approach to minimize errors from lateral heat conduction due to the presence of strong surface temperature gradients induced by the shock interactions. The general trends in the leading-edge heat transfer behavior were similar for the three shock-shock interactions, respectively, between the test articles with varying leading-edge radius. The dimensional peak heat transfer coefficient augmentation increased with decreasing leading-edge radius. The dimensional peak heat transfer output from the two-dimensional code was about 20% higher than the value from a standard, semi-infinite one-dimensional method.
Quantum networking protocols relying on interference and precise time-of-flight measurements require high precision clock synchronization. This study describes the design, implementation, and characterization of two optical time transfer methods in a metropolitan-scale quantum networking research testbed. The first technique called active electronic stabilization (ELSTAB), was able to achieve sub-picosecond time deviation (TDEV) at integration times between 1 s and 105s over multiple different fiber lengths of both underground and aerial deployed fiber. The second technique, White Rabbit-Precision Time Protocol (WR-PTP), was able to achieve 10-picosecond TDEV’s over similar fiber lengths at comparable integration times. Environmental fluctuations affected the stability of clock synchronization over deployed fiber. These fluctuations resulted in path delay gradients, chromatic dispersion, polarization drift, and optical power variations. The results from this study will inform future work in the development of compensation methods essential for enabling experimental research in developing practical quantum networking protocols.
Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.
A chromatographic method for the quantitation of promethazine (PMZ) and its three metabolites in urine employing on-line solid-phase extraction and column-switching has been developed. The column-switching system described here uses an extraction column for the purification of PMZ and its metabolites from a urine matrix. The extraneous matrix interference was removed by flushing the extraction column with a gradient elution. The analytes of interest were then eluted onto an analytical column for further chromatographic separation using a mobile phase of greater solvent strength. This method is specific and sensitive with a range of 3.75-1400 ng/ml for PMZ and 2.5-1400 ng/ml for the metabolites promethazine sulfoxide, monodesmethyl promethazine sulfoxide and monodesmethyl promethazine. The lower limits of quantitation (LLOQ) were 3.75 ng/ml with less than 6.2% C.V. for PMZ and 2.50 ng/ml with less than 11.5% C.V. for metabolites based on a signal-to-noise ratio of 10:1 or greater. The accuracy and precision were within +/- 11.8% in bias and not greater than 5.5% C.V. in intra- and inter-assay precision for PMZ and metabolites. Method robustness was investigated using a Plackett-Burman experimental design. The applicability of the analytical method for pharmacokinetic studies in humans is illustrated.
Since 1986, the NASA-Langley Research Center has sponsored the NASA-UVa Light Alloy and Structures Technology (LA2ST) Program at the University of Virginia (UVa). The fundamental objective of the LA2ST program is to conduct interdisciplinary graduate student research on the performance of next generation, light-weight aerospace alloys, composites and thermal gradient structures. The LA2ST program has aimed to product relevant data and basic understanding of material mechanical response, environmental/corrosion behavior, and microstructure; new monolithic and composite alloys; advanced processing methods; measurement and modeling advances; and a pool of educated graduate students for aerospace technologies. The scope of the LA2ST Program is broad. Research areas include: (1) Mechanical and Environmental Degradation Mechanisms in Advanced Light Metals and Composites, (2) Aerospace Materials Science, (3) Mechanics of materials for Aerospace Structures, and (4) Thermal Gradient Structures. A substantial series of semi-annual progress reports issued since 1987 documents the technical objectives, experimental or analytical procedures, and detailed results of graduate student research in these topical areas.
Approximate analytic expressions are derived for resonance-line wing diagnostics, accounting for frequency redistribution effects, for homogeneous slabs, and slabs with a constant Planck function gradient. Resonance-line emission profiles from a simplified conceptual standpoint are described in order to elucidate the basic physical parameters of the line-forming layers prior to the performance of detailed numerical calculations. An approximate analytic expression is derived for the dependence on stellar surface gravity of the location of the Ca II and Mg II resonance-line profile peaks. An approximate radiative transfer equation using generalized second-order escape probabilities, applicable even in the presence of nearly coherent scattering in the damping wings of resonance lines, is derived. Approximate analytic solutions that can be applied in special regimes and achieve good agreement with accurate numerical results are found.