Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive state estimation”

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 451 records · Page 25

CFD flowfield simulation of Delta Launch Vehicles in a power-on configuration

This paper summarizes recent work at McDonnell Douglas Aerospace (MDA) to develop and validate computational fluid dynamic (CFD) simulations of under expanded rocket plume external flowfields for multibody expendable launch vehicles (ELVs). Multi engine reacting gas flowfield predictions of ELV base pressures are needed to define vehicle base drag and base heating rates for sizing external nozzle and base region insulation thicknesses. Previous ELV design programs used expensive multibody power-on wind tunnel tests that employed chamber/nozzle injected high pressure cold or hot-air. Base heating and pressure measurements were belatedly made during the first flights of past ELV's to correct estimates from semi-empirical engineering models or scale model tests. Presently, CFD methods for use in ELV design are being jointly developed at the Space Transportation Division (MDA-STD) and New Aircraft Missiles Division (MDA-NAMD). An explicit three dimensional, zonal, finite-volume, full Navier-Stokes (FNS) solver with finite rate hydrocarbon/air and aluminum combustion kinetics was developed to accurately compute ELV power-on flowfields. Mississippi State University's GENIE++ general purpose interactive grid generation code was chosen to create zonal, finite volume viscous grids. Axisymmetric, time dependent, turbulent CFD simulations of a Delta DSV-2A vehicle with a MB-3 liquid main engine burning RJ-1/LOX were first completed. Hydrocarbon chemical kinetics and a k-epsilon turbulence model were employed and predictions were validated with flight measurements of base pressure and temperature. Zonal internal/external grids were created for a Delta DSV-2C vehicle with a MB-3 and three Castor-1 solid motors burning and a Delta-2 with an RS-27 main engine (LOX/RP-1) and 9 GEM's attached/6 burning. Cold air, time dependent FNS calculations were performed for DSV-2C during 1992. Single phase simulations that employ finite rate hydrocarbon and aluminum (solid fuel) combustion chemistry are currently in progress. Reliable and efficient Eulerian algorithms are needed to model two phase (solid-gas) momentum and energy transfer mechanisms for solid motor fuel combustion products.

Pavish, D. L.↗

Some fundamental aspects of the fatigue and fracture response of composite materials

A simple analytical model, based on the concept that the constraining plies surrounding cracks through the thickness of a given ply or plies control the spacing of those cracks by controlling the surrounding stress state, appears to agree well with experimental data. Such a model can be used to evaluate critically various constraint situations including the prediction of the ability of a given constraint layer to arrest flaw development. The analytical model, which uses an equilibrium element approach in the present case, is sensitive to lamina and laminate elastic properties, flaw size, stacking sequence, and an interlaminar stress transfer distance which is estimated from experimental observations. The analysis is easily extended to study flaw growth and plastic region formation.

Reifsnider, K. L.↗

A study of fiber volume fraction effects in notched unidirectional SCS-6/Ti-15V-3Cr-3Al-3Sn composite

Notched unidirectional SCS-6/Ti-15-3 composite of three different fiber volume fractions (vf = 0.15, 0.37, and 0.41) was investigated for various room temperature microstructural and material properties including: fatigue crack initiation, fatigue crack growth, and fracture toughness. While the matrix hardness is similar for all fiber volume fractions, the fiber/matrix interfacial shear strength and matrix residual stress increases with fiber volume fraction. The composite fatigue crack initiation stress is shown to be matrix controlled and occurs when the net maximum matrix stress approaches the endurance limit stress of the matrix. A model is presented which includes residual stresses and presents the composite initiation stress as a function of fiber volume fraction. This model predicts a maximum composite initiation stress at vf approximately 0.15 which agrees with the experimental data. The applied composite stress levels were increased as necessary for continued crack growth. The applied Delta(K) values at crack arrest increase with fiber volume fraction by an amount better approximated using an energy based formulation rather than when scaled linear with modulus. After crack arrest, the crack growth rate exponents for vf37 and vf41 were much lower and toughness much higher, when compared to the unreinforced matrix, because of the bridged region which parades with the propagating fatigue crack. However, the vf15 material exhibited a higher crack growth rate exponent and lower toughness than the unreinforced matrix because once the bridged fibers nearest the crack mouth broke, the stress redistribution broke all bridged fibers, leaving an unbridged crack. Degraded, unbridged behavior is modeled using the residual stress state in the matrix ahead of the crack tip. Plastic zone sizes were directly measured using a metallographic technique and allow prediction of an effective matrix stress intensity which agrees with the fiber pressure model if residual stresses are considered. The sophisticated macro/micro finite element models of the 0.15 and 0.37 fiber volume fractions presented show good agreement with experimental data and the fiber pressure model when an estimated effective fiber/matrix debond length is used.

