Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic strength”

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.

1,077 records · Page 60

Multidisciplinary Model Transformation Through Simplified Intermediate Representations

There has long been a challenge of making engineering tools from multiple disciplines interoperate. This problem extends to system modeling practices. This challenge has been confronted with a wide variety of techniques. These techniques include attempting to interface tools together into combined suites, attempting to find underlying commonalities in mathematics, supporting connections through semantic encoding, various graph mappings and transformations, and code wrappers. All of these approaches have strengths and weaknesses. These are measured in multiple areas: relative freedom of action of individual domain engineers in developing their own tools, speed of execution, ease of creation, traceability, fidelity of information transfer, and degree of alignment between the concepts of different domains. This paper presents an approach to this interoperation problem currently being used in the World-Wide Web. The approach is to develop easy-to-parse formats that allow flexibility to both the file author and file interpreter. Many of the formats that are currently deployed sacrifice runtime performance for the ability of third parties to easily understand what to do with the data. XML became popular earlier as a de-facto standard format for many web applications, but is now being replaced by JSON to enhance human readability and provide a simpler data model. This is the basis for work in this paper. Our approach, which provides the key to interoperation, is a simplified “shrapnel” intermediate collection of objects and relationships that is the result of a breakdown of the system model into minimal pieces. It is then reassembled on the destination side, forming a two-step transformation. Previous efforts with single-step transformations have proven too difficult to create efficiently. In contrast, the use of this approach leads to an almost automatic procedure for transformation development. The Europa project is a large engineering project that must coordinate the efforts of many different teams with different specialties. The traditional form of exchanging engineering information has been documentation. The vision of model-based systems engineering is to make this information exchange much more digital. This paper presents the application of our simplified format to connecting two different engineering tools to the system model, with a focus on a dynamic mission simulation encoded in Modelica.

Cole, Bjorn↗

Navier-Stokes flowfield computation of wing/rotor interaction for a tilt rotor aircraft in hover

The download on the wing produced by the rotor-induced downwash of a tilt rotor aircraft in hover is of major concern because of its severe impact on payload-carrying capability. A method has been developed to help gain a better understanding of the fundamental fluid dynamics that causes this download, and to help find ways to reduce it. In particular, the method is employed in this work to analyze the effect of a tangential leading edge circulation-control jet on download reduction. Because of the complexities associated with modeling the complete configuration, this work focuses specifically on the wing/rotor interaction of a tilt rotor aircraft in hover. The three-dimensional, unsteady, thin-layer compressible Navier-Stokes equations are solved using a time-accurate, implicit, finite difference scheme that employs LU-ADI factorization. The rotor is modeled as an actuator disk which imparts both a radical and an azimuthal distribution of pressure rise and swirl to the flowfield. A momentum theory blade element analysis of the rotor is incorporated into the Navier-Stokes solution method. Solution blanking at interior points of the mesh has been shown here to be an effective technique in introducing the effects of the rotor and tangential leading edge jet. Results are presented both for a rotor alone and for wing/rotor interaction. The overall mean characteristics of the rotor flowfield are computed including the flow acceleration through the rotor disk, the axial and swirl velocities in the rotor downwash, and the slipstream contraction. Many of the complex tilt rotor flow features are captured including the highly three-dimensional flow over the wing, the recirculation fountain at the plane of symmetry, wing leading and trailing edge separation, and the large region of separated flow beneath the wing. Mean wing surface pressures compare fairly well with available experimental data, but the time-averaged download/thrust ratio is 20-30 percent higher than the measured value. The discrepancy is due to a combination of factors that are discussed. Leading edge tangential blowing, of constant strength along the wing span, is shown to be effective in reducing download. The jet serves primarily to reduce the pressure on the wing upper surface. The computation clearly shows that, because of the three-dimensionality of the flowfield, optimum blowing would involve a spanwise variation in blowing strength.

Fejtek, Ian G.↗

Biomechanical Modeling of the Deadlift Exercise to Improve the Efficacy of Resistive Exercise Microgravity Countermeasures

