SEARCH · Engineering Papers

Results for “Local time stepping”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Using Total Lightning Observations to Enhance Lightning Safety

Lightning is often the underrated threat faced by the public when it comes to dangerous weather phenomena. Typically, larger scale events such as floods, hurricanes, and tornadoes receive the vast majority of attention by both the general population and the media. This comes from the fact that these phenomena are large, longer lasting, can impact a large swath of society at one time, and are dangerous events. The threat of lightning is far more isolated on a case by case basis, although millions of cloud-to-ground lightning strikes hit this United States each year. While attention is given to larger meteorological events, lightning is the second leading cause of weather related deaths in the United States. This information raises the question of what steps can be taken to improve lightning safety. Already, the meteorological community s understanding of lightning has increased over the last 20 years. Lightning safety is now better addressed with the National Weather Service s access to the National Lightning Detection Network data and enhanced wording in their severe weather warnings. Also, local groups and organizations are working to improve public awareness of lightning safety with easy phrases to remember, such as "When Thunder Roars, Go Indoors." The impacts can be seen in the greater array of contingency plans, from airports to sports stadiums, addressing the threat of lightning. Improvements can still be made and newer technologies may offer new tools as we look towards the future. One of these tools is a network of sensors called a lightning mapping array (LMA). Several of these networks exist across the United States. NASA s Short-term Prediction Research and Transition Center (SPoRT), part of the Marshall Spaceflight Center, has access to three of these networks from Huntsville, Alabama, the Kennedy Space Center, and Washington D.C. The SPoRT program s mission is to help transition unique products and observations into the operational forecast environment. SPoRT has been collaborating with the Huntsville National Weather Service (NWS) Office since 2003 and has since included several other offices to better implement LMA observations into real-time applications. Much of that work has focused on the LMA s ability to detect intra-cloud lightning in addition to cloud-to-ground lightning strikes. Combined, these observations are called total lightning. With total lightning observations, NWS offices can enhance their situational awareness and improve severe weather warnings. Just as importantly, the observed intra-cloud flashes often precede the first cloud-to-ground strike by a few minutes. SPoRT and its partner NWS offices are working to develop visualizations and applications to better utilize these data. However, there is a drawback. The LMAs have a short range of no more than 200 km. This is being addressed with the next generation geostationary satellite, GOES-R, which will boast the Geostationary Lightning Mapper (GLM). SPoRT, in conjunction with NOAA s GOES-R Proving Ground, is working to prepare the end user community for the GLM era using the LMA observations as a demonstration tool. Working collaboratively with our NWS partners, SPoRT is working to determine how best to integrate these future observations to improve both severe storm warnings and lightning safety.

Stano, Geoffrey T.↗

NASA NTRS1992

MHD flow past an obstacle - Large-scale flow in the magnetosheath

As a step to the study of the large-scale flow in the magnetosheath, the MHD flow past an obstacle is investigated. A time asymptotic method is used to obtain 3D steady-state solutions. The results indicate the formation of a depletion layer near the obstacle due to the increase of the magnetic field. Along the earth-sun line the plasma density increases first and then decreases from the post bow shock to the magnetopause. The local density maximum in front of the magnetopause may correspond to what was recently observed. When the interplanetary magnetic field direction is tilted from the solar wind flow, the IMF influences the shape of the bow shock as well as the location of the stagnation point at the magnetopause in a way consistent with observation. In addition, the results show the existence of a magnetosheath current in the post parallel-shock region.

Wu, C. C.↗

NASA NTRSDate not supplied

Assessing the Viability of Using GEOS-Forecast Product for Landslides Forecasting: A Step Toward Early Warning System

