Search NASA⌕ Search

SEARCH · Search NASA

Results for “deformation function”

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.

At least 613 records · Page 34

Damage-mitigating control of a reusable rocket engine for high performance and extended life

The goal of damage mitigating control in reusable rocket engines is to achieve high performance with increased durability of mechanical structures such that functional lives of the critical components are increased. The major benefit is an increase in structural durability with no significant loss of performance. This report investigates the feasibility of damage mitigating control of reusable rocket engines. Phenomenological models of creep and thermo-mechanical fatigue damage have been formulated in the state-variable setting such that these models can be combined with the plant model of a reusable rocket engine, such as the Space Shuttle Main Engine (SSME), for synthesizing an optimal control policy. Specifically, a creep damage model of the main thrust chamber wall is analytically derived based on the theories of sandwich beam and viscoplasticity. This model characterizes progressive bulging-out and incremental thinning of the coolant channel ligament leading to its eventual failure by tensile rupture. The objective is to generate a closed form solution of the wall thin-out phenomenon in real time where the ligament geometry is continuously updated to account for the resulting deformation. The results are in agreement with those obtained from the finite element analyses and experimental observation for both Oxygen Free High Conductivity (OFHC) copper and a copper-zerconium-silver alloy called NARloy-Z. Due to its computational efficiency, this damage model is suitable for on-line applications of life prediction and damage mitigating control, and also permits parametric studies for off-line synthesis of damage mitigating control systems. The results are presented to demonstrate the potential of life extension of reusable rocket engines via damage mitigating control. The control system has also been simulated on a testbed to observe how the damage at different critical points can be traded off without any significant loss of engine performance. The research work reported here is built upon concepts derived from the disciplines of Controls, Thermo-fluids, Structures, and Materials. The concept of damage mitigation, as presented in this report, is not restricted to control of rocket engines. It can be applied to any system where structural durability is an important issue.

Ray, Asok↗

A Probabilistic Mass Estimation Algorithm for a Novel 7- Channel Capacitive Sample Verification Sensor

A document describes an algorithm created to estimate the mass placed on a sample verification sensor (SVS) designed for lunar or planetary robotic sample return missions. A novel SVS measures the capacitance between a rigid bottom plate and an elastic top membrane in seven locations. As additional sample material (soil and/or small rocks) is placed on the top membrane, the deformation of the membrane increases the capacitance. The mass estimation algorithm addresses both the calibration of each SVS channel, and also addresses how to combine the capacitances read from each of the seven channels into a single mass estimate. The probabilistic approach combines the channels according to the variance observed during the training phase, and provides not only the mass estimate, but also a value for the certainty of the estimate. SVS capacitance data is collected for known masses under a wide variety of possible loading scenarios, though in all cases, the distribution of sample within the canister is expected to be approximately uniform. A capacitance-vs-mass curve is fitted to this data, and is subsequently used to determine the mass estimate for the single channel s capacitance reading during the measurement phase. This results in seven different mass estimates, one for each SVS channel. Moreover, the variance of the calibration data is used to place a Gaussian probability distribution function (pdf) around this mass estimate. To blend these seven estimates, the seven pdfs are combined into a single Gaussian distribution function, providing the final mean and variance of the estimate. This blending technique essentially takes the final estimate as an average of the estimates of the seven channels, weighted by the inverse of the channel s variance.

Wolf, Michael↗

The Extensional Rheology of Non-Newtonian Materials