During long-duration spaceflight missions, astronauts exposure to microgravity without adequate countermeasures can result in losses of muscular strength and endurance, as well as loss of bone mass. As a countermeasure to this challenge, astronauts engage in resistive exercise during spaceflight to maintain their musculoskeletal function. The Hybrid Ultimate Lifting Kit (HULK) has been designed as a prototype exercise device for an exploration-class vehicle; the HULK features a much smaller footprint than previous devices such as the Advanced Resistive Exercise Device (ARED) on the International Space Station (ISS), which makes the HULK suitable for extended spaceflight missions in vehicles with limited volume. As current ISS exercise countermeasure equipment represents an improvement over previous generations of such devices, the ARED is being employed as a benchmark of functional performance. This project involves the development of a biomechanical model of the deadlift exercise, and is novel in that it is the first exercise analyzed in this context to include the upper limbs in the loading path, in contrast to the squat, single-leg squat, and heel raise exercises also being modeled by our team. OpenSim software is employed to develop these biomechanical models of humans performing resistive exercises to assess and improve the new exercise device designs. Analyses include determining differences in joint and muscle forces when using different loading strategies with the device, comparing and contrasting with the ARED benchmark, and determining whether the loading is sufficient to maintain musculoskeletal health. During data collection, the number of repetitions, load, cadence, stance, and grip width are controlled in order to facilitate comparisons between loading configurations. To date, data have been collected for two human subjects performing the deadlift exercise on the HULK device using two different loading conditions. Recorded data include motion capture, electromyography (EMG), ground reaction forces, device load cell data, photos and videos, and anthropometric data. Work is ongoing to perform biomechanical analyses including inverse kinematics and inverse dynamics to compare different versions of the deadlift model in order to determine which provides an appropriate level of detail to study this exercise. This work is supported by the National Space Biomedical Research Institute through NCC 9-58.

physical exercise↗

Systematic Differences between Satellite-Based Presipitation Climatologies over the Tropical Oceans

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. It is clear however that there still remain significant uncertainties with satellite precipitation retrievals which limit their usefulness for many purposes. Systematic differences i'A tropical precipitation estimates have been brought to light in comparison activities such as the GPCP Algorithm Intercomparison Project and more recent Wetnet Precipitation Intercomparison Project 3. These uncertainties are assuming more importance because of the demands for validation associated with global climate modeling and data assimilation methodologies. The objective of the present study is to determine the physical basis for systematic differences in spatial structure of tropical precipitation as portrayed by several different satellite-based data sets. The study is limited to oceanic regions only and deals primarily with aspects of spatial variability. We are specifically interested in why MSU channel 1 and GPI precipitation differences are so striking over the Eastern Pacific ITCZ and why they both differ from other microwave emission-based precipitation estimates from SSM/I and a scattering-based deep convective ice index from MSU channel 2. Our results to date have shown that MSU channel I precipitation estimates are biased high over the Eastern Pacific ITCZ because of two factors: (1) the hypersensitivity of this frequency to cloud water in contrast to falling rain drops, and (2) unaccounted for scattering effects by precipitation-size ice which depresses the signal of the liquid water emission. Likewise, cold cloud top climatologies such as the GPI show an excess (a deficit) in estimated rainfall over the E. Pacific ITCZ (Warm Pool region). We show that these algorithms need to account for regionally varying heights (or temperatures) at which tropical convection detrains to form cirrus shields. A second objective we pursue is to identify variations in the macroscale cloud physical and thermodynamic properties of precipitation regimes" and relate these differences to tropical dynamical mechanisms of tropical heat and moisture balance. Finally, we interpret the algorithm differences and their associations with tropical dynamics in terms of WCRP GPCP goals for constructing precipitation climatologies.

Robertson, Frankin R.↗

The Broad-Band X-Ray Spectrum of IC 4329A from a Joint NuSTAR/Suzaku Observation

We have obtained a deep, simultaneous observation of the bright, nearby Seyfert galaxy IC 4329A with Suzaku and NuSTAR (Nuclear Spectroscopic Telescope Array). Through a detailed spectral analysis, we are able to robustly separate the continuum, absorption, and distant reflection components in the spectrum. The absorbing column is found to be modest (approximately 6 times 10 (sup 21) per square centimeter), and does not introduce any significant curvature in the Fe K band. We are able to place a strong constraint on the presence of a broadened Fe K alpha line (E (sub rest) equals 6.46 (from plus 0.08 to minus 0.07) kiloelectronvolts with sigma equals 0.33 (from plus 0.08 to minus 0.07) kiloelectronvolts and EW equals 34 (from plus 8 to minus 7) electronvolts), though we are not able to constrain any of the parameters of a relativistic reflection model. These results highlight the range in broad Fe K line strengths observed in nearby, bright, active galactic nuclei (roughly an order of magnitude), and imply a corresponding range in the physical properties of the inner accretion disk in these sources. We have also updated our previously reported measurement of the high-energy cutoff of the hard X-ray emission using both observatories rather than just NuSTAR alone: E (sub cut) equals 186 plus or minus 14 kiloelectronvolts. This high-energy cutoff acts as a proxy for the temperature of the coronal electron plasma, enabling us to further separate this parameter from the plasmas optical depth and to update our results for these parameters as well. We derive kT 50 (from plus 6 to minus 3) kiloelectronvolts with tau equals 2.34 (from plus +0.16 to minus 0.11) using a spherical geometry, kT equals 61 plus or minus1 kiloelectronvolt with tau equals 0.68 plus or minus 0.02 for a slab geometry, with both having an equivalent goodness-of-fit.