Landslides across the globe are mostly triggered by extreme rainfall events affecting infrastructure, transportation and livelihoods. The risks are rarely quantified due to lack of data, analytical skills and limited modeling techniques. Knowledge of local to global scale landslide risks provides communities and national agencies the ability to adapt disaster management practices to mitigate and recover from these hazards. In order to minimize the risks and improve characterization of community resilience to landslides, it is vital to have reliable information about the factors triggering landslides such as rainfall, well ahead in time. Forecasting potential landslide activity and impacts can be achieved through reliable precipitation forecast models. However, it is challenging because of the temporal and spatial variability of precipitation, an important factor in triggering landslides. Evaluation of the precipitation field, associated errors, and sampling uncertainties is integral for development of efficient and reliable landslide forecasting and early warning system. This study develops a methodology to assess the viability of using a precipitation field provided by a global model and its potential integration in the landslide forecasting system. The study focuses on the comparison between the IMERG (Integrated Multi-satellitE Retrievals for Global Precipitation Mission) and GEOS (NASA Goddard Earth Observing System)-Forecast product over contiguous United States (CONUS). GEOS model assimilates new observations every 6 hours, at 00, 06, 12, and 18 UTC. The framework is tested on the GEOS-Forecast Model initialized at 00 UTC using daily IMERG early product as reference using both categorical and continuous statistics. The categorical statistics includes the probability of detection (POD), success ratio (SR), critical success index (CSI), and the hit bias. Continuous statistics such as correlation, normalized standard deviation, and root-mean-square error are also evaluated. Overall, GEOS-Forecast precipitation field over the analysis period (~1 year) show underestimation with respect to IMERG early for the daily accumulated rainfall. However, the probability distribution function and cumulative distribution function of both show similar patterns. In terms of correlations, POD, SR, CSI, hit bias, the performance varies with respect to the rainfall threshold used.

Sana Khan↗

NASA NTRSDate not supplied

Inferred Sea Level Prediction in the NASA GMAO Seasonal Forecasting System

Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.

Richard I Cullather↗

NASA NTRS2020

Demonstration of Autonomous Nested Search for Local Maxima using an Unmanned Underwater Vehicle, ICRA Submission

Ocean Worlds represent one of the best chances for extra-terrestrial life in our solar system. A new mission concept must be developed to explore these oceans. This mission would require traversing the 10s of km thick icy shell and releasing a submersible into the ocean below. During the transit of the icy shell and the exploration of the ocean, the vehicle(s) would be out of contact with Earth for weeks or potentially months at a time. During this time the vehicle must have sufficient autonomy to locate and study scientific targets of interest. One such target of interest is hydrothermal venting. We have previously developed an autonomous nested search method to locate and investigate sources of hydrothermal venting by locating local maxima in hydrothermal vent emissions. In this work we demonstrate this approach on board an OceanServer Iver2 AUV in Chesapeake Bay, MD using simulated sensor data from a hydrothermal plume model. This represents the first step towards the deployment of this approach in conditions analogous to those that we might expect on an Ocean World.

Branch, Andrew↗

NASA NTRS2000

Development of a Miniature, Two-Axis, Triple-Helmholtz-Driven Gimbal

This paper details the development of a Helmholtz-driven, 2-axis gimbal to position a flat mirror within 50 microradian (fine positioning) in a space environment. The gimbal is intended to travel on a deep space mission mounted on a miniature "rover" vehicle. The gimbal will perform both pointing and scanning functions. The goal for total mass of the gimbal was 25 grams. The primary challenge was to design and build a bearing system that would achieve the required accuracy in addition to supporting the relatively large mass of the mirror and the outer gimbal. The mechanism is subjected to 100-G loading without the aid of any additional caging mechanism. Additionally, it was desired to have the same level of accuracy during Earth-bound, 1-G testing. Due to the inherent lack of damping in a zero-G, vacuum environment; the ability of the gimbal to respond to very small amounts of input energy is paramount. Initial testing of the first prototype revealed exceedingly long damping times required even while exposed to the damping effects of air and 1-G friction. It is envisioned that fine positioning of the gimbal will be accomplished in very small steps to avoid large disturbances to the mirror. Various bearing designs, including materials, lubrication options and bearing geometry will be discussed. In addition various options for the Helmholtz coil design will be explored with specific test data given. Ground testing in the presence of 1-G was compounded by the local magnetic fields due to the "compass" effect on the gimbal. The test data will be presented and discussed. Additionally, rationale for estimating gimbal performance in a zero-G environment will be presented and discussed.

Sharif, Boz↗

NASA NTRS1999

Integration of a Decentralized Linear-Quadratic-Gaussian Control into GSFC's Universal 3-D Autonomous Formation Flying Algorithm

