Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS COMPUTER PROGRAM”

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 685 records · Page 38

Intracellular stress tomography reveals stress focusing and structural anisotropy in cytoskeleton of living cells

We describe a novel synchronous detection approach to map the transmission of mechanical stresses within the cytoplasm of an adherent cell. Using fluorescent protein-labeled mitochondria or cytoskeletal components as fiducial markers, we measured displacements and computed stresses in the cytoskeleton of a living cell plated on extracellular matrix molecules that arise in response to a small, external localized oscillatory load applied to transmembrane receptors on the apical cell surface. Induced synchronous displacements, stresses, and phase lags were found to be concentrated at sites quite remote from the localized load and were modulated by the preexisting tensile stress (prestress) in the cytoskeleton. Stresses applied at the apical surface also resulted in displacements of focal adhesion sites at the cell base. Cytoskeletal anisotropy was revealed by differential phase lags in X vs. Y directions. Displacements and stresses in the cytoskeleton of a cell plated on poly-L-lysine decayed quickly and were not concentrated at remote sites. These data indicate that mechanical forces are transferred across discrete cytoskeletal elements over long distances through the cytoplasm in the living adherent cell.

NASA Discipline Cell Biology↗

Rescue Shuttle Flight Re-Entry: Controlling Astronaut Thermal Exposure

A rescue mission for the STS-125 Hubble Telescope Repair Mission requires reentry from space with 11 crew members aboard, exceeding past cabin thermal load experience and risking crew thermal stress potentially causing cognitive performance and physiological decrements. The space shuttle crew cabin air revitalization system (ARS) was designed to support a nominal crew complement of 4 to 7 crew and 10 persons in emergencies, all in a shirt-sleeve environment. Subsequent to the addition of full pressure suits with individual cooling units, the ARS cannot maintain a stable temperature in the crew cabin during reentry thermal loads. Bulk cabin thermal models, used for rescue mission planning and analysis of crew cabin air, were unable to accurately represent crew workstation values of air flow, carbon dioxide, and heat content for the middeck. Crew temperature models suggested significantly elevated core temperatures. Planning for an STS-400 potential rescue of seven stranded crew utilized computational fluid dynamics (CFD) models to demonstrate inhomogeneous cabin thermal properties and improve analysis compared to bulk models. In the absence of monitoring of crew temperature, heart rate, metabolic rate and incomplete engineering data on the performance of the integrated cooling garment/cooling unit (ICG/CU) at cabin temperatures above 75 degrees F, related systems & models were reevaluated and tests conducted with humans in the loop. Changes to the cabin ventilation, ICU placement, crew reentry suit-donning procedures, Orbiter Program wave-off policy and post-landing power down and crew extraction were adopted. A second CFD and core temperature model incorporated the proposed changes and confirmed satisfactory cabin temperature, improved air distribution, and estimated core temperatures within safe limits. CONCLUSIONS: These changes in equipment, in-flight and post-landing procedures, and policy were implemented for the STS-400 rescue shuttle & will be implemented in any future rescue flights from the International Space Station of stranded shuttle crews.

Gillis, David B.↗

On Shaft Data Acquisition System (OSDAS)