The evolution of the transient extensional stresses in dilute and semi-dilute viscoelastic polymer solutions are measured with a filament stretching rheometer of a design similar to that first introduced by Sridhar, et al. The solutions are polystyrene-based (PS) Boger fluids that are stretched at constant strain rates ranging from 0.6 less than or equal to epsilon(0) less than or equal to 4s(exp -1) and to Hencky strains of epsilon greater than 4. The test fluids all strain harden and Trouton ratios exceeding 1000 are obtained at high strains. The experimental data strain hardens at lower strain levels than predicted by bead-spring FENE models. In addition to measuring the transient tensile stress growth, we also monitor the decay of the tensile viscoelastic stress difference in the fluid column following cessation of uniaxial elongation as a function of the total imposed Hencky strain and the strain rate. The extensional stresses initially decay very rapidly upon cessation of uniaxial elongation followed by a slower viscoelastic relaxation, and deviate significantly from FENE relaxation predictions. The relaxation at long times t is greater than or equal to 5 s, is compromised by gravitational draining leading to non-uniform filament profiles. For the most elastic fluids, partial decohension of the fluid filament from the endplates of the rheometer is observed in tests conducted at high strain rates. This elastic instability is initiated near the rigid endplate fixtures of the device and it results in the progressive breakup of the fluid column into individual threads or 'fibrils' with a regular azimuthal spacing. These fibrils elongate and bifurcate as the fluid sample is elongated further. Flow visualization experiments using a modified stretching device show that the instability develops as a consequence of an axisymmetry-breaking meniscus instability in the nonhomogeneous region of highly deformed fluid near the rigid endplate.

Spiegelberg, Stephen H.↗

Buckling and postbuckling behavior of square compression-loaded graphite-epoxy plates with circular cutouts

Results are presented for unidirectional (0, 10)(sub s) and (90,10)(sub s) plates, ((0/90)(sub 5)(sub s)) plates, and for aluminum plates. Results are also presented for ((+/- theta)(sub 6)(sub s)) angle-ply plates for values of theta = 30, 45, and 60 degrees. The results indicate that the change in axial stiffness of a plate at buckling is strongly dependent upon cutout size and plate orthotropy. The presence of a cutout gives rise to an internal load distribution that changes, sometimes dramatically, as a function of cutout size coupled with the plate orthotropy. In the buckled state, the role of orthotropy becomes more significant since bending in addition to membrane orthotropy is present. Most of the plates with cutouts exhibited less postbuckling stiffness than the corresponding plate without a cutout, and the postbuckling stiffness decreased with increasing cutout size. However, some of the highly orthotropic plates with cutouts exhibited more postbuckling stiffness than the corresponding plate without a cutout. These results suggest the possibility of tailoring the cutout size and the stacking sequence of a composite plate to optimize postbuckling stiffness. It was found that plates with large radius cutouts do exhibit some postbuckling strength. The results also indicate that a cutout can influence modal interaction in a plate. Specifically, results are presented that show a plate with a relatively small cutout buckling at a higher load than the corresponding plate without a cutout, due to modal interaction. Other results are presented that indicate the presence of nonlinear prebuckling deformations, due to material nonlinearity, in the angle-ply plates with theta = 45 and 60 degrees. The nonlinear prebuckling deformations are more pronounced in the plates with theta = 45 degrees and become even more pronounced as the cutout size increases. Results are also presented that show how load-path eccentricity due to improper machining of the test specimens affects the buckling behavior. Some of the plates with cutouts and eccentricity exhibited a snap-through type of buckling behavior.

Nemeth, Michael P.↗

Embossed Teflon AF Laminate Membrane Microfluidic Diaphragm Valves