A decentralized control is investigated for applicability to the autonomous formation flying control algorithm developed by GSFC for the New Millenium Program Earth Observer-1 (EO-1) mission. This decentralized framework has the following characteristics: The approach is non-hierarchical, and coordination by a central supervisor is not required; Detected failures degrade the system performance gracefully; Each node in the decentralized network processes only its own measurement data, in parallel with the other nodes; Although the total computational burden over the entire network is greater than it would be for a single, centralized controller, fewer computations are required locally at each node; Requirements for data transmission between nodes are limited to only the dimension of the control vector, at the cost of maintaining a local additional data vector. The data vector compresses all past measurement history from all the nodes into a single vector of the dimension of the state; and The approach is optimal with respect to standard cost functions. The current approach is valid for linear time-invariant systems only. Similar to the GSFC formation flying algorithm, the extension to linear LQG time-varying systems requires that each node propagate its filter covariance forward (navigation) and controller Riccati matrix backward (guidance) at each time step. Extension of the GSFC algorithm to non-linear systems can also be accomplished via linearization about a reference trajectory in the standard fashion, or linearization about the current state estimate as with the extended Kalman filter. To investigate the feasibility of the decentralized integration with the GSFC algorithm, an existing centralized LQG design for a single spacecraft orbit control problem is adapted to the decentralized framework while using the GSFC algorithm's state transition matrices and framework. The existing GSFC design uses both reference trajectories of each spacecraft in formation and by appropriate choice of coordinates and simplified measurement modeling is formulated as a linear time-invariant system. Results for improvements to the GSFC algorithm and a multiple satellite formation will be addressed. The goal of this investigation is to progressively relax the assumptions that result in linear time-invariance, ultimately to the point of linearization of the non-linear dynamics about the current state estimate as in the extended Kalman filter. An assessment will then be made about the feasibility of the decentralized approach to the realistic formation flying application of the EO-1/Landsat 7 formation flying experiment.

Folta, David C.↗

NASA NTRS2001

On the Application of Hilbert Spectral Analysis for Climate Studies

The Hilbert spectral analysis (Huang et al, 1998, Proceedings of the Royal Society of London, A 454, pp 903-995) consisted of two steps: First, the data has to be reduced into a finite number of Intrinsic Mode Function by the Empirical Mode Decomposition method, then the resulting Intrinsic Mode Functions are converted to time-frequency-energy distribution through Hilbert transform. In this approach, the Empirical Mode Functions served as the basis functions with which the data is expanded. This basis function is adaptive, and the decomposition is nonlinear. Furthermore, as the Hilbert transform is a singular transform, it retains a high degree of local information. The instantaneous frequency is determined by differentiation of the phase function; therefore, there is no restriction of the 'uncertainty principle' for all the time-frequency analysis resulting from a priori basis approach. With the adaptive basis and the instantaneous frequency, the Hilbert Spectral analysis can represent data from nonlinear and nonstationary processes without resorting to the harmonics. Another advantage of using instantaneous frequency is the ability to find out frequency from limited length of data, which is a critical problem in climate studies. As the processes driving the climate changes could be both nonlinear and nonstationary, the Hilbert Spectral Analysis could be of great use in examining the underlying mechanisms. A preliminary study based on the length of day data will be presented as example for the application of the Hilbert Spectral Analysis for climate study.

Huang, Norden E.↗

NASA NTRS2005

Control Algorithms and Simulated Environment Developed and Tested for Multiagent Robotics for Autonomous Inspection of Propulsion Systems

The NASA Glenn Research Center and academic partners are developing advanced multiagent robotic control algorithms that will enable the autonomous inspection and repair of future propulsion systems. In this application, on-wing engine inspections will be performed autonomously by large groups of cooperative miniature robots that will traverse the surfaces of engine components to search for damage. The eventual goal is to replace manual engine inspections that require expensive and time-consuming full engine teardowns and allow the early detection of problems that would otherwise result in catastrophic component failures. As a preliminary step toward the long-term realization of a practical working system, researchers are developing the technology to implement a proof-of-concept testbed demonstration. In a multiagent system, the individual agents are generally programmed with relatively simple controllers that define a limited set of behaviors. However, these behaviors are designed in such a way that, through the localized interaction among individual agents and between the agents and the environment, they result in self-organized, emergent group behavior that can solve a given complex problem, such as cooperative inspection. One advantage to the multiagent approach is that it allows for robustness and fault tolerance through redundancy in task handling. In addition, the relatively simple agent controllers demand minimal computational capability, which in turn allows for greater miniaturization of the robotic agents.

Wong, Edmond↗

NASA NTRS2011

Tracking Object Existence From an Autonomous Patrol Vehicle