On Shaft Data Acquisition System (OSDAS) is a rugged, compact, multiple-channel data acquisition computer system that is designed to record data from instrumentation while operating under extreme rotational centrifugal or gravitational acceleration forces. This system, which was developed for the Heritage Fuel Air Turbine Test (HFATT) program, addresses the problem of recording multiple channels of high-sample-rate data on most any rotating test article by mounting the entire acquisition computer onboard with the turbine test article. With the limited availability of slip ring wires for power and communication, OSDAS utilizes its own resources to provide independent power and amplification for each instrument. Since OSDAS utilizes standard PC technology as well as shared code interfaces with the next-generation, real-time health monitoring system (SPARTAA Scalable Parallel Architecture for Real Time Analysis and Acquisition), this system could be expanded beyond its current capabilities, such as providing advanced health monitoring capabilities for the test article. High-conductor-count slip rings are expensive to purchase and maintain, yet only provide a limited number of conductors for routing instrumentation off the article and to a stationary data acquisition system. In addition to being limited to a small number of instruments, slip rings are prone to wear quickly, and introduce noise and other undesirable characteristics to the signal data. This led to the development of a system capable of recording high-density instrumentation, at high sample rates, on the test article itself, all while under extreme rotational stress. OSDAS is a fully functional PC-based system with 48 channels of 24-bit, high-sample-rate input channels, phase synchronized, with an onboard storage capacity of over 1/2-terabyte of solid-state storage. This recording system takes a novel approach to the problem of recording multiple channels of instrumentation, integrated with the test article itself, packaged in a compact/rugged form factor, consuming limited power, all while rotating at high turbine speeds.

Pedings, Marc↗

Model-based parameter estimation using cardiovascular response to orthostatic stress

This paper presents a cardiovascular model that is capable of simulating the short-term (< or approximately equal to 3 min) transient hemodynamic response to gravitational stress and a gradient-based optimization method that allows for the automated estimation of model parameters from simulated or experimental data. We perform a sensitivity analysis of the transient heart rate response to determine which parameters of the model impact the heart rate dynamics significantly. We subsequently include only those parameters in the estimation routine that impact the transient heart rate dynamics substantially. We apply the estimation algorithm to both simulated and real data and showed that restriction to the 20 most important parameters does not impair our ability to match the data.

Non-NASA Center↗

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗

Relating a Spectral Index from MODIS and Tower-based Measurements to Ecosystem Light Use Efficiency for a Fluxnet-Canada Coniferous Forest