A microfluidic system has been designed to survive spaceflight and to function autonomously on the Martian surface. It manipulates microscopic quantities of liquid water and performs chemical analyses on these samples to assay for the presence of molecules associated with past or present living processes. This technology lies at the core of the Urey Instrument, which is scheduled for inclusion on the Pasteur Payload of the ESA ExoMars rover mission in 2013. Fabrication processes have been developed to make the microfabricated Teflon-AF microfluidic diaphragm pumps capable of surviving extreme temperature excursions before and after exposure to liquid water. Two glass wafers are etched with features and a continuous Teflon membrane is sandwiched between them (see figure). Single valves are constructed using this geometry. The microfabricated devices are then post processed by heating the assembled device while applying pneumatic pressure to force the Teflon diaphragm against the valve seat while it is softened. After cooling the device, the embossed membrane retains this new shape. This solves previous problems with bubble introduction into the fluid flow where deformations of the membrane at the valve seat occurred during device bonding at elevated temperatures (100-150 C). The use of laminated membranes containing commercial Teflon AF 2400 sheet sandwiched between spun Teflon AF 1600 layers performed best, and were less gas permeable than Teflon AF 1600 membranes on their own. Spinning Teflon AF 1600 solution (6 percent in FLOURINERT(Registered TradeMark) FC40 solvent, 3M Company) at 500 rpm for 1.5 seconds, followed by 1,000 rpm for 3 seconds onto Borofloat glass wafers, results in a 10-micron-thick film of extremely smooth Teflon AF. This spinning process is repeated several times on flat, blank, glass wafers in order to gradually build a thick, smooth membrane. After running this process at least five times, the wafer and Teflon coating are heated under vacuum at 220 C for one hour in order to drive off any residual solvent present in the composite film. After this, a second blank, glass wafer is brought down from above and the stack is held under vacuum at 3 atm mechanical pressure for ten 10 hours.

Willis, Peter↗

An Improved Plastically Dilatant Unified Viscoplastic Constitutive Formulation for Multiscale Analysis of Polymer Matrix Composites Under High Strain Rate Loading

Polymer matrix composites are commonly used to fabricate energy-absorbing structures expected to experience impact loading. As such, a detailed understanding of the dynamic response of the constituent materials is necessary. Since the rate, temperature, and pressure dependence of carbon fiber reinforced polymer matrix composites are primarily manifestations of the rate, temperature, and pressure dependence of the polymer matrix, it is crucial that the constitutive behavior of the matrix be accurately characterized. In this work, an existing unified viscoplastic constitutive formulation is extended to ensure thermodynamic consistency and to more accurately account for the tension-compression asymmetry observed in the response of polymeric materials. A new plastic potential function is proposed, and elementary loading conditions are utilized to determine relations between model constants to ensure nonnegative plastic dissipation, a necessary thermodynamic requirement. Expressions for plastic Poisson’s ratios are derived and are bounded by enforcing nonnegative plastic dissipation. The model is calibrated against available experimental data from tests conducted over a range of strain rates, temperatures, and loading cases on a representative thermoset epoxy; good correlation between simulations and experimental data is obtained. Temperature rises due to the conversion of plastic work to heat are computed via the adiabatic heat energy equation. The viscoplastic polymer model is then used as a constitutive model in the generalized method of cells micromechanics theory to investigate the effects of matrix adiabatic heating on the high strain rate response of a unidirectional composite. The thermodynamic consistency of the model ensures plastic dissipation can only cause an increase in temperature. Simulation results indicate that significant thermal softening due to the conversion of plastic work to heat is observed in the composite for matrix dominated deformation modes.

Impact loading↗

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LS-DYNA (Registered), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic material model with a non-associative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Plasticity↗

Detection of atmospheric pressure loading using very long baseline interferometry measurements

Loading of the Earth by the temporal redistribution of global atmospheric mass is likely to displace the positions of geodetic monuments by tens of millimeters both vertically and horizontally. Estimates of these displacements are determined by convolving National Meteorological Center (NMC) global values of atmospheric surface pressure with Farrell's elastic Green's functions. An analysis of the distances between radio telescopes determined by very long baseline interferometry (VLBI) between 1984 and 1992 reveals that in many of the cases studied there is a significant contribution to baseline length change due to atmospheric pressure loading. Our analysis covers intersite distances of between 1000 and 10,000 km and is restricted to those baselines measured more than 100 times. Accounting for the load effects (after first removing a best fit slope) reduces the weighted root-mean-square (WRMS) scatter of the baseline length residuals on 11 of the 22 baselines investigated. The slight degradation observed in the WRMS scatter on the remaining baselines is largely consistent with the expected statistical fluctuations when a small correction is applied to a data set having a much larger random noise. The results from all baselines are consistent with approximately 60% of the computed pressure contribution being present in the VLBI length determinations. Site dependent coefficients determined by fitting local pressure to the theoretical radial displacement are found to reproduce the deformation caused by the regional pressure to within 25% for most inland sites. The coefficients are less reliable at near coastal and island stations.