Covey, Steven J.↗

Verification of thermal analysis codes for modeling solid rocket nozzles

One of the objectives of the Solid Propulsion Integrity Program (SPIP) at Marshall Space Flight Center (MSFC) is development of thermal analysis codes capable of accurately predicting the temperature field, pore pressure field and the surface recession experienced by decomposing polymers which are used as thermal barriers in solid rocket nozzles. The objective of this study is to provide means for verifications of thermal analysis codes developed for modeling of flow and heat transfer in solid rocket nozzles. In order to meet the stated objective, a test facility was designed and constructed for measurement of the transient temperature field in a sample composite subjected to a constant heat flux boundary condition. The heating was provided via a steel thin-foil with a thickness of 0.025 mm. The designed electrical circuit can provide a heating rate of 1800 W. The heater was sandwiched between two identical samples, and thus ensure equal power distribution between them. The samples were fitted with Type K thermocouples, and the exact location of the thermocouples were determined via X-rays. The experiments were modeled via a one-dimensional code (UT1D) as a conduction and phase change heat transfer process. Since the pyrolysis gas flow was in the direction normal to the heat flow, the numerical model could not account for the convection cooling effect of the pyrolysis gas flow. Therefore, the predicted values in the decomposition zone are considered to be an upper estimate of the temperature. From the analysis of the experimental and the numerical results the following are concluded: (1) The virgin and char specific heat data for FM 5055 as reported by SoRI can not be used to obtain any reasonable agreement between the measured temperatures and the predictions. However, use of virgin and char specific heat data given in Acurex report produced good agreement for most of the measured temperatures. (2) Constant heat flux heating process can produce a much higher heating rate than the radiative heating process. The results show that heating rates of 125 C/s (225 F/s) can be achieved. (3) The electrical resistance of the FM 5055 samples were about 150 Omega in the virgin state, and decreased to about 35 Omega when charred. A reliable scheme must be developed to electrically insulate the composite from the heater foil in order to prevent any current leakage into the sample which can result in volumetric heating of the char zone of the composite.

Keyhani, M.↗

Investigating the potential of a global precipitation forecast to inform landslide prediction

Extreme rainfall events within landslide-prone areas can be catastrophic, resulting in loss of property, infrastructure, and life. A global Landslide Hazard Assessment for Situational Awareness (LHASA) model provides routine near-real time estimates of landslide hazard using Integrated Multi-Satellite Precipitation Retrievals for the Global Precipitation Mission (IMERG). However, it does not provide information on potential landslide hazard in the future. Forecasting potential landslide events at a global scale presents an area of open research. This study compares a global precipitation forecast provided by NASA's Goddard Earth Observing System (GEOS) with near-real time satellite precipitation estimates. The Multi-Radar Multi-Sensor gauge corrected (MRMS-GC) reference is used to assess the performance of both satellite and model-based precipitation products over the contiguous United States (CONUS). The forecast lead time of 24hrs is considered, with a focus on extreme precipitation events. The performance of IMERG and GEOS-Forecast products is assessed in terms of the probability of detection, success ratio, critical success index and hit bias as well as continuous statistics. The results show that seasonality influences the performance of both satellite and model-based precipitation products. Comparison of IMERG and GEOS-Forecast globally as well as in several event case studies (Colombia, southeast Asia, and Tajikistan) reveals that GEOS-Forecast detects extreme rainfall more frequently relative to IMERG for these specific analyses. For recent landslide points across the globe, the 24hr accumulated precipitation forecast >100 mm corresponds well with near-real time daily accumulated IMERG precipitation estimates. GEOS-Forecast and IMERG precipitation match more closely for tropical cyclones than for other types of storms. The main intention of this study is to assess the viability of using a global forecast for landslide predictions and understand the extent of the variability between these products to inform where we would expect the landslide modeling results to most prominently diverge. Results of this study will be used to inform how forecasted precipitation estimates can be incorporated into the LHASA model to provide the first global predictive view of landslide hazards.