An autonomous vehicle patrols a large region, during which an algorithm receives measurements of detected potential objects within its sensor range. The goal of the algorithm is to track all objects in the region over time. This problem differs from traditional multi-target tracking scenarios because the region of interest is much larger than the sensor range and relies on the movement of the sensor through this region for coverage. The goal is to know whether anything has changed between visits to the same location. In particular, two kinds of alert conditions must be detected: (1) a previously detected object has disappeared and (2) a new object has appeared in a location already checked. For the time an object is within sensor range, the object can be assumed to remain stationary, changing position only between visits. The problem is difficult because the upstream object detection processing is likely to make many errors, resulting in heavy clutter (false positives) and missed detections (false negatives), and because only noisy, bearings-only measurements are available. This work has three main goals: (1) Associate incoming measurements with known objects or mark them as new objects or false positives, as appropriate. For this, a multiple hypothesis tracker was adapted to this scenario. (2) Localize the objects using multiple bearings-only measurements to provide estimates of global position (e.g., latitude and longitude). A nonlinear Kalman filter extension provides these 2D position estimates using the 1D measurements. (3) Calculate the probability that a suspected object truly exists (in the estimated position), and determine whether alert conditions have been triggered (for new objects or disappeared objects). The concept of a probability of existence was created, and a new Bayesian method for updating this probability at each time step was developed. A probabilistic multiple hypothesis approach is chosen because of its superiority in handling the uncertainty arising from errors in sensors and upstream processes. However, traditional target tracking methods typically assume a stationary detection volume of interest, whereas in this case, one must make adjustments for being able to see only a small portion of the region of interest and understand when an alert situation has occurred. To track object existence inside and outside the vehicle's sensor range, a probability of existence was defined for each hypothesized object, and this value was updated at every time step in a Bayesian manner based on expected characteristics of the sensor and object and whether that object has been detected in the most recent time step. Then, this value feeds into a sequential probability ratio test (SPRT) to determine the status of the object (suspected, confirmed, or deleted). Alerts are sent upon selected status transitions. Additionally, in order to track objects that move in and out of sensor range and update the probability of existence appropriately a variable probability detection has been defined and the hypothesis probability equations have been re-derived to accommodate this change. Unsupervised object tracking is a pervasive issue in automated perception systems. This work could apply to any mobile platform (ground vehicle, sea vessel, air vehicle, or orbiter) that intermittently revisits regions of interest and needs to determine whether anything interesting has changed.

Wolf, Michael↗

NASA NTRS2022

A preliminary assessment of olivine phenocrysts from the monogenetic basalt of the McCartys Flow, Zuni-Bandera Volcanic Field, New Mexico

Monogenetic small-volume basaltic volcanoes are the most abundant subaerial volcanic landforms on Earth but are some of the most poorly understood systems. Their short durations, small volumes, and lack of recurrence make monitoring and hazard assessment difficult. The Zuni-Bandera volcanic field in western New Mexico contains small-volume basaltic centers erupting tholeiitic to alkalic basalts. Evidence shows no correlation of magma composition with eruption age, location, or volumetric output, prompting questions about the influence of magma ascent rates, magma storage conditions, and mantle source characteristics on lava compositions. Here, we present olivine major and minor element mineral chemistry from the 3200-year-old McCartys Flow, the youngest tholeiite basalt in the volcanic field. Olivine displays four phenocryst types with unique textures and major and minor element compositions. Multiple olivine types co-exist at the thin section scale. Major and minor element diffusion at frozen melt–phenocryst interfaces was modeled, revealing magma residence times ranging from 3–9 months. Type 3 olivine phenocrysts require step function initial conditions and record diffusion re-equilibration followed by magma mixing. These profiles indicate the magma resided in the reservoir for 10–15 years and accumulated from multiple batches of mixed magmas less than 10 days before the eruption. Our results show that primitive magmas in small-volume monogenetic volcanoes have complex lithospheric magmatic histories and stored in magma bodies influenced by an open system to develop different local chemical environments.

Gary S. Michelfelder↗

NASA NTRS2018

Data, Meet Compute: NASA's Cumulus Ingest Architecture