Vandam, T. M.↗

Reflector Surface Error Compensation in Dual-Reflector Antennas

By probing the field on a small subreflector at a minimal number of points, the main reflector surface errors can be obtained and subsequently used to design a phase-correction subreflector that can compensate for main reflector errors. The compensating phase-error profile across the subreflector can be achieved either by a surface deformation or by the use of an array of elements such as patch antennas that can cause a phase shift between the incoming and outgoing fields. The second option is of primary interest here, but the methodology can be applied to either case. The patch array is most easily implemented on a planar surface. Therefore, the example of a flat subreflector and a parabolic main reflector (a Newtonian dual reflector system) is considered in this work. The subreflector is assumed to be a reflector array covered with patch elements. The phase variation on a subreflector can be detected by a small number of receiving patch elements (probes). By probing the phase change at these few selected positions on the subreflector, the phase error over the entire surface can be recovered and used to change the phase of all the patch elements covering the subreflector plane to compensate for main reflector errors. This is accomplished by using a version of sampling theorem on the circular aperture. The sampling is performed on the phase-error function on the circular aperture of the main reflector by a method developed using Zernike polynomials. This method is based upon and extended from a theory previously proposed and applied to reflector aperture integration. This sampling method provides for an exact retrieval of the coefficients of up to certain orders in the expansion of the phase function, from values on a specifically calculated set of points in radial and azimuthal directions in the polar coordinate system, on the circular reflector aperture. The corresponding points on the subreflector are then obtained and, by probing the fields at these points, a set of phase values is determined that is then transferred back to the main reflector aperture for recovering the phase function. Once this function is recovered, the corresponding phase function on the subreflector is calculated and used to compensate for main reflector surface errors. In going back and forth between sub and main reflectors, geometrical (ray) optics is employed, which even though it ignores edge diffraction and other effects, is shown to be accurate for phase recovery.

Jamnejad, Vahraz↗

Functionality Improvements to Overaero

The functionality of the overset, static aeroelasticity, Navier-Stokes flow solver OVERAERO was increased by adding capability to the flow solver and enhancing code performance. Improvements were made to the fluids/structure interface, an MLP version of the parallel OVERAERO code was developed, and the OVERAERO-MPI code was ported to the Cray T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested successfully on the IPG testbed and a means of reducing communication overhead within OVERFLOW-MPI was investigated. To solve an aeroelastic problem computationally, a structures grid surface definition and a fluids grid surface definition are required. Typically, the structures grid surface has a lower fidelity than the fluids grid surface. Thus, the methods developed to transfer data between the two grid systems are vital to the accuracy and efficiency of the aeroelasticity code. The fluids/structures interface developed for the OVERAERO code was improved to more accurately treat fluids surfaces that bridge between two different structural surfaces. For example, the method allowed the forward portion of a flap track fairing to deform with the wing and the aft end of the fairing to deform with the flap. A tightly-coupled version of the code based on OVERFLOW-MLP was developed to improve code performance on the SGI Origin 2000. This required a new parallelization strategy to couple the fluids and structures codes. The OVERAERO-MPI code was ported to the Cray T3E to extend the usability of the code. The port required extensive use of dynamic memory management techniques to fit large problems within the memory limitations of the T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested on the IPG testbed being developed within NASA. For small problems with minimal data transfer between grids, there was little to no performance penalty spreading the computation across two machines. For very large problems, methods were developed to minimize intermachine communication via the grid partitioning scheme. By minimizing the intermachine communication requirements of the problem, it may still be beneficial to run a tightly-coupled flow solver across two machines within the IPG.

Gee, Ken↗