As part of the North American Carbon Program effort to quantify the terrestrial carbon budget of North America, we have been examining the possibility of retrieving ecosystem light use efficiency (LUE, the carbon sequestered per unit photosynthetically active radiation) directly from satellite observations. Our novel approach has been to compare LUE derived from tower fluxes with LUE estimated using spectral indices computed from MODIS satellite observations over forests in the Fluxnet-Canada Research Network, using the MODIS narrow ocean bands acquired over land. We matched carbon flux data collected around the time of the MODIS mid-day overpass for over one hundred relatively clear days in five years (2001-2006) from a mature Douglas fir forest in British Columbia. We also examined hyperspectral reflectance data collected diurnally from the tower in conjunction with the eddy correlation fluxes and meteorological measurements made throughout the 2006 growing season at this site. The tower-based flux data provided an opportunity to examine diurnal and seasonal LUE processes and their relationship to spectral indices at the scale of the forest stand. We evaluated LUE in conjunction with the Photochemical Reflectance Index (PRI), a normalized difference spectral index that uses 531 nm and a reference band to capture responses to high light induced stress afforded by the xanthophyll cycle. Canopy structure information, retrieved from airborne laser scanning radar (LiDAR) observations, was used to partition the forest canopy into sunlit and shaded fractions throughout the day, on numerous days during 2006. At each observation period throughout a day, the PRI was examined for the sunlit, shaded, and intermediate canopy segments defined by their instantaneous position relative to the solar principal plane (SPP). The sunlit sector was associated with the illumination "hotspot" (the reflectance backscatter maximum), the shaded sector with the "cold or dark spot" (the reflectance forward scatter minimum), while the intermediate, mixed sunlit/shade sector was located in the cross-plane to the SPP. The PRI indices clearly captured the differences in leaf groups, with sunlit foliage exhibiting the lowest values on sunny days throughout the 2006 season. When tower-based canopy-level LUE was recalculated to estimate foliage-based values (LUE(sub foilage) for the three foliage groups under their incident light environments, a strong linear relationship for PRI:LUE(sub foilage) was demonstrated (0.6 less than or equal to r(sup 2) less than or equal to 0.8, n=822, P<0.0001). The MODIS data represent relatively large areas when acquired at nadir (approx.1 sq km) or at variable off-nadir view angles (greater than or equal to 1 sq km) looking forward or aft. Nevertheless, a similar relationship between MODIS PRI and tower-based LUE was obtained from satellite observations (r(sup 2) = 0.76, n=105, P= 0.026) when the azimuth offsets from the SPP for off-nadir observations were considered. At this relatively high latitude of 50 degrees, the MODIS directional observations were offset from the SPP by approximately 50 degrees, but still represented backscatter or forward scatter sectors of the bidirectional reflectance distribution function (BRDF). The backscatter observations sampled the sunlit forest and provided lower PRI values, in general, than the forward scatter observations from the shaded forest. Since the hotspot and darkspot were not typically directly observed, the dynamic range for MODIS PRI was less than that observed in the SPP at the canopy level; therefore, MODIS PRI values were more similar to those observed in sifu in the BRDF cross-plane. While not ideal in terms of spatial resolution or optimal viewing configuration, the MODIS observations nevertheless provide a means to monitor forest under stress using narrow spectral band indices and off-nadir observations. This research has stimulated several spin-off studies for remote sensinf LUE, and demonstrates the importance of the connection between ecosystem structure and physiological function.

Middleton, Elizabeth M.↗

Development of a finite element based delamination analysis for laminates subject to extension, bending, and torsion

Delamination is a common failure mode of laminated composite materials. This type of failure frequently occurs at the free edges of laminates where singular interlaminar stresses are developed due to the difference in Poisson's ratios between adjacent plies. Typically the delaminations develop between 90 degree plies and adjacent angle plies. Edge delamination has been studied by several investigators using a variety of techniques. Recently, Chan and Ochoa applied the quasi-three-dimensional finite element model to the analysis of a laminate subject to bending, extension, and torsion. This problem is of particular significance relative to the structural integrity of composite helicopter rotors. The task undertaken was to incorporate Chan and Ochoa's formulation into a Raju Q3DG program. The resulting program is capable of modeling extension, bending, and torsional mechanical loadings as well as thermal and hygroscopic loadings. The addition of the torsional and bending loading capability will provide the capability to perform a delamination analysis of a general unsymmetric laminate containing four cracks, each of a different length. The solutions obtained using this program are evaluated by comparing them with solutions from a full three-dimensional finite element solution. This comparison facilitates the assessment of three dimensional affects such as the warping constraint imposed by the load frame grips. It wlso facilitates the evaluation of the external load representation employed in the Q3D formulation. Finally, strain energy release rates computed from the three-dimensional results are compared with those predicted using the quasi-three-dimensional formulation.

Hooper, Steven J.↗

Design and aerodynamic performance evaluation of a high-work mixed flow turbine stage

As axial and radial turbine designs have been pushed to their aerothermodynamic and mechanical limits, the mixed-flow turbine (MFT) concept has been projected to offer performance and durability improvements, especially when ceramic materials are considered. The objective of this NASA/U.S. Army sponsored mixed-flow turbine (AMFT) program was to determine the level of performance attainable with MFT technology within the mechanical constraints of 1997 projected ceramic material properties. The MFT geometry is similar to a radial turbine, exhibiting a large radius change from inlet to exit, but differing in that the inlet flowpath is not purely radial, nor axial, but mixed; it is the inlet geometry that gives rise to the name 'mixed-flow'. The 'mixed' orientation of the turbine inlet offers several advantages over radial designs by allowing a nonzero inlet blade angle yet maintaining radial-element blades. The oblique inlet not only improves the particle-impact survivability of the design, but improves the aerodynamic performance by reducing the incidence at the blade inlet. The difficulty, however, of using mixed-flow geometry lies in the scarcity of detailed data and documented design experience. This paper reports the design of a MFT stage designed with the intent to maximize aerodynamic performance by optimizing design parameters such as stage reaction, rotor incidence, flowpath shape, blade shape, vane geometry, and airfoil counts using 2-D, 3-D inviscid, and 3-D viscous computational fluid dynamics code. The aerodynamic optimization was accomplished while maintaining mechanical integrity with respect to vibration and stress levels in the rotor. A full-scale cold-flow rig test was performed with metallic hardware fabricated to the specifications of the hot ceramic geometry to evaluate the stage performance.

Neri, Remo N.↗

NASGRO 3.0: A Software for Analyzing Aging Aircraft

Structural integrity analysis of aging aircraft is a critical necessity in view of the increasing numbers of such aircraft in general aviation, the airlines and the military. Efforts are in progress by NASA, the FAA and the DoD to focus attention on aging aircraft safety. The present paper describes the NASGRO software which is well-suited for effectively analyzing the behavior of defects that may be found in aging aircraft. The newly revised Version 3.0 has many features specifically implemented to suit the needs of the aircraft community. The fatigue crack growth computer program NASA/FLAGRO 2.0 was originally developed to analyze space hardware such as the Space Shuttle, the International Space Station and the associated payloads. Due to popular demand, the software was enhanced to suit the needs of the aircraft industry. Major improvements in Version 3.0 are the incorporation of the ability to read aircraft spectra of unlimited size, generation of common aircraft fatigue load blocks, and the incorporation of crack-growth models which include load-interaction effects such as retardation due to overloads and acceleration due to underloads. Five new crack-growth models, viz., generalized Willenborg, modified generalized Willenborg, constant closure model, Walker-Chang model and the deKoning-Newman strip-yield model, have been implemented. To facilitate easier input of geometry, material properties and load spectra, a Windows-style graphical user interface has been developed. Features to quickly change the input and rerun the problem as well as examine the output are incorporated. NASGRO has been organized into three modules, the crack-growth module being the primary one. The other two modules are the boundary element module and the material properties module. The boundary-element module provides the ability to model and analyze complex two-dimensional problems to obtain stresses and stress-intensity factors. The material properties module allows users to store and curve-fit fatigue-crack growth data. On-line help and documentation are provided for each of the modules. In addition to the popular PC windows version, a unix-based X-windows version of NASGRO is also available. A portable C++ class library called WxWindows was used to facilitate cross-platform availability of the software.

Mettu, S. R.↗

Small Engine Technology (Set) Task 8 Aeroelastic Prediction Methods

AlliedSignal Engines, in cooperation with NASA LeRC, completed an evaluation of recently developed aeroelastic computer codes using test cases from the AlliedSignal Engines fan blisk database. Test data for this task includes strain gage, light probe, performance, and steady-state pressure information obtained for conditions where synchronous or flutter vibratory conditions were found to occur. Aeroelastic codes evaluated include the quasi 3-D UNSFLO (developed at MIT and modified to include blade motion by AlliedSignal), the 2-D FREPS (developed by NASA LeRC), and the 3-D TURBO-AE (under development at NASA LeRC). Six test cases each where flutter and synchronous vibrations were found to occur were used for evaluation of UNSFLO and FREPS. In addition, one of the flutter cases was evaluated using TURBO-AE. The UNSFLO flutter evaluations were completed for 75 percent radial span and provided good agreement with the experimental test data. Synchronous evaluations were completed for UNSFLO but further enhancement needs to be added to the code before the unsteady pressures can be used to predict forced response vibratory stresses. The FREPS evaluations were hindered as the steady flow solver (SFLOW) was unable to converge to a solution for the transonic flow conditions in the fan blisk. This situation resulted in all FREPS test cases being attempted but no results were obtained during the present program. Currently, AlliedSignal is evaluating integrating FREPS with our existing steady flow solvers to bypass the SFLOW difficulties. ne TURBO-AE steady flow solution provided an excellent match with the AlliedSignal Engines calibrated DAWES 3-D viscous solver. Finally, the TURBO-AE unsteady analyses also matched experimental observations by predicting flutter for the single test case evaluated.

Eick, Chris D.↗

Veggie Hardware Validation Test Preliminary Results and Lessons Learned

The Veggie hardware validation test, VEG-01, was conducted on the International Space Station during Expeditions 39 and 40 from May through June of 2014. The Veggie hardware and the VEG-01 experiment payload were launched to station aboard the SpaceX-3 resupply mission in April, 2014. Veggie was installed in an Expedite-the-Processing-of-Experiments-to-Space-Station (ExPRESS) rack in the Columbus module, and the VEG-01 validation test was initiated. Veggie installation was successful, and power was supplied to the unit. The hardware was programmed and the root mat reservoir and plant pillows were installed without issue. As expected, a small amount of growth media was observed in the sealed bags which enclosed the plant pillows when they were destowed. Astronaut Steve Swanson used the wet/dry vacuum to clean up the escaped particles. Water insertion or priming the first plant pillow was unsuccessful as an issue prevented water movement through the quick disconnect. All subsequent pillows were successfully primed, and the initial pillow was replaced with a backup pillow and successfully primed. Six pillows were primed, but only five pillows had plants which germinated. After about a week and a half it was observed that plants were not growing well and that pillow wicks were dry. This indicated that the reservoir was not supplying sufficient water to the pillows via wicking, and so the team reverted to an operational fix which added water directly to the plant pillows. Direct watering of the pillows led to a recovery in several of the stressed plants; a couple of which did not recover. An important lesson learned involved Veggie's bellows. The bellows tended to float and interfere with operations when opened, so Steve secured them to the baseplate during plant tending operations. Due to the perceived intensity of the LED lights, the crew found it challenging to both work under the lights and read crew procedures on their computer. Although the lights are not a safety hazard, for visual comfort crewmembers were advised to wear sunglasses when working with the plants and then they can lift glasses to read procedures. Steve Swanson had already trail-blazed this procedure when he initiated VEG-01. The temperature and humidity data logger was relocated mid-experiment to provide measurements on both sides of the unit. Images of the plants were downlinked weekly, and videos of installation and harvest were recorded. This imaging frequency was not sufficient to monitor and respond to changes in plant growth. Plants, samples, and data loggers will be returned on SpaceX-4, scheduled to return the fall of 2014. Lessons learned will be translated into hardware and operational modifications for future Veggie payloads.

International Space Station↗

Slow Crack Growth and Fatigue Life Prediction of Ceramic Components Subjected to Variable Load History

Present capabilities of the NASA CARES/Life (Ceramic Analysis and Reliability Evaluation of Structures/Life) code include probabilistic life prediction of ceramic components subjected to fast fracture, slow crack growth (stress corrosion), and cyclic fatigue failure modes. Currently, this code has the capability to compute the time-dependent reliability of ceramic structures subjected to simple time-dependent loading. For example, in slow crack growth (SCG) type failure conditions CARES/Life can handle the cases of sustained and linearly increasing time-dependent loads, while for cyclic fatigue applications various types of repetitive constant amplitude loads can be accounted for. In real applications applied loads are rarely that simple, but rather vary with time in more complex ways such as, for example, engine start up, shut down, and dynamic and vibrational loads. In addition, when a given component is subjected to transient environmental and or thermal conditions, the material properties also vary with time. The objective of this paper is to demonstrate a methodology capable of predicting the time-dependent reliability of components subjected to transient thermomechanical loads that takes into account the change in material response with time. In this paper, the dominant delayed failure mechanism is assumed to be SCG. This capability has been added to the NASA CARES/Life (Ceramic Analysis and Reliability Evaluation of Structures/Life) code, which has also been modified to have the ability of interfacing with commercially available FEA codes executed for transient load histories. An example involving a ceramic exhaust valve subjected to combustion cycle loads is presented to demonstrate the viability of this methodology and the CARES/Life program.

Jadaan, Osama↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

Working toward a three-dimensional fatigue closure model for surface cracks

The first reliable elastic fracture mechanics solutions for a surface crack in a plate were obtained by Newman and Raju. The authors, both from the Mechanics of Materials Branch at NASA-Langley, used a highly detailed finite element solution requiring substantial computational resources. Computers have since become more powerful and available; however, many important related problems remain computationally expensive. The problem of three-dimensional fatigue crack growth taking into account plasticity-induced crack closure is one such problem. It is the goal of this research to provide an efficient method to account for three-dimensional crack closure in fatigue. Newman developed a two-dimensional plasticity-induced crack closure model for center cracked specimens. This model requires iterations to determine both the contact solution at each growth step and the extent of the plastic zone at the crack tip. A three-dimensional version of this model would require obtaining these nonlinear variables all along the crack front. This model must be efficient so that repeated calculations can be performed during crack growth simulations. The highly versatile line spring model (LSM) with contact, fatigue, and plasticity will form the basis of the closure model. There are several required additions to past work to address the three-dimensional crack closure problem. Initially, these additions will include (1) an improved LSM to more accurately obtain the crack opening displacement, stress intensity factors, and elastic T-stress near the ends of the surface crack; (2) a method to determine the extent of the plastic zone all along the crack front; (3) a method to determine the contact zone given a perfectly plastic layer of material on the crack surfaces; (4) a method to determine the magnitude of the compressive contact stress; and (5) a way to implement the degree of constraint along the curved crack front. During the summer ASEE program an enhanced LSM was developed. A method similar to that of 'strip synthesis' first introduced by Fujimoto was used. Briefly, the crack opening displacements of 'slices' of the surface crack in a direction parallel to the plate surface are considered in addition to the standard LSM approach that makes use of springs obtained from slices perpendicular to the plate surface. This enhancement is necessary so that an accurate three-dimensional representation of quantities such as contact zone size, plastic zone size, stress intensity factors, T-stress, and crack opening displacement can be determined. By combining results of previous investigations with the LSM, the problem of three-dimensional crack closure will be addressed. In addition to a closure model, the enhanced LSM can be used for many other problems including interacting surface cracks and fatigue crack growth of a through crack with a curved crack front.

Joseph, Paul F.↗

Development and Implementation of Non-Newtonian Rheology Into the Generalized Fluid System Simulation Program (GFSSP)

One of the goals of this program was to develop the experimental and analytical/computational tools required to predict the flow of non-Newtonian fluids through the various system components of a propulsion system: pipes, valves, pumps etc. To achieve this goal we selected to augment the capabilities of NASA's Generalized Fluid System Simulation Program (GFSSP) software. GFSSP is a general-purpose computer program designed to calculate steady state and transient pressure and flow distributions in a complex fluid network. While the current version of the GFSSP code is able to handle various systems components the implicit assumption in the code is that the fluids in the system are Newtonian. To extend the capability of the code to non-Newtonian fluids, such as silica gelled fuels and oxidizers, modifications to the momentum equations of the code have been performed. We have successfully implemented in GFSSP flow equations for fluids with power law behavior. The implementation of the power law fluid behavior into the GFSSP code depends on knowledge of the two fluid coefficients, n and K. The determination of these parameters for the silica gels used in this program was performed experimentally. The n and K parameters for silica water gels were determined experimentally at CFDRC's Special Projects Laboratory, with a constant shear rate capillary viscometer. Batches of 8:1 (by weight) water-silica gel were mixed using CFDRC s 10-gallon gelled propellant mixer. Prior to testing the gel was allowed to rest in the rheometer tank for at least twelve hours to ensure that the delicate structure of the gel had sufficient time to reform. During the tests silica gel was pressure fed and discharged through stainless steel pipes ranging from 1", to 36", in length and three diameters; 0.0237", 0.032", and 0.047". The data collected in these tests included pressure at tube entrance and volumetric flowrate. From these data the uncorrected shear rate, shear stress, residence time, and viscosity were evaluated using formulae for non-Newtonian, power law fluids. The maximum shear rates (corrected for entrance effects) obtained in the rheometer with the current setup were in the 150,000 to 170,000sec- range. GFSSP simulations were performed with a flow circuit simulating the capillary rheometer and using Power Law gel viscosity coefficients from the experimental data. The agreement between the experimental data and the simulated flow curves was within +/-4% given quality entrance effect data.

DiSalvo, Roberto↗

High-lift flow-physics flight experiments on a subsonic civil transport aircraft (B737-100)

As part of the subsonic transport high-lift program, flight experiments are being conducted using NASA Langley's B737-100 to measure the flow characteristics of the multi-element high-lift system at full-scale high-Reynolds-number conditions. The instrumentation consists of hot-film anemometers to measure boundary-layer states, an infra-red camera to detect transition from laminar to turbulent flow, Preston tubes to measure wall shear stress, boundary-layer rakes to measure off-surface velocity profiles, and pressure orifices to measure surface pressure distributions. The initial phase of this research project was recently concluded with two flights on July 14. This phase consisted of a total of twenty flights over a period of about ten weeks. In the coming months the data obtained in this initial set of flight experiments will be analyzed and the results will be used to finalize the instrumentation layout for the next set of flight experiments scheduled for Winter and Spring of 1995. The main goal of these upcoming flights will be: (1) to measure more detailed surface pressure distributions across the wing for a range of flight conditions and flap settings; (2) to visualize the surface flows across the multi-element wing at high-lift conditions using fluorescent mini tufts; and (3) to measure in more detail the changes in boundary-layer state on the various flap elements as a result of changes in flight condition and flap deflection. These flight measured results are being correlated with experimental data measured in ground-based facilities as well as with computational data calculated with methods based on the Navier-Stokes equations or a reduced set of these equations. Also these results provide insight into the extent of laminar flow that exists on actual multi-element lifting surfaces at full-scale high-life conditions. Preliminary results indicate that depending on the deflection angle, the slat and flap elements have significant regions of laminar flow over a wide range of angles of attack. Boundary-layer transition mechanisms that were observed include attachment-line contamination on the slat and inflectional instability on the slat and fore flap. Also, the results agree fairly well with the predictions reported in a paper presented at last year's AIAA Fluid Dynamics Conference. The fact that extended regions of laminar flow are shown to exist on the various elements of the high-lift system raises the question what the effect is of loss of laminar flow as a result of insect contamiantion, rain or ice accumulation on high-life performance.

Vandam, Cornelis P.↗

Development of a knoweledge-based system for validating finite element models

The finite element modeling of an airframe structure requires knowledge of general purpose programs such as NASTRAN as well as a detailed understanding of the airframe structure. Due to the sophistication of general purpose programs such as NASTRAN, a substantial investment in time and effort is required to gain expertise in using them effectively. This paper describes the development of an expert system used in the validation of NASTRAN based finite element models. Experts in the NASTRAN based finite element modeling of airframe structures were interviewed to document, understand, and represent their knowledge and reasoning in the expert system. Finite element stress analysis and internal loads reports generated by the experts were reviewed to determine expert resolution of problem areas. As a result, areas requiring expert assistance in the modeling of the airframe structures were identified. The finite element input data is represented as a set of 'facts'. A rule based representation is used to code the expert knowledge. A hierarchical set of rules are applied. The expert system first acts as an intelligent front end to insure that all the prerequisites needed to perform the analysis are present. This includes material properties, boundary conditions, and connectivity information. The next step examines if incompatible sets of elements are connected. The succeeding step examines if the specific airframe component is modeled by the appropriate set of elements. The next and final step examines if the airframe members used are adequate to represent the anticipated state of stress. The expert system is designed to inform the user the severity of the error, the likely consequence and possible remedial action. The shell used in the development of the expert system is CLIPS. CLIPS contains a forward chaining inference engine based on the Rete algorithm. CLIPS may be implemented on most personal computers as well as mini computers and mainframes. The approach taken is general and can be implemented using commercially available expert system shells.

Munir, N. I.↗