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.
Generalized method, including computer program, to solve problem of heat transfer and material response for stagnation region of blunt bodies experiencing ablation
This paper presents a path planning method for sea surface vehicles that prevents capsizing and bow-diving in a high sea-state. A key idea is to use response amplitude operators (RAOs) or, in control terminology, the transfer functions from a sea state to a vessel's motion, in order to find a set of speeds and headings that results in excessive pitch and roll oscillations. This information is translated to arithmetic constraints on the ship's velocity, which are passed to a model predictive control (MPC)-based path planner to find a safe and optimal path that achieves specified goals. An obstacle avoidance capability is also added to the path planner. The proposed method is demonstrated by simulations.
The present paper reviews the basic concepts of nonlinear response of panels to surface flow and acoustic pressures, simulation of random processes, time domain solutions and the Monte Carlo Method. Applications of this procedure to the orbit-on-demand space vehicles, acoustic fatigue and composite materials are discussed. Numerical examples are included for a variety of nonlinear problems to illustrate the applicability of this method.
The measurement conditions are described for an intensive field campaign at White Sands Missile Range for the calibration of the AVHRRs on NOAA-9, NOAA-10 and NOAA-11, LANDSAT-4 TM and SPOT. Three different methods for calibration of AVHRRs by reference to a ground surface site are reported, and results from these methods are compared. Significant degradations in NOAA-9 and NOAA-10 AVHRR responsivities occurred since prelaunch calibrations were completed. As of February 1988, degradations in NOAA-9 AVHRR responsivities were on the order of 37 percent in channel and 41 percent in channel 2, and for the NOAA-10 AVHRR these degradations were 42 and 59 percent in channels 1 and 2, respectively.
An objective assessment of the impact of a new cloud scheme, called Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme (McRAS) (together with its radiation modules), on the finite volume general circulation model (fvGCM) was made with a set of ensemble forecasts that invoke performance evaluation over both weather and climate timescales. The performance of McRAS (and its radiation modules) was compared with that of the National Center for Atmospheric Research Community Climate Model (NCAR CCM3) cloud scheme (with its NCAR physics radiation). We specifically chose the boreal summer months of May and June 2003, which were characterized by an anomalously wet eastern half of the continental United States as well as northern regions of Amazonia. The evaluation employed an ensemble of 70 daily 10-day forecasts covering the 61 days of the study period. Each forecast was started from the analyzed initial state of the atmosphere and spun-up soil moisture from the first-day forecasts with the model. Monthly statistics of these forecasts with up to 10-day lead time provided a robust estimate of the behavior of the simulated monthly rainfall anomalies. Patterns of simulated versus observed rainfall, 500-hPa heights, and top-of-the-atmosphere net radiation were recast into regional anomaly correlations. The correlations were compared among the simulations with each of the schemes. The results show that fvGCM with McRAS and its radiation package performed discernibly better than the original fvGCM with CCM3 cloud physics plus its radiation package. The McRAS cloud scheme also showed a reasonably positive response to the observed sea surface temperature on mean monthly rainfall fields at different time leads. This analysis represents a method for helpful systematic evaluation prior to selection of a new scheme in a global model.
The surface responses due to impact load on an infinite uniaxial graphite/epoxy plate with and without delamination cracks are investigated both in time and frequency domain by using a hybrid method combining the finite element discretization of the near-field with boundary integral representation of the field outside a contour enclosing completely the crack. This combined method leads to a set of linear unsymmetric complex matrix equations which are solved to obtain the response in the frequency domain by biconjugate gradient method. The time domain response is then obtained by using an FFT. In order to capture the time-domain characteristics accurately, high order finite elements have been used. Also, both the six node singular elements and eight node transition elements are used around the crack tips to model the crack-tip singularity. It is shown that from the numerical results for surface responses both depth and length of this crack can be identified.
The report describes a method for performing optimization of a system whose analysis is so expensive that it is impractical to let the optimization code invoke it directly because excessive computational cost and elapsed time might result. In such situation it is imperative to have user control the number of times the analysis is invoked. The reported method achieves that by two techniques in the Design of Experiment category: a uniform dispersal of the trial design points over a n-dimensional hypersphere and a response surface fitting, and the technique of krigging. Analyses of all the trial designs whose number may be set by the user are performed before activation of the optimization code and the results are stored as a data base. That code is then executed and referred to the above data base. Two applications, one of the airborne laser system, and one of an aircraft optimization illustrate the method application.
The 3m Spherical Primary Optical Telescope (SPOT) will utilize a single ring of 0.86111 point-to-point hexagonal mirror segments. The f2.85 spherical mirror blanks will be fabricated by the same replication process used for mass-produced commercial telescope mirrors. Diffraction-limited phasing will require segment-to-segment radius of curvature (ROC) variation of approx.1 micron. Low-cost, replicated segment ROC variations are estimated to be almost 1 mm, necessitating a method for segment ROC adjustment & matching. A mechanical architecture has been designed that allows segment ROC to be adjusted up to 400 microns while introducing a minimum figure error, allowing segment-to-segment ROC matching. A key feature of the architecture is the unique back profile of the mirror segments. The back profile of the mirror was developed with shape optimization in MSC.Nastran(TradeMark) using optical performance response equations written with SigFit. A candidate back profile was generated which minimized ROC-adjustment-induced surface error while meeting the constraints imposed by the fabrication method. Keywords: optimization, radius of curvature, Pyrex spherical mirror, Sigfit
In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code, are used with a blowing boundary condition in the DPLR hypersonic CFD code. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver, is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.
The need for a computer program to perform kinematic and dynamic analyses of large truss structures while deploying from a packaged configuration in space led to the evaluation of several existing programs. ADAMS (automatic dynamic analysis of mechanical systems), a generalized program from performing the dynamic simulation of mechanical systems undergoing large displacements, is applied to two concepts of deployable space antenna units. One concept is a one cube folding unit of Martin Marietta's Box Truss Antenna and the other is a tetrahedral truss unit of a Tetrahedral Truss Antenna. Adequate evaluation of dynamic forces during member latch-up into the deployed configuration is not yet available from the present version of ADAMS since it is limited to the assembly of rigid bodies. Included is a method for estimating the maximum bending stress in a surface member at latch-up. Results include member displacement and velocity responses during extension and an example of member bending stresses at latch-up.
An adaptive control method is presented for the minimization of drag during flight for transport aircraft. The minimization of drag is achieved by taking advantage of the redundant control capability available in the pitch axis, with the horizontal tail used as the primary surface and symmetric deflection of the ailerons and cruise flaps used as additional controls. The additional control surfaces are excited with sinusoidal signals, while the altitude and velocity loops are closed with guidance and control laws. A model of the throttle response as a function of the additional control surfaces is formulated and the parameters in the model are estimated from the sensor measurements using a least squares estimation method. The estimated model is used to determine the minimum drag positions of the control surfaces. The method is presented for the optimization of one and two additional control surfaces. The adaptive control method is extended to optimize rate of climb with the throttle fixed. Simulations that include realistic disturbances are presented, as well as the results of a Monte Carlo simulation analysis that shows the effects of changing the disturbance environment and the excitation signal parameters.
A document discusses the use of ultraviolet spectroscopy and imaging for the non-destructive evaluation of the degree of cure, aging, and other properties of resin-based composite materials. This method can be used in air, and is portable for field use. This method operates in reflectance, absorbance, and luminescence modes. The ultraviolet source is used to illuminate a composite surface of interest. In reflectance mode, the reflected response is acquired via the imaging system or via the spectrometer. The spectra are analyzed for organic compounds (conjugated organics) and inorganic compounds (semiconducting band-edge states; luminescing defect states such as silicates, used as adhesives for composite aerospace applications; and metal oxides commonly used as thermal coating paints on a wide range of spacecraft). The spectra are compared with a database for variation in conjugation, substitution, or length of molecule (in the case of organics) or band edge position (in the case of inorganics). This approach is useful in the understanding of material quality. It lacks the precision in defining the exact chemical structure that is found in other materials analysis techniques, but it is advantageous over methods such as nuclear magnetic resonance, infrared spectroscopy, and chromatography in that it can be used in the field to assess significant changes in chemical structure that may be linked to concerns associated with weaknesses or variations in structural integrity, without disassembly of or destruction to the structure of interest.
STUDY DESIGN: Brief (30-second) isometric trunk extensions at 5%, 20%, 40%, 60%, and 80% of maximal voluntary contraction (MVC) and 3 minutes of prolonged trunk extension (20% MVC) in erect position were studied in nine healthy male subjects. OBJECTIVES: To investigate the intercorrelation between intramuscular pressure and tissue oxygenation of the paravertebral muscles during submaximal isometric contractions and further, to evaluate paravertebral electromyogram and intramuscular pressure as indicators of force development. SUMMARY OF BACKGROUND DATA: Local physiologic responses to muscle contraction are incompletely understood. METHODS: Relative oxygenation was monitored with noninvasive near-infrared spectroscopy, intramuscular pressure was measured with a transducer-tipped catheter, and surface electromyogram was monitored at three recording sites. RESULTS: The root mean square amplitudes of the paravertebral electromyogram (L4, left and right; T12, right) and intramuscular pressure measured in the lumbar multifidus muscle at L4 increased with greater force development in a curvilinear manner. A significant decrease in the oxygenation of the lumbar paravertebral muscle in response to muscle contraction was found at an initial contraction level of 20% MVC. This corresponded to a paravertebral intramuscular pressure of 30-40 mm Hg. However, during prolonged trunk extension, no further decrease in tissue oxygenation was found compared with the tissue oxygenation level at the end of the brief contractions, indicating that homeostatic adjustments (mean blood pressure and heart rate) over time were sufficient to maintain paravertebral muscle oxygen levels. CONCLUSION: At a threshold intramuscular pressure of 30-40 mm Hg during muscle contraction, oxygenation in the paravertebral muscles is significantly reduced. The effect of further increase in intramuscular pressure on tissue oxygenation over time may be compensated for by an increase in blood pressure and heart rate. Surface electromyogram amplitudes and intramuscular pressure can be used as indicators of paravertebral muscle force.
It is now possible to grow silicon nanowires at chosen positions and orientations by a method that involves a combination of standard microfabrication processes. Because their positions and orientations can be chosen with unprecedented precision, the nanowires can be utilized as integral parts of individually electronically addressable devices in dense arrays. Nanowires made from silicon and perhaps other semiconductors hold substantial promise for integration into highly miniaturized sensors, field-effect transistors, optoelectronic devices, and other electronic devices. Like bulk semiconductors, inorganic semiconducting nanowires are characterized by electronic energy bandgaps that render them suitable as means of modulating or controlling electronic signals through electrostatic gating, in response to incident light, or in response to molecules of interest close to their surfaces. There is now potential for fabricating arrays of uniform, individually electronically addressable nanowires tailored to specific applications. The method involves formation of metal catalytic particles at the desired positions on a substrate, followed by heating the substrate in the presence of silane gas. The figure illustrates an example in which a substrate includes a silicon dioxide surface layer that has been etched into an array of pillars and the catalytic (in this case, gold) particles have been placed on the right-facing sides of the pillars. The catalytic thermal decomposition of the silane to silicon and hydrogen causes silicon columns (the desired nanowires) to grow outward from the originally catalyzed spots on the substrate, carrying the catalytic particles at their tips. Thus, the position and orientation of each silicon nanowire is determined by the position of its originally catalyzed spot on the substrate surface, and the orientation of the nanowire is perpendicular to the substrate surface at the originally catalyzed spot.
Thermal Infrared (TIR) remote sensing data can provide important measurements of surface energy fluxes and temperatures, which are integral to understanding landscape processes and responses. One example of this is the successful application of TIR remote sensing data to estimate evapotranspiration and soil moisture, where results from a number of studies suggest that satellite-based measurements from TIR remote sensing data can lead to more accurate regional-scale estimates of daily evapotranspiration. With further refinement in analytical techniques and models, the use of TIR data from airborne and satellite sensors could be very useful for parameterizing surface moisture conditions and developing better simulations of landscape energy exchange over a variety of conditions and space and time scales. Thus, TIR remote sensing data can significantly contribute to the observation, measurement, and analysis of energy balance characteristics (i.e., the fluxes and redistribution of thermal energy within and across the land surface) as an implicit and important aspect of landscape dynamics and landscape functioning. The application of TIR remote sensing data in landscape ecological studies has been limited, however, for several fundamental reasons that relate primarily to the perceived difficulty in use and availability of these data by the landscape ecology community, and from the fragmentation of references on TIR remote sensing throughout the scientific literature. It is our purpose here to provide evidence from work that has employed TIR remote sensing for analysis of landscape characteristics to illustrate how these data can provide important data for the improved measurement of landscape energy response and energy flux relationships. We examine the direct or indirect use of TIR remote sensing data to analyze landscape biophysical characteristics, thereby offering some insight on how these data can be used more robustly to further the understanding and modeling of landscape ecological processes.
In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases [1]. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code [2], are used with a blowing boundary condition in the DPLR hypersonic CFD code [3]. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver [4], is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.
The hardware configuration used in the present study of turbine blade planform friction damping, by means of the lumped parameter analysis, is the first turbine stage of the Space Shuttle Main Engine's High Pressure Fuel Turbopump. The analysis procedure solves the nonlinear equations of motion for a turbine blade that is acted on by a platform friction damper, using an iterative matrix method. Attention is given to the effects on blade deflection response of variations in friction coefficient, the normal force on the friction surface interface, blade hysteretic damping, the blade-to-blade phase angle of the harmonic forcing function, and the amplitude of the forcing function.
A pilot data center to compute ocean atmosphere heat exchange over the tropical ocean is prposed at the Jet Propulsion Laboratory (JPL) in response to the scientific needs of the Tropical Ocean and Global Atmosphere (TOGA) Program. Optimal methods will be used to estimate sea surface temperature (SET), surface wind speed, and humidity from spaceborne observations. A monthly summary of these parameters will be used to compute ocean atmosphere latent heat exchanges. Monthly fields of surface heat flux over tropical oceans will be constructed using estimations of latent heat exchanges and short wave radiation from satellite data. Verification of all satellite data sets with in situ measurements at a few locations will be provided. The data center will be an experimental active archive where the quality and quantity of data required for TOGA flux computation are managed. The center is essential to facilitate the construction of composite data sets from global measurements taken from different sensors on various satellites. It will provide efficient utilization and easy access to the large volume of satellite data available for studies of ocean atmosphere energy exchanges.