Using Digital Image Correlation for Material Testing at Cryogenic Temperatures

Recent efforts were made to parametrize an aluminum alloy at cryogenic temperatures to calibrate finite element models simulating material deformation. Digital image correlation (DIC) was used to measure material properties using front/back strains of specimens during tensile tests. Data acquired from DIC was compared to physical extensometer measurements. Two camera pairs (front and back) were mounted via aluminum extrusion to the exterior of an environmental chamber of a hydraulic test stand. Images were taken through single-pane optical glass windows of the environment chamber. Positive pressure bags encased the cameras, attaching to the perimeter of the glass to keep out moisture and prevent windows from frosting. Figure 1 shows the environmental setup and the camera setup during cryogenic testing. Some of the challenges during testing at cryogenic temperatures include the gaseous flow of liquid nitrogen in the optical path of the mounted cameras, frost buildup on the windows as seen in Figure 1b, and frost on the test coupon. The temperature gradient present during the introduction of liquid nitrogen into the chamber can increase the projection error of images captured of the specimen. Figure 2 shows the estimated projection error from DIC as a function of temperature during the cool down process. This presentation addresses some of the concerns and observations when using DIC for testing at cryogenic temperatures using methods of isolated camera system and techniques for maintaining the glass panel clear from frost.

Joseph Cochran↗

Acceleration effects in solid propellant rocket motors

The performance variations due to acceleration loads imposed on spinning solid propellant rocket motors are investigated. The four potentially most significant modes of acceleration-induced phenomena are identified from a study of the literature and modeled. The four modes are a mechanical mode which deals with deformations of the propellant and case: a thermodynamic mode which covers acceleration-induced combustion phenomena; a stress mode which covers the stressed propellant's effect on burn rate; and a gas dynamic mode which deals with changes in gas flow in the chamber and through the nozzle. Simplified models of each mode are developed or taken from the literature and are added to an internal ballistics evaluation computer program. The resulting analysis is the first to include all of the modes. In order to do this an original analysis of the mechanical and stress modes was necessary. However, the analysis shows that the stress mode is not important for the circular perforated grains studied. The other effects are shown to have a significant influence on solid rocket motor performance. The magnitude of the different mode effects are such that one may not be ignored over the others as has been done in the past. The results of the analysis are compared to published rocket motor data. The comparisons indicate an erosive burning effect that is a function of spin rate. A qualitative explanation of the erosive effect is presented.

Langhenry, M. T.↗

Coronagraph Design Optimization for Segmented Aperture Telescopes

The goal of directly imaging Earth-like planets in the habitable zone of other stars has motivated the design of coronagraphs for use with large segmented aperture space telescopes. In order to achieve an optimal trade-o between planet light throughput and di racted starlight suppression, we consider coronagraphs comprised of a stage of phase control implemented with deformable mirrors (or other optical elements), pupil plane apodization masks (gray scale or complex valued), and focal plane masks (either amplitude only or complex-valued, including phase only such as the vector vortex coronagraph). The optimization of these optical elements, with the goal of achieving 10 or more orders of magnitude in the suppression of on-axis (starlight) di racted light, represents a challenging non-convex optimization problem with a nonlinear dependence on control degrees of freedom. We develop a new algorithmic approach to the design optimization problem, which we call the "Auxiliary Field Optimization" (AFO) algorithm. The central idea of the algorithm is to embed the original optimization problem, for either phase or amplitude (apodization) in various planes of the coronagraph, into a problem containing additional degrees of freedom, speci cally ctitious "auxiliary" electric elds which serve as targets to inform the variation of our phase or amplitude parameters leading to good feasible designs. We present the algorithm, discuss details of its numerical implementation, and prove convergence to local minima of the objective function (here taken to be the intensity of the on-axis source in a "dark hole" region in the science focal plane). Finally, we present results showing application of the algorithm to both unobscured o -axis and obscured on-axis segmented telescope aperture designs. The application of the AFO algorithm to the coronagraph design problem has produced solutions which are capable of directly imaging planets in the habitable zone, provided end-to-end telescope system stability requirements can be met. Ongoing work includes advances of the AFO algorithm reported here to design in additional robustness to a resolved star, and other phase or amplitude aberrations to be encountered in a real segmented aperture space telescope.