S. Khan↗

Mass Inferencing Model Creation and Deployment to the RASSOR Lunar Excavation Robot

The Regolith Advanced Surface Systems Operations Robot (RASSOR) Excavator is a mobile robotic bucket-drum excavator platform with a unique space regolith excavation capability. The Intelligent Capabilities Enhanced RASSOR research project developed functionality for estimating the quantity of regolith mass ingested during RASSOR operation, enhancing RASSOR’s ability to successfully complete In-Situ Resource Utilization (ISRU) missions. To teleoperate or run autonomously, it is crucial for the amount of regolith mass ingested to be available as a system state for efficient operation. For example, during autonomous operation, RASSOR should navigate and move to a processing plant to offload the collected regolith when the drums are full; without knowledge of the total mass in the drums, this type of high-level planning is not possible. Three distinct modeling approaches were employed in developing a mass inferencing approach that could work on RASSOR, none of which require modification to the hardware. All take in system states, such as arm/drum motor positions, velocities, currents, voltages, and robot pose, and output a mass prediction for each set of the robot’s bucket drums. The developed models run in real-time, outputting predictions for the front drum mass, rear drum mass, timestamp of the last prediction, and total drum mass (sum of front and rear) in RASSOR’s drums. Models deployed to the hardware have low error (<7.5% mean error over the mass range, and <2.6% mean error when drums are more than half full) when making predictions in real-time. Our modeling approach can be adapted to use lunar excavation data to create models that are reflective of RASSOR’s dynamics when operating on the lunar surface. The results of this work are promising and show that models can be developed to accurately estimate excavated regolith mass.

Bucket Drum Excavators↗

Mass Inferencing Model Creation and Deployment to the RASSOR Lunar Excavation Robot

The Regolith Advanced Surface Systems Operations Robot (RASSOR) Excavator is a mobile robotic bucket-drum excavator platform with a unique space regolith excavation capability. The Intelligent Capabilities Enhanced RASSOR research project developed functionality for estimating the quantity of regolith mass ingested during RASSOR operation, enhancing RASSOR’s ability to successfully complete In-Situ Resource Utilization (ISRU) missions. To teleoperate or run autonomously, it is crucial for the amount of regolith mass ingested to be available as a system state for efficient operation. For example, during autonomous operation, RASSOR should navigate and move to a processing plant to offload the collected regolith when the drums are full; without knowledge of the total mass in the drums, this type of high-level planning is not possible. Three distinct modeling approaches were employed in developing a mass inferencing approach that could work on RASSOR, none of which require modification to the hardware. All take in system states, such as arm/drum motor positions, velocities, currents, voltages, and robot pose, and output a mass prediction for each set of the robot’s bucket drums. The developed models run in real-time, outputting predictions for the front drum mass, rear drum mass, timestamp of the last prediction, and total drum mass (sum of front and rear) in RASSOR’s drums. Models deployed to the hardware have low error (<7.5% mean error over the mass range, and <2.6% mean error when drums are more than half full) when making predictions in real-time. Our modeling approach can be adapted to use lunar excavation data to create models that are reflective of RASSOR’s dynamics when operating on the lunar surface. The results of this work are promising and show that models can be developed to accurately estimate excavated regolith mass.

bucket drum excavators↗

Investigating the Potential of A Global Precipitation Forecast to Inform Landslide Prediction

