Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chaotic dynamics”

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.

193 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Atmospheric Dynamics of the Outer Planets

The giant planets-Jupiter, Saturn, Uranus, and Neptune-are fluid objects. The winds are powered by absorbed sunlight, as on earth, and by internal heat left over from planetary formation. The main constituents of the atmospheres are hydrogen and helium. The clouds are made of ammonia, hydrogen sulphide, and water. All four giant planets are banded, with multiple zonal jet streams. Even Uranus, whose spin axis is tipped by 98deg relative to the orbit axis, shows latitudinal banding and zonal jets. Equator-to-pole temperature differences are close to zero. Wind speeds are larger than on earth and do not decrease with distance from the sun. Although the power/area at Neptune is only 1/20 that at Jupiter, the winds at Neptune are three times stronger. Stable vortices like the Great Red Spot of Jupiter and similar spots on Neptune come in all size ranges and exhibit a variety of behaviours including merging, orbiting, filament ejection, and oscillating in both shape and position. At least at cloud-top levels, 90% of the long-lived vortices are anticyclonic and sit in anticyclonic shear zones. Features in the cyclonic zones tend to be chaotic, with lifetimes of several days or less. These mesoscale eddies tend to have lightning in them, which suggests that they get their energy from moist convection. The rate of conversion of eddy kinetic energy into kinetic energy of the zonal jets is more than 10% of the power/area radiated by Jupiter. This fraction is more than an order of magnitude larger than on earth. Several lines of evidence now indicate that the winds at cloud-top levels are the surface manifestation of deep-rooted motions that extend into the interior and are presumably driven by internal heat.

Ingersoll, Andrew P.↗
NASA NTRS2022

Computing Sensitivities in Evolutionary Systems: A Real-time Reduced Order Modeling Strategy