Redding, Dave↗

Observations on the long-period variability of the Gulf Stream downstream of Cape Hatteras

To examine the long-period variability of the Gulf Stream, sea level residuals relative to a 2-year mean sea level in the Gulf Stream downstream of Cape Hatteras (between 75 deg W and 60 deg W longitude) are used. Residuals, as derived from Geosat altimetry between November 1986 and December 1988, were gridded in space and time at a temporal resolution of 10 days and spatial resolution of 1/4 deg. Complex empirical orthogonal function (CEOF) analysis was applied to the data set to extract the spatially correlated signal with the original data subsampled to 1/2 deg. In addition to determining the space-time scales and propagation characterisitics of the different modes, wavenumber-frequency spectral techniques were used to separate the variability into propagating and stationary components. The CEOF technique applied to the data set indicated that the first four CEOF modes accounted for 60% of the variability and were found to be above the noise leve 99% of the time. CEOF 1 was associated with westward propagation at 5 km/d at a wavelength of 2000 km and eastward propagation at 1-2 km/d centered at a 500-km wavelength. This first CEOF is in good agreement with thin-jet equivalent barotropic models which predict westward propagation for wavelengths greater than 1130 km. A deflection of the wavelike pattern at 65 deg W also indicates a possible topographic effect. A simple scaling of the effect of topography indicates that for length scales longer than the internal Rossby radius of deformation, the topographic term is at least of the same order of magnitude as the beta effect. The second CEOF was more broadbanded in wavenumber space, with eastward propagation occurring in a wavenumber-frequency band between 300 and 1400 km and 0.5 and 2.0 cycles/yr. The third CEOF is similar in structure to the first, but with less energy. CEOF 4 was clearly identifiable with higher frequencies than the first three with westward propagation at 4 km/d. The spatial location of this mode along with the westward propagation indicates possible influences from eddy-stream interactions. Thus topography, Rossby wave dynamics and eddy-stream interactions all appear to have a significant role in determining the space-time scales and propagation properties of the long-period response of sea level in the Gulf Stream.

Vazquez, Jorge↗

MOD Tool (Microwave Optics Design Tool)

The Jet Propulsion Laboratory (JPL) is currently designing and building a number of instruments that operate in the microwave and millimeter-wave bands. These include MIRO (Microwave Instrument for the Rosetta Orbiter), MLS (Microwave Limb Sounder), and IMAS (Integrated Multispectral Atmospheric Sounder). These instruments must be designed and built to meet key design criteria (e.g., beamwidth, gain, pointing) obtained from the scientific goals for the instrument. These criteria are frequently functions of the operating environment (both thermal and mechanical). To design and build instruments which meet these criteria, it is essential to be able to model the instrument in its environments. Currently, a number of modeling tools exist. Commonly used tools at JPL include: FEMAP (meshing), NASTRAN (structural modeling), TRASYS and SINDA (thermal modeling), MACOS/IMOS (optical modeling), and POPO (physical optics modeling). Each of these tools is used by an analyst, who models the instrument in one discipline. The analyst then provides the results of this modeling to another analyst, who continues the overall modeling in another discipline. There is a large reengineering task in place at JPL to automate and speed-up the structural and thermal modeling disciplines, which does not include MOD Tool. The focus of MOD Tool (and of this paper) is in the fields unique to microwave and millimeter-wave instrument design. These include initial design and analysis of the instrument without thermal or structural loads, the automation of the transfer of this design to a high-end CAD tool, and the analysis of the structurally deformed instrument (due to structural and/or thermal loads). MOD Tool is a distributed tool, with a database of design information residing on a server, physical optics analysis being performed on a variety of supercomputer platforms, and a graphical user interface (GUI) residing on the user's desktop computer. The MOD Tool client is being developed using Tcl/Tk, which allows the user to work on a choice of platforms (PC, Mac, or Unix) after downloading the Tcl/Tk binary, which is readily available on the web. The MOD Tool server is written using Expect, and it resides on a Sun workstation. Client/server communications are performed over a socket, where upon a connection from a client to the server, the server spawns a child which is be dedicated to communicating with that client. The server communicates with other machines, such as supercomputers using Expect with the username and password being provided by the user on the client.