Extreme rainfall events within landslide-prone areas can be catastrophic, resulting in loss of property, infra-structure, and life. A global Landslide Hazard Assessment for Situational Awareness (LHASA) model provides routine near-real time estimates of landslide hazard using Integrated Multi-Satellite Precipitation Retrievals for the Global Precipitation Mission (IMERG). However, it does not provide information on potential landslide hazard in the future. Forecasting potential landslide events at a global scale presents an area of open research. This study compares a global precipitation forecast provided by NASA’s Goddard Earth Observing System (GEOS) with near-real time satellite precipitation estimates. The Multi-Radar Multi-Sensor gauge corrected (MRMS-GC) reference is used to assess the performance of both satellite and model-based precipitation products over the contiguous United States (CONUS). The forecast lead time of 24hrs is considered, with a focus on extreme precipitation events. The performance of IMERG and GEOS-Forecast products is assessed in terms of the probability of detection, success ratio, critical success index and hit bias as well as continuous statistics. The results show that seasonality influences the performance of both satellite and model-based precipitation products. Comparison of IMERG and GEOS-Forecast globally as well as in several event case studies (Colombia, southeast Asia, and Tajikistan) reveals that GEOS-Forecast detects extreme rainfall more frequently relative to IMERG for these specific analyses. For recent landslide points across the globe, the 24hr accumulated precipitation forecast >100 mm corresponds well with near-real time daily accumulated IMERG precipitation estimates. GEOS-Forecast and IMERG precipitation match more closely for tropical cyclones than for other types of storms. The main intention of this study is to assess the viability of using a global forecast for landslide predictions and understand the extent of the variability between these products to inform where we would expect the landslide modeling results to most prominently diverge. Results of this study will be used to inform how forecasted precipitation estimates can be incorporated into the LHASA model to provide the first global predictive view of landslide hazards.

S. Khan↗

Theoretical spectrum of AlN

Theoretical spectroscopic constants and radiative lifetimes are presented for the low-lying electronic states of AlN and AlN(+). The ground state of AlN is predicted to be X 3Pi, in agreement with earlier theoretical calculations, but the A 3Pi- state is extremely low-lying. The near degeneracy of these states is analogous to the situation for Al2 and implies that sigma and pi bonds are of nearly equal strength. The 3Pi - 3Pi AlN band system observed in emission is assigned as the C 3Pi - X 3Pi transition. The C 3Pi radiative lifetime is estimated as 151, 163, and 175 ns for v-prime = 0, 1, 2 respectively. These lifetimes may be significantly shortened by dissociation by the (1) 5Pi stated. The calculated D(e) of 2.35 eV for the ground state is about half that of the isovalent BN molecule. Thus the bonding in AlN is significantly weaker than for BN owing to the smaller overlaps of the valence orbitals.

Langhoff, Stephen R.↗

The Interaction of a Reflected Shock Wave with the Boundary Layer in a Shock Tube

Ideally, the reflection of a shock from the closed end of a shock tube provides, for laboratory study, a quantity of stationary gas at extremely high temperature. Because of the action of viscosity, however, the flow in the real case is not one-dimensional, and a boundary layer grows in the fluid following the initial shock wave. In this paper simplifying assumptions are made to allow an analysis of the interaction of the shock reflected from the closed end with the boundary layer of the initial shock afterflow. The analysis predicts that interactions of several different types will exist in different ranges of initial shock Mach number. It is shown that the cooling effect of the wall on the afterflow boundary layer accounts for the change in interaction type. An experiment is carried out which verifies the existence of the several interaction regions and shows that they are satisfactorily predicted by the theory. Along with these results, sufficient information is obtained from the experiments to make possible a model for the interaction in the most complicated case. This model is further verified by measurements made during the experiment. The case of interaction with a turbulent boundary layer is also considered. Identifying the type of interaction with the state of turbulence of the interacting boundary layer allows for an estimate of the state of turbulence of the boundary layer based on an experimental investigation of the type of interaction. A method is proposed whereby the effect of the boundary-layer interaction on the strength of the reflected shock may be calculated. The calculation indicates that the reflected shock is rapidly attenuated for a short distance after reflection, and this result compares favorably with available experimental results.

Mark, Herman↗

Curie isotherm surfaces inferred from high-altitude magnetic anomaly data