Dynamics - stars↗

Global Carbon Budget 2022

Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO 2 emissions (E FOS ) are based on energy statistics and cement production data, while emissions from land-use change (E LUC ), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO 2 concentration is measured directly, and its growth rate (G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink (S OCEAN ) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO 2 sink (S LAND ) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, E FOS increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr −1 (9.9 ± 0.5 GtC yr −1 when the cement carbonation sink is included), and E LUC was 1.1 ± 0.7 GtC yr −1 , for a total anthropogenic CO 2 emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr −1 (40.0 ± 2.9 GtCO 2 ). Also, for 2021, G ATM was 5.2 ± 0.2 GtC yr −1 (2.5 ± 0.1 ppm yr −1 ), S OCEAN was 2.9 ± 0.4 GtC yr −1 , and S LAND was 3.5 ± 0.9 GtC yr −1 , with a B IM of −0.6 GtC yr −1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO 2 concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in E FOS relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO 2 concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr −1 persist for the representation of annual to semi-decadal variability in CO 2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO 2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

Pierre Friedlingstein↗

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

Actuated Hybrid Mirror Telescope

The figure depicts the planned Actuated Hybrid Mirror Telescope (AHMT), which is intended to demonstrate a new approach to the design and construction of wide-aperture spaceborne telescopes for astronomy and Earth science. This technology is also appropriate for Earth-based telescopes. The new approach can be broadly summarized as using advanced lightweight mirrors that can be manufactured rapidly at relatively low cost. More specifically, it is planned to use precise replicated metallic nanolaminate mirrors to obtain the required high-quality optical finishes. Lightweight, dimensionally stable silicon carbide (SiC) structures will support the nanolaminate mirrors in the required surface figures. To enable diffraction- limited telescope performance, errors in surface figures will be corrected by use of mirror-shape-control actuators that will be energized, as needed, by a wave-front-sensing and control system. The concepts of nanolaminate materials and mirrors made from nanolaminate materials were discussed in several previous NASA Tech Briefs articles. Nanolaminates constitute a relatively new class of materials that can approach theoretical limits of stiffness and strength. Nanolaminate mirrors are synthesized by magnetron sputter deposition of metallic alloys and/or compounds on optically precise master surfaces to obtain optical-quality reflector surfaces backed by thin shell structures. As an integral part of the deposition process, a layer of gold that will constitute the reflective surface layer is deposited first, eliminating the need for a subsequent and separate reflective-coating process. The crystallographic textures of the nanolaminate will be controlled to optimize the performance of the mirror. The entire deposition process for making a nanolaminate mirror takes less than 100 hours, regardless of the mirror diameter. Each nanolaminate mirror will be bonded to its lightweight SiC supporting structure. The lightweight nanolaminate mirrors and SiC supporting structures will be fabricated from reusable master molds. The mirror-shape-control actuators will be low-power, high-capacitance lead magnesium niobate electrostrictive actuators that will be embedded in the SiC structures. The mode of operation of these actuators will be such that once power was applied, they will change in length and once power was removed, they will maintain dimensional stability to nanometer precision. This mode of operation will enable the use of low-power, minimally complex electronic control circuitry. The wave-front-sensing and control system will be designed and built according to a two-stage architecture. The first stage will be implemented by a Shack- Hartmann (SH) sensor subsystem, which will provide a large capture range. The second, higher-performance stage will be implemented by an image-based wave-front-sensing subsystem that will include a phase-retrieval camera (PRC), and will utilize phase retrieval and other techniques to measure wavefront error directly. Phase retrieval is a process in which multiple images of an unresolved object are iterated to estimate the phase of the optical system that acquired the images. The combination of SH and phase-retrieval sensors will afford the virtues of both a dynamic range of 105 and an accuracy of <10 nm.

Hickey, Gregory↗

Structural load control during construction

In the absence of gravitational pull, the major design considerations for large space structures are stiffness for controllability, and transient dynamic loadings (as opposed to the traditional static load associated with earth-based structures). Because of the absence of gravitational loading, space structures can be designed to be significantly lighter than their counterparts on Earth. For example, the Space Shuttle manipulator arm is capable of moving and positioning a 60,000 lb payload, yet weighs less than 1,000 lbs. A recent design for the Space Station which had a total weight of about 500,000 lbs. used a primary loadcarrying keel beam which weighed less than 10,000 lbs. For many large space structures designs it is quite common for the load-carrying structure to have a mass fraction on the order of one or two percent of the total spacecraft mass. This significant weight reduction for large space structures is commonly accompanied by very low natural frequencies. These low frequencies cause an unprecedented level of operational complexity for mission applications which require a high level of positioning and control accuracy. This control problem is currently the subject of considerable research directed towards reducing the flexibility problem. In addition, however, the small mass fraction typically results in structures which are quite unforgiving to inadvertent high loadings. In other words, the structures are 'fragile.' In order to deal with the fragility issue CSC developed a load-limiting concept for space truss structures. This concept is aimed at limiting the levels of load which can occur in a large space structure during the construction process as well as during subsequent operations. Currently, the approach for dealing with large loadings is to make the structure larger. The impact this has on construction is significant. The larger structures are more difficult to package in the launch vehicle, and in fact in some instances the concept must be changed from a deployable truss to an erectable truss to permit packaging. The new load-limiting concept is aimed at permitting the use in large space structures of smaller trusses with a high level of strength robustness, in order to simplify the construction process. To date several analyses conducted on the concept have demonstrated its feasibility, and an experiment is currently being designed to demonstrate its operation.

Mikulas, Martin M., Jr.↗

On the role of exchange of power and information signals in control and stability of the human-robot interaction

A human's ability to perform physical tasks is limited, not only by his intelligence, but by his physical strength. If, in an appropriate environment, a machine's mechanical power is closely integrated with a human arm's mechanical power under the control of the human intellect, the resulting system will be superior to a loosely integrated combination of a human and a fully automated robot. Therefore, we must develop a fundamental solution to the problem of 'extending' human mechanical power. The work presented here defines 'extenders' as a class of robot manipulators worn by humans to increase human mechanical strength, while the wearer's intellect remains the central control system for manipulating the extender. The human, in physical contact with the extender, exchanges power and information signals with the extender. The aim is to determine the fundamental building blocks of an intelligent controller, a controller which allows interaction between humans and a broad class of computer-controlled machines via simultaneous exchange of both power and information signals. The prevalent trend in automation has been to physically separate the human from the machine so the human must always send information signals via an intermediary device (e.g., joystick, pushbutton, light switch). Extenders, however are perfect examples of self-powered machines that are built and controlled for the optimal exchange of power and information signals with humans. The human wearing the extender is in physical contact with the machine, so power transfer is unavoidable and information signals from the human help to control the machine. Commands are transferred to the extender via the contact forces and the EMG signals between the wearer and the extender. The extender augments human motor ability without accepting any explicit commands: it accepts the EMG signals and the contact force between the person's arm and the extender, and the extender 'translates' them into a desired position. In this unique configuration, mechanical power transfer between the human and the extender occurs because the human is pushing against the extender. The extender transfers to the human's hand, in feedback fashion, a scaled-down version of the actual external load which the extender is manipulating. This natural feedback force on the human's hand allows him to 'feel' a modified version of the external forces on the extender. The information signals from the human (e.g., EMG signals) to the computer reflect human cognitive ability, and the power transfer between the human and the machine (e.g., physical interaction) reflects human physical ability. Thus the information transfer to the machine augments cognitive ability, and the power transfer augments motor ability. These two actions are coupled through the human cognitive/motor dynamic behavior. The goal is to derive the control rules for a class of computer-controlled machines that augment human physical and cognitive abilities in certain manipulative tasks.

Kazerooni, H.↗

Shock-Wave/Boundary-Layer Interaction (SWBLI) Experiments in the Presence of Transition-to-Turbulence on a Flat Plate Model in the NASA LaRC 31-Inch Mach 10 Air Tunnel

Historically, there have been a limited number of studies involving Transitional Shock-Wave/Boundary-Layer Interactions (XSWBLI) in hypersonic flows, the majority of which have been performed at Mach Numbers of 7.5 and below. Additionally, there have been even fewer such experiments which incorporate non-intrusive and optical diagnostics. NASA Langley Research Center (NASA LaRC), in collaboration with The University of Texas at San Antonio (UTSA) and The University of Tennessee Space Institute (UTSI), is currently planning and preparing for a series of experimental aerodynamic tests using the NASA LaRC 31-inch Mach 10 Air Tunnel. The experiments will primarily focus on studying the dynamics of Shock-Wave/Boundary-Layer Interactions (SWBLI) in the presence of transition (XSWBLI) and turbulence (SWTBLI) on a large, flat plate model at a freestream Mach number of 10. Besides continued interest in understanding the flow at high Mach numbers, executing such tests at Mach 10 avoids some known concerns and will aid in solving new problems. XSWBLI are a highly unsteady phenomena and the generation of XSWBLI at lower Mach numbers can prove to be a significant challenge. This concern will be somewhat mitigated in the Mach 10 flow, as boundary layers become increasingly stabilized, but some difficulty in achieving transition is still expected. Modelling SWBLI in the presence of transitional and turbulent boundary layers has also proven to be difficult, so a high-Mach number experiment which applies non-intrusive and optical diagnostics will aid in solving a unique problem as well as advancing the understanding and characterization of an aerodynamic surface at Mach 10.The NASA LaRC 31-inch Wind Tunnel incorporates a test section having a cross-section of 31-inches × 31-inches and provides optical access to the test section via three (3) ultra-violet (UV) transmitting windows on the top, side, and bottom. The Mach 10 operating conditions will consist of pressures (P0) ranging from 2.4 to 10MPa (348 to 1,450 psi) at a temperature (T0) of 1,000 K (1,800 °R). The NASA LaRC 31-inch wind tunnel will deliver 1-minute blow-down runs, of which 30-45 seconds will be consumed to reach the aforementioned test conditions. The pump down time between runs will be approximately 45-60 minutes, which will provide for 6-10 test runs per day over the course of a weeklong entry. The proposed flat plate model will be constructed of stainless steel with geometry consisting of a10-inch × 30-inch upper surface and a 2-inch thickness. The model will be designed to support multiple leading-edge inserts in order to examine blunt, round, and sharp leading edges at Mach 10. Adjustments will be made to manage boundary layer thickness as well as the strength and size of the shock interaction region for each test. Measurements and diagnostics will be performed though high-speed Schlieren, on-body high-speed pressure transducers, IR thermography, and oil flow visualization. Hypersonic vehicle applications have and will continue to emerge at the forefront of aerospace. The NASA-UTSA-UTSI team realizes the value of characterizing SWBLI, XSWBLI, and SWTBLI behaviors well beyond single-digit hypersonic Mach numbers. As hypersonic vehicle speed capabilities continue to increase, measurement and diagnostic methods for Mach numbers of 10 and beyond will provide the foundation for a firm understanding of flow field behavior which will directly influence the advancement of technologies towards the design and manufacturing of high-speed aerodynamic surfaces, controls, thermal protection systems, acoustic treatments, and structural components. The use of non-intrusive and optical diagnostic methods in such experimentation is pivotal in developing the visualization and empirical data necessary to advance the aforementioned technology areas.

Shockwave/Boundary-Layer Interactions↗

Free-to-Roll Testing of Airplane Models in Wind Tunnels

A free-to-roll (FTR) test technique and test rig make it possible to evaluate both the transonic performance and the wingdrop/ rock behavior of a high-strength airplane model in a single wind-tunnel entry. The free-to-roll test technique is a single degree-of-motion method in which the model is free to roll about the longitudinal axis. The rolling motion is observed, recorded, and analyzed to gain insight into wing-drop/rock behavior. Wing-drop/rock is one of several phenomena symptomatic of abrupt wing stall. FTR testing was developed as part of the NASA/Navy Abrupt Wing Stall Program, which was established for the purposes of understanding and preventing significant unexpected and uncommanded (thus, highly undesirable) lateral-directional motions associated with wing-drop/rock, which have been observed mostly in fighter airplanes under high-subsonic and transonic maneuvering conditions. Before FTR testing became available, wingrock/ drop behavior of high-performance airplanes undergoing development was not recognized until flight testing. FTR testing is a reliable means of detecting, and evaluating design modifications for reducing or preventing, very complex abrupt wing stall phenomena in a ground facility prior to flight testing. The FTR test rig was designed to replace an older sting attachment butt, such that a model with its force balance and support sting could freely rotate about the longitudinal axis. The rig (see figure) includes a rotary head supported in a stationary head with a forward spherical roller bearing and an aft needle bearing. Rotation is amplified by a set of gears and measured by a shaft-angle resolver; the roll angle can be resolved to within 0.067 degrees at a rotational speed up to 1,000 degrees/s. An assembly of electrically actuated brakes between the rotary and stationary heads can be used to hold the model against a rolling torque at a commanded roll angle. When static testing is required, a locking bar is used to fix the rotating head rigidly to the stationary head. Switching between the static and FTR test modes takes only about 30 minutes. The FTR test rig was originally mounted in a 16-ft (approximately 4.0-m) transonic wind tunnel, but could just as well be adapted to use in any large wind tunnel. In one series of tests on the FTR rig, static and dynamic characteristics of models of four different fighter airplanes were measured. Two of the models exhibited uncommanded lateral motions; the other two did not. A figure of merit was developed to discern the severity of lateral motions. Using this figure of merit, it was shown that the FTR test technique enabled identification of conditions under which the uncommanded lateral motions occurred. The wind-tunnel conditions thus identified were found to be correlated with flight conditions under which the corresponding full-size airplanes exhibited uncommanded lateral motions.

Capone, Francis J.↗

Muscle pain perception and sympathetic nerve activity to exercise during opioid modulation

The purpose of this experiment was to examine the effects of the endogenous opioid system on forearm muscle pain and muscle sympathetic nerve activity (MSNA) during dynamic fatiguing exercise. Twelve college-age men (24 +/- 4 yr) performed graded (1-min stages; 30 contractions/min) handgrip to fatigue 1 h after the ingestion of either 60 mg codeine, 50 mg naltrexone, or placebo. Pain (0-10 scale) and exertion (0-10 and 6-20 scales) intensities were measured during the last 15 s of each minute of exercise and every 15 s during recovery. MSNA was measured continuously from the peroneal nerve in the left leg. Pain threshold occurred earlier [1.8 +/- 1, 2. 2 +/- 1, 2.2 +/- 1 J: codeine, naltrexone, and placebo, respectively] and was associated with a lower rating of perceived exertion (RPE) (2.7 +/- 2, 3.6 +/- 2, 3.8 +/- 2: codeine, naltrexone, and placebo, respectively) in the codeine condition compared with either the naltrexone or placebo conditions. There were no main effects (i.e., drugs) or interaction (i.e., drugs x time) for either forearm muscle pain or RPE during exercise [pain: F (2, 22) = 0.69, P = 0.51]. There was no effect of drug on MSNA, heart rate, or blood pressure during baseline, exercise, or recovery. Peak exercise MSNA responses were 21 +/- 1, 21 +/- 2.0, and 21 +/- 2.0 bursts/30 s for codeine, naltrexone, and placebo conditions, respectively. Peak mean arterial pressure responses were 135 +/- 4, 131 +/- 3, and 132 +/- 4 mmHg for codeine, naltrexone, and placebo conditions, respectively. It is concluded that neither 60 mg codeine nor 50 mg naltrexone has an effect on forearm muscle pain, exertion, or MSNA during high- intensity handgrip to fatigue.

Clinical Trial↗

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗

Development of a Numerical Model of Hypervelocity Impact into a Pressurized Composite Overwrapped Pressure Vessel

As the outlook for space exploration becomes more ambitious and spacecraft travel deeper into space than ever before, it is increasingly important that propulsion systems perform reliably within the space environment. The increased reliability compels designers to increase design margin at the expense of system mass, which contrasts with the need to limit vehicle mass to maximize payload. Such are the factors that motivate the integration of high specific strength composite materials in the construction of pressure vessels commonly referred to as composite overwrapped pressure vessels (COPV). The COPV consists of a metallic liner for the inner shell of the COPV that is stiff, negates fluid permeation and serves as the anchor for composite laminates or filaments, but the liner itself cannot contain the stresses from the pressurant it contains. The compo-site-fiber reinforced polymer (CFRP) is wound around the liner using a combination of hoop (circumferential) and helical orientations. Careful consideration of wrap orientation allows the composite to evenly bear structural loading and creates the COPV's characteristic high strength to weight ratio. As the CFRP overwrap carries most of the stresses induced by pressurization, damage to the overwrap can affect mission duration, mission success and potentially cause loss-of-vehicle/loss-of-crew. For this reason, it is critical to establish a fundamental understanding of the mechanisms involved in the failure of a stressed composite such as that of the COPV. One of the greatest external threats to the integrity of a spacecraft's COPV is an impact from the meteoroid and orbital debris environments (MMOD). These impacts, even from submillimeter particles, generate extremely high stress states in the CFRP that can damage numerous fibers. As a result of this possibility, initial assumptions in survivability analysis for some human-rated NASA space-craft have assumed that any alteration of the vessel due to impact is considered a catastrophic failure. This assumption is conservative and made due to lack of knowledge on the level of allow-able damage to the composite overwrap that can be sustained and still allow successful completion of the mission. To quantify the allowable damage level to the composite overwrap involves assessing stress redistribution following damage as well as evaluating possible time-dependent mechanisms involved in the COPV response to an impact event. Limited published work in this subject has shown that COPV can withstand at least some level of damage due to high energy impacts. These observations have been confirmed and expanded upon in recent experimental research performed by NASA. This research has demonstrated that there is not only robustness in a COPV to compensate for CFRP damage, but has also identified two significant failure modes for pressurized COPV. The lowest threshold failure mode involves the perforation of the vessel, and the highest threshold failure mode is the catastrophic rupture. While both of these failure modes mean a loss of the COPV, system robustness affords some tolerance to the venting as opposed to the more catastrophic rupture. As a consequence, it is necessary to understand the conditions that result in the transition between these failure modes. The aforementioned experimental research has been performed in both the unpressurized and pressurized condition to identify the damage level that triggered the failure thresh-old. This COPV test program was sponsored by the NASA Engineering and Safety Center (NESC), and tests were performed at NASA White Sands Test Facility (WSTF). Planning and coordination were provided by NASA JSC Hypervelocity Impact Technology (HVIT) group, and the COPVs were provided by the ISS Program. Unpressurized testing has been conducted at the pressure of the vacuum test chamber, while, the pressurized testing has been conducted at 290 +/- 10 bar (4,200  100 psi) using nitrogen as the pressurizing gas, which corresponds to the design pressure for the target COPV. In this research, spherical aluminum projectiles with varying diameter has been chosen as the impactor. For the unpressurized COPV, the dependence of penetration up to the dependence of hole size in the liner has been obtained as a function of impact conditions. For the pressurized research, the dependence of penetration up to rupture has been obtained as a function of im-pact conditions. Two representative post-test photographs of the failed COPV's from a nor-mal impact into the COPV surface are shown in Fig. 1. These images display the dramatic difference between failure modes, venting (Fig. 1a) and rupture (Fig. 1b). For venting, liner perforation, severed composite fibers/tows and ply delamination are commonly observed damage characteristics of this COPV failure mode. In the case of rupture, the COPV typically experienced a separation of its domed regions and severe break-up of the cylindrical region. Fully understanding the transition from venting to rupture experimentally is costly and potentially unachievable for conditions that cannot be generated in the laboratory. These shortcomings have motivated the performance of three-dimensional numerical simulations to expand the existing experimental database. These simulations have been carried out with the nonlinear-structural-dynamics, analysis-tool, CTH. A typical pressure contour plot from an impact simulation of an entire COPV is shown in Fig. 2. To generate the COPV stress state without initiating a shock wave, the pressure in the simulated COPV is ramped up to the final pressure over a millisecond prior to impact of the projectile with nitrogen gas. Figure 2a shows the system in this initial condition. After one millisecond, a projectile is initiated into the simulation and impacts the COPV. Figure 2b shows the system after this impact. In the figure, the onset of venting is represented as the change in pressure (μbar), red to green, at the perforation site. Also seen in the figure is the eroded projectile that had passed into the COPV vessel with the generated shock wave in the pressurant propagating just ahead of the material. In this paper, pertinent experimental details and the development of the material constitutive models necessary for this work along with the efforts to validate their use are dis-cussed. The simulation results are presented and compared with the NASA experimental observations. While work is on-going from this effort, early observations pertinent to the failure threshold are presented.

Garcia, M. A.↗