Katz, Daniel S.↗

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗

Reduction of Effective Acceleration to Microgravity Levels

Acceleration due to earth's gravity causes buoyancy driven convection and sedimentation in solutions. In addition. pressure gradients occur as a function of the height within a liquid column. Hence gravity effects both equilbria conditions and phase transitions as a result of hydrostatic pressure gradients. The affect of gravity on the rate of heat and man transfer in solutal processes can be particularly important in polymer processing due to the high sensitivity of polymeric materials to processing conditions. The term microgravity has been coined to describe an environment in which the affects of gravitational acceleration am greatly reduced. It may seem odd to talk in term of reducing the effects of gravitational acceleration since gravitational attraction is a basic property of matter. However, die presence of gravity on in situ processing or measurements can be negated by achieving conditions in which the laboratory, or more specifically the container of the experimental materials, a subjected to the same acceleration as the materials themselves. With regard to the laboratory reference frame, there is virtually no force on the experimental solutions. This is difficult to achieve but can be done. A short review of Newtonian physics provides an explanation on both how processes we affected by gravity and how microgravity conditions are achieved. The fact that fluids deform when subject to a force bid solids do not indicates that solids have a structure able to exert an opposing force that negates an externally applied force. Liquids deform when a force is applied, indicating that a liquid structure cannot completely negate an applied force. Just how easily a liquid resists deformation is related to its viscosity. Spaceflight provides an environment in which the laboratory reference frame i.e. the spacecraft and all the equipment therein an experiencing virtually identical forces. There is no solid foundation underneath such a laboratory, so the laboratory accelerates according to the force of gravity as do the experimental fluids within the lab. Hence, the magnitude of the form excited by the laboratory on the experimental solutions within are greatly reduced. When compared with a laboratory on the ground and averaged over time, the fluids in a spaceflight laboratory experience approximately a 10 (sup -6)decrease in acceleration relative to their laboratory reference frame hence the term microgravity.

Downey, James P.↗

Animal models and their importance to human physiological responses in microgravity

Two prominent theories to explain the physiological effects of microgravity relate to the cascade of changes associated with the cephalic shifts of fluids and the absence of tissue deformation forces. One-g experiments for humans used bed rest and the head-down tilt (HDT) method, while animal experiments have been conducted using the tail-suspended, head-down, and hindlimbs non-weightbearing model. Because of the success of the HDT approach with rats to simulate the gravitational effects on the musculoskeletal system exhibited by humans, the same model has been used to study the effects of gravity on the cardiopulmonary systems of humans and other vertebrates. Results to date indicate the model is effective in producing comparable changes associated with blood volume, erythropoiesis, cardiac mass, baroreceptor responsiveness, carbohydrate metabolism, post-flight VO2max, and post-flight cardiac output during exercise. Inherent with these results is the potential of the model to be useful in investigating responsible mechanisms. The suspension model has promise in understanding the capillary blood PO2 changes in space as well as the arterial PO2 changes in subjects participating in a HDT experiment. However, whether the model can provide insights on the up-or-down regulation of adrenoreceptors remains to be determined, and many investigators believe the HDT approach should not be followed to study gravitational influences on pulmonary function in either humans or animals. It was concluded that the tail-suspended animal model had sufficient merit to study in-flight and post-flight human physiological responses and mechanisms.

NASA Discipline Cardiopulmonary↗