We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first order optimality condition of the variational principle leads to a system of closed form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by the forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any input/output load and can be used in applications in which real-time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rössler system, (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto--Sivashinsky equation, and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.

Reduced order model↗
NASA NTRS2006

Designing capture trajectories to unstable periodic orbits around Europa

The hostile environment of third body perturbations restricts a mission designer's ability to find well-behaved reproducible capture trajectories when dealing with limited control authority as is typical with low-thrust missions. The approach outlined in this paper confronts this shortcoming by utilizing dynamical systems theory and an extensive preexisting database of Restricted Three Body Problem (RTBP) periodic orbits. The stable manifolds of unstable periodic orbits are utilized to attract a spacecraft towards Europa. By selecting an appropriate periodic orbit, a mission designer can control important characteristics of the captured state including stability, minimum altitudes, characteristic inclinations, and characteristic radii among others. Several free parameters are optimized in the non-trivial mapping from the RTBP to a more realistic model. Although the ephemeris capture orbit is ballistic by design, low-thrust is used to target the state that leads to the capture orbit, control the spacecraft after arriving on the unstable quasi-periodic orbit, and begin the spiral down towards the science orbit. The approach allows a mission designer to directly target fuel efficient captures at Europa in an ephemeris model. Furthermore, it provides structure and controllability to the design of capture trajectories that reside in a chaotic environment.

low-thrust↗
NASA NTRS1996

Cascade Error Projection: A Learning Algorithm for Hardware Implementation

In this paper, we workout a detailed mathematical analysis for a new learning algorithm termed Cascade Error Projection (CEP) and a general learning frame work. This frame work can be used to obtain the cascade correlation learning algorithm by choosing a particular set of parameters. Furthermore, CEP learning algorithm is operated only on one layer, whereas the other set of weights can be calculated deterministically. In association with the dynamical stepsize change concept to convert the weight update from infinite space into a finite space, the relation between the current stepsize and the previous energy level is also given and the estimation procedure for optimal stepsize is used for validation of our proposed technique. The weight values of zero are used for starting the learning for every layer, and a single hidden unit is applied instead of using a pool of candidate hidden units similar to cascade correlation scheme. Therefore, simplicity in hardware implementation is also obtained. Furthermore, this analysis allows us to select from other methods (such as the conjugate gradient descent or the Newton's second order) one of which will be a good candidate for the learning technique. The choice of learning technique depends on the constraints of the problem (e.g., speed, performance, and hardware implementation); one technique may be more suitable than others. Moreover, for a discrete weight space, the theoretical analysis presents the capability of learning with limited weight quantization. Finally, 5- to 8-bit parity and chaotic time series prediction problems are investigated; the simulation results demonstrate that 4-bit or more weight quantization is sufficient for learning neural network using CEP. In addition, it is demonstrated that this technique is able to compensate for less bit weight resolution by incorporating additional hidden units. However, generation result may suffer somewhat with lower bit weight quantization.

Duong, Tuan A.↗
NASA NTRS1992

Long-term orbital evolution of short-period comets found in Project Cosmo-DICE

Orbital evolutions of about 160 short-period (SP) comets are numerically integrated for 4400 years in the framework of a realistic dynamical model. By the round-trip error in closure test, a reliable time space of the integrated orbits is estimated for each comet. Majority of the SP comets with their Tisserand's constant(J) between 2.8 and 3.1 are found to evolve within the past 1000-2000 years from the orbits whose perihelia are near the Jovian orbit to the orbits with perihelia of 1-2 AU. This evolution is much more rapid than that expected from Monte Carlo simulations based on symmetric distribution of planetary perturbations, thus suggesting that asymmetry of perturbation distribution play an important role in cometary evolution. Several comets are shown to evolve from the near-Saturn orbits and then to be handed over under the control of Jupiter. We also find that a few comets were captured from long-period orbits (a = 75-125 AU) via only a few close encounters with Jupiter. It is confirmed that the captured SP comets of low-inclination with 2.7 less than J less than 3.1 show more or less strong chaotic behavior. On the other hand, comets with longer orbital period and/or of high inclination reveal slow or quasi-periodic orbital evolution.

Nakamura, Tsuko↗
NASA NTRS1994

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗
NASA NTRS1988

Dynamic analysis of nonlinear rotor-housing systems

Nonlinear analysis methods are developed which will enable the reliable prediction of the dynamic behavior of the space shuttle main engine (SSME) turbopumps in the presence of bearing clearances and other local nonlinearities. A computationally efficient convolution method, based on discretized Duhamel and transition matrix integral formulations, is developed for the transient analysis. In the formulation, the coupling forces due to the nonlinearities are treated as external forces acting on the coupled subsystems. Iteration is utilized to determine their magnitudes at each time increment. The method is applied to a nonlinear generic model of the high pressure oxygen turbopump (HPOTP). As compared to the fourth order Runge-Kutta numerical integration methods, the convolution approach proved to be more accurate and more highly efficient. For determining the nonlinear, steady-state periodic responses, an incremental harmonic balance method was also developed. The method was successfully used to determine dominantly harmonic and subharmonic responses fo the HPOTP generic model with bearing clearances. A reduction method similar to the impedance formulation utilized with linear systems is used to reduce the housing-rotor models to their coordinates at the bearing clearances. Recommendations are included for further development of the method, for extending the analysis to aperiodic and chaotic regimes and for conducting critical parameteric studies of the nonlinear response of the current SSME turbopumps.

Noah, Sherif T.↗
NASA NTRS2001

Seismoball: A Small Europa Orbiter Drop-Off Probe for Early Exploration of the Europan Surface

Recent magnetometry data received from Galileo indicate that the most likely explanation for the magnetic signature there is indeed a global conducting layer below the surface. This conducting layer is well matched by a salty, mineral rich strata beneath the Europan ice crust or a salt water ocean. Galileo imaging results show a variety of terrain types thought to contain young material; for example, lineaments, chaotic terrain, and eruption features. Additionally, Galileo images have shown indications of areas of up-welling where subsurface material periodically gets pushed to the surface due to the forces of fracturing, butting, and refreezing of the ice sheet. While Europa Orbiter will provide close-flyby high resolution images, as well as magnetometry, spectroscopy and other remote sensing data of the surface, it will not be able to provide essential engineering data like surface hardness and surface ice structure needed to support eventual landed missions. Additionally, ice chemical composition at microscopic scales can only be studied in detail through in situ instrumentation. Seismoball is a small probe designed to be injected into a surface intersect orbit around Europa. Using small reverse thrusters, the probe will be capable of nulling the high horizontal injection velocity as it approaches the 2 km surface injection altitude, thus allowing it to fall to the surface at an impact velocity of < 100m/sec (much less than the DS-2 impact velocities). The external breakaway thruster structure and crushable exterior shell absorb the impact energy while allowing the science instrument suite to remain intact. JPL has already started analyzing the entry dynamics and designing/building a small, low mass probe which will withstand the impact g-forces and fit as a 'carry-on' on board the Europa Orbiter. The probe will carry a suite of 5-6 micro-instruments for imaging the surface (both microscopic and far-field), surface and shallow subsurface ice temperatures, surface hardness, crustal dynamics and periodicity, and compositional chemistry. If selected, this flight development activity will provide a unique science opportunity and adjunct to the primary Orbiter science mission. The final flight system will be designed to accommodate orbiter mass, volume, and power interface constraints, as well as entry dynamics, g-load mitigation, and arbitrary landing orientation.

Tamppari, L.↗
NASA NTRS1994

Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties

A semi-implicit fully spectral collocation method for the simulation of three-dimensional mantle convection with depth-dependent thermo-dynamic and transport properties is presented. The variable property Navier-Stokes equation expressed in terms of the primitive variable velocity and pressure is solved with the mass continuity and temperature equations. The periodic horizontal boundary conditions allow a Fourier expansion for the two horizontal directions. The stress-free, impermeable isothermal boundary conditions along with the depth dependent coefficients are handled with a Chebyshev expansion in the vertical direction. In the limit of an infinite Prandtl number appropriate to mantle convection, the inertial terms in the momentum equation are unimportant. In this case an explicit solution of a Poisson equation for pressure can be avoided; instead a fourth-order equation for vertical velocity can be solved. Simultaneous imposition of both impermeable and continuity boundary conditions during the vertical velocity evaluation is discussed. The pressure distributions on the top and bottom bounding planes were determined by means of an influence matrix technique. The numerical method employed here avoids time-splitting errors and enforces velocity boundary conditions and continuity over the entire domain, including the boundaries, to machine accuracy. Strongly time-dependent three-dimensional solutions up to a surface Rayleigh number of 1 x 10(exp 7) have been obtained. Strong upwellings, pulsating chaotically, are formed by the collective merging of cylindrical plumes.

Balachandar, S.↗
NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗
NASA NTRS1992

Hydrogen leak detection in the Space Shuttle

This study focuses on a helium gas jet flowing into room air. Measurements of helium concentration and velocity in the jet-air mixture are reported. The objective is to learn about jet characteristics so that dynamically similar hydrogen leaks may be located in the Space Shuttle. The hazardous gas detection system (HGDS) in the mobile launch pad uses mass spectrometers to monitor the shuttle environment for leaks. The mass spectrometers are fed by long sample tubes which draw gas from the payload bay, mid body, aft engine compartment and external tank. The overall purpose of this study is to improve the HGDS especially in its potential for locating hydrogen leaks. A rapid-response leak detection experiment was designed, built, and tested, following on the work done in this program last summer. The apparatus included a Perkin Elmer MGA-1200 mass spectrometer and air velocity transducer, both monitored by a Macintosh IIFX computer using LabVIEW software. A jet of helium flowing into the lab air simulated a gas leak. Steady helium or hydrogen-nitrogen jets were logged for concentration and velocity, and the power spectral density of each was computed. Last year, large eddies and vortices were visually seen with Schlieren imaging, and they were detected in the time plots of the various instruments. The response time of the MGA-1200 was found in the range of 0.05 to 0.1 sec. Pulsed concentration waves were clearly detected at 25 cycles per sec by spectral analysis of MGA data. No peaks were detected in the power spectrum, so in the present study, 10 Hz bandwidth-averaged power levels were examined at regular frequency intervals. The practical consequences of last year's study are as follows: sampling frequency should be increased above the present rate of 1 sample per second so that transients could be observed and analyzed with frequency response methods. Many more experiments and conditions were observed in this second summer, including the effects of orifice diameter, jet velocity, sample tube design, radial effects, vertical flow, and low hydrogen concentrations (1 percent). A frequent observation was that the power spectrum, calculated from the Fourier transform of concentration fluctuations, gives a separate piece of information from concentration. Many of the tests suggest that power is high where mixing occurs at the helium-air interface. This fact is apparently independent of the concentration level, which could be high or low, but depends on the sample location relative to the jet (leak) origin, whereas high concentration may be due to a strong leak far away or a small leak close to the sample tube. If the power is low for any concentration level, this would signify helium is arriving at the sample tube by diffusion, not chaotic mixing caused by the jet interaction with air. The practical result is to propose a modification of the HGDL mass spectrometer data sampling and software so that sampling rates could be capable of observing at least 25 Hz fluctuations.

Barile, Ronald G↗
NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗
NASA NTRS1998

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