Two-dimensional Curie depth models along two sections across the western United States are derived using an equivalent layer magnetization model derived from Magsat data, and the results are used to constrain finite element thermal models developed by Mayhew and Majer (1980). Regional heat flow variations predicted by the models compare favorably with those inferred from direct measurements. The methodology is applied to upward-continued aeromagnetic data, and the results are found to be in good agreement with previous Curie depth estimates and average measured heat flow.

Mayhew, M. A.↗

Assessment of the 1997-1998 Asian Monsoon Anomalies

Using State-of-the-art satellite-gauge monthly rainfall estimate and optimally interpolated sea surface temperature (SST) data, we have assessed the 1997-98 Asian monsoon anomalies in terms of three basic causal factors: basin-scale SST, regional coupling, and internal variability. Singular Value Decomposition analysis of rainfall and SST are carried out globally over the entire tropics and regionally over the Asian monsoon domain. Contributions to monsoon rainfall predictability by various factors are evaluated from cumulative anomaly correlation with dominant regional SVD modes. Results reveal a dominant, large-scale monsoon-El Nino coupled mode with well-defined centers of action in the near-equatorial monsoon regions. it is noted that some subcontinental regions such as all-India, or arbitrarily chosen land regions over East Asia, while important socio-economically, are not near the centers of influence from El Nino, hence are not necessarily representative of the response of the entire monsoon region to El Nino. The observed 1997-98 Asian monsoon anomalies are found to be very complex with approximately 34% of the anomalies attributable to basin- scale SST influence associated with El Nino. Regional coupled processes contribute an additional 19%, leaving about 47% due to internal dynamics. Also noted is that the highest monsoon predictability is not necessary associated with major El Nino events (e.g. 1997, 1982) but rather in non-El Nino years (e.g. 1980, 1988) when contributions from the regional coupled modes far exceed those from the basin-scale SST. The results suggest that in order to improve monsoon seasonal-to-interannual predictability, there is a need to exploit not only monsoon-El Nino relationship, but also monsoon regional coupled processes and their modulation by long-term climate change.

Lau, William K.-M.↗

Mars InSight Entry, Descent, and Landing Trajectory and Atmosphere Reconstruction

The InSight mission landed on the surface of Mars on November 26th, 2018. The InSight system performance met all design requirements, although several performance metrics fell near the boundaries of the predictions. The peak deceleration was high, the overall timeline was short, and the landing site was uprange and crossrange from the target. This paper describes the reconstruction of the entry, descent, and landing trajectory and atmosphere. The approach utilizes a Kalman filter to blend sensor data to obtain the vehicle trajectory. The aerodynamic database is used in combination with the sensed accelerations to obtain estimates of the atmosphere-relative state, which in turn is used to derive the free-stream atmospheric conditions during entry, until the time of parachute deployment. The results indicate that the reconstructed atmosphere was approximately 1σbelow the preflight atmosphere. Analysis of the reconstructed vehicle attitude angles indicate that the aerodynamic lift was oriented downward at entry. The vehicle developed a roll rate during entry, which directed a component of the lift to the north. The low density and aerodynamic lift direction are determined to be the primary causes of the high deceleration, short timeline, and location of the landing site relative to the target.

Christopher D Karlgaard↗

Mars InSight Entry, Descent, and Landing Trajectory and Atmosphere Reconstruction

The InSight mission landed on the surface of Mars on November 26th, 2018. The InSight system performance met all design requirements, although several performance metrics fell near the boundaries of the predictions. The peak deceleration was high, the overall timeline was short, and the landing site was uprange and crossrange from the target. This paper describes the reconstruction of the entry, descent, and landing trajectory and atmosphere. The approach utilizes a Kalman filter to blend sensor data to obtain the vehicle trajectory. The aerodynamic database is used in combination with the sensed accelerations to obtain estimates of the atmosphere-relative state, which in turn is used to derive the free-stream atmospheric conditions during entry, until the time of parachute deployment. The results indicate that the reconstructed atmosphere was approximately 1σbelow the preflight atmosphere. Analysis of the reconstructed vehicle attitude angles indicate that the aerodynamic lift was oriented downward at entry. The vehicle developed a roll rate during entry, which directed a component of the lift to the north. The low density and aerodynamic lift direction are determined to be the primary causes of the high deceleration, short timeline, and location of the landing site relative to the target.