NASA's Earth Observing System Data and Information System (EOSDIS) houses nearly 30PBs of critical Earth Science data and with upcoming missions is expected to balloon to between 200PBs-300PBs over the next seven years. In addition to the massive increase in data collected, researchers and application developers want more and faster access - enabling complex visualizations, long time-series analysis, and cross dataset research without needing to copy and manage massive amounts of data locally. NASA has looked to the cloud to address these needs, building its Cumulus system to manage the ingest of diverse data in a wide variety of formats into the cloud. In this talk, we look at what Cumulus is from a high level and then take a deep dive into how it manages complexity and versioning associated with multiple AWS Lambda and ECS microservices communicating through AWS Step Functions across several disparate installations

EOSDIS↗

NASA NTRS1985

A 3-dimensional Navier-Stokes-Euler code for blunt-body flow computations

The shock-layer flowfield is obtained with or without viscous and heat-conducting dissipations from the conservative laws of fluid dynamics equations using a shock-fitting implicity finite-difference technique. The governing equations are cast in curvilinear-orthogonal coordinates and transformed to the domain between the shock and the body. Another set of equations is used for the singular coordinate axis, which, together with a cone generator away from the stagnation point, encloses the computation domain. A time-dependent alternating direction implicit factorization technique is applied to integrate the equations with local-time increment until a steady solution is reached. The shock location is updated after the flowfield computation, but the wall conditions are implemented into the implicit procedure. Innovative procedures are introduced to define the initial flowfield, to treat both perfect and equilibrium gases, to advance the solution on a coarse-to-fine grid sequence, and to start viscous flow computations from their corresponding inviscid solutions. The results are obtained from a grid no greater than 28 by 18 by 7 and converged within 300 integration steps. They are of sufficient accuracy to start parabolized Navier-Stokes or Euler calculations beyond the nose region, to compare with flight and wind-tunnel data, and to evaluate conceptual designs of reentry spacecraft.

Li, C. P.↗

NASA NTRS2009

Tracking Debris Shed by a Space-Shuttle Launch Vehicle

The DEBRIS software predicts the trajectories of debris particles shed by a space-shuttle launch vehicle during ascent, to aid in assessing potential harm to the space-shuttle orbiter and crew. The user specifies the location of release and other initial conditions for a debris particle. DEBRIS tracks the particle within an overset grid system by means of a computational fluid dynamics (CFD) simulation of the local flow field and a ballistic simulation that takes account of the mass of the particle and its aerodynamic properties in the flow field. The computed particle trajectory is stored in a file to be post-processed by other software for viewing and analyzing the trajectory. DEBRIS supplants a prior debris tracking code that took .15 minutes to calculate a single particle trajectory: DEBRIS can calculate 1,000 trajectories in .20 seconds on a desktop computer. Other improvements over the prior code include adaptive time-stepping to ensure accuracy, forcing at least one step per grid cell to ensure resolution of all CFD-resolved flow features, ability to simulate rebound of debris from surfaces, extensive error checking, a builtin suite of test cases, and dynamic allocation of memory.

Stuart, Phillip C.↗

NASA NTRSDate not supplied

The Impact of Drifting Orbits on the Monthly Regional TOA Flux Assuming Constant Meteorology

The NASA Clouds and the Earth's Radiant Energy System (CERES) gridded Single Scanner Footprint (SSF1deg) product provides TOA SW and LW monthly 1° regional all-sky fluxes, which are used to monitor the Earth’s energy balance. The CERES long-term climate data record relies on Terra and Aqua satellite sun-synchronous orbits that are maintained at 10:30 and 1:30 local equator crossing times (LECT), respectively. The Terra and Aqua satellites are expected to drift outside of their respective LECT during mid 2022. Both Terra and Aqua will drift over several years towards sunrise and sunset, respectively, and eventually will be deorbited. The CERES SSF1deg product monthly regional fluxes are based on the well calibrated and stable CERES observed fluxes and are temporally interpolated assuming constant meteorology between measurements to resolve the regional diurnal flux cycle to obtain a daily averaged flux. The drifting orbits may impact the monthly regional TOA flux over regions with systematic diurnal cycles, because the observations will shift in local time. The CERES project would like to determine the maximum Terra and Aqua LECT time shift before the monthly regional fluxes are diurnally impacted and become unreliable for long-term climate monitoring. To determine the impact of the drifting orbits on the SSF1deg monthly regional fluxes,15-minute Geostationary Earth Radiation Budget (GERB) broadband observed fluxes over the Meteosat geostationary satellite domain (±60° in longitude and latitude) are used as proxy CERES observations. The drifting orbit sampling pattern is achieved by simply incrementing the observation time by steps of 15 minutes from the CERES footprint time. For each 15-minute time interval, the CERES observed fluxes are replaced by the GERB observed fluxes. The-15 minute incremented monthly regional fluxes based on constant meteorology are compared to the reference 10:30 and 1:30 LECT fluxes. Regions with systematic diurnal cycles, include morning maritime stratus, where the clouds dissipate during the morning, and land afternoon convection, where clouds increase in the afternoon will impact the regional flux differences. Based on January and July 2010 GERB data, even a 15-minute LECT change caused regional monthly flux differences that would impact long term regional trend analysis. Results will be shown at the conference.