Christopher D Karlgaard↗

A comparison of SIR-B directional ocean wave spectra with aircraft scanning radar spectra

Directional ocean wave spectra derived from Shuttle Imaging Radar-B (SIR-B) L-band imagery collected off the coast of Southern Chile on 11 and 12 October 1984 were compared with independent spectral estimates from two airborne scanning radars. In sea states with significant wave heights ranging from 3 to 5 meters, the SIR-B-derived sspectra at 18 deg and 25 deg off nadir yielded reasonable estimates of wavelengths, directions, and spectral shapes for all wave systems encountered, including a purely azimuth-traveling system. A SIR-B image intensity variance spectrum containing predominantly range-traveling waves closely resembles an independent aircraft estimate of the slope variance spectrum. The prediction of a U.S. Navy global spectral ocean wave model on 11 October 1984 exhibited no significant bias in dominant wave number but contained a directional bias of about 30 deg with respect to the mean of the aircraft and spacecraft estimates.

Beal, R. C.↗

Subsonic flight test evaluation of a propulsion system parameter estimation process for the F100 engine

An adaptive-performance-seeking control system which optimizes the quasi-steady-state performance of the F-15 propulsion system is discussed. This paper presents flight- and ground-test evaluations of the propulsion system parameter-estimation process used by the performance seeking control system. The estimator consists of a compact propulsion system model and an extended Kalman filter. The extended Kalman filter estimates five engine component deviation parameters from measured inputs. The compact model uses measurements and Kalman-filter estimates as inputs to predict unmeasured propulsion parameters such as net propulsive force and fan stall margin. The ability to track trends and estimate absolute values of propulsion system parameters was demonstrated. For example, thrust stand results show a good correlation especially in trends between the performance seeking control estimated and measured thrust.

Orme, John S.↗

Turbine Blade Research

Under contract with the New York State Energy Research and Development Authority, GE's Energy Systems Programs Department used a COSMIC program in assessing the problem of blade erosion in a PFB (pressurized fluid bed) environment. Data provided by this program and an associated program enabled the company engineers to determine gas velocities and the velocities of the particles striking the blades, calculations necessary for predicting blade erosion and potential damage. The assessment resulted in a new estimate for the allowable dust load for a modern heavy duty gas turbine.

Source record↗

The Predicted Growth of the Low Earth Orbit Space Debris Environment: An Assessment of Future Risk for Spacecraft

Space debris is a worldwide-recognized issue concerning the safety of commercial, military, and exploration spacecraft. The space debris environment includes both naturally occuring meteoroids and objects in Earth orbit that are generated by human activity, termed orbital debris. Space agencies around the world are addressing the dangers of debris collisions to both crewed and robotic spacecraft. In the United States, the Orbital Debris Program Office at the NASA Johnson Space Center leads the effort to categorize debris, predict its growth, and formulate mitigation policy for the environment from low Earth orbit (LEO) through geosynchronous orbit (GEO). This paper presents recent results derived from the NASA long-term debris environment model, LEGEND. It includes the revised NASA sodium potassium droplet model, newly corrected for a factor of two over-estimation of the droplet population. The study indicates a LEO environment that is already highly collisionally active among orbital debris larger than 1 cm in size. Most of the modeled collision events are non-catastrophic (i.e., They lead to a cratering of the target, but no large scale fragmentation.). But they are potentially mission-ending, and take place between impactors smaller than 10 cm and targets larger than 10 cm. Given the small size of the impactor these events would likely be undetectable by present-day measurement means. The activity continues into the future as would be expected. Impact rates of about four per year are predicted by the current study within the next 30 years, with the majority of targets being abandoned intacts (spent upper stages and spacecraft). Still, operational spacecraft do show a small collisional activity, one that increases over time as the small fragment population increases.

Krisko, Paula H.↗