David Doelling↗

NASA NTRSDate not supplied

The Orbital Drift Impact on the Monthly Regional TOA Flux Assuming Constant Meteorology

D R Doelling↗

NASA NTRS1989

Application of unstructured grid methods to steady and unsteady aerodynamic problems

The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.

Batina, John T.↗

NASA NTRSDate not supplied

Design and Analysis of Shape Memory Spring Tires for Martian and Lunar Rover Vehicles

Shape memory alloys (SMAs) have played an important role in various innovative engineering and medical applications, such as aerospace actuators, vibration damping devices, and coronary stents. In applications, shape memory alloys are commonly utilized in two fundamentally different ways: (i) making use of the superelasticity/ pseudoelasticity (SE/PE) phenomena, as in applications in biomedical engineering, and (ii) taking advantages of the shape memory effect (SME), as is used for actuators. Their ability to act in such vastly different capacities is mainly due to their unique capability to recover large amounts of deformation produced by either applied stresses or temperature changes. One recent emerging application in use of SMAs has been in the area of non-pneumatic tire designs for Martian or Lunar roving vehicles. These vehicles require tires that are capable of traversing rugged terrain while withstanding extreme temperatures and atmospheric conditions. Inspired by the flexible wire mesh tires used on three Lunar Roving Vehicle (LRV) missions to the Moon on Apollo 15, 16, and 17, a new compliant tire technology was developed by the NASA Glenn Research Center (GRC) and Goodyear Tire & Rubber, known as the Spring Tire. The Spring Tire consists of several hundred coiled springs woven into a flexible mesh and formed into the shape of a tire. Like the LRV wire mesh tires, the original Spring Tires were made from spring steel and were prone to permanent deformation when undergoing high localized loads. Later, a new iteration of this technology was invented, known as the ‘Superelastic Tire’. This new technology incorporated the use of superelastic SMA springs, which could effectively undergo approximately 30 times more reversible deformation than the steel spring. It also provided even greater durability and allowed for more flexibility in design, such as the use of other structural forms to reduce mass or increase load carrying capacity. Because of the unique nature of both the SMA material and the complex interactions between the springs, designing Spring Tires for a specific application requires extensive effort. Historically, design decisions have relied on full-scale empirical testing; however, this is very expensive and time consuming, especially when multiple iterations are needed. Therefore, developing a large-scale, robust, and predictive numerical model entailing complex spring interactions and the shape memory material behavior within a tire construct is the first essential step towards a successful design program. The current work focuses on implementation of the user-defined Shape Memory Alloy (SMA) model, otherwise known as SMA-GVIPs, in the Finite Element analysis (FEA) program ABAQUS for large-scale simulations of the GRC-developed Spring Tires made of SMA. The novelty of this work lies in the thorough, detail-oriented, and computationally efficient finite element analysis of full-scale SMA tires. A well-thought material characterization plan followed by model validation and a systemic sensitivity study on SMA tires has never been reported in the previous literature. The main objective of this study is to help the team improve and optimize the structural design of the SMA tires through in-depth numerical analysis and sensitivity studies. Various design variables (wire diameter, coil diameter, pitch, bead angle, and number of springs) were varied to study their influence on the global load-displacement response of the tire construct. A detailed investigation of the three-dimensional stress states was also carried out to enhance our understanding of the local changes as the tire goes through global deformation. It was concluded that a robust numerical model with a good predictive capability, together with a thoughtfully crafted sensitivity study can result in improved design iterations required to reach a desired tire performance while, significantly reducing manufacturing, labor and testing expenses. A summary of the Finite Element (FE) model construction will be presented together with a description of the user-defined SMA model, characterization process, experimental results, model validation, and numerical sensitivity study results.

shape memory alloys↗